Abstract

The anomalous microwave emission (AME) has been proved to be an important component of the galactic diffuse emission in the range from 20 to 60 GHz. To discriminate between different models of AME, low frequency microwave data from 10 to 20 GHz are needed. We present here a reanalysis of published and unpublished Tenerife data from 10 to 33 GHz at large angular scales (from 5 to 15 degrees). We cross-correlate the Tenerife data to templates of the main galactic diffuse emissions: synchrotron, free-free, and thermal dust. We find evidence of dust-correlated emission in the Tenerife data that could be explained as spinning dust grain emission.

1. Introduction

The anomalous microwave emission (AME) is an important contributor of the galactic diffuse emission in the range from 20 to 60 GHz. It was first identified by [1, 2] as free-free emission from electrons with temperature,  K. Draine and Lazarian [3] argued that AME may result from electric dipole radiation due to small rotating grains, the so-called spinning dust. Models of the spinning dust emission, Draine and Lazarian [4] show that an emissivity spectrum peaking at around 20–50 GHz is able to reproduce the observations [512]. The initial spinning dust model has been refined regarding the shape and rotational properties of the dust grains [1315]. An alternative explanation of AME was proposed by Draine and Lazarian [16] based on magnetic dipole radiation arising from hot ferromagnetic grains. Observations have placed limits of a few percent on the fractional polarization towards AME targets [9, 1720]. This excludes perfectly aligned single-domain magnetic grains; however, other alignments and grain compositions produce similarly low levels of polarization [21].

A correlation between microwave and infrared maps, mainly dominated by dust thermal emission [22], was observed for various experiments, for example, on COBE/DMR [23, 24], OVRO [1, 2], Saskatoon [25, 26], survey at 19 GHz [27, 28], and Tenerife [29, 30]. A similar signal was found in compact regions by [5] and in some molecular clouds based on data from COSMOSOMAS [7, 31], AMI (Ami-Consortium: [32, 33]), CBI [9, 34], VSA [12], and Planck [35]. A recent study of the Small Magellanic Cloud also claims a detection of AME [36].

Independently, Bennett et al. [41] proposed an alternative explanation of AME based on flat-spectrum synchrotron emission associated to star-forming regions to explain part of the WMAP first-year observations. This hypothesis seems to be in disagreement with results from de Oliveira-Costa et al. [6]; Fernández-Cerezo et al. [42]; Hildebrandt et al. [43]; Ysard et al. [44] which showed that spinning dust best explained the excess below 20 GHz. Furthermore, Davies et al. [45] showed the existence of important correlation between microwave and infrared emission in regions outside star-forming areas. More recently, Kogut et al. [46] discussed the fact that spinning dust fits better to ARCADE data (3.8 and 10 GHz) than a flat-spectrum synchrotron.

To discriminate between the different AME models and the alternative explanations such as those discussed earlier low frequency microwave data in the range, from 10 to 20 GHz at different angular resolutions are needed. Indeed in this frequency range, we expect the AME spectrum to be significantly distinct from magnetic dust and flat-spectrum synchrotron. At this respect, the Tenerife data set, from 10 to 33 GHz and at large angular scales (from 5 to 15 degrees), is unique. We present in the following a re-analysis of these data including previously published data by Gutierrez et al. [47] and un-published data since January 1998 to December 2000. The paper is structured as follows.Section 2 presents the Tenerife data and discusses the reprocessing of these data. In Section 3 we describe the main galactic emission mechanisms and the associated templates used in the analysis. Section 4 discusses the point source contribution to the Tenerife data. In Section 5, we present the cross correlation between the Tenerife data and the galactic templates. The main results are discussed in Section 6. Finally, we draw conclusions in Section 7.

2. The Tenerife Data

The Tenerife experiment observed unidimensional scans at constant declination at 10, 15, and 33 GHz using a triple beam pattern of FWHM 5° and of 8.1° beam spacing. The region of the sky centered at declination was chosen for observations because it corresponds to the largest area of the sky at high latitudes where contamination from foregrounds is a minimum. To reconstruct 2D maps of the sky, consecutive declinations separated by half the beamwidth (2.5°) are observed. Each single declination is repeatedly observed until sufficient sensitivity and full RA coverage is achieved. The scheduling of the observations takes into account the position of the sun so that its contribution to the data is a minimum. This requires observations of the same declination at different times of the year for full RA coverage. Observations were performed day and night. Day-time observations present an increase in total power and noise with respect to night-time observations due to receiver gain changes. In extreme cases, data observed at day time have to be removed.

The atmospheric contribution to the data depends on frequency, being severe at 33 GHz and relatively small at 10 GHz. The observing efficiency depends mainly on the weather conditions. Data strongly affected by atmospheric contamination cannot be used. In addition, technical problems, such as power cutoff, warmed and oscillating HEMT amplifiers, RF interference, and failure in the electronic systems, can also affect the data. At 10 and 15 GHz, more than of the observed data are useful. However, at 33 GHz only about of data are kept, due mainly to atmospheric effects.

We use in the following the full Tenerife data set which includes new data (taken from January 1998 to December 2000) with respect to previously published releases (see, e.g., [47, 48]). These data have been completely reprocessed including hand editing to remove clearly systematic contaminated regions and accurate removal of the atmospheric emission using an MEM-based baseline removal technique as described in Gutierrez et al. [47]. The latter has been improved to enlarge the sky region for which atmospheric residuals are negligible. In particular, this allows us to consider in the following analysis regions at low galactic latitudes. Figures 1, 2, and 3 present the cleaned Tenerife data at 10, 15, and 33 GHz, respectively. The main properties of these clean data: central declination of the scan, mean temperature, mean noise per beam, and r.m.s per 1° pixel in the RA interval 150° to 250°, are presented in Tables 1, 2 and 3 for the 10, 15, and 33 GHz channels, respectively.

3. Galactic Diffuse Emission

We present in this section the main known galactic diffuse mechanisms: synchrotron, free-free, and vibrational and rotational dust.

3.1. Synchrotron

Synchrotron emission results from cosmic-ray electrons accelerated in magnetic fields, and thus, depends on the energy spectrum of the electrons and the intensity of the magnetic field [49, 50]. The local energy spectrum of the electrons has been measured to be, for energies contributing to the observed radio synchrotron emission, a power law to good approximation with index from about −2.7 to −3.3 over this energy range [51]. Such an increase of energy spectrum slope is expected, as the energy loss mechanism for electrons increases as the square of the energy.

Radio surveys at frequencies less than 2 GHz are dominated by synchrotron emission [52]. The only all-sky survey available at these frequencies is the 408 MHz map [53]. This survey has a resolution of 0.85° and was produced using the Parkes 64 m telescope in Australia for the southern sky and the Bonn 100 m and Jodrell Bank MK1A telescopes for the northern sky. The scanning strategy with the Bonn telescope was to fix the azimuth at the local meridian and scan up and down in elevation at a rate of about 6°/min. This technique reduced the atmospheric contribution to the map but led to a set of vertical stripes (constant RA) separated by 7°. The quoted errors in the temperature scale are of the order of and K in the absolute brightness temperature levels.

In addition, there is the 1420 MHz survey [54] which covers the declination range −19° to 90° and has a FWHM of 0.58°. Stripes are also present in this map due to the scanning strategy which consisted of azimuthal scans at constant elevation. The errors in the temperature scale are of the order of and K in the absolute brightness temperature.

The 408 MHz and 1420 MHz maps have been used to determine the synchrotron spectral index at radio frequencies [54, 55]. Assuming that and after destriping and correction for zero levels, spectral indexes, , in the range to were found. The spatial angular power spectrum of the synchrotron emission is poorly understood but is believed to be [56]. At high galactic latitudes in the region observed by the Tenerife experiment, the synchrotron spatial angular power spectrum is slightly flatter than [57]. In the following sections, we will use the destriped version of the 408 MHz map [55] as a template of the synchrotron emission. This map is shown in Figure 4.

3.2. Free-Free

When a charged particle is accelerated in a Coulomb field it will emit radiation which is called braking radiation or Bremsstrahlung. The galactic free-free emission is the thermal bremsstrahlung from hot electrons (~104 K) produced in the interstellar gas by the UV radiation field [51]. This emission is not easily identified at radio frequencies, except near the galactic plane. At higher latitudes, it must be separated from synchrotron emission by virtue of their different spectral indices, since the spectral index of optically thin free-free emission is .

The diffuse galactic recombination radiation, is a good tracer of free-free emission since both are emitted in the same Warm Ionised Medium (WIM) and both have intensities proportional to the emission measure, The ratio of free-free brightness temperature to the surface brightness in (Rayleigh) is [58] where is the observing frequency in GHz, is the temperature of the electrons in units of 104 K, and the last factor corresponds to helium which is assumed completely ionised and creates free-free-emission like hydrogen but does not emit light.

Recently, a full sky survey of light has been released by the WHAM (Wisconsin Mapper) collaboration [59]. The WHAM instrument consists of a 6 inch dual-etalon Fabry-Perot spectrometer with a narrow band filter of FWHM ~20 Å which images onto a cryogenically cooled CCD. The resolution of the survey is ~1°. The spatial power spectrum of free-free emission has not yet been determined from the WHAM map; however, analysis of images of the North Celestial Pole area [60] suggests that .

In Figure 4, we present a map of the northern sky produced by the WHAM survey. In the following sections, we will use this map as a template for the free-free emission. At intermediate Galactic latitudes (say |b| > 30°), about 10 of the light is absorbed by dust and therefore estimates of free-free emission from will be systematically lower than the true value. At latitudes below 10°, this correction becomes increasingly uncertain.

3.3. Vibrational Dust

At the higher frequency range (≥100 GHz) of the microwave background experiments, dust emission starts to become dominant. Dust grains are heated by interstellar radiation, absorbing optical and UV photons, and emitting energy in the far infrared.

The intensity of the emission from an ensemble of dust grains is given by where is the emissivity at frequency , and the integral is along the line of sight. In the Rayleigh-Jeans regime and assuming a constant line-of-sight density of dust, where is the brightness temperature of the dust emission.

The spectrum of the dust emission has been measured at millimeter and submillimeter wavelengths by the Far Infrared Absolute Spectrophotometer (FIRAS) and can be fitted by a single greybody spectrum of temperature K and emissivity [61] at high galactic latitudes. From IRAS observations of dust emission [62], it was found that the spatial power spectrum of the dust fluctuations is . This has also been confirmed at larger angular scales by the COBE-DIRBE satellite.

In the following sections, we will use the combined IRAS-DIRBE map at 100 μm [63] as a template for the dust emission. This map has a resolution of FWHM ~6 arcmin and covers the full sky. Zodiacal light and point sources have been removed from the map, and the regions of the sky which were not observed by the IRAS satellite have been replaced by DIRBE data. The combined map preserves the DIRBE zero point and calibration. This map is in units of MJy/sr. In Figure 4, we present the northern sky part of the combined IRAS-DIRBE map.

3.4. Rotational Dust

Small spinning interstellar dust grains containing atoms can produce detectable rotational emission in the 10–100 GHz range. This emission depends on the component of the electric dipole moment perpendicular to the angular velocity of the grain and on the physical properties of the interstellar medium [3]. For these small grains, rotational excitation is dominated by direct collisions with ions and plasma drag. The very smallest grains () have their rotation damped primarily by electric dipole emission; for plasma drag dominates.

In the following sections, we show models for the spectrum of the spinning dust emission provided by Dr. Draine (private communication). These models depend on a large set of parameters such as the distribution of grain sizes, the charge of the grains, the composition of the grains, and the physical properties of the interstellar medium which were fixed by the authors. However, the normalization of the model can be assumed as a free parameter although it is related to the hydrogen column density in the interstellar medium which is of the order of a few times .

The spatial distribution of small dust grains is not well-known although it seems reasonable to believe that it is not different from that of larger grains but for dense regions where dust coagulation may deplete small grains. For the purpose of this work, we will use the IRAS-DIRBE 100 m map as a template for the spinning dust grains.

4. Galactic and Extragalactic Point Sources

The contribution from resolved point sources to the Tenerife data at 10 and 15 GHz was extrapolated from data of the Michigan monitoring program (M. Aller and H. Aller 2000, private communication). The Michigan program regularly monitors point sources with fluxes above 0.5 Jy at 4.8, 8.0, and 14.5 GHz. The Michigan catalog is neither complete in flux nor time domain. Well-known strong variable sources such as 3C345 are regularly observed, and data at all observing frequencies are available. Weak sources are poorly observed, and often data are available at a single frequency of the three possible. Further, data are available up to June 1999 while the Tenerife experiment operated until September 2000. For the 33 GHz Tenerife data, we used the Metsahovi catalog (Metsahovi group, private communication) which regularly monitors sources above 1 Jy at 22 and 37 GHz. This catalog is neither complete in flux nor time, although observations are available up to January 2000.

We have developed a software to produce time and frequency complete point source catalogs. The flux of the sources for each frequency was interpolated in time by fitting Fourier series to the data. If the number of independent observations per frequency was smaller than 10, the sources were assumed not variable, and the data was fitted to a constant with time. The extrapolation of the source fluxes into the future time was performed assuming no variability since last observed data point. Spectral indices were calculated for each source from the available data. If only observations at a single frequency were available, we assumed a flat spectrum. We constructed four point source catalogs at 10, 13, 15, and 33 GHz covering the time range from January 1984 to December 2000 with one Julian day sampling and a flux limit of about 0.5 Jy.

For each single day of observations with the Tenerife radiometers, we estimated the total contribution from point sources to the data. This contribution was calculated from the extrapolated catalog at 10, 15, and 33 GHz and the flux, , was converted into antenna temperature, , using [64] where is the beam solid angle, and is the wavelength. The daily point source contributions were stacked in the same manner as the Tenerife data and subtracted from the Tenerife stacks. In the left column of Figure 5, we plot the Tenerife stacked data at 15 GHz for the 10 declinations presented in this paper and overplot in green the contribution from point sources to the data. In the right column, we plot the stacks after subtraction of the point sources.

4.1. Unresolved Sources

A study of the contribution of unresolved point sources (i.e., weak point sources not detected individually) to CMB experiments has been produced by Franceschini et al. [65]. They used numerous surveys, including VLA and IRAS data, to put limits on the contribution to single beam CMB experiments by a random distribution of point sources. The contribution from unresolved point sources to the Tenerife data deduced from the previous analysis is presented in Table 4. We consider sources with flux under 1 Jy at 5 GHz and a main beam of 5°.

5. Assessment of the Galactic Contribution to the Tenerife Data

To assess galactic contribution to the Tenerife data, we correlated observations with the galactic synchrotron, dust, and emission templates described in the previous sections. We convolved the foreground maps with the Tenerife beam at each of the Tenerife frequencies before the correlation was performed. The Tenerife data used in the correlation are the final stacked data presented in Section 2. At 15 GHz we used ten declination stacks in the range 25°–47.5°; seven declination stacks at 10 GHz in the range 30°–45°; and only two declination stacks at 33 GHz covering declinations 40° and 42.5°. The data were processed so that data at were excluded from the baseline fit and therefore were neither stacked or reconstructed. The discrete point sources were subtracted from the Tenerife stacks as discussed in the previous section.

5.1. The Method

To simultaneously correlate the Tenerife data to the three galactic templates, we use a method which was first applied to this problem by Górski et al. [66] to fit galactic and extragalactic templates in Fourier space to the COBE-DMR data. This method was applied to the Tenerife data first by de Oliveira-Costa et al. [29, 30] and then by Mukherjee et al. [39] to study the possible emission of spinning dust at the Tenerife frequencies.

Assuming that the microwave data is a superposition of CMB, noise, and galactic components, we can write where is a Tenerife data vector of pixels; is an element matrix containing foreground templates convolved with the Tenerife beam; is a vector of size that represents the levels at which these foreground templates are present in the Tenerife data—correlation coefficients for each foreground template; is the instrumental noise in the data; and is the CMB component of the data. For this analysis, we assume the noise and CMB to be uncorrelated.

The minimum variance estimate of is given by with errors given by , where is given as In the before mentioned, is the total covariance matrix, the sum of the noise covariance matrix, and the CMB covariance matrix. The noise covariance matrix of the Tenerife data is taken to be diagonal—no correlation between different pixels. The CMB covariance matrix was obtained analytically following Zaldarriaga [67].

This correlation method produces minimum variance and unbiased estimates of if the following holds. (i)The error in the Tenerife data is Gaussian and with zero mean. To a very good approximation the instrumental noise in the Tenerife data is uncorrelated.(ii)The templates perfectly trace the foreground emissions at microwave frequencies. There may be components of emission present in the Tenerife data, apart from those we have identified, for which there are not obvious counterparts at other frequencies. Moreover, the templates we use do not perfectly trace the microwave emission of the galactic foreground. The latter is specially important in the case of vibrational and rotational dust for which the template comes from data at much higher frequencies. For example, Draine and Lazarian [3] proposed that a 12 μm map would be a much better tracer of rotational dust emission because small grains will also emit at this frequency. Also, the emission is absorbed at low galactic latitudes by interstellar dust and therefore does not fully trace the free-free emission at those latitudes. Moreover, the 408 MHz map has residual striping. (iii)The templates are perfect—errors bars equal to zero. The error estimated for the templates are about 5%–10 % of the signal. They are not corrected for in this analysis because they are small compared to the errors in the Tenerife data. (iv)The templates are not correlated. It there is a correlation between the different templates, the set of minimum variance correlation coefficients is degenerate and, therefore, we would not be able to discriminate among the different foreground components. At low galactic latitudes, the galactic plane emission dominates and Mukherjee et al. [39] have found that the foreground templates are correlated in this region. (v)The correlation coefficient is the same throughout the area of sky for which the correlation is performed. If it is not, the error associated with is systematically underestimated. As a double check, we divided the sky observed into independent areas and performed the same correlation test in each, calculating a mean correlation coefficient and the dispersion of the individual values which in most cases was in good agreement with the error of calculated for the total area.

5.2. Cross-Correlation Results

The cross-correlation results are presented in Table 5. The correlation was performed for three different galactic cuts , , and . represents the r.m.s. of the galactic maps after convolution with the Tenerife beam. and are the correlation coefficient and the error associated with it. They have units of K/K, K/(MJy sr) and K/R for the correlation with the 408 MHz, 100 μm, and WHAM maps respectively. is the r.m.s. contribution from the galactic foregrounds to the Tenerife data, which was obtained as . This analysis is improved with respect to previous analyses by de Oliveira-Costa et al. [29, 30] and Mukherjee et al. [39], First, because we present Tenerife data at 10 and 15 GHz for a much larger area of the sky and with improved sensitivity; and second, because we also include in the cross-correlation a template for the free-free emission. Moreover, data at 33 GHz have also been included in the analysis although the area of the sky observed is significantly smaller than at 10 and 15 GHz and consequently the error bars in the estimated cross-correlation coefficients are much larger. We have pictorially summarised the cross-correlation results in Figure 6. In the top row of the figure, we plot the r.m.s. contribution from synchrotron to the Tenerife data at the three galactic cuts analysed. In black, we overplot the expected synchrotron contribution for a spectral index of and derived from the r.m.s. level in the 408 MHz map. The contributions are systematically lower than expected at 10 GHz, and the spectral index is much flatter such that, for example, the contribution at 15 GHz is consistent with that at 10 GHz for b > 40°. We have derived from the Tenerife data spectral indexes of , and at , and respectively, and the r.m.s contributions based on these spectral indices is overplotted in red.

In the middle row, we plot the r.m.s. contribution from free-free emission to the Tenerife data. In black, we overplot the expected free-free emission at microwave frequencies derived from (2) which are actually in very good agreement with the observations. In red, we plot the free-free contribution for spectral indexes of , , and at , , and , respectively.

In the bottom row of Figure 6, we plot the r.m.s. contribution from dust emission to the Tenerife data which is significantly larger (few orders of magnitude) than expected from vibrational dust. In solid black, we overplot an r.m.s brightness temperature expected from the CNM model of rotational dust proposed by Draine and Lazarian [3]. The models were rescaled to fit the data. The intensity spectrum for this model peaks at 50 GHz which corresponds to a peak in brightness temperature around 20 GHz of about 300 μK. From the best fit model to the data, we have estimated fluctuations in the temperature of the rotational dust of ~8% at the angular scales of the Tenerife experiment.

The r.m.s. contributions from dust to the 15 GHz data at and , as well as the contribution from free-free at seem to be underestimated. This could be caused by correlations between the synchrotron, dust, and free-free templates which would reduce the validity of the minimum variance solution and could bias the estimates of the correlation. In Figure 7, we plot the following correlations, from top to bottom synchrotron versus dust; free-free versus dust, and free-free versus synchrotron for , , and . We observe no correlation between free-free and synchrotron, moderate correlation between free-free and dust and a quite strong correlation between synchrotron and dust at low galactic latitudes weakening down at high latitudes. Note that at , the correlation plots show negative versus negative points which do not follow the correlation pattern of the main body of points. The former points correspond to the lowest galactic latitude data which run into a negative beam of the triple beam pattern for the galactic plane crossing. The observed correlations between templates can not justify the observed lack of dust correlation at 15 GHz data at and . In the following, we concentrate on the region for which we detect significant dust correlated signal at all Tenerife frequencies.

6. Dust-Correlated Emission

Hinshaw et al. [37] cross-correlated the COBE differential microwave radiometer (DMR) maps with DIRBE far infrared maps and discovered that statistically significant correlation did exist at each DMR frequency, which was inconsistent with vibrational dust alone. This extra correlation was explained at the time as free-free emission. Following this publication, other authors have cross-correlated CMB data sets with dust templates and found an excess of correlation which has been interpreted as free-free emission, flat-spectrum synchrotron, or emission from spinning dust. In Table 6, we present an up-to-date list of CMB data for which the correlation with dust has been performed. In all cases, a multitemplate method has been used to perform the correlation but no free-free template has been used. We have also added in this table previous analyses of the Tenerife data by de Oliveira-Costa et al. [29, 30] and Mukherjee et al. [39] which did not include the 33 GHz data and covered a much smaller area of the sky at 10 and 15 GHz.

The correlation coefficients we deduced are significantly smaller at than those calculated by de Oliveira-Costa et al. [29, 30] and Mukherjee et al. [39]. This is probably due to 50 % larger area (localized signals are diluted) and a more careful subtraction of baselines at low galactic latitudes in our data. The inclusion of an extra free-free template seems to play a minor role on this at 10 GHz and none at 15 GHz.

The cross-correlation results for the Tenerife experiment presented in this paper confirm the existence of extra dust-correlated emission at microwave frequencies and suggest that it is not due to free-free emission. The moderate correlation found between the dust and free-free templates could confuse the results obtained but clearly cannot account for all the observed dust-correlated component. A more detailed study is needed to take into account correlations between templates in the calculation of the correlation coefficients.

In Figure 8, we plot the dust cross-correlation coefficients () for the Tenerife (diamonds) and COBE-DMR (triangles) data. The rest of the experiments presented in Table 6 were not included in this plot because they observe at different angular scales and different areas of the sky (see table for details) and therefore no direct comparison with the Tenerife and COBE data is possible. The Tenerife data suggest a peak in the spectrum at about 30 GHz although this is mainly based on the data at 33 GHz which is significantly much noisier and cover a much smaller part of the sky. However, DMR data combined with the 10 and 15 GHz Tenerife data points prefer a peak in the range 15–20 GHz. In color, we overplot the six spinning dust models proposed by Draine and Lazarian [3] combined with the standard vibrational dust model described in the previous sections, for which the amplitude has been taken as a free parameter and fitted to the COBE and Tenerife 10 and 15 GHz data. We observe that the shape of the spectrum given by the data is similar to some of the predicted spectra although the models do not indicate a sharp rise in the range 10–15 GHz.

7. Conclusions

We have presented in this paper a reanalysis of the Tenerife data including previously published data and new data from January 1998 to December 2000. This analysis leads to evidences for an excess of dust-correlated emission at microwave frequencies in the range 10–33 GHz and at large angular scales, from 5 to 15 degrees. This correlation can not all be associated with free-free emission. A combination of the Tenerife and COBE data suggests spinning dust emission could account for the extra correlation. However, the scatter observed in the data and the discrepancy in the spectrum shape indicate other components may also be responsible for the extra correlation. Furthermore, the analysis does not take into account correlations between the galactic emission templates and this can bias the estimate of the correlation coefficients. To correct this, a more detailed analysis, which will account for the cross-correlation between templates, is needed.

Acknowledgment

The authors would like to thank the editors of this paper.