Dust in the diffuse emission of the galactic plane - The Herschel/Spitzer SED fitting

The first Herschel Hi-Gal images of the galactic plane unveil the far-infrared diffuse emission of the interstellar medium with an unprecedented angular resolution and sensitivity. In this paper, we present the first analysis of these data in combination with that of Spitzer Glimpse & Mipsgal. We selected a relatively diffuse and low excitation region of the l~59\,^{\circ} Hi-Gal Science Demonstration Phase field to perform a pixel by pixel fitting of the 8 to 500 microns SED using the DustEM dust emission model. We derived maps of the Very Small Grains (VSG) and PAH abundances from the model. Our analysis allows us to illustrate that the Aromatic Infrared Bands (AIB) intensity does not trace necessarily the PAH abundance but rather the product of "abundance x column density x intensity of the exciting radiation field". We show that the spatial structure of PACS70microns map resembles the shorter wavelengths (e.g. IRAC8microns) maps, because they trace both the intensity of exciting radiation field and column density. We also show that the modeled VSG contribution to PACS70microns (PACS160microns) band intensity can be up to 50% (7%). The interpretation of diffuse emission spectra at these wavelengths must take stochastically heated particles into account. Finally, this preliminary study emphasizes the potential of analyzing the full dust SED sampled by Herschel and Spitzer data, with a physical dust model (DustEM) to reach the properties of the dust at simultaneously large and small scales.Comment: accepted for publication in ApJ

Spitzer characterisation of dust in an anomalous emission region: the Perseus cloud

Anomalous microwave emission is known to exist in the Perseus cloud. One of the most promising candidates to explain this excess of emission is electric dipole radiation from rapidly rotating very small dust grains, commonly referred to as spinning dust. Photometric data obtained with the Spitzer Space Telescope have been reprocessed and used in conjunction with the dust emission model DUSTEM to characterise the properties of the dust within the cloud. This analysis has allowed us to constrain spatial variations in the strength of the interstellar radiation field ($\chi_\mathrm{ISRF}$), the mass abundances of the PAHs and VSGs relative to the BGs (Y$_\mathrm{PAH}$ and Y$_\mathrm{VSG}$), the column density of hydrogen (N$_\mathrm{H}$) and the equilibrium dust temperature (T$_\mathrm{dust}$). The parameter maps of Y$_\mathrm{PAH}$, Y$_\mathrm{VSG}$ and $\chi_\mathrm{ISRF}$ are the first of their kind to be produced for the Perseus cloud, and we used these maps to investigate the physical conditions in which anomalous emission is observed. We find that in regions of anomalous emission the strength of the ISRF, and consequently the equilibrium temperature of the dust, is enhanced while there is no significant variation in the abundances of the PAHs and the VSGs or the column density of hydrogen. We interpret these results as an indication that the enhancement in $\chi_\mathrm{ISRF}$ might be affecting the properties of the small stochastically heated dust grains resulting in an increase in the spinning dust emission observed at 33 GHz. This is the first time that such an investigation has been performed, and we believe that this type of analysis creates a new perspective in the field of anomalous emission studies, and represents a powerful new tool for constraining spinning dust models.Comment: 13 pages, 9 figures, accepted for publication in MNRA

Radiative transfer on hierarchial grids

We present new methods for radiative transfer on hierarchial grids. We develop a new method for calculating the scattered flux that employs the grid structure to speed up the computation. We describe a novel subiteration algorithm that can be used to accelerate calculations with strong dust temperature self-coupling. We compute two test models, a molecular cloud and a circumstellar disc, and compare the accuracy and speed of the new algorithms against existing methods. An adaptive model of the molecular cloud with less than 8 % of the cells in the uniform grid produced results in good agreement with the full resolution model. The relative RMS error of the surface brightness <4 % at all wavelengths, and in regions of high column density the relative RMS error was only 10^{-4}. Computation with the adaptive model was faster by a factor of ~5. The new method for calculating the scattered flux is faster by a factor of ~4 in large models with a deep hierarchy structure, when images of the scattered light are computed towards several observing directions. The efficiency of the subiteration algorithm is highly dependent on the details of the model. In the circumstellar disc test the speed-up was a factor of two, but much larger gains are possible. The algorithm is expected to be most beneficial in models where a large number of small, dense regions are embedded in an environment with a lower mean density.Comment: Accepted to A&A; 13 pages, 8 figures; (v2: minor typos corrected

Modelling the spinning dust emission from dense interstellar clouds

Electric dipole emission arising from PAHs is often invoked to explain the anomalous microwave emission (AME). This assignation is based on an observed tight correlation between the mid-IR emission of PAHs and the AME; and a good agreement between models of spinning dust and the broadband AME spectrum. So far often detected at large scale in the diffuse interstellar medium, the AME has recently been studied in detail in well-known dense molecular clouds with the help of Planck data. While much attention has been given to the physics of spinning dust emission, the impact of varying local physical conditions has not yet been considered in detail. Our aim is to study the emerging spinning dust emission from interstellar clouds with realistic physical conditions and radiative transfer. We use the DustEM code from Compiegne et al. to describe the extinction and IR emission of all dust populations. The spinning dust emission is obtained with SpDust, as described by Silsbee et al., that we have coupled to DustEM. We carry out full radiative transfer simulations and carefully estimate the local gas state as a function of position within interstellar clouds. We show that the spinning dust emission is sensitive to the abundances of the major ions and we propose a simple scheme to estimate these abundances. We also investigate the effect of changing the cosmic-ray rate. In dense media, where radiative transfer is mandatory, we show that the relationship between the spinning and mid-IR emissivities of PAHs is no longer linear and that the spinning dust emission may actually be strong at the centre of clouds where the mid-IR PAH emission is weak. These results provide new ways to trace grain growth from diffuse to dense medium and will be useful for the analysis of AME at the scale of interstellar clouds.Comment: 7 pages, 10 figures, accepted by A&

Evolution of dust in the Orion Bar with Herschel: I. Radiative transfer modelling

Interstellar dust is a key element in our understanding of the interstellar medium and star formation. The manner in which dust populations evolve with the excitation and the physical conditions is a first step in the comprehension of the evolution of inter- stellar dust. Within the framework of the Evolution of interstellar dust Herschel key program, we have acquired PACS and SPIRE spec- trophotometric observations of various photodissociation regions, to characterise this evolution. The aim of this paper is to trace the evolution of dust grains in the Orion Bar photodissociation region. We use Herschel/PACS (70 and 160 mic) and SPIRE (250, 350 and 500 mic) together with Spitzer/IRAC observations to map the spatial distribution of the dust populations across the Bar. Brightness profiles are modelled using the DustEM model coupled with a radiative transfer code. Thanks to Herschel, we are able to probe finely the dust emission of the densest parts of the Orion Bar with a resolution from 5.6" to 35.1". These new observations allow us to infer the temperature of the biggest grains at different positions in the Bar, which reveals a gradient from \sim 80 K to 40 K coupled with an increase of the spectral emissivity index from the ionization front to the densest regions. Combining Spitzer/IRAC observations, which are sensitive to the dust emission from the surface, with Herschel maps, we have been able to measure the Orion Bar emission from 3.6 to 500 mic. We find a stratification in the different dust components which can be re- produced quantitatively by a simple radiative transfer model without dust evolution. However including dust evolution is needed to explain the brightness in each band. PAH abundance variations, or a combination of PAH abundance variations with an emissivity enhancement of the biggest grains due to coagulation give good results.Comment: 11 pages, 12 figure

A Novel Approach to Constrain the Escape Fraction and Dust Content at High Redshift Using the Cosmic Infrared Background Fractional Anisotropy

The Cosmic Infrared Background (CIB) provides an opportunity to constrain many properties of the high redshift (z>6) stellar population as a whole. This background, specifically, from 1 to 200 microns, will contain any information about the era of reionization and the stars responsible for producing these ionizing photons. In this paper, we look at the fractional anisotropy delta I/I of this high redshift population, which is the ratio of the magnitude of the fluctuations (delta I) and the mean intensity (I). We show that this can be used to constrain the escape fraction of the population as a whole. The magnitude of the fluctuations of the CIB depend on the escape fraction, while the mean intensity does not. This results in lower values of the escape fraction producing higher values of the fractional anisotropy. This difference is predicted to be larger at the longer wavelengths bands (above 10 microns), albeit it is also much harder to observe in that range. We show that the fractional anisotropy can also be used to separate a dusty from a dust-free population. Finally, we discuss the constraints provided by current observations on the CIB fractional anisotropy.Comment: 8 pages, 4 figures, accepted to ApJ, some clarifications added, matches accepted versio

Dust processing in photodissociation regions - Mid-IR emission modelling

Mid-infrared spectroscopy of dense illuminated ridges (or photodissociation regions, PDRs) suggests dust evolution. Such evolution must be reflected in the gas physical properties through processes like photo-electric heating or H_2 formation. With Spitzer Infrared Spectrograph (IRS) and ISOCAM data, we study the mid-IR emission of closeby, well known PDRs. Focusing on the band and continuum dust emissions, we follow their relative contributions and analyze their variations in terms of abundance of dust populations. In order to disentangle dust evolution and excitation effects, we use a dust emission model that we couple to radiative transfer. Our dust model reproduces extinction and emission of the standard interstellar medium that we represent with diffuse high galactic latitude clouds called Cirrus. We take the properties of dust in Cirrus as a reference to which we compare the dust emission from more excited regions, namely the Horsehead and the reflection nebula NGC 2023 North. We show that in both regions, radiative transfer effects cannot account for the observed spectral variations. We interpret these variations in term of changes of the relative abundance between polycyclic aromatic hydrocarbons (PAHs, mid-IR band carriers) and very small grains (VSGs, mid-IR continuum carriers). We conclude that the PAH/VSG abundance ratio is 2.4 times smaller at the peak emission of the Horsehead nebula than in the Cirrus case. For NGC2023 North where spectral evolution is observed across the northern PDR, we conclude that this ratio is ~5 times lower in the dense, cold zones of the PDR than in its diffuse illuminated part where dust properties seem to be the same as in Cirrus. We conclude that dust in PDRs seems to evolve from "dense" to "diffuse" properties at the small spatial scale of the dense illuminated ridge.Comment: 11 pages, 11 figures, accepted for publication in A&

90GHz and 150GHz observations of the Orion M42 region. A sub-millimeter to radio analysis

We have used the new 90GHz MUSTANG camera on the Robert C. Byrd Green Bank Telescope (GBT) to map the bright Huygens region of the star-forming region M42 with a resolution of 9" and a sensitivity of 2.8mJy/beam. 90GHz is an interesting transition frequency, as MUSTANG detects both the free-free emission characteristic of the HII region created by the Trapezium stars, normally seen at lower frequencies, and thermal dust emission from the background OMC1 molecular cloud, normally mapped at higher frequencies. We also present similar data from the 150GHz GISMO camera taken on the IRAM telescope. This map has 15" resolution. By combining the MUSTANG data with 1.4, 8, and 21GHz radio data from the VLA and GBT, we derive a new estimate of the emission measure (EM) averaged electron temperature of Te = 11376K by an original method relating free-free emission intensities at optically thin and optically thick frequencies. Combining ISO-LWS data with our data, we derive a new estimate of the dust temperature and spectral emissivity index within the 80" ISO-LWS beam toward OrionKL/BN, Td = 42K and Beta=1.3. We show that both Td and Beta decrease when going from the HII region and excited OMC1 interface to the denser UV shielded part of OMC1 (OrionKL/BN, Orion S). With a model consisting of only free-free and thermal dust emission we are able to fit data taken at frequencies from 1.5GHz to 854GHz.Comment: 18 pages, 8 figures, submitted to the Astrophysical Journa

UV-driven chemistry in simulations of the interstellar medium. I. Post-processed chemistry with the Meudon PDR code

Our main purpose is to estimate the effect of assuming uniform density on the line-of-sight in PDR chemistry models, compared to a more realistic distribution for which total gas densities may well vary by several orders of magnitude. A secondary goal of this paper is to estimate the amount of molecular hydrogen which is not properly traced by the CO (J = 1 -> 0) line, the so-called "dark molecular gas". We use results from a magnetohydrodynamical (MHD) simulation as a model for the density structures found in a turbulent diffuse ISM with no star-formation activity. The Meudon PDR code is then applied to a number of lines of sight through this model, to derive their chemical structures. It is found that, compared to the uniform density assumption, maximal chemical abundances for H2, CO, CH and CN are increased by a factor 2 to 4 when taking into account density fluctuations on the line of sight. The correlations between column densities of CO, CH and CN with respect to those of H2 are also found to be in better overall agreement with observations. For instance, at N(H2) > 2.10^{20} cm-2, while observations suggest that d[log N(CO)]=d[log N(H2)] = 3.07 +/- 0.73, we find d[log N(CO)]=d[log N(H2)] =14 when assuming uniform density, and d[log N(CO)]=d[log N(H2)] = 5.2 when including density fluctuations.Comment: 14 pages, 16 figures, accepted for publication in Astronomy & Astrophysic

Physical structure of the photodissociation regions in NGC 7023: Observations of gas and dust emission with <i>Herschel</i>

The determination of the physical conditions in molecular clouds is a key step towards our understanding of their formation and evolution of associated star formation. We investigate the density, temperature, and column density of both dust and gas in the photodissociation regions (PDRs) located at the interface between the atomic and cold molecular gas of the NGC 7023 reflection nebula. We study how young stars affect the gas and dust in their environment. Our approach combining both dust and gas delivers strong constraints on the physical conditions of the PDRs. We find dense and warm molecular gas of high column density in the PDRs