# On the Limitations of the Anomalous Microwave Emission Emissivity

Authors: Christopher T. Tibbs; Roberta Paladini; Clive Dickinson
Journal: Advances in Astronomy
 (2012)

---

## Abstract

Many studies of anomalous microwave emission (AME) have computed an AME emissivity to compare the strength of the AME detected in different regions. Such a value is usually defined as the ratio between the intensity of the AME at 1 cm and the thermal dust emission at 100μm. However, as studies of Galactic dust emission have shown, the intensity of the thermal dust emission at 100 μm is strongly dependent on the dust temperature, which has severe implications for the AME emissivity defined in this way. In this work, we illustrate and quantify this effect and find that the AME emissivity decreases by a factor of 11.1 between dust temperatures of 20 and 30 K. We, therefore, conclude that computing the AME emissivity relative to the 100 μm emission does not allow for accurate comparisons between the AME observed in different environments. With this in mind, we investigate the use of other tracers of the dust emission with which to compute the AME emissivity and we ultimately conclude that, despite the difficulty in deriving its value, the column density of the dust would be the most suitable quantity with which to compute the AME emissivity.


---

## Body

## 1. Introduction

In recent years there has been growing evidence for the existence of a new component of microwave emission present in the interstellar medium (ISM). This emission appears to be spatially correlated with the dust in the ISM, although it produces an excess of emission with respect to the predicted thermal, vibrational dust emission at these frequencies (e.g., [1]), and as such this component is often described as anomalous. There have been only a handful of detections of this anomalous microwave emission (AME) component originating from both the diffuse ISM at mid-to-high latitudes (e.g., [2–4]) and specific Galactic sources such as the Perseus and ρ Oph molecular clouds [5–9], the dark clouds LDN1622 [10, 11], LDN1621 [12], and LDN1111 [13], and a variety of Hii regions [14–16]. Additionally, the first detection of AME in a star forming region in the external galaxy, NGC6946, was reported [17, 18], confirming that this emission mechanism is truly ubiquitous.Observations have shown that the AME occurs in the frequency range ~10–100 GHz and is highly correlated with the mid-infrared (IR) dust emission. It is this correlation with the mid-IR emission that led Draine and Lazarian [19] to postulate their model of spinning dust emission. In this model, the observed excess emission is produced by the very smallest dust grains (very small grains, VSGs or polycyclic aromatic hydrocarbons, PAHs), characterised by an electric dipole moment, rotating rapidly, resulting in the production of electric dipole radiation. The model predicts a well-defined peaked spectrum, rising below ~20–30 GHz and then falling off at higher frequencies as expected from a Boltzmann cutoff in grain rotation frequencies. The Draine and Lazarian [19] spinning dust model has been refined and updated [20–24], and current spinning dust models incorporate a variety of grain rotational excitation and damping processes: collisions with neutral and ionised gas particles, plasma drag (the interaction between the electric field of ions and the electric dipole moment of the dust grains), absorption and emission of a photon, the photoelectric effect, microwave emission, and the formation of H2 on the grain surface. Some of the more recent models incorporate additional features including dust grains that are not only spinning about their axis of greatest inertia [21, 24] and dust grains of irregular shape [23].Understanding how the AME varies between different phases of the ISM is extremely important in improving our understanding of the AME. To help compare the strength of AME detected in one region to another, some authors have computed an AME emissivity. Given the strong association between the AME and the dust grains in the ISM, the AME emissivity is usually defined as the ratio between the antenna temperature of the AME at wavelengths of 1 cm and the surface brightness of the thermal dust emission at 100μm. This quantity, first calculated as a cross-correlation coefficient over large areas of sky (e.g., [28, 29]), has now been computed for a variety of individual regions in which AME has been detected (see Table 1). However, the AME emissivity defined in this manner is a highly biased method for making these comparisons because the surface brightness at 100 μm is significantly dependent on the dust temperature.Table 1
Sample of AME emissivity values for AME detections in a variety of different environments computed relative to the intensity of the 100μm emission.
Source
AME emissivity (μK (MJy sr-1)-1)
Reference
Hii regions
6 southern Hii  regions
3.3
±
1.7
Dickinson et al. [14]
9 northern Hii  regions
3.9
±
0.8
Todorović et al. [15]
Pleiades
2.01
±
0.09
Génova-Santos et al. [25]
RCW49
13.6
±
4.2
Dickinson et al. [14]
RCW175
14.2
±
2.7
Tibbs et al. [16]
High latitudes
15 regions WMAP
11.2
±
1.5
Davies et al. [3]
All-sky WMAP
10.9
±
1.1
Davies et al. [3]
Molecular clouds
Perseus
15.7
±
0.3
Watson et al. [5]
Perseus A1
2.8
±
0.7
Tibbs et al. [7]
Perseus A2
16.4
±
4.1
Tibbs et al. [7]
Perseus A3
12.8
±
6.1
Tibbs et al. [7]
Perseus B
13.2
±
3.6
Tibbs et al. [7]
Perseus C
13.0
±
3.2
Tibbs et al. [7]
Dark clouds
LDN1621
18.1
±
4.4
Dickinson et al. [12]
LDN1622
21.3
±
0.6
Casassus et al. [11]In this work, we aim to (1) illustrate the weakness of such a definition and (2) present and discuss possible alternatives. In Section2, we point out why the commonly used AME emissivity definition is flawed and estimate the amplitude of the bias it introduces when comparing regions with different dust temperatures. In Section 3, we discuss more suitable definitions of the AME emissivity including the use of gas and dust column densities. Finally, we present our conclusions in Section 4.


## 2. The AME Emissivity Defined Relative to the 100μm Emission

The concept of an AME emissivity is particularly useful as it provides a normalisation of the strength of the AME detected in a variety of regions, and hence it facilitates a comparison between such AME detections. In Table1, we list AME emissivities defined relative to the 100 μm emission from the literature for a variety of different phases of the ISM, including Hii regions, molecular clouds, dark clouds, and diffuse emission at intermediate latitudes. Looking at the values in Table 1, it is possible to identify the observed range of values of the AME emissivity. A typical value for the AME emissivity is of the order 10 μK (MJy sr−1)−1 with a range from ~3 to 25 μK (MJy sr−1)−1 [2, 3]. In terms of flux density, this corresponds to approximately 1 Jy of AME at 1 cm for every 3000 Jy at 100 μm.Observations at far-IR and submillimeter wavelengths (e.g., [30, 31]) have revealed that the spectrum of the thermal dust emission is well represented by a modified black body function of the form
(1)Sν∝2hν3+βc21exp(hν/kTdust)-1,
where β is the dust emissivity spectral index and Tdust is the equilibrium temperature of the big dust grains. The thermal dust emission observed at such far-IR and submillimeter wavelengths is produced by big dust grains absorbing ultraviolet photons from the exciting radiation field of the environment and reradiating this emission thermally. It is this balance between absorption and emission that produces the well-defined spectrum of thermal dust emission (see Figure 1) and keeps the big grains in thermal equilibrium. In the ISM, the dust temperature is known to vary from one environment to another [32] with dust temperatures ranging from ~15 K in dense molecular clouds to 25–30 K in Hii regions. Additionally, the dust temperature has been observed to vary on small angular scales and within regions in which AME has been detected (e.g., between 15 and 25 K in the Perseus molecular cloud [33]). Since β is dependent on the physical properties of the dust grains, it is also expected to vary between different environments; Dupac et al. [34] found β to vary from 0.8 to 2.4 in a range of Galactic environments with dust temperatures between 11 and 80 K. Therefore, given that the emission at 100 μm is an observed quantity that depends strongly on Tdust and β, and since both Tdust and β are known to vary in the different phases of the ISM, the 100 μm emission is not a suitable quantity with which to estimate the AME emissivity.Figure 1
The modified black body spectrum representing the thermal dust emission for a range of dust temperatures from 15 to 45 K with a fixed dust emissivity spectral index of 1.8. This illustrates how the 100μm emission (vertical dashed line) increases with increasing dust temperature, and hence how the AME emissivity defined relative to the 100 μm emission is biased by changes in dust temperature.To quantify the effect of dust temperature variations on the AME emissivity, we consider a dust emissivity spectral index of 1.8, which is the value derived for the Galactic solar neighbourhood [32], and two dust temperatures (20 and 30 K). We then calculate that for a given intensity at 1 cm, the AME emissivity at 20 K is 11.1 times higher than at 30 K. This effect is substantial and could by itself, without any intrinsic variation of the AME intensity, explain the variations of the AME emissivities listed in Table 1. This clearly demonstrates the bias in the AME emissivity when computed relative to the 100 μm emission. With this in mind, we can now revisit the AME emissivity values listed in Table 1. For example, we focus on the five regions of AME detected in the Perseus molecular cloud (A1, A2, A3, B, and C) by Tibbs et al. [7]. It is apparent that region A1 is much less emissive (by a factor of ~4-5) than the other four regions. However, as discussed by Tibbs et al. [7, 33], the physical environment of region A1, which corresponds to the open cluster IC348, is very different from the other four regions. In fact, simply because of the difference in dust temperature between IC348 (~22 K) and the other regions (~18 K), we estimate that the AME emissivity is reduced by a factor of 4.3 in region A1 relative to the other regions. Therefore, the inconsistency between region A1 and the other regions can be solely accounted for by dust temperature variations, illustrating the impact of the dust temperature in biasing the AME emissivity computed using the 100 μm emission.


## 3. Redefining the AME Emissivity

In Section2, we have shown that the commonly used definition of the AME emissivity computed relative to the 100 μm emission is flawed. Therefore, an unbiased definition of the AME emissivity is required to help quantify the AME variations.As part of their investigation of the diffuse emission at intermediate latitudes, Davies et al. [3] computed the AME emissivity relative to the dust emission at 94 GHz rather than 100 μm. Since 94 GHz is further from the peak of the thermal dust emission (see Figure 1), the emission at 94 GHz is much less dependent on the dust temperature than the emission at 100 μm. However, even using the intensity of the 94 GHz emission, there is still an effect due to the dust temperature, as the AME emissivity decreases by a factor of 1.6 when the dust temperature increases from 20 to 30 K.Rather than using the intensity of the dust emission to compute an AME emissivity, other authors (e.g., [10, 26]) have computed an AME emissivity relative to the hydrogen column density, NH. This definition has several advantages. First, it results in an AME emissivity that is more directly suited to comparing with the theoretical values as the current spinning dust theories compute an emissivity in units of Jy sr-1 cm2 per H atom. Using NH also mitigates the effects of dust temperature, hence allowing a much more impartial comparison between measurements of AME in different environments. To illustrate this, we use the NH map of the Perseus molecular cloud from Tibbs et al. [33], which is based on a visual extinction map, and compute the AME emissivity for the five regions of AME in the Perseus molecular cloud (see Table 2). When computing the AME emissivity relative to NH, we find that region A1 is now consistent (within 2-3σ) with the four other regions in the Perseus molecular cloud, which leads to a completely different interpretation. We now want to explore how this new definition can help us to perform a less biased comparison between AME detections in different environments. Vidal et al. [26] compared the AME emissivity of six regions covering a range of column densities from Galactic cirrus to dense dark clouds (see Table 2). Adding these six measurements to the five we computed for the Perseus molecular cloud gives a sample of 11 AME emissivities based on NH. It is apparent that the AME emissivities for the Perseus molecular cloud are consistent with the values computed for other AME detections. It is also noticeable that there is less scatter in the AME emissivities listed in Table 2 compared to those listed in Table 1.Table 2
AME emissivity values for a variety of phases of the ISM computed relative to the hydrogen column density.
Source
N
H  (1022 H cm-2)
AME emissivity  (10-18 Jy sr-1 cm2/H)
Reference
Perseus A1
0.82
±
0.19
2.2
±
0.5
This work
Perseus A2
1.19
±
0.27
2.7
±
0.6
This work
Perseus A3
1.28
±
0.29
1.2
±
0.3
This work
Perseus B
0.79
±
0.18
3.4
±
0.8
This work
Perseus C
0.74
±
0.17
3.6
±
0.8
This work
Cirrus
0.15
±
0.07
4.6
±
2.0
Leitch et al. [1]
ζ Oph
0.22
±
0.02
4.1
±
0.6
Vidal et al. [26]
LDN1780
0.45
±
0.04
3.5
±
0.4
Vidal et al. [26]
LDN1622
1.50
±
0.15
2.0
±
0.2
Casassus et al. [11]
ρ Oph
5.00
±
0.50
3.2
±
0.5
Casassus et al. [6]
M78
22.80
±
0.23
0.9
±
0.1
Castellanos et al. [27]Following the example of Vidal et al. [26], we plotted the AME emissivities listed in Table 2 as a function of NH. This plot is displayed in Figure 2. We fitted the data with a power-law and found that the best fitting spectral index is -0.31±0.03. Although there is some scatter in the data, the results hint that the AME emissivity is decreasing with increasing NH, suggesting that the AME becomes less emissive as the column density increases. This is consistent with the AME being due to the smallest dust grains, as it is known that as NH increases, the abundance of the smallest dust grains decreases due to dust grain coagulation (e.g., [35, 36]). The idea of the AME being associated with the smallest dust grains is additional support for the AME being due to spinning dust.Figure 2
AME emissivity computed usingNH for the AME detections listed in Table 2 plotted as a function of NH. Also plotted is the fit to the data (dashed line), which we estimate to have a spectral index of -0.31±0.03.Although we have shown that defining the AME emissivity usingNH removes the bias introduced by the dust temperature and allows for a much more impartial comparison of AME between different environments, estimating NH, including molecular, atomic and ionised forms, is technically challenging. There are three methods used to estimate NH: near-IR extinction mapping; far-IR thermal dust emission, and the use of gas tracers. However, each method has its limitations. For example, near-IR extinction is only reliable if there are a sufficient number of background sources available, thermal dust emission estimates are affected by variations in the dust opacity, while observations of gas tracers, although they do provide kinematical information, are limited by the critical density of the observed molecule. A detailed analysis performed by Goodman et al. [37] clearly illustrates the various uncertainties and limitations in estimating NH using the three separate methods.In addition to the difficulties in measuringNH, the AME is believed to be associated with interstellar dust grains, and hence it would make more sense to normalise the strength of the AME with a physical property of the dust. For example, the column density of dust, Ndust, would be a useful quantity as it has the same advantages as using NH with the addition that it is a direct tracer of the dust. However, Ndust is even more difficult to estimate than NH. It is possible to assume a given Mdust/MH and convert from NH to Ndust; however, this assumes that we know Mdust/MH. Globally, within our own Galaxy, it is known that Mdust/MH is ~0.01 (e.g., [38]), but AME has recently been observed in external galaxies (e.g., [17, 39]) where it is known that the metallicity, and hence Mdust/MH, differs from our own Galaxy (e.g., [40]). To mitigate the effects introduced by variations in Mdust/MH, Ndust may be computed by fitting the far-IR SED with a dust model which incorporates the physical properties of the dust grains. However, to date, no such modelling has yet been performed systematically for the environments in which AME has been detected.


## 4. Conclusions

Having discussed various definitions of the AME emissivity, it is clear that using the 100μm emission introduces a bias due to the effect of dust temperature, and that using NH allows for more accurate comparisons. However, given the association between the AME and the interstellar dust, an intrinsic property of the dust such as Ndust, would represent the best quantity with which to compute the AME emissivity. Therefore, until such a time as a new AME emissivity is defined, we stress that care should be taken before using it to directly compare the strength of AME in different environments.



---

*Source: 124931-2012-12-12.xml*