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

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

Satellite Radiation Products for Ocean Biology and Biogeochemistry: Needs, State-of-the-Art, Gaps, Development Priorities, and Opportunities

Knowing the spatial and temporal distribution of the underwater light field, i.e., the spectral and angular structure of the radiant intensity at any point in the water column, is essential to understanding the biogeochemical processes that control the composition and evolution of aquatic ecosystems and their impact on climate and reaction to climate change. At present, only a few properties are reliably retrieved from space, either directly or via water-leaving radiance. Existing satellite products are limited to planar photosynthetically available radiation (PAR) and ultraviolet (UV) irradiance above the surface and diffuse attenuation coefficient. Examples of operational products are provided, and their advantages and drawbacks are examined. The usefulness and convenience of these products notwithstanding, there is a need, as expressed by the user community, for other products, i.e., sub-surface planar and scalar fluxes, average cosine, spectral fluxes (UV to visible), diurnal fluxes, absorbed fraction of PAR by live algae (APAR), surface albedo, vertical attenuation, and heating rate, and for associating uncertainties to any product on a pixel-by-pixel basis. Methodologies to obtain the new products are qualitatively discussed in view of most recent scientific knowledge and current and future satellite missions, and specific algorithms are presented for some new products, namely sub-surface fluxes and average cosine. A strategy and roadmap (short, medium, and long term) for usage and development priorities is provided, taking into account needs and readiness level. Combining observations from satellites overpassing at different times and geostationary satellites should be pursued to improve the quality of daily-integrated radiation fields, and products should be generated without gaps to provide boundary conditions for general circulation and biogeochemical models. Examples of new products, i.e., daily scalar PAR below the surface, daily average cosine for PAR, and sub-surface spectral scalar fluxes are presented. A procedure to estimate algorithm uncertainties in the total uncertainty budget for above-surface daily PAR, based on radiative simulations for expected situations, is described. In the future, space-borne lidars with ocean profiling capability offer the best hope for improving our knowledge of sub-surface fields. To maximize temporal coverage, space agencies should consider placing ocean-color instruments in L1 orbit, where the sunlit part of the Earth can be frequently observed

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&

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&

The global dust SED: Tracing the nature and evolution of dust with DustEM

The Planck and Herschel missions are currently measuring the farIR-mm emission of dust, which combined with existing IR data, will for the first time provide the full SED of the galactic ISM dust emission with an unprecedented sensitivity and angular resolution. It will allow a systematic study of the dust evolution processes that affect the SED. Here we present a versatile numerical tool, DustEM, that predicts the emission and extinction of dust given their size distribution and their optical and thermal properties. In order to model dust evolution, DustEM has been designed to deal with a variety of grain types, structures and size distributions and to be able to easily include new dust physics. We use DustEM to model the dust SED and extinction in the diffuse interstellar medium at high-galactic latitude (DHGL), a natural reference SED. We present a coherent set of observations for the DHGL SED. The dust components in our DHGL model are (i) PAHs, (ii) amorphous carbon and (iii) amorphous silicates. We use amorphous carbon dust, rather than graphite, because it better explains the observed high abundances of gas-phase carbon in shocked regions of the interstellar medium. Using the DustEM model, we illustrate how, in the optically thin limit, the IRAS/Planck HFI (and likewise Spitzer/Herschel for smaller spatial scales) photometric band ratios of the dust SED can disentangle the influence of the exciting radiation field intensity and constrain the abundance of small grains relative to the larger grains. We also discuss the contributions of the different grain populations to the IRAS, Planck and Herschel channels. Such information is required to enable a study of the evolution of dust as well as to systematically extract the dust thermal emission from CMB data and to analyze the emission in the Planck polarized channels. The DustEM code described in this paper is publically available.Comment: accepted for publication in A&

Excitation of H2_2 in photodissociation regions as seen by Spitzer

We present spectroscopic observations obtained with the infrared Spitzer Space Telescope, which provide insight into the H$_2$ physics and gas energetics in photodissociation Regions (PDRs) of low to moderate far-ultraviolet (FUV) fields and densities. We analyze data on six well known Galactic PDRs (L1721, California, N7023E, Horsehead, rho Oph, N2023N), sampling a poorly explored range of excitation conditions ($\chi \sim 5-10^3$), relevant to the bulk of molecular clouds in galaxies. Spitzer observations of H$_2$ rotational lines are complemented with H$_2$ data, including ro-vibrational line measurements, obtained with ground-based telescopes and ISO, to constrain the relative contributions of ultraviolet pumping and collisions to the H$_2$ excitation. The data analysis is supported by model calculations with the Meudon PDR code. The observed column densities of rotationally excited H$_2$ are observed to be much higher than PDR model predictions. In the lowest excitation PDRs, the discrepancy between the model and the data is about one order of magnitude for rotational levels $J \ge$3. We discuss whether an enhancement in the H$_2$ formation rate or a local increase in photoelectric heating, as proposed for brighter PDRs in former ISO studies, may improve the data-model comparison. We find that an enhancement in the H$_2$ formation rates reduces the discrepancy, but the models still fall short of the data. This large disagreement suggests that our understanding of the formation and excitation of H$_2$ and/or of PDRs energetics is still incomplete. We discuss several explanations, which could be further tested using the Herschel Space TelescopeComment: A&A in pres

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