Planck early results
Free Access
Issue
A&A
Volume 536, December 2011
Planck early results
Article Number A21
Number of page(s) 18
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201116455
Published online 01 December 2011

© ESO, 2011

1. Introduction

In order to understand our own Galaxy, it is necessary to explore Galactic plane (|b| ≲ 10°), where we are able to observe the emission coming from a large range of distances. However the observed emission is the sum of a large number of line-of-sight components, often probing very different environments.

Several previous studies (Bloemen et al. 1986; Bloemen et al. 1990; Giard et al. 1994; Sodroski et al. 1997; Paladini et al. 2007) have separated observed integrated emission into a number of Galactocentric radii in order to study its properties as a function of Galactic position and in different phases of the interstellar gas (e.g., atomic, molecular and ionised). The radial velocity of the gas is used to separate the Galactic gas emission into a number of Galactocentric rings and then the spectral energy distribution (SED) of each ring/gas phase is fitted with a physical model of dust and gas emissions. These methods have been used successfully in previous studies to map out basic properties of the ISM throughout the Galaxy.

For instance, Giard et al. (1994) demonstrated that polycyclic aromatic hydrocarbons (PAHs) are a ubiquitous component of the ISM, and Bloemen et al. (1990) showed that the dust temperature decreases with distance to the Galactic Centre in a way which is fully consistent with an exponential decrease of the interstellar radiation field (ISRF) in the stellar disk. Sodroski et al. (1997) suggested that the abundance of large dust grains within each gas phase exhibits a gradient that is equivalent, within the uncertainties, to the metallicity gradient in the Galactic disk. Paladini et al. (2007) showed that the dust in molecular clouds appears to be heated in a significant way by young massive stars still embedded in their parent clouds.

With the advent of the Planck satellite1, it is now possible to perform an inversion on the emission arising from the entire infrared, millimetre and centimetre range, making it possible to determine the dust and gas properties in many different environments in the Milky Way. The High Frequency Instrument (HFI) channels (100–857 GHz) allow us to properly constrain the big grain temperature and emissivity, while using the Low Frequency Instrument (LFI) channels (30–70 GHz), it is possible to constrain non-thermal emission mechanisms such as free-free, synchrotron or anomalous dust. We have used these data along with other ancillary data to perform a large-scale, low-resolution analysis of the Milky Way ISM emission.

Knowledge of the dust emission and physical conditions throughout the Milky Way provides a self-consistent context for other Planck studies: for example, the environmental effect on the properties of cold cores (Planck Collaboration 2011r,s) can be evaluated, the modelling of individual anomalous dust regions can be placed in a Galaxy-wide context (Planck Collaboration 2011p), or the Milky Way values can simply be compared to other galaxies (Planck Collaboration 2011l,m).

In Sect. 2 we describe the data used in the study, including both those to be inverted and those used as templates to represent a particular phase of the ISM. In Sect. 3, we describe how we optimise the separation of the gas observation into a number of Galactocentric rings, as well as how we report uncertainties on our best solution. We present the results in Sect. 4, and discuss potential sources of uncertainty and bias in Sect. 5.1. Finally we conclude in Sect. 6.

2. Data

There are two categories of data used in this study, the data to be inverted, and data used to create templates which will be used to perform the inversion.

2.1. Data to be inverted

We use data from 12μm (25 THz) to 1.4 GHz in order to study the dust and gas properties in a self consistent manner throughout the Milky Way. Here we briefly present these data in order of decreasing (increasing) frequency (wavelength).

The Infrared Astronomical Satellite (IRAS), a joint project of the US, UK, and the Netherlands (Neugebauer et al. 1984), performed a survey of 98% of the sky at four wavelengths: 12, 25, 60, and 100μm. We use a reprocessed version of this data set (IRIS, Miville-Deschênes & Lagache 2005) which benefits from a better zodiacal light subtraction and from a calibration and zero level compatible with DIRBE, as well as from a better destriping.

The DIRBE instrument (Diffuse Infrared Background Experiment, Hauser et al. 1998) on board the COBE satellite imaged the full sky in 10 broad photometric bands from 1 to 240μm with a beam of 0.7°. The bands at 140 and 240μm allow us to bridge the gap between IRAS and the highest frequencies of HFI, thus covering the peak of the thermal dust SED.

The Planck data we used for this analysis are the DR2 LFI and HFI maps, with the cosmic microwave background (CMB) removed, as described in detail in Zacchei et al. (2011) and Planck HFI Core Team (2011b), respectively. It consists of a set of 9 frequency maps with central frequencies of 28.5, 44.1, 70.3, 100, 143, 217, 353, 545 and 857GHz), smoothed to a common resolution of 1° FWHM, at a common HEALPix (Górski et al. 2005) pixelization resolution of Nside = 256 (~15′ pixels). For simplicity, in the remainder of the text we will refer to the first three bands as 30, 44 and 70 GHz, respectively. The CMB component that was removed from these maps was found by the needlet internal linear combination (NILC) CMB extraction method presented in Planck HFI Core Team (2011b).

Finally, to provide constraints at radio frequencies we use the full-sky 1.4 GHz data from (Reich 1982; Reich & Reich 1986; Reich et al. 2001).

All of the above data sets have been smoothed with a Gaussian function to a common resolution of 1° FWHM and gridded to a common HEALPix Nside of 256.

2.2. Templates

2.2.1. Atomic phase

We aim to characterise the different emission mechanisms in the different phases of the ISM (atomic, molecular and ionised). For the atomic phase we use the Leiden/Argentine/Bonn Galactic Hi Survey (Kalberla et al. 2005). The survey results are expressed in antenna temperature, TA. In order to convert this into column density, the effect of optical depth should be taken into account, for example (Binney & Merrifield 1998) (1)for a given spin temperature (Ts). Recent studies in the plane of the Milky Way (Dickey et al. 2003, 2009) advocate the use of a higher spin temperature of around 250 K. This higher value is justified on the grounds that along any line of sight there is a mix of cold and warm Hi gas, and so an effective spin temperature is needed to correct for the opacity. This will most certainly vary as a function of Galactic longitude. However, constraining this variation is beyond the scope of this study. We have therefore assumed a constant hydrogen spin temperature of 250 K for the entire Galactic plane. Varying this value in the range 150–400 K has only a small effect on our results and they remain consistant within the quoted uncertainties.

2.2.2. Molecular phase

We assume that the bulk of the molecular gas mass is at first order well traced by CO emission (Tielens 2005). The CO data we use come from the Composite CO Survey of Dame et al. (2001). The column density of molecular hydrogen can be expressed as a function of antenna temperature (TCO) (2)where XCO = N(H2)/WCO is the ratio of molecular hydrogen column density to the velocity-integrated intensity of the 12CO line. We use a value for XCO which is compatible with the recent values from the Fermi collaboration (Abdo et al. 2010; Ackermann et al. 2011), XCO = 1.8 × 1020, which results from a study on diffuse gamma-ray emission in the Galactic plane.

Both the atomic and molecular data contain line of sight velocity information. It is therefore possible to use this information to separate the emission into a set of Galactocentric rings, as a given line of sight velocity measurement uniquely identifies a Galactocentric distance, assuming circular orbits. This should not be confused with heliocentric distances for which there is a distance ambiguity corresponding to the near or far side of the derived Galactocentric orbit. The velocity information alone does not allow one to determine if the observed gas is on the near side of the determined orbit or on the far side. We separate the observed Hi and CO observations each into a set of 5 Galactocentric rings. The reasons for choosing this decomposition is explained in Sect. 3.2

2.2.3. Dark gas

It has been shown (e.g., Reach et al. 1994; Grenier et al. 2005) that the CO and Hi do not account for all the atomic and molecular gas in interstellar clouds: there is in addition a dark gas phase. This phase is also extensively studied in Planck Collaboration (2011o), where the total gas column density, as traced by Hi and CO observations, is correlated with the optical depth of the dust calculated using IRAS 100μm along with Planck HFI data at 350 and 500μm (857 and 545 GHz). The correlation is linear at low and high column densities, but it departs from linearity at column densities of 7.4 × 1020 < NH < 5.0 × 1021. This excess is attributed to thermal emission by dust associated with a dark gas phase, undetected by Hi and CO measurements.

The dark gas column density used here is defined in Planck Collaboration (2011o) as (3)where τD is the thermal dust optical depth, τM = (τD/NH)ref  [NHi + 2XCOWCO] ) is the modelled dust opacity and (τD/NH)ref is the reference dust emissivity measured in low NH regions. We use their dark gas estimate at 857 GHz as an additional template.

The dark gas phase represents a significant fraction of the total gas column density, as shown in Fig. 2. The Hi dominates over almost the entire longitude range, and the CO is strongest in the molecular ring 40° > l >  − 40°. Outside of the molecular ring the dark gas is stronger than the CO, except towards the anti-centre direction where the dark gas and CO are roughly comparable. The proportion of dark gas to total hydrogen column density ranges from roughly 10% in the anti-centre direction up to nearly 60% towards the Galactic centre.

2.2.4. Synchrotron emission

An additional emission mechanism which is important in the frequency range we are studying and is not traced by usual gas tracers is synchrotron emission, which arises from the acceleration of cosmic-ray electrons in magnetic fields. This emission follows approximately a power law but the spectral index most certainly varies on the sky and as a function of frequency. More than 90% of the observed synchrotron emission arises from a diffuse component (Bennett et al. 2003) and in the Galactic plane recent studies have shown that there is no hardening of the synchrotron spectrum (Finkbeiner 2004; Boughn & Pober 2007), but that emission from spinning dust is likely to dominate between 8 and 60 GHz. For our synchrotron template we have chosen the 408 MHz all-sky survey from Haslam et al. (1982), which traces the diffuse, soft synchrotron component.

thumbnail Fig.1

Templates used in the inversion. The top five images are the Hi templates, the following five are the molecular templates, all of which are presented in increasing Galactocentric radius (the radii defining each ring are described in Sect. 3.2). The final three templates are the free-free, synchrotron, and dark gas templates, respectively. The units are all expressed in 1020 atoms, except for synchrotron, expressed in MJy sr-1. The centre and anti-centre regions have been masked as no distance information can be obtained from the radial velocity of the gas.

2.2.5. Ionised phase

To represent the ionised component we have used the maximum entropy method (MEM) free-free template from the WMAP seven year data (Jarosik et al. 2011). The intensity map (Iν) has been converted to ionised hydrogen column density as described by Sodroski et al. (1997): (4)for effective electron density neff and electron temperature Te. Following Sodroski et al. (1989), we use neff = 10 cm-3 and Te = 8000 K. In the future, the availability of radio recombination line (RRL) surveys will provide a more direct measurement of the free-free component in the Galaxy and its distribution separated into distinct Galactocentric bins using the radial velocity information included in such surveys.

As has been done for the data, the templates have all been smoothed to 1° FWHM and gridded to an Nside of 256. These templates are all displayed in Fig. 1.

3. The Galactic inversion

Component separation techniques have been widely tested in the framework of multi-frequency observations of the CMB, especially in the context of Planck (see Leach et al. 2008, , for a review of the methods developed within the Planck Collaboration). These methods exploit spectral and spatial correlations and independence in observations made at separate frequencies, as well as external constraints and physical modelling, in order to distinguish between different physical sources of emission. Most of them are designed to efficiently separate CMB from foregrounds and possibly foregrounds into “components”: Galactic synchrotron, free-free and dust emissions; extra-galactic and far-IR point sources; Sunyaev-Zeldovich effect, etc.

The approach presented here is different from these methods as it tries to quantify the separate components responsible for diffuse Galactic emission. Nevertheless it is not strictly speaking a component separation method, in the sense presented above, as it fits one map of the sky with spatial templates at different Galactocentric distances and is not exploiting the spectral signature of each component to perform the separation.

thumbnail Fig.2

Longitude profile of Hi, CO and dark gas. The black solid line shows the total contribution from Hi & CO. The dark gas represents a significant fraction of the gas column density, and dominates the CO outside of the molecular ring.

3.1. Inversion method

We assume that each of the frequency maps (Iν) can be expressed as a linear combination of the spatial templates described in Sect. 2.2 and presented in Fig. 1. Then, for each frequency map, we are solving for the mean emissivity of each template (in units of MJy sr-1/1020   atoms) which best reproduces the observed data. Note that we are not taking into account any spatial variations of the emissivity within a given template, but are assuming a constant value for each of our templates. The modelled emission can be written, for any frequency map Iν(x,y), as: (5)where the sum is over a number (n) of Galactocentric rings (i) for the atomic and molecular phases, denoted NHI and NH2 respectively. NHII,Ns and Nd denote the ionised, synchrotron and dark gas components, respectively, and Cν is a constant. Note that the last three templates do not have any rings, as we have no way of separating the information into different Galactocentric distances and so they are a function of frequency only. As such, the quoted ϵ are a function of ν only.

In order to find the values for the different emissivities, ϵ, we minimise the chi-squared value defined by (6)where σν is the noise in the map Iν, estimated as the median value of the intensity for |b| > 70°.

In our case, the difference between model and observations is not at all dominated by the low noise present in the Planck bands, but by the simplistic model that we are adopting for the inversion. As such we suppose that the noise does not vary spatially.

The solution which minimises Eq. (6) is found via matrix inversion. By writing the problem as (7)where b is the observed intensity map, A is an N × M matrix of the different templates and a is the vector of the parameters. The solution is simply (Press et al. 1992) (8)

The optimal mathematical solution which solves Eqs. (5)–(6) may not be physical; in particular, it may include negative emissivities. However, we do not have any cases where the emissivity is more than one sigma below zero, so remaining consistent with zero. As such we do not place any constraints on the parameters of the inversion.

thumbnail Fig.3

Results of the inversion for the 857 GHz band. From top to bottom are: observed emission; modelled emission; atomic contribution; molecular contribution; ionised contribution; dark gas contribution; synchrotron contribution; residual; and relative residual. All images except for the residual maps are at the same intensity scale. The centre and anti-centre regions have been masked, as the kinematic distance method is inapplicable to these regions. Areas where no CO observations are available have also been masked.

thumbnail Fig.4

Same as Fig. 3 for the 30 GHz band.

thumbnail Fig.5

Same as Fig. 3 for the 1.4 GHz band.

3.2. Galactocentric ring optimisation

The observed radial velocities of Hi and CO are used to separate the observed emission into a number of Galactocentric rings, assuming circular orbits for the gas and using the rotation curve from Clemens (1985). The choice of the distance intervals for the rings is not arbitrary, but is subject to a number of constraints. Choosing too many would create a set of rings that would be highly correlated and result in a degenerate solution. An optimised set of rings has been defined by minimising the correlation between rings, and thus minimising the degeneracy between the various components of the inversion. The total number of rings and their thickness are also driven by the requirement of a minimal number of data falling into each ring in order to constrain it and separate its contribution from the others. Further, we have chosen to fix the solar circle as one of the rings, including all velocity channels that contain high latitude, low velocity clouds.

The ring optimisation process consists first in dividing the Hi and CO templates into a large number of small rings, and then building the covariance matrix (ATA, see Eq. (8)) associated with this highly discrete basis of annuli. An iterative process is then performed: neighbouring rings with the highest correlation are merged and the new covariance matrix is updated with a reduced number of annuli. The process continues, merging the rings which have strongest correlation between them, until a reasonable number of rings is obtained. This leads to the following decomposition: RGAL =  [1.5,4.2,5.4,7.1,9.8,25.0]  kpc, where the values define the ring edges.

3.3. Estimating uncertainty

The result of the inversion is an average gas emissivity at each frequency, in each gas phase and in each ring (when available). However, within a given gas phase the dust and gas properties may vary significantly and we wish to report this uncertainty. In other words, the residual differences between the model and the data are not due to noise in the data, meaning that the uncertainty on the resulting emissivities cannot be computed analytically (i.e., the residuals are not well-modelled by the uniform Gaussian noise implicit in the solution, Eq. (8)). The residuals are rather due to unknown differences between the templates and the data. An alternate method is then needed to estimate the uncertainties on the results.

The dust properties in the northern Galaxy may differ from those in the southern part, and of course the variation may be even more localised than that. To quantify this uncertainty we perform jackknife tests by masking approximately 30% of the Galaxy in 30 random blocks, each 3° wide, and performing the inversion on the unmasked zones. This is repeated 500 times, and the rms value of the jackknife tests is used to provide the total variation of the emissivity within each ring and gas phase. We therefore provide robust estimates of the dust emission and provide their variation for each ring and gas phase.

It should be stressed that this estimation does not take into account the errors on either the Hi opacity correction nor the uncertainty on the assumed XCO value for each phase globally.

4. Properties of the ISM in the Galactic plane

The above inversion method has been applied to 16 different frequency bands from 12μm (25 THz) to 1.4 GHz. The results are presented in table form in the appendix (Table A.1).

Figures 35 compare our best-fit model to the data for three frequencies (two Planck frequencies, 857 and 30 GHz, as well as the 1.4 GHz radio frequency) and the contributions of the different gas components to the model. The model reproduces the data rather well, especially when one notices the scale of the residual.

The main contributors to the 857 GHz map (see Fig. 3) are the atomic and molecular phases (including the dark gas) with a modest contribution from the ionised component. At 30 GHz (Fig. 4) the emission is a combination of the molecular, ionised, dark gas and synchrotron, with very little hydrogen. As for the solution at 1.4 GHz (Fig. 5) it is the synchrotron component that dominates the solution, but with a significant contribution from the ionised component.

4.1. Dust in thermal equilibrium

thumbnail Fig.6

SEDs for the atomic component, each sorted in increasing Galactocentric radius. The colour of the symbols refer to the mission the data came from: IRAS (blue), DIRBE (yellow), HFI (red), LFI (green) and 1.4 GHz (black). The SEDs also show various fitted laws: the dotted line is the thermal dust SED, the short dashed line spinning dust in an atomic medium, and the solid line the sum of all contributions. The spinning dust contribution is not a fit, but simply the result of using typical values in the model (see Sect. 4.2.3 for details). The contribution of the atomic component in the outer Galaxy to the 1.4 GHz signal is lower than the noise level derived for the 1.4 GHz map and so should not be viewed as significant.

We have constructed SEDs for each of our templates across all frequencies that we have analysed. Where the emissivity is compatible with zero, we plot the 2σ upper limit as a downward arrow.

We have fit the SEDs using a sum of thermal dust (big grains), spinning PAH molecules, free-free and synchrotron emission. The large dust grains in thermal equilibrium are well described by a modified blackbody of the form (9)where T is the dust temperature and β is the emissivity spectral index. When fitting using both T and β as free parameters, we find a mean value of β = 1.8 ± 0.1 throughout the Galaxy. However, to increase the reliability of the fit, we follow Planck Collaboration (2011t,u,o) and fix β = 1.8. We use all data points down to 143 GHz in the fit, but exclude 217 GHz in the molecular and dark gas phases due to contamination from CO line emission. All points are colour corrected assuming νIν = constant as described in Planck HFI Core Team (2011b). The SEDs for the atomic, molecular and ionised phases are shown in Figs. 68, respectively. We are able to constrain the dust temperature for each component and in the various Galactocentric radii. The resulting temperature profile is shown in Fig. 9. The temperature in the atomic phase is seen to decline steadily with increasing Galactocentric distance. This decrease in temperature follows the trend of the ambient ISRF (Mathis et al. 1983), shown by the solid line in Fig. 9, indicating that the dust in the atomic phase is predominantly heated by it.

The dust in the molecular component, however, is seen to be more steady, peaking in the molecular ring (3 ≲ RG ≲ 6 kpc) and once again outside the solar circle. This phase seems to be heated by the presence of star formation. The ionised phase (dotted line) is seen to be dominated by grains that are on average much warmer than either of the other two phases. These results are in agreement with previous Galactic analyses (Sodroski et al. 1997; Paladini et al. 2007). The temperature within the solar circle (T = 17.6 ± 0.1   K) is completely in agreement with the average value at high Galactic latitude (Planck Collaboration 2011t).

The dust opacity at a given wavelength is defined as τλ/NH where τλ = ϵλ/B(ν,T) is the dust optical depth. The dust opacity in the atomic phase for the solar circle is equal to 0.92 ± 0.05 × 10-25   cm2 at 250 μm, slightly below, but in agreement with, the same value found by Planck Collaboration (2011t) at high Galactic latitude. The variation of the opacity as a function of Galactic radius is shown in Fig. 10. We do not detect any significant variation of the dust opacity as a function of Galactic radius, in either the atomic or molecular phases. Note, however, that the opacity of the molecular component is completely degenerate with the chosen value for XCO. By increasing XCO, we decrease the opacity. Further, the appropriate spin temperature for the Hi opacity correction may differ from the constant 250 K we have assumed.

Finally, after fitting the thermal dust in the molecular phase, there is an obvious excess at 100 and 217 GHz due to the presence of CO line emission in these Planck bands (Planck HFI Core Team 2011b). This is a clear indication that the inversion is working as it should, as this excess is present in the molecular phase. The excess with respect to the fit is 16.2 ± 1.0 and 20.2 ± 4.5   μKCMB/KRJkms-1 at 100 and 217 GHz, respectively. Planck HFI Core Team (2011b) find 14.2 ± 1.0 and 44.2 ± 6.0 at 100 and 217 GHz, respectively. At 100 GHz our value is slightly higher, which is to be expected as it is simply an excess with respect to 12CO and does not take into account the other CO isotopes present along the line of sight. The difference in the values at 217 GHz is not as easy to understand, but may come from our simple assumption of a fixed β.

4.2. Analysis of non-thermal SED components

The SEDs for all components at frequencies lower than 100 GHz show non-thermal dust emission features. The emission mechanisms at work in this spectral range are thermal bremsstrahlung (free-free), synchrotron emission, as well as a third commonly referred to as anomalous dust emission.

4.2.1. Synchrotron

The SED for the synchrotron component resulting from our inversion is shown in Fig. 11. A power-law has been recovered in the two lowest frequency LFI points and the 1.4 GHz radio point.

For electrons with a power-law distribution of energies N(E) ∝ E − p the frequency dependence of the emission is characterised by antenna temperature T(ν) ∝ ν   β with spectral index β =  − (p + 3)/2, with typically β ~  −3. As we are not working in antenna temperature, but flux density, the power law takes the form Iν ∝ ν   α with α = β − 2.

We have performed a fit of the rising synchrotron spectrum with a power law and obtained a spectral index α =  −1.0 ± 0.1. This value is slightly flatter than the average value found between 408 MHz and 23 GHz by WMAP (Dunkley et al. 2009; Davies et al. 2006). However, this value falls within the range of values found in the Galaxy (Bennett et al. 2003,  −0.5 < α <  −1.1), and most likely represents the stronger power injection in the Galactic plane. There is no evidence in the data for any hardening of the synchrotron spectrum below 70 GHz.

thumbnail Fig.7

SEDs for the CO component, each sorted in increasing Galactocentric radius. The different lines and colours are as described in Fig. 6, except for the long dashed line which represents spinning dust in a molecular medium.

thumbnail Fig.8

The resulting SEDs for the free-free component. The different lines and colours are as described in Fig. 6, except for the dash-dot line and the dash triple dot line which represent free-free emission and spinning dust in an ionised medium, respectively.

thumbnail Fig.9

Average dust temperature as a function of Galactocentric radius for the atomic (blue) and molecular (red) components. The dust temperature in the ionised phase and the dark gas phase are shown by the dashed line and dot-dashed line, respectively. The solid line shows the behaviour of the ISRF with Galactic radius. Horizontal bars show the range in distance over which the mean temperature applies. The temperature of dust mixed with atomic Hi decreases with radius, consistent with the fact that this component is heated by the ambient ISRF. Dust in the molecular component is heated by embedded star formation.

thumbnail Fig.10

Dust opacity at 250μm as a function of Galactocentric radii for Hi (blue, diamonds) and CO (red, triangles). The dashed line and the dash dot line show the dust opacity for the ionised phase and the dark gas phase, respectively.

thumbnail Fig.11

The resulting SEDs for synchrotron component. The different lines and colours are as described in Fig. 6. The point at 408 MHz is not part of the inversion, but has been plotted here to show consistency.

4.2.2. Free-free

Free-free emission dominates the ionised hydrogen template (Fig. 8) at frequencies below 143 GHz. This emission follows a power-law with a spectral index of α =  −0.1 for 70 GHz < ν < 1.4 GHz. As we have used the WMAP free-free template to trace the ionised medium, it is interesting to compare the total free-free emission that we attribute to this phase with the WMAP estimate. Using the power-law fit shown in Fig. 8, we can extrapolate to the frequency of our choice and, using our H+ template, produce an estimate of the free-free emission in any of the WMAP bands.

Performing this at 23 GHz, we recover  ~80% of the emission that the WMAP team find in their MEM free-free component separation. Using radio recombination lines (RRLs) from the Hi Parkes zone of avoidance (ZOA) survey, Alves et al. (2010) recently showed that diffuse ionised emission can be recovered in the Galactic plane using RRL surveys. They then compared the level of free-free derived with the WMAP MEM estimate. As in this work, they showed it to be too high in the Galactic plane by at least 10%, albeit in a much smaller (8° × 8°) region.

Looking at the molecular SEDs (Fig. 7), and the dark gas SED (Fig. 13), the LFI points depart from the thermal dust and seem to describe a free-free spectrum. It is not possible, however, to describe the three LFI points and the radio point using a free-free power law, except for the first and last molecular templates. Even so, all the radio points in the molecular templates are upper limits and the predicted free-free spectrum would be above the 4-σ level for three of the templates, showing that the data are not compatible with this model. In addition, if this were indeed free-free emission, our estimate of the total Galactic free-free emission would surpass that of the WMAP MEM free-free estimate by over 30%. This seems unlikely in light of the Alves et al. (2010) finding discussed above. In that case, the LFI-radio points in the atomic, molecular and the dark gas phases are anomalous microwave emission. Indeed, Boughn & Pober (2007) recently showed that the dominant emission mechanism at 20 GHz is spinning dust, and so is identified as anomalous microwave emission.

4.2.3. Anomalous microwave emission

Anomalous microwave emission, arising between 10 and 100 GHz and first detected in the 1990s (Kogut et al. 1996), cannot be explained using classical emission mechanisms known in this frequency range (Banday et al. 2003). Indeed it is too bright to be free-free emission, and it is not polarised as thermal dust or synchrotron emissions would be. However, it is correlated with dust IR emission and especially with the interstellar PAH and very-small-grain (VSG) emission in the mid-IR (Casassus et al. 2006; Ysard et al. 2010; Scaife et al. 2010). Anomalous emission was first proposed by Draine & Lazarian (1998) to come from rapidly rotating PAHs (Rouan et al. 1992), the spinning dust grains. Since, it has been observed in various interstellar environments: dark clouds (Watson et al. 2005; Casassus et al. 2006; Scaife et al. 2009; Dickinson et al. 2010); Hii regions (Dickinson et al. 2007; Scaife et al. 2008); planetary nebulae (Casassus et al. 2007) and diffuse ISM (Miville-Deschênes et al. 2008). Several processes can excite or damp the grains’ rotation: photon emission (IR and radio); gas/grain interactions; formation of H2 molecules on the grain surface; or photoelectric emission.

thumbnail Fig.12

Depiction of the spinning dust emission at 30 GHz in the Galactic plane in MJy sr-1 (top) and as a fraction of the total 30 GHz signal (bottom). In making this image, it is assumed that the 30 GHz emission associated with the atomic, molecular and dark gas phases are solely due to spinning dust. Some features may be simply due to residuals in our inversion; however, this illustrative view highlights regions where spinning dust may be particularly strong compared to other emission mechanisms.

Here we investigate the anomalous component in CO-correlated emission. All the rings exhibit a flattening of the SED in the three LFI bands (Fig. 7), as well as upper limits on the emission at 1.4 GHz. This clear excess in the microwave range cannot be free-free emission due to the reasons discussed in Sect. 4.2.2. The peak frequency of this excess lies approximately from 20 to 60 GHz.

We attempt to fit the excess with the spinning dust model described in Silsbee et al. (2011). Models show that spinning dust emission is sensitive to the grain size distribution, to the intensity of the ISRF, and to the density of the medium in which the grains are embedded (Ali-Haïmoud et al. 2009; Ysard & Verstraete 2010). The smallest grains, namely the PAHs, are responsible for most of the spinning dust emissivity. Fitting the thermal dust emission with a modified blackbody, we estimate the radiation field to be close to the standard ISRF (Mathis et al. 1983; Porter et al. 2008). To describe the small dust grains, we use two grain populations: PAHs and VSGs. We adopt a size distribution that is the sum of two log-normal components (Weingartner & Draine 2001). The first one is peaked at 3 nm (VSG), and the second one at 0.6 nm (PAHs), in line with Compiègne et al. (2011), with the width of both log-normal components set to 0.4. Most of the spinning dust emission comes from the smallest grains, the second log-normal component, for which we assume an abundance equal to the value measured in the solar neighbourhood. Finally we describe the grains’ environment as molecular gas with typical parameters (nH = 350 cm-3 and Tgas = 20 K).

The fit to the data with this simple model is surprisingly good, and shows that not only is anomalous emission present in the molecular phase throughout the Milky Way, but that spinning dust provides a compelling model with which to describe it. The predominance of spinning dust in the molecular phase is in agreement with more detailed modelling of individual anomalous emission regions (Planck Collaboration 2011p).

The ISRF could in theory also have an effect on this emission mechanism, and the big grain temperature gives us an idea of its fluctuations: χ = (T/17.5   K)β + 4, where χ is dimensionless and equal to one in the solar neighbourhood. However this range of ISRF fluctuations has little effect on the spinning dust spectrum (Ali-Haïmoud et al. 2009; Ysard & Verstraete 2010), so we assume a solar value for all SEDs. The average density of the molecular gas may vary with radius, which would have enough impact to significantly alter the fit. However, we have not proceeded with a detailed fit of this component, as there are many parameters and the solution would be degenerate (e.g. grain size distribution, gas density, ISRF intensity, as well as radiative transfer considerations). We have kept the analysis simple, which allows us to say that we have definitely detected anomalous microwave emission, that it is associated with the molecular phase, and that spinning dust provides a good model with which to understand the spectrum of this excess emission. Some of this emission could be explained by free-free emission, but it cannot be the dominant emission mechanism in this phase.

To study the possibility that spinning dust emission is present in either the atomic or ionised phase, we have also estimated its contribution to these. The different environments will have an impact on the spinning dust emission, so we have adapted the gas parameters in the modelling. For the atomic phase we have used the results of Heiles & Troland (2003) and modelled the atomic ISM as composed of 40% cold neutral medium (CNM) and 60% warm neutral medium (WNM). The parameters for these phases are nH = 30 cm-3 and Tgas = 100 K for the CNM, and nH = 0.4 cm-3 and Tgas = 6000 K for the WNM. As can be seen in Fig. 6, the SEDs show that spinning dust emission is able to reproduce the departure from the thermal dust spectrum for ν < 70 GHz. The 1.4 GHz point in the last atomic template is almost certainly an artifact. The atomic and synchrotron templates both describe diffuse smooth emission in the Galaxy, leading to some “crosstalk” between the two components. However it should be noted that when the emissivity found for the atomic phase in the outer Galaxy is multiplied by the column density of hydrogen in the template, the total emission is very low, lower than the noise level derived for the 1.4 GHz map. As such, it should not be viewed as significant.

For the ionised phase, we use the same gas parameters as in Sect. 2.2, namely nH = 10 cm-3 and Tgas = 8000 K. The SED in Fig. 8 shows that spinning dust may also be present in this phase, but that it is completely dominated by free-free emission. In fact, the situation may be even worse, as we have assumed solar value PAH abundance, whereas recently Dobler et al. (2009) showed that the PAH abundance in this type of region is about a factor of three lower.

Assuming that the atomic, molecular and dark gas contributions at 30 GHz are due solely to spinning dust, and that the ionised phase has a negligible contribution, we have produced a Galactic-plane image of spinning dust emission. This provides a complementary image to that presented in Planck Collaboration (2011p), as it has been obtained using a different approach. The intensity map is shown at the top of Fig. 12; at the bottom we show the fraction of the total signal at 30 GHz due to spinning dust. Several bright regions show up in this map, highlighting areas where spinning dust emission is the dominant source of emission at these frequencies. Overall, spinning dust accounts for 25 ± 5% (statistical) of the total emission at 30 GHz. This error estimate does not take into account any systematic uncertainty due to our simple model, which is difficult to quantify, and so will likely be higher.

4.3. Dark gas

Dark gas has been detected in numerous studies (Reach et al. 1994; Grenier et al. 2005) as an infrared, and a correlated infrared/gamma-ray, excess over and above what is expected given the gas column density derived from Hi and CO data. The inversion performed without this component results in positive and negative residuals, as the other components compete to account for this missing gas phase. Not only does adding this component lower our reduced χ2 fits across all bands, but it also allows us to examine the nature of this dark component.

thumbnail Fig.13

The resulting SED for the dark gas component. The different lines and colours are as described in Fig. 6. The dashed and short-dashed lines shows the spinning dust contribution from a molecular phase, and from an atomic phase, respectively.

Looking at Fig. 13, the resulting SED for this component looks quite similar to both the atomic and molecular SEDs. Its temperature does not differentiate it immediately from either of the two phases, and it shows a low-frequency tail beyond the thermal dust SED. These points look similar to the anomalous dust in the molecular phase (dashed line); however, the spinning dust model lies above the LFI points. Assuming that the spinning dust is in an atomic phase (dotted line), it can be seen that it does not predict the observed emission. It thus seems to be more of a mix of the two phases, with a larger influence from the molecular phase. The different parameters discussed in Sect. 4.2.3 could also be altered, but we do not have enough data to constrain such a fit.

The interpretation of these results is somewhat complicated for several reasons. Firstly, the amount of dark gas is not expected to be independent of Galactic position, but using a 2D template we are forced to assume just that. Secondly, the estimate of this gas phase was obtained at high galactic latitude and we are using an extrapolation into the Galactic plane. It would be preferable to use an independent template, derived from gamma-ray observations for example, in future inversion analyses. Reconstructing the line-of-sight distribution of the dark gas would enable a more exact determination of the properties of the dark gas, by coupling the analysis with a 3D extinction map (e.g. Drimmel et al. 2003; Marshall et al. 2006; Sale et al. 2009). This type of analysis will be considered in a future study of the Galactic inversion.

In summary, the dark gas SED is consistent with molecular gas, most likely diffuse H2 and optically thick CO. The best fit at longer wavelengths comes from the spinning dust model, in a molecular medium.

thumbnail Fig.14

Correlation matrix of the emissivities for each ring at 30 GHz. Overplotted in red are the results with corresponding error bars. One-dimensional histograms replace the contour plots where only one component is concerned.

5. Possible sources of bias

Our inversion algorithm determines the flux density per column density of each template, and hence the result depends on the different templates themselves. Care should therefore be taken to use the most physically representative input for each phase of the ISM that we want to include in the inversion. Nonetheless, the jackknife tests described in Sect. 3.3 show that the method is self-consistent. Moreover, many different combinations for the Hi+CO Galactocentric rings has been attempted, and the main results presented here remain coherent with those. However, the results are notably sensitive to one particular input, the synchrotron template. Without it, the Hi templates trace the Galactic synchrotron emission, as the Hi observations trace the large-scale smooth structure of the Milky Way, similar to the diffuse soft synchrotron emission. This artificial component in the Hi templates then masks any contribution for the anomalous microwave emission. We have attempted to trace the synchrotron emission with two different synchrotron templates: the WMAP MCMC synchrotron template (Gold et al. 2009), and the 408 MHz full sky survey (Haslam et al. 1982). We have decided to use the latter as it more successfully removes the synchrotron contribution from the Hi templates and allows us to detect anomalous emission in the atomic phase.

5.1. Correlation between parameters and uniqueness of the solution

In Fig. 14 we show the correlation between the different emissivities of our best solution at 30 GHz, where many emission mechanisms are present and the inversion is more challenging. These were obtained from the result of the jackknife tests presented in Sect. 3.3. It can be seen that for almost all the templates the correlation is low. Three elements are to be highlighted. The first one is that the emissivities for the outer Galaxy molecular template, as well as the ionised template, both have pronounced departures from a Gaussian distribution. Looking at Fig. 1, both templates are patchy; the masking of several regions simultaneously will change the morphology of each map significantly and therefore have a large impact on the derived emissivity.

Secondly, rings 1–3 of the atomic and molecular phase are slightly anti-correlated (the contours are elliptical). This region of the Galaxy is indeed difficult to disentangle, and due to the low filling factor of these rings the algorithm has less to work with than in the solar circle, for example. Nevertheless the error bars deduced from the jackknife analysis account for this slight degeneracy. Higher resolution maps will help reduce this crosstalk between templates in future analyses.

Thirdly a small bias can be seen, as the results of the inversion are not always centred in the middle of the jackknife distributions. We calculated this bias by looking at the difference between the mean deviation of each result with respect to the jackknife realisations, in units of result uncertainty. That is the bias B per template i and per band ν can be written: (10)where Nr is the number of jackknife realisations, ϵ(i,ν) and σ(i,ν) are the result and uncertainty of the inversion, and  − ϵJ(r,i,ν) is the result of one jackknife realisation. We have found this effect to be small, with a median bias of 0.005σ over all templates and frequencies, and with a mean absolute deviation of 0.08σ. In all cases, the bias is always lower than 0.25σ.

5.2. Effect of CMB subtraction

The CMB removal in the Planck maps we have used may create significant residuals in the Galactic plane. Therefore, to test the impact of the CMB removal on our results, a different CMB estimate was reconstructed through a classical internal linear combination (ILC) by means of Lagrange multipliers (Eriksen et al. 2004). A  ± 15° strip in Galactic latitude was extracted from the HFI CMB frequency channels maps (100, 143, 217 and 353 GHz) reduced to a common resolution (10 arcmin), in units of KCMB. This CMB component was subtracted from all the Planck LFI and HFI channels maps, resulting in a new and independent set of CMB-removed Planck maps.

Even though the methods are similar, the resulting CMB in the Galactic plane is quite different. In constructing the the standard CMB subtraction (Needlet ILC, “NILC”), an intensity mask was applied (Planck HFI Core Team 2011b) which results in CMB residuals remaining in a substantial part of the Galactic plane. No such mask was used in the Galactic plane ILC. Outside of this mask, the differences between the two CMB estimates are of the same order as the residuals of our fit at low frequencies, 44 GHz    < ν < 143 GHz.

thumbnail Fig.15

Comparison of emissivities derived after application of two different CMB subtractions, the official Needlet ILC and the Galactic plane ILC. The atomic templates are the blue points, molecular are black, ionised is orange, dark gas is red and synchrotron is green.

Using these maps to perform the same inversion, we have concluded that the impact is small, and almost always within the error bars. In Fig. 15 we plot the derived emissivities using the NILC and the Galactic plane centred CMB removal (“ILC”).

To get a better idea of where the solution might be off by more than 1σ, we show the difference in units of σ per band in Fig. 16. The value plotted is, for each component: (11)The agreement is generally very good, except for some points concentrated in the bands 143, 100 and 70 GHz. In these cases the CO templates are higher in the Galactic plane ILC whereas the Hi and H+ templates are lower. Nevertheless, there are no points that are different by more than 3σ. Therefore, these differences do not change our main conclusions, but would slightly modify the amount of anomalous emission, for example. As we are not modelling the specifics of any emission mechanisms at the lower frequencies, we are not sensitive to the impact of the CMB removal for this study.

5.3. Impact of using dark gas template

thumbnail Fig.16

Comparison of emissivities derived after application of two different CMB subtractions. The value plotted is shown in Eq. (11). The symbols represent Hi (plus), H2 (diamond), H+ (triangle), Synchrotron (square) and dark gas (asterisk). Differences are generally small, but become larger in the range 143–44 GHz.

The dark gas template used in this study has been calculated by Planck Collaboration (2011o). It is the result of finding at what column densities the correlation between total hydrogen column density (Hi+2) versus optical depth of the dust departs from linearity. The frequencies used to obtain the map at 857 GHz that we use as a template are 857 and 545 GHz, along with the IRAS 100μm band. There is some circularity in using a template based on a residual in the Planck data to fit the same data. However, the data used is restricted to the two highest HFI frequencies. Here we present the impact at 30 GHz, which in no respect suffers from the same circularity, as the data are far from the frequencies used to derive the dark gas template, and they come from LFI, a completely separate instrument from HFI.

We have performed the inversion with and without the dark gas template at 30 GHz, and looked at the correlation between the residuals and the column density of the dark gas. The results are shown in Fig. 17; the inversion performed without the dark gas template is on the left, the same inversion with the dark gas template is on the right. A distinct component of the residuals can be seen to correlate very well with the dark gas column density when the dark gas template is not used. This correlation completely disappears when it is used. In addition, the reduced χ2 for the fit shows a significant improvement of approximately 50%.

The choice of both the opacity correction for the Hi and the XCO value chosen will have an impact on our results. We have chosen to keep the analysis simple for this study, as the spin temperature we assume, together with the dark gas template, should account for most of the optically thick gas in the plane. The ability of the dark gas template used to adequately account for all the optically thick gas remains to be studied, so future work could involve using more sophisticated methods to account for this effect.

6. Conclusions

The analysis performed in this study provides a realistic description of the dust and gas properties as a function of Galactic radius in the Milky Way. The dust temperature in the Hi is seen to decrease as a function of Galactocentric distance from 24.0 to 13.9 K, and the temperature in the molecular phase is heated by star formation associated to the Galactic spiral arms. The opacity found for grains in the solar circle is τ/NH = 0.92 ± 0.05 × 10-25   cm2 at 250 μm, and with no significant variation with Galactic radius, even though the dust temperature is seen to drop by over 10 K from the centre to the outer Galaxy. The dust temperature in the Hi gas in the solar circle, 17.6 ± 0.1 K is also compatible with the recent Planck value of high-latitude cirrus (17.8 ± 0.9 K).

The extension of our analysis to lower frequencies than previously used in inversion techniques has allowed us to place constraints on the free-free, synchrotron and anomalous microwave emission:

  • anomalous dust emission is clearly seen in the atomic,molecular and dark gas phases. It is well fit by a very simplemodel consisting of spinning PAH molecules embedded ineach of the gas phases. We highlight regions where spinningdust emission may be strong with respect to other emissionmechanisms. According to our simple model, in the Galacticplane spinning dust accounts for 25 ± 5%n (statistical) of the total emission at 30 GHz. However, systematic uncertainties linked to our model may make this value more uncertain;

    thumbnail Fig.17

    Dark gas intensity as a function of the residuals of the inversion at 30 GHz without (Left) and with (Right) the dark gas template. A distinct component of the residual (shown with the solid line) is seen to be correlated with the dark gas template. This correlated residual disappears when the dark gas is used in the inversion.

  • the dark gas phase has been explored spectrally across the Galactic plane for the first time. Its SED is similar to the molecular phase and therefore this phase seems to be tracing diffuse H 2 where CO is disassociated by energetic photons, as well as tracing optically thick CO and Hi;

  • the free-free emission is completely traced by our ionised component and indicates that the free-free estimate from WMAP component separation methods may be too high. In the ionised phase, spinning dust may be present but it is dominated by free-free emission;

  • synchrotron emission is well characterised by a power law with a spectral index of  −1.0.

This first analysis of the large-scale Galactic emission as seen by Planck has constrained the large scale ISM properties in the Galaxy. It could be improved in the future via several modifications. Individual regions could be studied using higher angular resolution data, which would probe smaller scale fluctuations that have been missed by the present study. The radial separation could be done differently to explore the difference between arm and inter-arm ISM, in order to focus more on the role of ISM conditions on the star formation process. It should also be possible to obtain an independent estimate of the dark gas in the plane, and ideally some radial information about its behaviour. In the near future, Planck component separation products could be used in the place of WMAP ones. Use of radio recombination line surveys currently underway would be very useful to constrain the free-free emission in the plane, and as a function of Galactic radius.


1

Planck (http://www.esa.int/Planck) is a project of the European Space Agency (ESA) with instruments provided by two scientific consortia funded by ESA member states (in particular the lead countries: France and Italy) with contributions from NASA (USA), and telescope reflectors provided in a collaboration between ESA and a scientific consortium led and funded by Denmark.

Acknowledgments

A description of the Planck Collaboration and a list of its members can be found at http://www.rssd.esa.int/index.php?project=PLANCK&page=Planck_Collaboration

References

  1. Abdo, A. A., Ackermann, M., Ajello, M., & Fermi/LAT Collaboration 2010, ApJ, 710, 133 [NASA ADS] [CrossRef] [Google Scholar]
  2. Ackermann, M., Ajello, M., Baldini, L., et al. 2011, ApJ, 726, 81 [NASA ADS] [CrossRef] [Google Scholar]
  3. Ali-Haïmoud, Y., Hirata, C. M., & Dickinson, C. 2009, MNRAS, 395, 1055 [NASA ADS] [CrossRef] [Google Scholar]
  4. Alves, M. I. R., Davies, R. D., Dickinson, C., et al. 2010, MNRAS, 405, 1654 [NASA ADS] [Google Scholar]
  5. Banday, A. J., Dickinson, C., Davies, R. D., Davis, R. J., & Górski, K. M. 2003, MNRAS, 345, 897 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bennett, C. L., Hill, R. S., Hinshaw, G., et al. 2003, ApJS, 148, 97 [NASA ADS] [CrossRef] [Google Scholar]
  7. Binney, J., & Merrifield, M. 1998, Galactic astronomy / James Binney and Michael Merrifield. Princeton [Google Scholar]
  8. Bloemen, J. B. G. M., Strong, A. W., Mayer-Hasselwander, H. A., et al. 1986, A&A, 154, 25 [NASA ADS] [Google Scholar]
  9. Bloemen, J. B. G. M., Deul, E. R., & Thaddeus, P. 1990, A&A, 233, 437 [NASA ADS] [Google Scholar]
  10. Boughn, S. P., & Pober, J. C. 2007, ApJ, 661, 938 [NASA ADS] [CrossRef] [Google Scholar]
  11. Casassus, S., Cabrera, G. F., Förster, F., et al. 2006, ApJ, 639, 951 [NASA ADS] [CrossRef] [Google Scholar]
  12. Casassus, S., Nyman, L.-Å., Dickinson, C., & Pearson, T. J. 2007, MNRAS, 382, 1607 [NASA ADS] [CrossRef] [Google Scholar]
  13. Clemens, D. P. 1985, Ph.D. Thesis Massachusetts Univ., 3 [Google Scholar]
  14. Compiègne, M., Verstraete, L., Jones, A., et al. 2011, A&A, 525, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Dame, T. M., Hartmann, D., & Thaddeus, P. 2001, ApJ, 547, 792 [NASA ADS] [CrossRef] [Google Scholar]
  16. Davies, R. D., Dickinson, C., Banday, A. J., et al. 2006, MNRAS, 370, 1125 [NASA ADS] [CrossRef] [Google Scholar]
  17. Dickey, J. M., McClure-Griffiths, N. M., Gaensler, B. M., & Green, A. J. 2003, ApJ, 585, 801 [NASA ADS] [CrossRef] [Google Scholar]
  18. Dickey, J. M., Strasser, S., Gaensler, B. M., et al. 2009, ApJ, 693, 1250 [NASA ADS] [CrossRef] [Google Scholar]
  19. Dickinson, C., Davies, R. D., Bronfman, L., et al. 2007, MNRAS, 379, 297 [NASA ADS] [CrossRef] [Google Scholar]
  20. Dickinson, C., Casassus, S., Davies, R. D., et al. 2010, MNRAS, 407, 2223 [NASA ADS] [CrossRef] [Google Scholar]
  21. Dobler, G., Draine, B., & Finkbeiner, D. P. 2009, ApJ, 699, 1374 [NASA ADS] [CrossRef] [Google Scholar]
  22. Draine, B. T., & Lazarian, A. 1998, ApJ, 494, L19 [NASA ADS] [CrossRef] [Google Scholar]
  23. Drimmel, R., Cabrera-Lavers, A., & López-Corredoira, M. 2003, A&A, 409, 205 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Dunkley, J., Komatsu, E., Nolta, M. R., et al. 2009, ApJS, 180, 306 [NASA ADS] [CrossRef] [Google Scholar]
  25. Eriksen, H. K., Banday, A. J., Górski, K. M., & Lilje, P. B. 2004, ApJ, 612, 633 [NASA ADS] [CrossRef] [Google Scholar]
  26. Finkbeiner, D. P. 2004, ApJ, 614, 186 [NASA ADS] [CrossRef] [Google Scholar]
  27. Giard, M., Lamarre, J. M., Pajot, F., & Serra, G. 1994, A&A, 286, 203 [NASA ADS] [Google Scholar]
  28. Gold, B., Bennett, C. L., Hill, R. S., et al. 2009, ApJS, 180, 265 [NASA ADS] [CrossRef] [Google Scholar]
  29. Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [NASA ADS] [CrossRef] [Google Scholar]
  30. Grenier, I. A., Casandjian, J.-M., & Terrier, R. 2005, Science, 307, 1292 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  31. Haslam, C. G. T., Salter, C. J., Stoffel, H., & Wilson, W. E. 1982, A&AS, 47, 1 [NASA ADS] [Google Scholar]
  32. Hauser, M. G., Arendt, R. G., Kelsall, T., et al. 1998, ApJ, 508, 25 [NASA ADS] [CrossRef] [Google Scholar]
  33. Heiles, C., & Troland, T. H. 2003, ApJ, 586, 1067 [NASA ADS] [CrossRef] [Google Scholar]
  34. Jarosik, N., Bennett, C. L., Dunkley, J., et al. 2011, ApJS, 192, 14 [NASA ADS] [CrossRef] [Google Scholar]
  35. Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Kogut, A., Banday, A. J., Bennett, C. L., et al. 1996, ApJ, 460, 1 [NASA ADS] [CrossRef] [Google Scholar]
  37. Leach, S. M., Cardoso, J., Baccigalupi, C., et al. 2008, A&A, 491, 597 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Mennella, A., Butler, R. C., Curto, A., et al. 2011, A&A, 536, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Marshall, D. J., Robin, A. C., Reylé, C., Schultheis, M., & Picaud, S. 2006, A&A, 453, 635 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Mathis, J. S., Mezger, P. G., & Panagia, N. 1983, A&A, 128, 212 [NASA ADS] [Google Scholar]
  41. Miville-Deschênes, M.-A., & Lagache, G. 2005, ApJS, 157, 302 [NASA ADS] [CrossRef] [Google Scholar]
  42. Miville-Deschênes, M.-A., Ysard, N., Lavabre, A., et al. 2008, A&A, 490, 1093 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Neugebauer, G., Habing, H. J., van Duinen, R., et al. 1984, ApJ, 278, L1 [Google Scholar]
  44. Paladini, R., Montier, L., Giard, M., et al. 2007, A&A, 465, 839 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Planck Collaboration 2011a, A&A, 536, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Planck Collaboration 2011b, A&A, 536, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Planck Collaboration 2011c, A&A, 536, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Planck Collaboration 2011d, A&A, 536, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Planck Collaboration 2011e, A&A, 536, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Planck Collaboration 2011f, A&A, 536, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Planck Collaboration 2011g, A&A, 536, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Planck Collaboration 2011h, A&A, 536, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Planck Collaboration 2011i, A&A, 536, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Planck Collaboration 2011j, A&A, 536, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Planck Collaboration 2011k, A&A, 536, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Planck Collaboration 2011l, A&A, 536, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Planck Collaboration 2011m, A&A, 536, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Planck Collaboration 2011n, A&A, 536, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Planck Collaboration 2011o, A&A, 536, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Planck Collaboration 2011p, A&A, 536, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Planck Collaboration 2011q, A&A, 536, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Planck Collaboration 2011r, A&A, 536, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Planck Collaboration 2011s, A&A, 536, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Planck Collaboration 2011t, A&A, 536, A24 [Google Scholar]
  65. Planck Collaboration 2011u, A&A, 536, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Planck Collaboration 2011v, The Explanatory Supplement to the Planck Early Release Compact Source Catalogue (ESA) [Google Scholar]
  67. Planck Collaboration 2011w, A&A, 536, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Planck HFI Core Team 2011a, A&A, 536, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Planck HFI Core Team 2011b, A&A, 536, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Porter, T. A., Moskalenko, I. V., Strong, A. W., Orlando, E., & Bouchet, L. 2008, ApJ, 682, 400 [NASA ADS] [CrossRef] [Google Scholar]
  71. Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in FORTRAN. The art of scientific computing, ed. W. H. Press, S. A. Teukolsky, W. T. Vetterling, & B. P. Flannery [Google Scholar]
  72. Reach, W. T., Koo, B.-C., & Heiles, C. 1994, ApJ, 429, 672 [NASA ADS] [CrossRef] [Google Scholar]
  73. Reich, P., & Reich, W. 1986, A&AS, 63, 205 [NASA ADS] [Google Scholar]
  74. Reich, P., Testori, J. C., & Reich, W. 2001, A&A, 376, 861 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Reich, W. 1982, A&AS, 48, 219 [NASA ADS] [Google Scholar]
  76. Rouan, D., Leger, A., Omont, A., & Giard, M. 1992, A&A, 253, 498 [NASA ADS] [Google Scholar]
  77. Sale, S. E., Drew, J. E., Unruh, Y. C., et al. 2009, MNRAS, 392, 497 [NASA ADS] [CrossRef] [Google Scholar]
  78. Scaife, A. M. M., Hurley-Walker, N., Davies, M. L., et al. 2008, MNRAS, 385, 809 [NASA ADS] [CrossRef] [Google Scholar]
  79. Scaife, A. M. M., Hurley-Walker, N., Green, D. A., et al. 2009, MNRAS, 400, 1394 [NASA ADS] [CrossRef] [Google Scholar]
  80. Scaife, A. M. M., Nikolic, B., Green, D. A., et al. 2010, MNRAS, 406, L45 [NASA ADS] [Google Scholar]
  81. Silsbee, K., Ali-Haïmoud, Y., & Hirata, C. M. 2011, MNRAS, 411, 2750 [NASA ADS] [CrossRef] [Google Scholar]
  82. Sodroski, T. J., Dwek, E., Hauser, M. G., & Kerr, F. J. 1989, ApJ, 336, 762 [NASA ADS] [CrossRef] [Google Scholar]
  83. Sodroski, T. J., Odegard, N., Arendt, R. G., et al. 1997, ApJ, 480, 173 [NASA ADS] [CrossRef] [Google Scholar]
  84. Tielens, A. G. G. M. 2005, The Physics and Chemistry of the Interstellar Medium, ed. A. G. G. M. Tielens [Google Scholar]
  85. Watson, R. A., Rebolo, R., Rubiño-Martín, J. A., et al. 2005, ApJ, 624, L89 [NASA ADS] [CrossRef] [Google Scholar]
  86. Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296 [NASA ADS] [CrossRef] [Google Scholar]
  87. Ysard, N., & Verstraete, L. 2010, A&A, 509, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Ysard, N., Miville-Deschênes, M. A., & Verstraete, L. 2010, A&A, 509, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Zacchei, A., Maino, D., Baccigalupi, C., et al. 2011, A&A, 536, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

Appendix A: Table of results

Table A.1

Result of the inversion per band and per component in MJy sr-1/1020 H atoms, except for the synchrotron component in which case it is in MJy sr-1(ν) per MJy sr-1(408 MHz).

All Tables

Table A.1

Result of the inversion per band and per component in MJy sr-1/1020 H atoms, except for the synchrotron component in which case it is in MJy sr-1(ν) per MJy sr-1(408 MHz).

All Figures

thumbnail Fig.1

Templates used in the inversion. The top five images are the Hi templates, the following five are the molecular templates, all of which are presented in increasing Galactocentric radius (the radii defining each ring are described in Sect. 3.2). The final three templates are the free-free, synchrotron, and dark gas templates, respectively. The units are all expressed in 1020 atoms, except for synchrotron, expressed in MJy sr-1. The centre and anti-centre regions have been masked as no distance information can be obtained from the radial velocity of the gas.

In the text
thumbnail Fig.2

Longitude profile of Hi, CO and dark gas. The black solid line shows the total contribution from Hi & CO. The dark gas represents a significant fraction of the gas column density, and dominates the CO outside of the molecular ring.

In the text
thumbnail Fig.3

Results of the inversion for the 857 GHz band. From top to bottom are: observed emission; modelled emission; atomic contribution; molecular contribution; ionised contribution; dark gas contribution; synchrotron contribution; residual; and relative residual. All images except for the residual maps are at the same intensity scale. The centre and anti-centre regions have been masked, as the kinematic distance method is inapplicable to these regions. Areas where no CO observations are available have also been masked.

In the text
thumbnail Fig.4

Same as Fig. 3 for the 30 GHz band.

In the text
thumbnail Fig.5

Same as Fig. 3 for the 1.4 GHz band.

In the text
thumbnail Fig.6

SEDs for the atomic component, each sorted in increasing Galactocentric radius. The colour of the symbols refer to the mission the data came from: IRAS (blue), DIRBE (yellow), HFI (red), LFI (green) and 1.4 GHz (black). The SEDs also show various fitted laws: the dotted line is the thermal dust SED, the short dashed line spinning dust in an atomic medium, and the solid line the sum of all contributions. The spinning dust contribution is not a fit, but simply the result of using typical values in the model (see Sect. 4.2.3 for details). The contribution of the atomic component in the outer Galaxy to the 1.4 GHz signal is lower than the noise level derived for the 1.4 GHz map and so should not be viewed as significant.

In the text
thumbnail Fig.7

SEDs for the CO component, each sorted in increasing Galactocentric radius. The different lines and colours are as described in Fig. 6, except for the long dashed line which represents spinning dust in a molecular medium.

In the text
thumbnail Fig.8

The resulting SEDs for the free-free component. The different lines and colours are as described in Fig. 6, except for the dash-dot line and the dash triple dot line which represent free-free emission and spinning dust in an ionised medium, respectively.

In the text
thumbnail Fig.9

Average dust temperature as a function of Galactocentric radius for the atomic (blue) and molecular (red) components. The dust temperature in the ionised phase and the dark gas phase are shown by the dashed line and dot-dashed line, respectively. The solid line shows the behaviour of the ISRF with Galactic radius. Horizontal bars show the range in distance over which the mean temperature applies. The temperature of dust mixed with atomic Hi decreases with radius, consistent with the fact that this component is heated by the ambient ISRF. Dust in the molecular component is heated by embedded star formation.

In the text
thumbnail Fig.10

Dust opacity at 250μm as a function of Galactocentric radii for Hi (blue, diamonds) and CO (red, triangles). The dashed line and the dash dot line show the dust opacity for the ionised phase and the dark gas phase, respectively.

In the text
thumbnail Fig.11

The resulting SEDs for synchrotron component. The different lines and colours are as described in Fig. 6. The point at 408 MHz is not part of the inversion, but has been plotted here to show consistency.

In the text
thumbnail Fig.12

Depiction of the spinning dust emission at 30 GHz in the Galactic plane in MJy sr-1 (top) and as a fraction of the total 30 GHz signal (bottom). In making this image, it is assumed that the 30 GHz emission associated with the atomic, molecular and dark gas phases are solely due to spinning dust. Some features may be simply due to residuals in our inversion; however, this illustrative view highlights regions where spinning dust may be particularly strong compared to other emission mechanisms.

In the text
thumbnail Fig.13

The resulting SED for the dark gas component. The different lines and colours are as described in Fig. 6. The dashed and short-dashed lines shows the spinning dust contribution from a molecular phase, and from an atomic phase, respectively.

In the text
thumbnail Fig.14

Correlation matrix of the emissivities for each ring at 30 GHz. Overplotted in red are the results with corresponding error bars. One-dimensional histograms replace the contour plots where only one component is concerned.

In the text
thumbnail Fig.15

Comparison of emissivities derived after application of two different CMB subtractions, the official Needlet ILC and the Galactic plane ILC. The atomic templates are the blue points, molecular are black, ionised is orange, dark gas is red and synchrotron is green.

In the text
thumbnail Fig.16

Comparison of emissivities derived after application of two different CMB subtractions. The value plotted is shown in Eq. (11). The symbols represent Hi (plus), H2 (diamond), H+ (triangle), Synchrotron (square) and dark gas (asterisk). Differences are generally small, but become larger in the range 143–44 GHz.

In the text
thumbnail Fig.17

Dark gas intensity as a function of the residuals of the inversion at 30 GHz without (Left) and with (Right) the dark gas template. A distinct component of the residual (shown with the solid line) is seen to be correlated with the dark gas template. This correlated residual disappears when the dark gas is used in the inversion.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.