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.

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., [24]) and specific Galactic sources such as the Perseus and Oph molecular clouds [59], the dark clouds LDN1622 [10, 11], LDN1621 [12], and LDN1111 [13], and a variety of HII regions [1416]. 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 [2024], 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.

In this work, we aim to (1) illustrate the weakness of such a definition and (2) present and discuss possible alternatives. In Section 2, 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 Table 1, 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 where is the dust emissivity spectral index and 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 and , and since both 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.

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 Section 2, 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, . 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  cm2 per H atom. Using 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 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 , 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 . 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.

Following the example of Vidal et al. [26], we plotted the AME emissivities listed in Table 2 as a function of . This plot is displayed in Figure 2. We fitted the data with a power-law and found that the best fitting spectral index is . Although there is some scatter in the data, the results hint that the AME emissivity is decreasing with increasing , 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 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.

Although we have shown that defining the AME emissivity using removes the bias introduced by the dust temperature and allows for a much more impartial comparison of AME between different environments, estimating , including molecular, atomic and ionised forms, is technically challenging. There are three methods used to estimate : 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 using the three separate methods.

In addition to the difficulties in measuring , 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, , would be a useful quantity as it has the same advantages as using with the addition that it is a direct tracer of the dust. However, is even more difficult to estimate than . It is possible to assume a given and convert from to ; however, this assumes that we know . Globally, within our own Galaxy, it is known that 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 , differs from our own Galaxy (e.g., [40]). To mitigate the effects introduced by variations in , 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 allows for more accurate comparisons. However, given the association between the AME and the interstellar dust, an intrinsic property of the dust such as , 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.

Acknowledgments

The authors thank the referee for useful comments which helped improve the content of the paper. This work has been performed within the framework of a NASA/ADP ROSES-2009 Grant no. 09-ADP09-0059. C. Dickinson acknowledges the support from an SFTC Advanced Fellowship and an EU Marie-Curie IRG Grant under the FP7.