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.

Excitation of H$_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 Telescope