Anomalous Microwave Emission: Theory, Modeling, and Observations
View this Special IssueResearch Article  Open Access
Spectrum of the Anomalous Microwave Emission in the North Celestial Pole with WMAP 7Year Data
Abstract
We estimate the frequency spectrum of the diffuse anomalous microwave emission (AME) on the North Celestial Pole (NCP) region of the sky with the Correlated Component Analysis (CCA) component separation method applied to WMAP 7yr data. The NCP is a suitable region for this analysis because the AME is weakly contaminated by synchrotron and freefree emission. By modeling the AME component as a peaked spectrum we estimate the peak frequency to be GHz, in agreement with previous analyses which favored < 23 GHz. The ability of our method to correctly recover the position of the peak is verified through simulations. We compare the estimated AME spectrum with theoretical spinning dust models to constrain the hydrogen density . The best results are obtained with densities around 0.2–0.3 cm^{−3}, typical of warm ionised medium (WIM) to warm neutral medium (WNM) conditions. The degeneracy with the gas temperature prevents an accurate determination of , especially for low hydrogen ionization fractions, where densities of a few cm^{−3} are also allowed.
1. Introduction
The anomalous microwave emission (AME) component is highly correlated with the far infrared dust emission [1–10] and is believed to be the electric dipole radiation from small spinning dust grains [11]. Spinning dust models predict a peaked spectrum ranging 10–150 GHz depending on the local physical conditions [12, 13]. Probing the peak of the emission enables us to compare models with observations and is the best way to distinguish the AME from synchrotron emission and freefree emission, both having a powerlaw spectrum.
The AME has been studied in individual dust clouds associated with reflection nebulae, molecular clouds, photodissociation regions, and HII regions [14, 15]. The best examples of peaked AME spectra are probably the Perseus and Ophiuchi molecular clouds, which have been probed with high accuracy thanks to the availability of many different datasets. In these regions the derived peak is at ~30 GHz [14].
The study of the AME in more diffuse regions essentially becomes a component separation problem, which is often complicated by the lack of lowfrequency data covering large sky areas. Diffuse dustcorrelated emission has been detected with COBEDMR [2] and WMAP data [4, 7, 16–18]; however the details of the spectrum are still unclear. Templatefitting analyses show that the dustcorrelated emission between 20 and 60 GHz is well described by a powerlaw [2, 4, 7, 18]. Given the error bars, a peaked spectrum is not ruled out, but a low peak frequency (<23 GHz) is favored. Alternatively, the observed spectrum could result from the superposition of multiple peaked components. The component separation analysis done by [17] found that a peaked AME model gives better results in terms of CMB cleaning, however no estimation of the AME peak frequency was attempted.
In this work we address the estimation of the peak frequency of the diffuse AME on WMAP 7yr data complemented by ancillary data with the Correlated Component Analysis (CCA, [19]) component separation method. We consider the North Celestial Pole (NCP) region of the sky, centered on Galactic coordinates . As already noted by [4], this region of the sky is particularly suited for this analysis since there is significant thermal and anomalous dust emission, while synchrotron emission and freefree emission (traced by the 408 MHz from [20] and the H map from [21]) are faint.
2. Description of the Method
In the following we describe the principles of operation of the harmonicdomain CCA; for further details we refer the reader to [19]. The CCA [19, 22] exploits secondorder statistics to estimate the frequency scaling of the components from the statistics of data and noise. For each position on the sky (each pixel) we write our data model as where is a vector whose elements contain the different components (CMB, freefree emission, synchrotron emission, thermal dust emission, AME); and are vectors too, each element containing, respectively, the data and the instrumental noise for each frequency. The diagonal matrix contains the instrumental beams for all frequency channels and denotes convolution. The matrix is called the mixing matrix and it contains the frequency scaling of the components for all the considered frequencies.
By translating (1) in the harmonic domain, the convolution becomes a multiplication and the data model becomes a linear mixture. For each transformed mode () we can write where , and are the transforms of , and , respectively, and is the transform of the matrix . The crossspectra of the data , sources and noise, , are related by where the dagger superscript denotes the adjoint matrix.
The lefthand side of (3) can be estimated from the data for a suitable set of spectral bins ; this is used by CCA to estimate the mixing matrix and the source crossspectra on the righthand side of the equation. Due to the scaling ambiguity of the problem, the mixing matrix is normalized at a reference frequency. To reduce the number of unknowns, the mixing matrix is parametrized through a parameter vector (such that ), by adopting suitable fitting relations for the spectra of the astrophysical components (as detailed in the next section).
By using twodimensional discrete Fourier transforms it is possible to apply the harmonicdomain CCA to square sky patches. This approach is advisable as the mixing matrix varies on the sky. The HEALPix ([24]) data on the sphere are projected on the plane tangential to the center of the patch and regridded with a suitable number of bins in order to correctly sample the original resolution. The patch size is obtained as a tradeoff between the need of having uniform spectral properties of the foregrounds and enough statistics for a robust computation of the auto and crossspectra of the data. In the present analysis we use a patch size of . By using higher resolution data (e.g., Planck data) it is possible to reduce the patch size by carrying the analysis up to a larger multipole. We analyzed the sky patch centered on (shown in Figure 1) and verified the stability of the results for a sample of patches shifted in latitude and longitude up to 10°.
(a)
(b)
(c)
(d)
3. Description of the Analysis
We used the following datasets: (i)WMAP 7yr K, Ka, Q, V, and W bands [25]. All maps have been used at the original resolution except the K band, which has been smoothed to 1° resolution to reduce the effect of beam asymmetry; (ii)408 MHz map [20] to trace the synchrotron component; (iii)predicted freefree emission at 23 GHz based on the H map by [21] corrected for dust absorption with the E(BV) map by [26] assuming ; (iv)predicted dust emission at 94 GHz by [23].
The different resolution of the data maps is accounted for with the beam matrix in (3). In particular, the beams are deconvolved for the mixing matrix estimation. The noise properties for the WMAP maps have been computed by simulating different noise realizations starting from the input maps and values. For the maps used as foreground templates (408 MHz map, 23 GHz freefree map, and 94 GHz dust prediction) we assumed a Gaussian noise at the 10% level, which is much higher than the instrumental noise. This extra noise mimics the error on the template as a tracer of the true component, the 10% level being indicative. We verified, however, that the results are not sensitive to the exact value assumed. We modeled the data as a mixture of five components: CMB, synchrotron emission, thermal dust emission, freefree emission, and AME. Each frequency map contributes five elements to the mixing matrix, one for each component. The synchrotron, freefree and dust templates are used as frequency channels having only one mixing matrix entry, that of the corresponding component, different from zero.
It is worth noting that, compared with previous work, we are not assuming correlation between the AME and thermal dust emission and we are not exploiting any template for the AME. Though the AME and the thermal dust emission are significantly correlated, the correlation is not supposed to be perfect, as, according to spinning dust models, the AME traces the distribution of smaller dust grains (PAHs). Indeed, the intensity ratio of AME and thermal dust is found to vary in the sky by a factor of ~2 (e.g., [4]).
We estimated the frequency spectrum of the AME while assuming the spectra of the other components to be known. For CMB we used the usual blackbody law with temperature of 2.726 K ([27]); for freefree we used the model where is the Gaunt factor, which is responsible for the departure from a pure powerlaw behavior, and is the electron temperature in units of K. The thermal dust emission has been modeled as a greybody with temperature K and spectral index . This is consistent with the 1component dust model by [23]; Planck data, extending the observations up to 857 GHz, will allow a refinement of this model. For synchrotron we assumed a powerlaw scaling with spectral index .
For the AME we used the parametric relation proposed by [17] which is a parabola in the plane parametrized in terms of peak frequency and slope at 60 GHz . The CCA solves for the two parameters and , so the accuracy of the estimation of the AME peak frequency is related to the ability of the parametric relation to reproduce the true spectrum. We verified, however, that this parametric model can fit a wide range of spinning dust spectra. This is shown in Figure 2, comparing spinning dust spectra produced with the SpDust [12, 28] code with the results of the fit of (4) obtained by minimizing the for the set of frequencies considered in this work. The physical models that we consider are: warm neutral medium (WNM), cold neutral medium (CNM), warm ionised medium (WIM), and molecular cloud (MC). In general, the fits are accurate up to 61 GHz, while at 94 GHz the parametric relation may not be able to reproduce the input spectra in detail. This is a consequence of fitting complex spectra with only a few parameters, the fit being less accurate where the AME is weaker. The bestfit parameters that we obtain, reported in Table 1, vary significantly from one input model to another.

4. Test on Simulated Data
We simulated WMAP 7yr data by assuming monochromatic bandpasses centered at the central frequency of the realistic bands, Gaussian beams at the nominal values, Gaussian noise generated according to realistic (spatially varying) RMS. Our model of the sky consists of the following components: (i)CMB emission constrained by the bestfit power spectrum model to WMAP 7yr; (ii)synchrotron emission given by [20] scaled in frequency with a powerlaw model with a spatially varying synchrotron spectral index as modeled by [29]; (iii)freefree component modeled by the H map of [21] corrected for dust absorption with the E(BV) map from [26] with and scaled in frequency with a powerlaw with fixed spectral index of ; (iv)thermal dust emission modeled with the 100 m map from [26] scaled in frequency with the bestfit model by [23] which consists of two greybody laws having different temperatures and emissivity indices; (v)AME emission modeled by the E(BV) map from [26] with intensity at 23 GHz calibrated on the results by [4] for average intermediatelatitude conditions and scaled according to SpDust models.
It is worth noting that this sky model is more complex than the spectral model assumed in the component separation. The synchrotron spectral index is spatially varying and both the thermal dust AME spectra, though spatially constant, are not generated with the CCA parametric models. This has been done purposely, to reflect a more realistic situation. To test the ability of our pipeline to recover the AME spectrum in different situations, we adopted two SpDust models for the AME: one peaking at 26 GHz and the other peaking at 19 GHz.
We performed an estimate of the AME spectrum on the simulated data following the same procedure applied to the real data. The results are presented in Figure 3. For AME peaking at 26 GHz the recovery of the spectrum is very accurate; the result on the peak frequency is . The accuracy of the estimation is helped by the weak synchrotron and freefree contamination in this region of the sky (in a more general case the error on would be of a few GHz). With the AME peaking at 19 GHz the errors increase, but the recovery of the spectrum is still satisfactory; we obtain . Our pipeline is able to distinguish very clearly between the two input models and to correctly estimate the peak frequency of the AME in both cases.
(a)
(b)
5. Results
The results for the AME spectrum on the NCP are shown in Figure 4; the estimated peak frequency is . The diamonds with error bars are the AME mixing matrix elements at WMAP frequencies and related uncertainties derived from the results on and .
By comparing these data points with theoretical models produced with SpDust we can constrain the physical properties of the medium with the hypothesis of spinning dust emission. Some relevant parameters of the SpDust code are: hydrogen number density (cm^{−3}); gas temperature (K); intensity of the radiation field with respect to average Galactic conditions, , and hydrogen ionization fraction, . Advanced parameters describe the dust grains and define the abundance of other elements.
The hydrogen density is sensitive to the position of the peak and in general to the shape of the spectrum in the considered frequency range. It is also degenerate with other parameters, in particular the gas temperature . We investigated this degeneracy by sampling the 2dimensional  space with a grid approach while conditioning the remaining parameters. For each point in the parameter space we computed the likelihood as where is the standard chisquare. In computing the we did not consider the 94 GHz data point, which could be biased, as verified through simulations (see Section 4). We repeated the likelihood estimation for different choices of the parameters and (varying all those parameters simultaneously would be computationally too demanding). As the CCA outputs are in normalized units, we compared both the model and the data after scaling each of them by the sum of the intensities at 23, 33, 41 and 61 GHz.
The results of the  parameter estimation are shown in Figure 5. The colored lines show the likelihoods for for different values of as detailed in the legend, while the black lines show the likelihood for marginalized over in the range 100–10000 K. Different panels correspond to different choices for the other parameters: , , on the top left; , , on the top right; , , on the bottom left; and , on the bottom right. The marginalized statistics for are reported in Table 2 ( and being the lower and upper 1 error on ). The  degeneracy weakens the constraint on in particular on the highdensity side of the distribution. The degeneracy decreases with increasing hydrogen ionization fraction.

(a)
(b)
(c)
(d)
The best results are obtained for densities of 0.2–0.4 cm^{−3}, typical of WIM/WNM conditions. In Figure 6 we show the data compared to some SpDust models to exemplify the quality of the fit. The comparison of the red ( cm^{−3}, K, ) and magenta ( cm^{−3}, K, ) lines illustrates the effect of . The blue line corresponds to cm^{−3}, K, , ; compared to the magenta line, it illustrates the effect of increasing . Finally, the cyan ( cm^{−3}, K, ) and green ( cm^{−3}, K, ) lines show the effect of . The availability of data at frequencies of a few GHz will be very valuable to improve the constraints on such parameters.
To complement the physical description of the medium we computed the hydrogen column density (cm^{−2}) from the brightness at 100 m, , through cm^{−2}] [30]. We used the IRIS [31] band 4 (100 m) data and integrated the emission inside a beam centered in , corresponding to bright dust emission (see Figure 1), to compute a representative value for . This analysis yields cm^{−2}, which is a plausible value for the diffuse ISM environment at intermediate latitudes.
6. Conclusions
We have studied the spectrum of the diffuse AME with WMAP 7yr and ancillary data in the North Celestial Pole (NCP) region of the sky. In this region the AME dominates the lowfrequency emission, as both synchrotron and freefree are faint. Previous templatefitting analysis by [4] found that the AME spectrum in this region is consistent with a powerlaw above 23 GHz; the same would apply to most diffuse AME at intermediate latitudes [2, 4, 7, 16–18]. This favors a low peak frequency (GHz).
For our analysis we rely on the CCA component separation method, which exploits the data auto and crossspectra to estimate the frequency spectra of the components in terms of a set of spectral parameters. Our method models the AME as a peaked spectrum and fits for the peak frequency, , and the slope at 60 GHz, . We verified with simulations that we are able to correctly recover the AME spectrum and more specifically the peak frequency, even when it is below 23 GHz. We get GHz, which is both a confirmation and an improvement of the previous results. This result relies on the assumption that our parametric model for the AME is a good representation of the true spectrum, which has been verified for a wide range of theoretical spinning dust models.
Using the SpDust code, we linked the estimated spectrum to the local physical conditions with the hypothesis of spinning dust emission. We investigated in particular the hydrogen density ( [cm^{−3}]), which is sensitive to the position of the peak, and its degeneracy with the gas temperature [K]. The densities that we get are those typical of WIM/WNM conditions (0.2–0.4 cm^{−3}). For a low hydrogen ionization fraction (), densities up to a few cm^{−3} are allowed by the  degeneracy. Lower radiation fields () require slightly lower . By considering 100 m data for our region we obtain a hydrogen column density cm^{−2}. Overall, the recovered AME spectrum is found to be consistent with that predicted by spinning dust models for plausible physical conditions.
Acknowledgments
S. Ricciardi acknowledges support by MIUR through PRIN 2009 Grant no. 2009XZ54H2.
References
 G. Lagache, “The largescale anomalous microwave emission revisited by WMAP,” Astronomy and Astrophysics, vol. 405, no. 3, pp. 813–819, 2003. View at: Google Scholar
 A. J. Banday, C. Dickinson, R. D. Davies, R. J. Davis, and K. M. Górski, “Reappraising foreground contamination in the COBEDMR data,” Monthly Notices of the Royal Astronomical Society, vol. 345, no. 3, pp. 897–911, 2003. View at: Publisher Site  Google Scholar
 A. Kogut, “Spatial correlation between dust and diffuse ionized gas,” in American Astronomical Society Meeting Abstracts, vol. 28 of Bulletin of the American Astronomical Society, p. 1295, 1996. View at: Google Scholar
 R. D. Davies, C. Dickinson, A. J. Banday, T. R. Jaffe, K. M. Górski, and R. J. Davis, “A determination of the spectra of Galactic components observed by the Wilkinson Microwave Anisotropy Probe,” Monthly Notices of the Royal Astronomical Society, vol. 370, no. 3, pp. 1125–1139, 2006. View at: Publisher Site  Google Scholar
 N. Ysard, M. A. MivilleDeschênes, and L. Verstraete, “Probing the origin of the microwave anomalous foreground,” Astronomy and Astrophysics, vol. 509, no. 1, 4 pages, 2010. View at: Publisher Site  Google Scholar
 B. Gold, N. Odegard, J. L. Weiland et al., “Sevenyear wilkinson microwave anisotropy probe (WMAP) observations: galactic foreground emission,” Astrophysical Journal, vol. 192, no. 2, 2011. View at: Publisher Site  Google Scholar
 G. Dobler and D. P. Finkbeiner, “Extended anomalous foreground emission in the WMAP threeyear data,” Astrophysical Journal Letters, vol. 680, no. 2, pp. 1222–1234, 2008. View at: Publisher Site  Google Scholar
 A. De OliveiraCosta, M. Tegmark, R. D. Davies et al., “The quest for microwave foreground X,” Astrophysical Journal, vol. 606, no. 2, pp. L89–L92, 2004. View at: Publisher Site  Google Scholar
 D. P. Finkbeiner, G. I. Langston, and A. H. Minter, “Microwave interstellar medium emission in the green bank galactic plane survey: evidence for spinning dust,” Astrophysical Journal, vol. 617, no. 1, pp. 350–359, 2004. View at: Publisher Site  Google Scholar
 E. M. Leitch, A. C. S. Readhead, T. J. Pearson, and S. T. Myers, “An anomalous component of galactic emission,” Astrophysical Journal Letters, vol. 486, no. 1, pp. L23–L26, 1997. View at: Google Scholar
 B. T. Draine and A. Lazarian, “Electric dipole radiation from spinning dust grains,” Astrophysical Journal Letters, vol. 508, no. 1, pp. 157–179, 1998. View at: Publisher Site  Google Scholar
 Y. AliHaïmoud, C. M. Hirata, and C. Dickinson, “A refined model for spinning dust radiation,” Monthly Notices of the Royal Astronomical Society, vol. 395, pp. 1055–1078, 2009. View at: Google Scholar
 N. Ysard and L. Verstraete, “The longwavelength emission of interstellar PAHs: characterizing the spinning dust contribution,” Astronomy and Astrophysics, vol. 509, no. A12, 13 pages, 2010. View at: Publisher Site  Google Scholar
 Planck Collaboration, “Planck early results 20: new light on anomalous microwave emission from spinning dust grains,” Astronomy & Astrophysics, vol. 536, no. A20, 17 pages, 2011. View at: Publisher Site  Google Scholar
 R. GénovaSantos, R. Rebolo, J. A. RubiñoMartín, C. H. LópezCaraballo, and S. R. Hildebrandt, “Detection of anomalous microwave emission in the pleiades reflection nebula with wilkinson microwave anisotropy probe and the COSMOSOMAS experiment,” Astrophysical Journal Letters, vol. 743, p. 67, 2011. View at: Google Scholar
 M. A. MivilleDeschênes, N. Ysard, A. Lavabre et al., “Separation of anomalous and synchrotron emissions using WMAP polarization data,” Astronomy and Astrophysics, vol. 490, no. 3, pp. 1093–1102, 2008. View at: Publisher Site  Google Scholar
 A. Bonaldi, S. Ricciardi, S. Leach, F. Stivoli, C. Baccigalupi, and G. De Zotti, “WMAP 3yr data with correlated component analysis: anomalous emission and impact of component separation on the CMB power spectrum,” Monthly Notices of the Royal Astronomical Society, vol. 382, no. 4, pp. 1791–1803, 2007. View at: Publisher Site  Google Scholar
 T. Ghosh, A. J. Banday, T. Jaffe et al., “Foreground analysis using croscorrelation of external templates on the 7year Wilkinson Microwave Anisotropy Probe data,” Monthly Notices of the Royal Astronomical Society, vol. 422, no. 4, pp. 3617–3642, 2012. View at: Google Scholar
 S. Ricciardi, A. Bonaldi, P. Natoli et al., “Correlated component analysis for diffuse component separation with error estimation on simulated Planck polarization data,” Monthly Notices of the Royal Astronomical Society, vol. 406, no. 3, pp. 1644–1658, 2010. View at: Publisher Site  Google Scholar
 C. G. T. Haslam, C. J. Salter, H. Stoffel, and W. E. Wilson, “A 408 MHz allsky continuum survey. II—the atlas of contour maps,” Astronomy and Astrophysics, vol. 47, p. 1, 1982. View at: Google Scholar
 C. Dickinson, R. D. Davies, and R. J. Davis, “Towards a freefree template for CMB foregrounds,” Monthly Notices of the Royal Astronomical Society, vol. 341, no. 2, pp. 369–384, 2003. View at: Publisher Site  Google Scholar
 A. Bonaldi, L. Bedini, E. Salerno, C. Baccigalupi, and G. De Zotti, “Estimating the spectral indices of correlated astrophysical foregrounds by a secondorder statistical approach,” Monthly Notices of the Royal Astronomical Society, vol. 373, no. 1, pp. 271–279, 2006. View at: Publisher Site  Google Scholar
 D. P. Finkbeiner, M. Davis, and D. J. Schlegel, “Extrapolation of galactic dust emission at 100 microns to cosmic microwave background radiation frequencies using FIRAS,” Astrophysical Journal Letters, vol. 524, no. 2, pp. 867–886, 1999. View at: Google Scholar
 M. Górski, E. Hivon, A. J. Banday et al., “HEALPix: a framework for highresolution discretization and fast analysis of data distributed on the sphere,” Astrophysical Journal, vol. 622, no. 2, pp. 759–771, 2005. View at: Publisher Site  Google Scholar
 N. Jarosik, C. L. Bennett, J. Dunkley et al., “Sevenyear wilkinson microwave anisotropy probe (WMAP) observations: sky maps, systematic errors, and basic results,” Astrophysical Journal, vol. 192, no. 2, 2011. View at: Publisher Site  Google Scholar
 D. J. Schlegel, D. P. Finkbeiner, and M. Davis, “Maps of dust infrared emission for use in estimation of reddening and cosmic microwave background radiation foregrounds,” Astrophysical Journal Letters, vol. 500, no. 2, pp. 525–553, 1998. View at: Publisher Site  Google Scholar
 D. J. Fixsen, J. L. Weiland, S. Brodd et al., “Comparison of the COBE FIRAS and DIRBE calibrations,” Astrophysical Journal Letters, vol. 490, no. 2, pp. 482–487, 1997. View at: Publisher Site  Google Scholar
 K. Silsbee, Y. AliHaïmoud, and C. M. Hirata, “Spinning dust emission: the effect of rotation around a nonprincipal axis,” Monthly Notices of the Royal Astronomical Society, vol. 411, pp. 2750–2769, 2011. View at: Google Scholar
 G. Giardino, A. J. Banday, K. M. Górski, K. Bennett, J. L. Jonas, and J. Tauber, “Towards a model of fullsky galactic synchrotron intensity and linear polarisation: a reanalysis of the Parkes data,” Astronomy and Astrophysics, vol. 387, no. 1, pp. 82–97, 2002. View at: Google Scholar
 F. Boulanger, A. Abergel, J. P. Bernard et al., “The dust/gas correlation at high Galactic latitude,” Astronomy and Astrophysics, vol. 312, no. 1, pp. 256–262, 1996. View at: Google Scholar
 M. A. MivilleDeschênes and G. Lagache, “IRIS: a new generation of IRAS maps,” in Astronomical Society of the Pacific Conference Series, L. Armus and W. T. Reach, Eds., vol. 357 of Astronomical Society of the Pacific Conference Series, p. 167, 2006. View at: Google Scholar
Copyright
Copyright © 2012 Anna Bonaldi and Sara Ricciardi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.