Search

[Abridged] The ionization of polycyclic aromatic hydrocarbons (PAHs), by ultraviolet (UV) photons from massive stars is expected to account for a large fraction of the heating of neutral gas in galaxies. Evaluation of this proposal, however, has been limited by our ability to directly compare observational diagnostics to the results of a molecular model describing PAH ionization. The objective of this article is to take advantage of the most recent values of molecular parameters derived from laboratory experiments and quantum chemical calculations on PAHs and provide a detailed comparison between modeled values and observational diagnostics for the PAH charge state and the heating efficiency for PAHs. Despite the use of a simple analytical model, we obtain a good agreement between model results and observational diagnostics over a wide range of radiation fields and physical conditions, in environments such as star-forming regions, galaxies, and protoplanetary disks. In addition, we found that the modeled photoelectric heating rates by PAHs are close to the observed cooling rates given by the gas emission. These results show that PAH ionization is the main source of neutral gas heating in these environments. The results of our photoelectric heating model by PAHs can thus be used to assess the contribution of UV radiative heating in galaxies (vs shocks, for instance). We provide the empirical formulas fitted to the model results, and the full python code itself, to calculate the heating rates and heating efficiencies for PAHs., Comment: Accepted for publication in Astronomy and Astrophysics

Context. The [C II] 158 m fine-structure line is the dominant cooling line of moderate-density photodissociation regions (PDRs) illuminated by moderately bright far-ultraviolet (FUV) radiation fields. This makes this line a prime diagnostic for extended regions illuminated by massive stars. Aims. We aim to understand the origin of [C II] emission and its relation to other tracers of gas and dust in PDRs. One focus is a study of the heating efficiency of interstellar gas as traced by the [C II] line to test models of the photoelectric heating of neutral gas by polycyclic aromatic hydrocarbon (PAH) molecules and very small grains. Methods. We make use of a one-square-degree map of velocity-resolved [C II] line emission toward the Orion Nebula complex, and split this out into the individual spatial components, the expanding Veil Shell, the surface of OMC4, and the PDRs associated with the compact H II region of M43 and the reflection nebula NGC 1977. We employed Herschel far-infrared photometric images to determine dust properties. Moreover, we compared with Spitzer mid-infrared photometry to trace hot dust and large molecules, and velocity-resolved IRAM 30m CO(2-1) observations of the molecular gas. Results. The [C II] intensity is tightly correlated with PAH emission in the IRAC 8 m band and far-infrared emission from warm dust, with small variations between the four studied subregions (Veil Shell, OMC4, M43, and NGC 1977). The correlation between [C II] and CO(2-1) is very different in the four subregions and is very sensitive to the detailed geometry of the respective regions. Constant-density PDR models are able to reproduce the observed [C II], CO(2-1), and integrated far-infrared (FIR) intensities. The physical conditions in the Veil Shell of the Orion Nebula, M43, and NGC 1977 reveal a constant ratio of thermal pressure pth over incident FUV radiation field measured by G0. We observe strong variations in the photoelectric heating efficiency in the Veil Shell be

The James Webb Space Telescope (JWST) will be launched in December 2021, with four instruments to perform imaging and spectroscopy. This paper presents work which is part of the Early Release Science (ERS) program "PDRs4All" aimed at observing the Orion Bar. It focuses on the Near Infrared Camera (NIRCam) imaging which will be performed as part of this project. The aim of this paper is to illustrate a methodology to simulate observations of an extended source that is similar to the Orion Bar with NIRCam, and to run the pipeline on these simulated observations. These simulations provide us with a clear idea of the observations that will be obtained, based on the "Astronomer's proposal tool" settings. The analysis also provides an assessment of the risks of saturation. The methodology presented in this document can be applied for JWST observing programs of extended objects containing bright point sources, e.g. for observations of nebulae or nearby galaxies.

The [CII] 158$\mu$m fine-structure line is the dominant cooling line of moderate-density photodissociation regions (PDRs) illuminated by moderately bright far-ultraviolet (FUV) radiation fields. We aim to understand the origin of [CII] emission and its relation to other tracers of gas and dust in PDRs. One focus is a study of the heating efficiency of interstellar gas as traced by the [CII] line to test models of the photoelectric heating of neutral gas by polycyclic aromatic hydrocarbon (PAH) molecules and very small grains. We make use of a one-square-degree map of velocity-resolved [CII] line emission toward the Orion Nebula complex, and split this out into the individual spatial components, the expanding Veil Shell, the surface of OMC4, and the PDRs associated with the compact HII region of M43 and the reflection nebula NGC 1977. We employed Herschel far-infrared photometric images to determine dust properties. Moreover, we compared with Spitzer mid-infrared photometry to trace hot dust and large molecules, and velocity-resolved IRAM 30m CO(2-1) observations of the molecular gas. The [CII] intensity is tightly correlated with PAH emission in the IRAC 8$\mu$m band and far-infrared emission from warm dust. The correlation between [CII] and CO(2-1) is very different in the four subregions and is very sensitive to the detailed geometry. Constant-density PDR models are able to reproduce the observed [CII], CO(2-1), and integrated far-infrared (FIR) intensities. We observe strong variations in the photoelectric heating efficiency in the Veil Shell behind the Orion Bar and these variations are seemingly not related to the spectral properties of the PAHs. The [CII] emission from the Orion Nebula complex stems mainly from moderately illuminated PDR surfaces. Future observations with the James Webb Space Telescope can shine light on the PAH properties that may be linked to these variations.

The [CII] $158\,\mu\mathrm{m}$ fine-structure line is one of the dominant coolants of the neutral interstellar medium. It is hence one of the brightest far-infrared emission lines and can be observed not only in star-forming regions throughout the Galaxy, but also in the diffuse interstellar medium and in distant galaxies. [CII] line emission has been suggested to be a powerful tracer of star-formation. We aim to understand the origin of [CII] emission and its relation to other tracers of interstellar gas and dust. This includes a study of the heating efficiency of interstellar gas as traced by the [CII] line to test models of gas heating. We make use of a one-square-degree map of velocity-resolved [CII] line emission towards the Orion Nebula complex, including M43 and NGC 1977. The [CII] intensity is tightly correlated with PAH emission in the IRAC $8\,\mu\mathrm{m}$ band and far-infrared emission from warm dust. The correlation between [CII] and CO(2-1) is affected by the detailed geometry of the region. We find particularly low [CII]-over-FIR intensity ratios towards large columns of (warm and cold) dust, which suggest the interpretation of the "[CII] deficit" in terms of a "FIR excess". A slight decrease in the FIR line-over-continuum intensity ratio can be attributed to a decreased heating efficiency of the gas. We find that, at the mapped spatial scales, predictions of the star-formation rate from [CII] emission underestimate the star-formation rate calculated from YSO counts in the Orion Nebula complex by an order of magnitude. [CII] emission from the Orion Nebula complex arises dominantly in the cloud surfaces, many viewed in edge-on geometry. [CII] emission from extended faint cloud surfaces may contribute significantly to the total [CII] emission on galactic scales.

The James Webb Space Telescope (JWST) will be launched in December 2021, with four instruments to perform imaging and spectroscopy. This paper presents work which is part of the Early Release Science (ERS) program "PDRs4All" aimed at observing the Orion Bar. It focuses on the Near Infrared Camera (NIRCam) imaging which will be performed as part of this project. The aim of this paper is to illustrate a methodology to simulate observations of an extended source that is similar to the Orion Bar with NIRCam, and to run the pipeline on these simulated observations. These simulations provide us with a clear idea of the observations that will be obtained, based on the "Astronomer's proposal tool" settings. The analysis also provides an assessment of the risks of saturation. The methodology presented in this document can be applied for JWST observing programs of extended objects containing bright point sources, e.g. for observations of nebulae or nearby galaxies.

The [CII] 158$\mu$m fine-structure line is the dominant cooling line of moderate-density photodissociation regions (PDRs) illuminated by moderately bright far-ultraviolet (FUV) radiation fields. We aim to understand the origin of [CII] emission and its relation to other tracers of gas and dust in PDRs. One focus is a study of the heating efficiency of interstellar gas as traced by the [CII] line to test models of the photoelectric heating of neutral gas by polycyclic aromatic hydrocarbon (PAH) molecules and very small grains. We make use of a one-square-degree map of velocity-resolved [CII] line emission toward the Orion Nebula complex, and split this out into the individual spatial components, the expanding Veil Shell, the surface of OMC4, and the PDRs associated with the compact HII region of M43 and the reflection nebula NGC 1977. We employed Herschel far-infrared photometric images to determine dust properties. Moreover, we compared with Spitzer mid-infrared photometry to trace hot dust and large molecules, and velocity-resolved IRAM 30m CO(2-1) observations of the molecular gas. The [CII] intensity is tightly correlated with PAH emission in the IRAC 8$\mu$m band and far-infrared emission from warm dust. The correlation between [CII] and CO(2-1) is very different in the four subregions and is very sensitive to the detailed geometry. Constant-density PDR models are able to reproduce the observed [CII], CO(2-1), and integrated far-infrared (FIR) intensities. We observe strong variations in the photoelectric heating efficiency in the Veil Shell behind the Orion Bar and these variations are seemingly not related to the spectral properties of the PAHs. The [CII] emission from the Orion Nebula complex stems mainly from moderately illuminated PDR surfaces. Future observations with the James Webb Space Telescope can shine light on the PAH properties that may be linked to these variations.

The [CII] $158\,\mu\mathrm{m}$ fine-structure line is one of the dominant coolants of the neutral interstellar medium. It is hence one of the brightest far-infrared emission lines and can be observed not only in star-forming regions throughout the Galaxy, but also in the diffuse interstellar medium and in distant galaxies. [CII] line emission has been suggested to be a powerful tracer of star-formation. We aim to understand the origin of [CII] emission and its relation to other tracers of interstellar gas and dust. This includes a study of the heating efficiency of interstellar gas as traced by the [CII] line to test models of gas heating. We make use of a one-square-degree map of velocity-resolved [CII] line emission towards the Orion Nebula complex, including M43 and NGC 1977. The [CII] intensity is tightly correlated with PAH emission in the IRAC $8\,\mu\mathrm{m}$ band and far-infrared emission from warm dust. The correlation between [CII] and CO(2-1) is affected by the detailed geometry of the region. We find particularly low [CII]-over-FIR intensity ratios towards large columns of (warm and cold) dust, which suggest the interpretation of the "[CII] deficit" in terms of a "FIR excess". A slight decrease in the FIR line-over-continuum intensity ratio can be attributed to a decreased heating efficiency of the gas. We find that, at the mapped spatial scales, predictions of the star-formation rate from [CII] emission underestimate the star-formation rate calculated from YSO counts in the Orion Nebula complex by an order of magnitude. [CII] emission from the Orion Nebula complex arises dominantly in the cloud surfaces, many viewed in edge-on geometry. [CII] emission from extended faint cloud surfaces may contribute significantly to the total [CII] emission on galactic scales.

The Orion Molecular Cloud is the nearest massive-star forming region. Massive stars have profound effects on their environment due to their strong radiation fields and stellar winds. Velocity-resolved observations of the [CII] $158\,\mu\mathrm{m}$ fine-structure line allow us to study the kinematics of UV-illuminated gas. Here, we present a square-degree-sized map of [CII] emission from the Orion Nebula complex obtained by the upGREAT instrument onboard SOFIA, covering the entire Orion Nebula (M42) plus M43 and the nebulae NGC 1973, 1975, and 1977. We compare the stellar characteristics of these three regions with the kinematics of the expanding bubbles surrounding them. The bubble blown by the O7V star $\theta^1$ Ori C in the Orion Nebula expands rapidly, at $13\,\mathrm{km\,s^{-1}}$. Simple analytical models reproduce the characteristics of the hot interior gas and the neutral shell of this wind-blown bubble and give us an estimate of the expansion time of $0.2\,\mathrm{Myr}$. M43 with the B0.5V star NU Ori also exhibits an expanding bubble structure, with an expansion velocity of $6\,\mathrm{km\,s^{-1}}$. Comparison with analytical models for the pressure-driven expansion of H\,{\sc ii} regions gives an age estimate of $0.02\,\mathrm{Myr}$. The bubble surrounding NGC 1973, 1975, and 1977 with the central B1V star 42 Orionis expands at $1.5\,\mathrm{km\,s^{-1}}$, likely due to the over-pressurized ionized gas as in the case of M43. We derive an age of $0.4\,\mathrm{Myr}$ for this structure. We conclude that the bubble of the Orion Nebula is driven by the mechanical energy input by the strong stellar wind from $\theta^1$ Ori C, while the bubbles associated with M43 and NGC 1977 are caused by the thermal expansion of the gas ionized by their central later-type massive stars.

The Orion Molecular Cloud is the nearest massive-star forming region. Massive stars have profound effects on their environment due to their strong radiation fields and stellar winds. Velocity-resolved observations of the [CII] $158\,\mu\mathrm{m}$ fine-structure line allow us to study the kinematics of UV-illuminated gas. Here, we present a square-degree-sized map of [CII] emission from the Orion Nebula complex obtained by the upGREAT instrument onboard SOFIA, covering the entire Orion Nebula (M42) plus M43 and the nebulae NGC 1973, 1975, and 1977. We compare the stellar characteristics of these three regions with the kinematics of the expanding bubbles surrounding them. The bubble blown by the O7V star $\theta^1$ Ori C in the Orion Nebula expands rapidly, at $13\,\mathrm{km\,s^{-1}}$. Simple analytical models reproduce the characteristics of the hot interior gas and the neutral shell of this wind-blown bubble and give us an estimate of the expansion time of $0.2\,\mathrm{Myr}$. M43 with the B0.5V star NU Ori also exhibits an expanding bubble structure, with an expansion velocity of $6\,\mathrm{km\,s^{-1}}$. Comparison with analytical models for the pressure-driven expansion of H\,{\sc ii} regions gives an age estimate of $0.02\,\mathrm{Myr}$. The bubble surrounding NGC 1973, 1975, and 1977 with the central B1V star 42 Orionis expands at $1.5\,\mathrm{km\,s^{-1}}$, likely due to the over-pressurized ionized gas as in the case of M43. We derive an age of $0.4\,\mathrm{Myr}$ for this structure. We conclude that the bubble of the Orion Nebula is driven by the mechanical energy input by the strong stellar wind from $\theta^1$ Ori C, while the bubbles associated with M43 and NGC 1977 are caused by the thermal expansion of the gas ionized by their central later-type massive stars.

Observations towards L1630 in the Orion B molecular cloud, comprising the iconic Horsehead Nebula, allow us to study the interplay between stellar radiation and a molecular cloud under relatively benign conditions, that is, intermediate densities and an intermediate UV radiation field. Contrary to the well-studied Orion Molecular Cloud 1 (OMC1), which hosts much harsher conditions, L1630 has little star formation. We aim to relate the [CII] fine-structure line emission to the physical conditions predominant in L1630 and compare it to studies of OMC1. The [CII] $158\,\mu\mathrm{m}$ emission from an area of $12' \times 17'$ in L1630 was observed using the upGREAT instrument onboard SOFIA. Of the [CII] emission from the mapped area 95%, $13\,L_{\odot}$, originates from the molecular cloud; the adjacent HII region contributes only 5%, that is, $1\,L_{\odot}$. From comparison with other data (CO (1-0)-line emission, far-infrared (FIR) continuum studies, emission from polycyclic aromatic hydrocarbons (PAHs)), we infer a gas density of the molecular cloud of $n_{\mathrm{H}}\sim 3\cdot 10^3\,\mathrm{cm^{-3}}$, with surface layers, including the Horsehead Nebula, having a density of up to $n_{\mathrm{H}}\sim 4\cdot 10^4\,\mathrm{cm^{-3}}$. The temperature of the surface gas is $T\sim 100\,\mathrm{K}$. The average [CII] cooling efficiency within the molecular cloud is $1.3\cdot 10^{-2}$. The fraction of the mass of the molecular cloud within the studied area that is traced by [CII] is only $8\%$. Our PDR models are able to reproduce the FIR-[CII] correlations and also the CO (1-0)-[CII] correlations. Finally, we compare our results on the heating efficiency of the gas with theoretical studies of photoelectric heating by PAHs, clusters of PAHs, and very small grains, and find the heating efficiency to be lower than theoretically predicted, a continuation of the trend set by other observations., Comment: 23 pages, 30 figures, accepted for publication in Astronomy & Astrophysics

Emission of fullerenes in their infrared vibrational bands has been detected in space near hot stars. The proposed attribution of the diffuse interstellar bands at 9577 and 9632 \AA\ to electronic transitions of the buckminsterfullerene cation (i.e. C$_{60}^+$ ) was recently supported by new laboratory data, confirming the presence of this species in the diffuse interstellar medium (ISM). In this letter, we present the detection, also in the diffuse ISM, of the 17.4 and 18.9 $\mu$m emission bands commonly attributed to vibrational bands of neutral C$_{60}$. According to classical models that compute the charge state of large molecules in space, C$_{60}$ is expected to be mostly neutral in the diffuse ISM. This is in agreement with the abundances of diffuse C$_{60}$ we derive here from observations., Comment: Accepted for publication in A&A. v3. including language editing and correction of values for abundances

Observations towards L1630 in the Orion B molecular cloud, comprising the iconic Horsehead Nebula, allow us to study the interplay between stellar radiation and a molecular cloud under relatively benign conditions, that is, intermediate densities and an intermediate UV radiation field. Contrary to the well-studied Orion Molecular Cloud 1 (OMC1), which hosts much harsher conditions, L1630 has little star formation. We aim to relate the [CII] fine-structure line emission to the physical conditions predominant in L1630 and compare it to studies of OMC1. The [CII] $158\,\mu\mathrm{m}$ emission from an area of $12' \times 17'$ in L1630 was observed using the upGREAT instrument onboard SOFIA. Of the [CII] emission from the mapped area 95%, $13\,L_{\odot}$, originates from the molecular cloud; the adjacent HII region contributes only 5%, that is, $1\,L_{\odot}$. From comparison with other data (CO (1-0)-line emission, far-infrared (FIR) continuum studies, emission from polycyclic aromatic hydrocarbons (PAHs)), we infer a gas density of the molecular cloud of $n_{\mathrm{H}}\sim 3\cdot 10^3\,\mathrm{cm^{-3}}$, with surface layers, including the Horsehead Nebula, having a density of up to $n_{\mathrm{H}}\sim 4\cdot 10^4\,\mathrm{cm^{-3}}$. The temperature of the surface gas is $T\sim 100\,\mathrm{K}$. The average [CII] cooling efficiency within the molecular cloud is $1.3\cdot 10^{-2}$. The fraction of the mass of the molecular cloud within the studied area that is traced by [CII] is only $8\%$. Our PDR models are able to reproduce the FIR-[CII] correlations and also the CO (1-0)-[CII] correlations. Finally, we compare our results on the heating efficiency of the gas with theoretical studies of photoelectric heating by PAHs, clusters of PAHs, and very small grains, and find the heating efficiency to be lower than theoretically predicted, a continuation of the trend set by other observations., Comment: 23 pages, 30 figures, accepted for publication in Astronomy & Astrophysics

Emission of fullerenes in their infrared vibrational bands has been detected in space near hot stars. The proposed attribution of the diffuse interstellar bands at 9577 and 9632 \AA\ to electronic transitions of the buckminsterfullerene cation (i.e. C$_{60}^+$ ) was recently supported by new laboratory data, confirming the presence of this species in the diffuse interstellar medium (ISM). In this letter, we present the detection, also in the diffuse ISM, of the 17.4 and 18.9 $\mu$m emission bands commonly attributed to vibrational bands of neutral C$_{60}$. According to classical models that compute the charge state of large molecules in space, C$_{60}$ is expected to be mostly neutral in the diffuse ISM. This is in agreement with the abundances of diffuse C$_{60}$ we derive here from observations., Comment: Accepted for publication in A&A. v3. including language editing and correction of values for abundances

We present mid-infrared photometry of the Orion Bar obtained with FORCAST aboard SOFIA at 6.4, 6.6, 7.7, 19.7, 31.5 and 37.1 \um. By complementing this observations with archival FORCAST and \emph{Herschel}/PACS images we are able to construct a complete infrared spectral energy distribution of the Huygens region in the Orion nebula By comparing the infrared images with gas tracers, we find that PACS maps trace the molecular cloud, while the FORCAST data trace the photodissociation region (PDR) and HII region. Analysis of the energetics of the region reveal that the PDR extends for 0.28~pc along the line-of-sight and that the Bar is inclined at an angle of $4\degr$. The infrared and submillimeter images reveal that the Orion Bar represents a swept up shell with a thickness of 0.1~pc. The mass of the shell implies a shock velocity of $\simeq 3$ km/s and an age of $\simeq 10^5$ yr for the HII region. Our analysis shows that the UV and infrared dust opacities in the HII region and the PDR are a factor 5 to 10 lower than in the diffuse interstellar medium. In the ionized gas, Ly$\alpha$ photons are a major source of dust heating at distances larger than $\simeq0.06$~pc from \toc. Dust temperatures can be explained if the size of the grains is between 0.1 to 1~\um. We derive the photo-electric heating efficiency of the atomic gas in the Orion Bar. The results are in good qualitative agreement with models and The quantitative differences indicate a decreased PAH abundance in this region.

Transitional disks are structures of dust and gas around young stars with large inner cavities in which planet formation may occur. Lopsided dust distributions are observed in the dust continuum emission at millimeter wavelengths. These asymmetrical structures can be explained as the result of an enhanced gas density vortex where the dust is trapped potentially promoting the rapid growth to the planetesimal scale. AB Aur hosts a transitional disk with a clear horseshoe morphology which strongly suggests the presence of a dust trap. Our goal is to investigate its formation and the possible effects on the gas chemistry. We used the NOEMA interferometer to image the 1mm continuum dust emission and the 13CO J=2->1, C18O J=2->1, SO J=56->45 and H2CO J=303->202 rotational lines. Line integrated intensity ratio images are built to investigate the chemical changes within the disk. We have used a single point (n,T) chemical model to investigate the lifetime of gaseous CO, H2CO and SO in the dust trap. Our model shows that for densities >10^7 cm^-3, the SO molecules are depleted (directly frozen or converted into SO2 and then frozen out) in less than 0.1 Myr. The lower SO abundance towards the dust trap could indicate that a larger fraction of the gas is in a high density environment. Gas dynamics, grain growth and gas chemistry are coupled in the planet formation process. Because of the strong dependence of SO abundance on the gas density, the sulfur chemistry can be used as a chemical diagnostic to detect the birthsites of future planets. However, the large uncertainties inherent to chemical models and the limited knowledge of the disk physical structure and initial conditions are important drawbacks.

NGC 7023 is a well-studied reflection nebula, which shows strong emission from polycyclic aromatic hydrocarbon (PAH) molecules in the form of aromatic infrared bands (AIBs). The spectral variations of the AIBs in this region are connected to the chemical evolution of the PAH molecules which, in turn, depends on the local physical conditions. We use the capabilities of SOFIA to observe a 3.2' x 3.4' region of NGC 7023 at wavelengths that we observe with high spatial resolution (2.7") at 3.3 and 11.2 um. We compare the SOFIA images with existing images of the PAH emission at 8.0 um (Spitzer), emission from evaporating very small grains (eVSG) extracted from Spitzer-IRS spectral cubes, the ERE (HST and CFHT), and H_2 (2.12 um). We create maps of the 11.2/3.3 um ratio to probe the morphology of the PAH size distribution and the 8.0/11.2 um ratio to probe the PAH ionization. We make use of an emission model and of vibrational spectra from the NASA Ames PAHdb to translate the 11.2/3.3 um ratio to PAH sizes. The 11.2/3.3 um map shows the smallest PAH concentrate on the PDR surface (H_2 and extended red emission) in the NW and South PDR. We estimated that PAHs in the NW PDR bear, on average, a number of carbon atoms (N_c) of ~70 in the PDR cavity and ~50 at the PDR surface. In the entire nebula, the results reveal a factor of 2 variation in the size of the PAH. We relate these size variations to several models for the evolution of the PAH families when they traverse from the molecular cloud to the PDR. The PAH size map enables us to follow the photochemical evolution of PAHs in NGC 7023. Small PAHs result from the photo-evaporation of VSGs as they reach the PDR surface. Inside the PDR cavity, the PAH abundance drops as the smallest PAH are broken down. The average PAH size increases in the cavity where only the largest species survive or are converted into C_60 by photochemical processing., Comment: accepted for publication in A&A

We present mid-infrared photometry of the Orion Bar obtained with FORCAST aboard SOFIA at 6.4, 6.6, 7.7, 19.7, 31.5 and 37.1 \um. By complementing this observations with archival FORCAST and \emph{Herschel}/PACS images we are able to construct a complete infrared spectral energy distribution of the Huygens region in the Orion nebula By comparing the infrared images with gas tracers, we find that PACS maps trace the molecular cloud, while the FORCAST data trace the photodissociation region (PDR) and HII region. Analysis of the energetics of the region reveal that the PDR extends for 0.28~pc along the line-of-sight and that the Bar is inclined at an angle of $4\degr$. The infrared and submillimeter images reveal that the Orion Bar represents a swept up shell with a thickness of 0.1~pc. The mass of the shell implies a shock velocity of $\simeq 3$ km/s and an age of $\simeq 10^5$ yr for the HII region. Our analysis shows that the UV and infrared dust opacities in the HII region and the PDR are a factor 5 to 10 lower than in the diffuse interstellar medium. In the ionized gas, Ly$\alpha$ photons are a major source of dust heating at distances larger than $\simeq0.06$~pc from \toc. Dust temperatures can be explained if the size of the grains is between 0.1 to 1~\um. We derive the photo-electric heating efficiency of the atomic gas in the Orion Bar. The results are in good qualitative agreement with models and The quantitative differences indicate a decreased PAH abundance in this region.

NGC 7023 is a well-studied reflection nebula, which shows strong emission from polycyclic aromatic hydrocarbon (PAH) molecules in the form of aromatic infrared bands (AIBs). The spectral variations of the AIBs in this region are connected to the chemical evolution of the PAH molecules which, in turn, depends on the local physical conditions. We use the capabilities of SOFIA to observe a 3.2' x 3.4' region of NGC 7023 at wavelengths that we observe with high spatial resolution (2.7") at 3.3 and 11.2 um. We compare the SOFIA images with existing images of the PAH emission at 8.0 um (Spitzer), emission from evaporating very small grains (eVSG) extracted from Spitzer-IRS spectral cubes, the ERE (HST and CFHT), and H_2 (2.12 um). We create maps of the 11.2/3.3 um ratio to probe the morphology of the PAH size distribution and the 8.0/11.2 um ratio to probe the PAH ionization. We make use of an emission model and of vibrational spectra from the NASA Ames PAHdb to translate the 11.2/3.3 um ratio to PAH sizes. The 11.2/3.3 um map shows the smallest PAH concentrate on the PDR surface (H_2 and extended red emission) in the NW and South PDR. We estimated that PAHs in the NW PDR bear, on average, a number of carbon atoms (N_c) of ~70 in the PDR cavity and ~50 at the PDR surface. In the entire nebula, the results reveal a factor of 2 variation in the size of the PAH. We relate these size variations to several models for the evolution of the PAH families when they traverse from the molecular cloud to the PDR. The PAH size map enables us to follow the photochemical evolution of PAHs in NGC 7023. Small PAHs result from the photo-evaporation of VSGs as they reach the PDR surface. Inside the PDR cavity, the PAH abundance drops as the smallest PAH are broken down. The average PAH size increases in the cavity where only the largest species survive or are converted into C_60 by photochemical processing., Comment: accepted for publication in A&A

Transitional disks are structures of dust and gas around young stars with large inner cavities in which planet formation may occur. Lopsided dust distributions are observed in the dust continuum emission at millimeter wavelengths. These asymmetrical structures can be explained as the result of an enhanced gas density vortex where the dust is trapped potentially promoting the rapid growth to the planetesimal scale. AB Aur hosts a transitional disk with a clear horseshoe morphology which strongly suggests the presence of a dust trap. Our goal is to investigate its formation and the possible effects on the gas chemistry. We used the NOEMA interferometer to image the 1mm continuum dust emission and the 13CO J=2->1, C18O J=2->1, SO J=56->45 and H2CO J=303->202 rotational lines. Line integrated intensity ratio images are built to investigate the chemical changes within the disk. We have used a single point (n,T) chemical model to investigate the lifetime of gaseous CO, H2CO and SO in the dust trap. Our model shows that for densities >10^7 cm^-3, the SO molecules are depleted (directly frozen or converted into SO2 and then frozen out) in less than 0.1 Myr. The lower SO abundance towards the dust trap could indicate that a larger fraction of the gas is in a high density environment. Gas dynamics, grain growth and gas chemistry are coupled in the planet formation process. Because of the strong dependence of SO abundance on the gas density, the sulfur chemistry can be used as a chemical diagnostic to detect the birthsites of future planets. However, the large uncertainties inherent to chemical models and the limited knowledge of the disk physical structure and initial conditions are important drawbacks.

[Context] Transitional disks are structures of dust and gas around young stars with large inner cavities in which planet formation may occur. Lopsided dust distributions are observed in the dust continuum emission at millimeter wavelengths. These asymmetrical structures can be explained as being the result of an enhanced gas density vortex where the dust is trapped, potentially promoting the rapid growth to the planetesimal scale., [Aims] AB Aur hosts a transitional disk with a clear horseshoe morphology which strongly suggests the presence of a dust trap. Our goal is to investigate its formation and the possible effects on the gas chemistry., [Methods] We used the NOrthern Extended Millimeter Array (NOEMA) interferometer to image the 1 mm continuum dust emission and the 13CO J = 2 → 1, C18O J = 2 → 1, SO J = 56 → 45, and H2CO J = 303 → 202 rotational lines., [Results] Line integrated intensity ratio images are built to investigate the chemical changes within the disk. The I(H2CO J = 303 → 202)/I(C18O J = 2 → 1) ratio is fairly constant along the disk with values of ~0.15 ± 0.05. On the contrary, the I(SO J = 56 → 45)/I(C18O J = 2 → 1) and I(SO J = 56 → 45)/I(H2CO J = 303 → 202) ratios present a clear northeast-southwest gradient (a factor of 3−6) with the minimum towards the dust trap. This gradient cannot be explained by a local change in the excitation conditions but by a decrease in the SO abundance. Gas densities up to ~109 cm-3 are expected in the disk midplane and two-three times larger in the high pressure vortex. We have used a single point (n,T) chemical model to investigate the lifetime of gaseous CO, H2CO, and SO in the dust trap. Our model shows that for densities >107 cm-3, the SO molecules are depleted (directly frozen, or converted into SO2 and then frozen out) in less than 0.1 Myr. The lower SO abundance towards the dust trap could indicate that a larger fraction of the gas is in a high density environment., [Conclusions] Gas dynamics, grain growth and gas chemistry are coupled in the planet formation process. We detect a chemical signature of the presence of a dust trap in a transitional disk. Because of the strong dependence of SO abundance on the gas density, the sulfur chemistry can be used as a chemical diagnostic to detect the birthsites of future planets. However, the large uncertainties inherent to chemical models and the limited knowledge of the disk’s physical structure and initial conditions are important drawbacks.

The nearby cloud L1642 is one of only two known very high latitude (|b| > 30 deg) clouds actively forming stars. It is a rare example of star formation in isolated conditions, and can reveal important details of star formation in general, e.g., of the effect of magnetic fields. We compare Herschel dust emission structures and magnetic field orientation revealed by Planck polarization maps in L1642. The high-resolution ($\sim20"$) Herschel data reveal a complex structure including a dense, compressed central clump, and low density striations. The Planck polarization data (at 10$'$ resolution) reveal an ordered magnetic field pervading the cloud and aligned with the surrounding striations. There is a complex interplay between the cloud structure and large scale magnetic field. This suggests that the magnetic field is closely linked to the formation and evolution of the cloud. CO rotational emission confirms that the striations are connected with the main clumps and likely to contain material either falling into or flowing out of the clumps. There is a clear transition from aligned to perpendicular structures approximately at a column density of $N_{\rm{H}} = 1.6 \times 10^{21}\, {{\rm cm}}^{-2}$. Comparing the Herschel maps with the Planck polarization maps shows the close connection between the magnetic field and cloud structure even in the finest details of the cloud., Comment: accepted to MNRAS; 14 pages; changes: moderate text and figure editing

Aims. Our goal is to determine the molecular composition of the circumstellar disk around AB Aurigae (hereafter, AB Aur). AB Aur is a prototypical Herbig Ae star and the understanding of its disk chemistry is of paramount importance to understand the chemical evolution of the gas in warm disks. Methods. We used the IRAM 30-m telescope to perform a sensitive search for molecular lines in AB Aur as part of the IRAM Large program ASAI (A Chemical Survey of Sun-like Star-forming Regions). These data were complemented with interferometric observations of the HCO+ 1-0 and C17O 1-0 lines using the IRAM Plateau de Bure Interferometer (PdBI). Single-dish and interferometric data were used to constrain chemical models. Results. Throughout the survey, several lines of CO and its isotopologues, HCO+, H2CO, HCN, CN and CS, were detected. In addition, we detected the SO 54-33 and 56-45 lines, confirming the previous tentative detection. Comparing to other T Tauri's and Herbig Ae disks, AB Aur presents low HCN 3-2/HCO+ 3-2 and CN 2-1/HCN 3-2 line intensity ratios, similar to other transition disks. AB Aur is the only protoplanetary disk detected in SO thus far. Conclusions. We modeled the line profiles using a chemical model and a radiative transfer 3D code. Our model assumes a flared disk in hydrostatic equilibrium. The best agreement with observations was obtained for a disk with a mass of 0.01 Msun , Rin=110 AU, Rout=550 AU, a surface density radial index of 1.5 and an inclination of 27 deg. The intensities and line profiles were reproduced within a factor of 2 for most lines. This agreement is reasonable taking into account the simplicity of our model that neglects any structure within the disk. However, the HCN 3-2 and CN 2-1 line intensities were predicted more intense by a factor of >10. We discuss several scenarios to explain this discrepancy.

The nearby cloud L1642 is one of only two known very high latitude (|b| > 30 deg) clouds actively forming stars. It is a rare example of star formation in isolated conditions, and can reveal important details of star formation in general, e.g., of the effect of magnetic fields. We compare Herschel dust emission structures and magnetic field orientation revealed by Planck polarization maps in L1642. The high-resolution ($\sim20"$) Herschel data reveal a complex structure including a dense, compressed central clump, and low density striations. The Planck polarization data (at 10$'$ resolution) reveal an ordered magnetic field pervading the cloud and aligned with the surrounding striations. There is a complex interplay between the cloud structure and large scale magnetic field. This suggests that the magnetic field is closely linked to the formation and evolution of the cloud. CO rotational emission confirms that the striations are connected with the main clumps and likely to contain material either falling into or flowing out of the clumps. There is a clear transition from aligned to perpendicular structures approximately at a column density of $N_{\rm{H}} = 1.6 \times 10^{21}\, {{\rm cm}}^{-2}$. Comparing the Herschel maps with the Planck polarization maps shows the close connection between the magnetic field and cloud structure even in the finest details of the cloud., Comment: accepted to MNRAS; 14 pages; changes: moderate text and figure editing

Aims. Our goal is to determine the molecular composition of the circumstellar disk around AB Aurigae (hereafter, AB Aur). AB Aur is a prototypical Herbig Ae star and the understanding of its disk chemistry is of paramount importance to understand the chemical evolution of the gas in warm disks. Methods. We used the IRAM 30-m telescope to perform a sensitive search for molecular lines in AB Aur as part of the IRAM Large program ASAI (A Chemical Survey of Sun-like Star-forming Regions). These data were complemented with interferometric observations of the HCO+ 1-0 and C17O 1-0 lines using the IRAM Plateau de Bure Interferometer (PdBI). Single-dish and interferometric data were used to constrain chemical models. Results. Throughout the survey, several lines of CO and its isotopologues, HCO+, H2CO, HCN, CN and CS, were detected. In addition, we detected the SO 54-33 and 56-45 lines, confirming the previous tentative detection. Comparing to other T Tauri's and Herbig Ae disks, AB Aur presents low HCN 3-2/HCO+ 3-2 and CN 2-1/HCN 3-2 line intensity ratios, similar to other transition disks. AB Aur is the only protoplanetary disk detected in SO thus far. Conclusions. We modeled the line profiles using a chemical model and a radiative transfer 3D code. Our model assumes a flared disk in hydrostatic equilibrium. The best agreement with observations was obtained for a disk with a mass of 0.01 Msun , Rin=110 AU, Rout=550 AU, a surface density radial index of 1.5 and an inclination of 27 deg. The intensities and line profiles were reproduced within a factor of 2 for most lines. This agreement is reasonable taking into account the simplicity of our model that neglects any structure within the disk. However, the HCN 3-2 and CN 2-1 line intensities were predicted more intense by a factor of >10. We discuss several scenarios to explain this discrepancy.

Planck observations at 353GHz provide the first fully-sampled maps of the polarized dust emission towards interstellar filaments and their backgrounds. The polarization data provide insight on the structure of their magnetic field (B). We present the polarization maps of three nearby star forming filament of moderate column density (NH~10^22cm^-2): Musca, B211, and L1506. We use the spatial information to separate Stokes I, Q, and U of the filaments from those of their backgrounds, an essential step to measure the intrinsic polarization fraction (p) and angle (psi) of each emission component. We find that the polarization angles in the three filaments (psi_fil) are coherent along their lengths and not the same as in their backgrounds (psi_bg). The differences between psi_fil and psi_bg are 12deg, 6deg, and 54deg for Musca, B211, and L1506, respectively. These differences for Musca and L1506 are larger than the dispersions of psi, both along the filaments and in their backgrounds. The observed changes of psi are direct evidence for variations of the orientation of the plane of the sky (POS) projection of the B-field. As in previous studies, we find a decrease of several percent of p with NH. We show that the drop in p cannot be explained by random fluctuations of the orientation of B within the filaments because they are too small (sigma_psi<10deg). We recognize the degeneracy between dust alignment efficiency and the structure of B in causing variations in p, but we argue that the decrease of p from the backgrounds to the filaments results in part from depolarization associated with the 3D structure of B: both its orientation in the POS and with respect to the POS. We do not resolve the inner structure of the filaments, but at the smallest scales accessible with Planck (~0.2pc), the observed changes of psi and p hold information on the B-field structure within filaments., Comment: 16 pages, 14 figures, Submitted to Astronomy & Astrophysics

The massive star-forming region Mon R2 hosts the closest ultra-compact HII region that can be spatially resolved with current single-dish telescopes. We used the IRAM-30m telescope to carry out an unbiased spectral survey toward two important positions (namely IF and MP2), in order to studying the chemistry of deuterated molecules toward Mon R2. We found a rich chemistry of deuterated species at both positions, with detections of C2D, DCN, DNC, DCO+, D2CO, HDCO, NH2D, and N2D+ and their corresponding hydrogenated species and isotopologs. Our high spectral resolution observations allowed us to resolve three velocity components: the component at 10 km/s is detected at both positions and seems associated with the layer most exposed to the UV radiation from IRS 1; the component at 12 km/s is found toward the IF position and seems related to the molecular gas; finally, a component at 8.5 km/s is only detected toward the MP2 position, most likely related to a low-UV irradiated PDR. We derived the column density of all the species, and determined the deuterium fractions (Dfrac). The values of Dfrac are around 0.01 for all the observed species, except for HCO+ and N2H+ which have values 10 times lower. The values found in Mon R2 are well explained with pseudo-time-dependent gas-phase model in which deuteration occurs mainly via ion-molecule reactions with H2D+, CH2D+ and C2HD+. Finally, the [H13CN]/[HN13C] ratio is very high (~11) for the 10 km/s component, which also agree with our model predictions for an age of ~0.01-0.1 Myr. The deuterium chemistry is a good tool for studying star-forming regions. The low-mass star-forming regions seem well characterized with Dfrac(N2H+) or Dfrac(HCO+), but it is required a complete chemical modeling to date massive star-forming regions, because the higher gas temperature together with the rapid evolution of massive protostars., Comment: 14 pages of manuscript, 17 pages of apendix, 7 figures in the main text, accepted for publication in A&A

Massive stars drive the evolution of the interstellar medium through their radiative and mechanical energy input. After their birth, they form bubbles of hot gas surrounded by a dense shell. Traditionally, the formation of bubbles is explained through the input of a powerful stellar wind, even though direct evidence supporting this scenario is lacking. Here we explore the possibility that interstellar bubbles seen by the Spitzer- and Herschel space telescopes, blown by stars with log(L/L_sun) < 5.2, form and expand due to the thermal pressure accompanying ionization of the surrounding gas. We show that density gradients in the natal cloud or a puncture in the swept up shell lead to an ionized gas flow through the bubble into the general interstellar medium, which is traced by a dust wave near the star, demonstrating the importance of radiation pressure during this phase. Dust waves provide a natural explanation for the presence of dust inside H II bubbles, offer a novel method to study dust in H II regions and provide direct evidence that bubbles are relieving their pressure into the ISM through a champagne flow, acting as a probe of the radiative interaction of a massive star with its surroundings. We create a parameter space connecting the ambient density, the ionizing source luminosity, and the position of the dust wave, while using the well-studied H II bubbles RCW 120 and RCW 82 as benchmarks of our model. Finally, we briefly examine the implications of our study for the environments of super star clusters formed in UltraLuminous InfraRed Galaxies (ULIRGs), merging galaxies and the early Universe, which occur in very luminous and dense environments and where radiation pressure is expected to dominate the dynamical evolution., Comment: 9 pages, 6 figures: corrected typos

Observations with the Spitzer Space Telescope and the WISE satellite have revealed a prominent arc-like structure at 50" ($\simeq0.1$ pc) from the O9.5V/B0.5V system $\sigma$ Ori AB. We attribute this dust structure to the interaction of radiation pressure from the star with dust carried along by the IC 434 photo-evaporative flow of ionized gas from the dark cloud L1630. We have developed a quantitative model for the interaction of a dusty ionized flow with nearby (massive) stars where radiation pressure stalls dust, piling it up at an appreciable distance (> 0.1 pc), and force it to flow around the star. The model demonstrates that for the conditions in IC 434, the gas will decouple from the dust and will keep its original flow lines. We argue that this dust structure is the first example of a dust wave created by a massive star moving through the interstellar medium. Dust waves (and bow waves) stratify dust grains according to their radiation pressure opacity, which reflects the size distribution and composition of the grain material. Comparison of our model with observations implies that dust-gas coupling through Coulomb interaction is less important than previously thought, challenging our understanding of grain dynamics in hot, ionized regions of space. We describe the difference between dust (and bow) waves and classical bow shocks. We conclude that dust waves and bow waves should be common around stars showing the weak-wind phenomenon and that these structures are best observed at mid-IR to FIR wavelengths. In particular, dust waves and bow waves are most efficiently formed around weak-wind stars moving through a high density medium. Moreover, they provide a unique opportunity to study the direct interaction between a (massive) star and its immediate surroundings., Comment: 17 pages, 14 figures; corrected typos and added references in intro

Planck observations at 353GHz provide the first fully-sampled maps of the polarized dust emission towards interstellar filaments and their backgrounds. The polarization data provide insight on the structure of their magnetic field (B). We present the polarization maps of three nearby star forming filament of moderate column density (NH~10^22cm^-2): Musca, B211, and L1506. We use the spatial information to separate Stokes I, Q, and U of the filaments from those of their backgrounds, an essential step to measure the intrinsic polarization fraction (p) and angle (psi) of each emission component. We find that the polarization angles in the three filaments (psi_fil) are coherent along their lengths and not the same as in their backgrounds (psi_bg). The differences between psi_fil and psi_bg are 12deg, 6deg, and 54deg for Musca, B211, and L1506, respectively. These differences for Musca and L1506 are larger than the dispersions of psi, both along the filaments and in their backgrounds. The observed changes of psi are direct evidence for variations of the orientation of the plane of the sky (POS) projection of the B-field. As in previous studies, we find a decrease of several percent of p with NH. We show that the drop in p cannot be explained by random fluctuations of the orientation of B within the filaments because they are too small (sigma_psi<10deg). We recognize the degeneracy between dust alignment efficiency and the structure of B in causing variations in p, but we argue that the decrease of p from the backgrounds to the filaments results in part from depolarization associated with the 3D structure of B: both its orientation in the POS and with respect to the POS. We do not resolve the inner structure of the filaments, but at the smallest scales accessible with Planck (~0.2pc), the observed changes of psi and p hold information on the B-field structure within filaments., Comment: 16 pages, 14 figures, Submitted to Astronomy & Astrophysics

The massive star-forming region Mon R2 hosts the closest ultra-compact HII region that can be spatially resolved with current single-dish telescopes. We used the IRAM-30m telescope to carry out an unbiased spectral survey toward two important positions (namely IF and MP2), in order to studying the chemistry of deuterated molecules toward Mon R2. We found a rich chemistry of deuterated species at both positions, with detections of C2D, DCN, DNC, DCO+, D2CO, HDCO, NH2D, and N2D+ and their corresponding hydrogenated species and isotopologs. Our high spectral resolution observations allowed us to resolve three velocity components: the component at 10 km/s is detected at both positions and seems associated with the layer most exposed to the UV radiation from IRS 1; the component at 12 km/s is found toward the IF position and seems related to the molecular gas; finally, a component at 8.5 km/s is only detected toward the MP2 position, most likely related to a low-UV irradiated PDR. We derived the column density of all the species, and determined the deuterium fractions (Dfrac). The values of Dfrac are around 0.01 for all the observed species, except for HCO+ and N2H+ which have values 10 times lower. The values found in Mon R2 are well explained with pseudo-time-dependent gas-phase model in which deuteration occurs mainly via ion-molecule reactions with H2D+, CH2D+ and C2HD+. Finally, the [H13CN]/[HN13C] ratio is very high (~11) for the 10 km/s component, which also agree with our model predictions for an age of ~0.01-0.1 Myr. The deuterium chemistry is a good tool for studying star-forming regions. The low-mass star-forming regions seem well characterized with Dfrac(N2H+) or Dfrac(HCO+), but it is required a complete chemical modeling to date massive star-forming regions, because the higher gas temperature together with the rapid evolution of massive protostars., Comment: 14 pages of manuscript, 17 pages of apendix, 7 figures in the main text, accepted for publication in A&A

Massive stars drive the evolution of the interstellar medium through their radiative and mechanical energy input. After their birth, they form bubbles of hot gas surrounded by a dense shell. Traditionally, the formation of bubbles is explained through the input of a powerful stellar wind, even though direct evidence supporting this scenario is lacking. Here we explore the possibility that interstellar bubbles seen by the Spitzer- and Herschel space telescopes, blown by stars with log(L/L_sun) < 5.2, form and expand due to the thermal pressure accompanying ionization of the surrounding gas. We show that density gradients in the natal cloud or a puncture in the swept up shell lead to an ionized gas flow through the bubble into the general interstellar medium, which is traced by a dust wave near the star, demonstrating the importance of radiation pressure during this phase. Dust waves provide a natural explanation for the presence of dust inside H II bubbles, offer a novel method to study dust in H II regions and provide direct evidence that bubbles are relieving their pressure into the ISM through a champagne flow, acting as a probe of the radiative interaction of a massive star with its surroundings. We create a parameter space connecting the ambient density, the ionizing source luminosity, and the position of the dust wave, while using the well-studied H II bubbles RCW 120 and RCW 82 as benchmarks of our model. Finally, we briefly examine the implications of our study for the environments of super star clusters formed in UltraLuminous InfraRed Galaxies (ULIRGs), merging galaxies and the early Universe, which occur in very luminous and dense environments and where radiation pressure is expected to dominate the dynamical evolution., Comment: 9 pages, 6 figures: corrected typos

Observations with the Spitzer Space Telescope and the WISE satellite have revealed a prominent arc-like structure at 50" ($\simeq0.1$ pc) from the O9.5V/B0.5V system $\sigma$ Ori AB. We attribute this dust structure to the interaction of radiation pressure from the star with dust carried along by the IC 434 photo-evaporative flow of ionized gas from the dark cloud L1630. We have developed a quantitative model for the interaction of a dusty ionized flow with nearby (massive) stars where radiation pressure stalls dust, piling it up at an appreciable distance (> 0.1 pc), and force it to flow around the star. The model demonstrates that for the conditions in IC 434, the gas will decouple from the dust and will keep its original flow lines. We argue that this dust structure is the first example of a dust wave created by a massive star moving through the interstellar medium. Dust waves (and bow waves) stratify dust grains according to their radiation pressure opacity, which reflects the size distribution and composition of the grain material. Comparison of our model with observations implies that dust-gas coupling through Coulomb interaction is less important than previously thought, challenging our understanding of grain dynamics in hot, ionized regions of space. We describe the difference between dust (and bow) waves and classical bow shocks. We conclude that dust waves and bow waves should be common around stars showing the weak-wind phenomenon and that these structures are best observed at mid-IR to FIR wavelengths. In particular, dust waves and bow waves are most efficiently formed around weak-wind stars moving through a high density medium. Moreover, they provide a unique opportunity to study the direct interaction between a (massive) star and its immediate surroundings., Comment: 17 pages, 14 figures; corrected typos and added references in intro

Context. Monoceros R2 (Mon R2), at a distance of 830 pc, is the only ultra-compact H ii region (UC H ii) where its associated photon-dominated region (PDR) can be resolved with the Herschel Space Observatory. Aims. Our aim is to investigate observationally the kinematical patterns in the interface regions (i.e., the transition from atomic to molecular gas) associated with Mon R2. Methods. We used the HIFI instrument onboard Herschel to observe the line profiles of the reactive ions CH+, OH+ and H2O+ toward different positions in Mon R2. We derive the column density of these molecules and compare them with gas-phase chemistry models. Results. The reactive ion CH+ is detected both in emission (at central and red-shifted velocities) and in absorption (at blue-shifted velocities). OH+ is detected in absorption at both blue- and red-shifted velocities, with similar column densities. H2O+ is not detected at any of the positions, down to a rms of 40 mK toward the molecular peak. At this position, we find that the OH+ absorption originates in a mainly atomic medium, and therefore is associated with the most exposed layers of the PDR. These results are consistent with the predictions from photo-chemical models. The line profiles are consistent with the atomic gas being entrained in the ionized gas flow along the walls of the cavity of the H ii region. Based on this evidence, we are able to propose a new geometrical model for this region. Conclusions. The kinematical patterns of the OH+ and CH+ absorption indicate the existence of a layer of mainly atomic gas for which we have derived, for the first time, some physical parameters and its dynamics., Comment: 6 pages, 5 figures. Accepted for publication in A&A

We study the chemistry of small hydrocarbons in the photon-dominated regions (PDRs) associated with the ultra-compact HII region Mon R2. Our goal is to determine the variations of the abundance of small hydrocarbons in a high-UV irradiated PDR and investigate their chemistry. We present an observational study of CH, CCH and c-C$_3$H$_2$ in Mon R2 combining data obtained with the IRAM 30m telescope and Herschel. We determine the column densities of these species, and compare their spatial distributions with that of polycyclic aromatic hydrocarbon (PAH). We compare the observational results with different chemical models to explore the relative importance of gas-phase, grain-surface and time-dependent chemistry in these environments. The emission of the small hydrocarbons show different patterns. The CCH emission is extended while CH and c-C$_3$H$_2$ are concentrated towards the more illuminated layers of the PDR. The ratio of the column densities of c-C$_3$H$_2$ and CCH shows spatial variations up to a factor of a few, increasing from $N(c-C$_3$H$_2$)/N(CCH)\approx0.004$ in the envelope to a maximum of $\sim0.015-0.029$ towards the 8$\mu$m emission peak. Comparing these results with other galactic PDRs, we find that the abundance of CCH is quite constant over a wide range of G$_0$, whereas the abundance of c-C$_3$H$_2$ is higher in low-UV PDRs. In Mon R2, the gas-phase steady-state chemistry can account relatively well for the abundances of CH and CCH in the most exposed layers of the PDR, but falls short by a factor of 10 to reproduce c-C$_3$H$_2$. In the molecular envelope, time-dependent effects and grain surface chemistry play a dominant role in determining the hydrocarbons abundances. Our study shows that CCH and c-C$_3$H$_2$ present a complex chemistry in which UV photons, grain-surface chemistry and time dependent effects contribute to determine their abundances., Comment: 18 pages, 11 figures, 7 tables. Proposed for acceptance in A&A. Abstract abridged

Context. Monoceros R2 (Mon R2), at a distance of 830 pc, is the only ultra-compact H ii region (UC H ii) where its associated photon-dominated region (PDR) can be resolved with the Herschel Space Observatory. Aims. Our aim is to investigate observationally the kinematical patterns in the interface regions (i.e., the transition from atomic to molecular gas) associated with Mon R2. Methods. We used the HIFI instrument onboard Herschel to observe the line profiles of the reactive ions CH+, OH+ and H2O+ toward different positions in Mon R2. We derive the column density of these molecules and compare them with gas-phase chemistry models. Results. The reactive ion CH+ is detected both in emission (at central and red-shifted velocities) and in absorption (at blue-shifted velocities). OH+ is detected in absorption at both blue- and red-shifted velocities, with similar column densities. H2O+ is not detected at any of the positions, down to a rms of 40 mK toward the molecular peak. At this position, we find that the OH+ absorption originates in a mainly atomic medium, and therefore is associated with the most exposed layers of the PDR. These results are consistent with the predictions from photo-chemical models. The line profiles are consistent with the atomic gas being entrained in the ionized gas flow along the walls of the cavity of the H ii region. Based on this evidence, we are able to propose a new geometrical model for this region. Conclusions. The kinematical patterns of the OH+ and CH+ absorption indicate the existence of a layer of mainly atomic gas for which we have derived, for the first time, some physical parameters and its dynamics., Comment: 6 pages, 5 figures. Accepted for publication in A&A

We study the chemistry of small hydrocarbons in the photon-dominated regions (PDRs) associated with the ultra-compact HII region Mon R2. Our goal is to determine the variations of the abundance of small hydrocarbons in a high-UV irradiated PDR and investigate their chemistry. We present an observational study of CH, CCH and c-C$_3$H$_2$ in Mon R2 combining data obtained with the IRAM 30m telescope and Herschel. We determine the column densities of these species, and compare their spatial distributions with that of polycyclic aromatic hydrocarbon (PAH). We compare the observational results with different chemical models to explore the relative importance of gas-phase, grain-surface and time-dependent chemistry in these environments. The emission of the small hydrocarbons show different patterns. The CCH emission is extended while CH and c-C$_3$H$_2$ are concentrated towards the more illuminated layers of the PDR. The ratio of the column densities of c-C$_3$H$_2$ and CCH shows spatial variations up to a factor of a few, increasing from $N(c-C$_3$H$_2$)/N(CCH)\approx0.004$ in the envelope to a maximum of $\sim0.015-0.029$ towards the 8$\mu$m emission peak. Comparing these results with other galactic PDRs, we find that the abundance of CCH is quite constant over a wide range of G$_0$, whereas the abundance of c-C$_3$H$_2$ is higher in low-UV PDRs. In Mon R2, the gas-phase steady-state chemistry can account relatively well for the abundances of CH and CCH in the most exposed layers of the PDR, but falls short by a factor of 10 to reproduce c-C$_3$H$_2$. In the molecular envelope, time-dependent effects and grain surface chemistry play a dominant role in determining the hydrocarbons abundances. Our study shows that CCH and c-C$_3$H$_2$ present a complex chemistry in which UV photons, grain-surface chemistry and time dependent effects contribute to determine their abundances., Comment: 18 pages, 11 figures, 7 tables. Proposed for acceptance in A&A. Abstract abridged

Context. Mon R2 is the only ultracompact HII region (UCHII) where the associated photon-dominated region (PDR) can be resolved with Herschel. Due to its brightness and proximity, it is the best source to investigate the chemistry and physics of highly UV-irradiated PDRs. Aims. Our goal is to estimate the abundance of H2O and NH3 in this region and investigate their origin. Methods. We present new observations obtained with HIFI and the IRAM-30m telescope. Using a large velocity gradient approach, we model the line intensities and derive an average abundance of H2O and NH3 across the region. Finally, we model the line profiles with a non-local radiative transfer model and compare these results with the abundance predicted by the Meudon PDR code. Results. The variations of the line profiles and intensities indicate complex geometrical and kinematical patterns. The H2O lines present a strong absorption at the ambient velocity and emission in high velocity wings towards the HII region. The spatial distribution of the o-H2^18O line shows that the its emission arises in the PDR surrounding the HII region. By modeling the o-H2^18O emission we derive a mean abundance of o-H2O of ~10^-8 relative to H2. The ortho-H2O abundance is however larger, ~1x10^-7, in the high velocity wings. Possible explanations for this larger abundance include an expanding hot PDR and/or an outflow. Ammonia seems to be present only in the envelope with an average abundance of ~2x10^-9 relative to H2. Conclusions. The Meudon PDR code can account for the measured water abundance in the high velocity gas as long as we assume that it originates from a <1 mag hot expanding layer of the PDR, i.e. that the outflow has only a minor contribution to this emission. To explain the abundances in the rest of the cloud the molecular freeze out and grain surface chemistry would need to be included., Comment: 12 pages, 7 figures, 3 tables. Accepted for publication in A&A. Abstract shortened. Updated references, language editing applied in v2

Context. In photo-dissociation regions (PDRs), Polycyclic Aromatic Hydrocarbons (PAHs) could be produced by evaporation of Very Small Grains (VSGs) by the impinging UV radiation field from a nearby star. Aims. We investigate quantitatively the transition zone between evaporating Very Small Grains (eVSGs) and PAHs in several PDRs. Methods. We study the relative contribution of PAHs and eVSGs to the mid-IR emission in a wide range of excitation conditions. We fit the observed mid-IR emission of PDRs by using a set of template band emission spectra of PAHs, eVSGs and gas lines. The fitting tool PAHTAT (PAH Toulouse Astronomical Templates) is made available to the community as an IDL routine. From the results of the fit, we derive the fraction of carbon f_eVSG locked in eVSGs and compare it to the intensity of the local UV radiation field. Results. We show a clear decrease of f_eVSG with increasing intensity of the local UV radiation field, which supports the scenario of photo-destruction of eVSGs. Conversely, this dependence can be used to quantify the intensity of the UV radiation field for different PDRs, including non resolved ones. Conclusions. PAHTAT can be used to trace the intensity of the local UV radiation field in regions where eVSGs evaporate, which correspond to relatively dense (nH = [100, 10^5 ] cm-3) and UV irradiated PDRs (G0 = [100, 5x10^4]) where H2 emits in rotational lines., Comment: 13 pages, 11 figures. Accepted for publication in A&A. Typos corrected

Ultracompact (UC) HII regions constitute one of the earliest phases in the formation of a massive star and are characterized by extreme physical conditions (Go>10^5 Habing field and n>10^6 cm^-3). The UC HII Mon R2 is the closest one and therefore an excellent target to study the chemistry in these complex regions. We carried out a 3mm and 1mm spectral survey using the IRAM 30-m telescope towards three positions that represent different physical environments in Mon R2: (i) the ionization front (IF) at (0",0"); two peaks in the molecular cloud (ii) MP1 at the offset (+15",-15") and (iii) MP2 at the farther offset (0",40"). In addition, we carried out extensive modeling to explain the chemical differences between the three observed regions. We detected more than thirty different species. We detected SO+ and C4H suggesting that UV radiation plays an important role in the molecular chemistry of this region. We detected the typical PDR molecules CN, HCN, HCO, C2H, and c-C3H2. While the IF and the MP1 have a chemistry similar to that found in high UV field and dense PDRs like the Orion Bar, the MP2 is more similar to lower UV/density PDRs like the Horsehead nebula. We also detected complex molecules that are not usually found in PDRs (CH3CN, H2CO, HC3N, CH3OH and CH3C2H). Sulfur compounds CS, HCS+, C2S, H2CS, SO and SO2 and the deuterated species DCN and C2D were also identified. [DCN]/[HCN]=0.03 and [C2D]/[C2H]=0.05, are among the highest in warm regions. Our results show that the high UV/dense PDRs present a different chemistry from that of the low UV case. Abundance ratios like [CO+]/[HCO+] or [HCO]/[HCO+] are good diagnostics to differentiate between them. In Mon R2 we have the two classes of PDRs, a high UV PDR towards the IF and the adjacent molecular bar and a low-UV PDR which extends towards the north-west following the border of the cloud., Comment: 31 pages

Context. Mon R2 is the only ultracompact HII region (UCHII) where the associated photon-dominated region (PDR) can be resolved with Herschel. Due to its brightness and proximity, it is the best source to investigate the chemistry and physics of highly UV-irradiated PDRs. Aims. Our goal is to estimate the abundance of H2O and NH3 in this region and investigate their origin. Methods. We present new observations obtained with HIFI and the IRAM-30m telescope. Using a large velocity gradient approach, we model the line intensities and derive an average abundance of H2O and NH3 across the region. Finally, we model the line profiles with a non-local radiative transfer model and compare these results with the abundance predicted by the Meudon PDR code. Results. The variations of the line profiles and intensities indicate complex geometrical and kinematical patterns. The H2O lines present a strong absorption at the ambient velocity and emission in high velocity wings towards the HII region. The spatial distribution of the o-H2^18O line shows that the its emission arises in the PDR surrounding the HII region. By modeling the o-H2^18O emission we derive a mean abundance of o-H2O of ~10^-8 relative to H2. The ortho-H2O abundance is however larger, ~1x10^-7, in the high velocity wings. Possible explanations for this larger abundance include an expanding hot PDR and/or an outflow. Ammonia seems to be present only in the envelope with an average abundance of ~2x10^-9 relative to H2. Conclusions. The Meudon PDR code can account for the measured water abundance in the high velocity gas as long as we assume that it originates from a <1 mag hot expanding layer of the PDR, i.e. that the outflow has only a minor contribution to this emission. To explain the abundances in the rest of the cloud the molecular freeze out and grain surface chemistry would need to be included., Comment: 12 pages, 7 figures, 3 tables. Accepted for publication in A&A. Abstract shortened. Updated references, language editing applied in v2

Context. In photo-dissociation regions (PDRs), Polycyclic Aromatic Hydrocarbons (PAHs) could be produced by evaporation of Very Small Grains (VSGs) by the impinging UV radiation field from a nearby star. Aims. We investigate quantitatively the transition zone between evaporating Very Small Grains (eVSGs) and PAHs in several PDRs. Methods. We study the relative contribution of PAHs and eVSGs to the mid-IR emission in a wide range of excitation conditions. We fit the observed mid-IR emission of PDRs by using a set of template band emission spectra of PAHs, eVSGs and gas lines. The fitting tool PAHTAT (PAH Toulouse Astronomical Templates) is made available to the community as an IDL routine. From the results of the fit, we derive the fraction of carbon f_eVSG locked in eVSGs and compare it to the intensity of the local UV radiation field. Results. We show a clear decrease of f_eVSG with increasing intensity of the local UV radiation field, which supports the scenario of photo-destruction of eVSGs. Conversely, this dependence can be used to quantify the intensity of the UV radiation field for different PDRs, including non resolved ones. Conclusions. PAHTAT can be used to trace the intensity of the local UV radiation field in regions where eVSGs evaporate, which correspond to relatively dense (nH = [100, 10^5 ] cm-3) and UV irradiated PDRs (G0 = [100, 5x10^4]) where H2 emits in rotational lines., Comment: 13 pages, 11 figures. Accepted for publication in A&A. Typos corrected

Ultracompact (UC) HII regions constitute one of the earliest phases in the formation of a massive star and are characterized by extreme physical conditions (Go>10^5 Habing field and n>10^6 cm^-3). The UC HII Mon R2 is the closest one and therefore an excellent target to study the chemistry in these complex regions. We carried out a 3mm and 1mm spectral survey using the IRAM 30-m telescope towards three positions that represent different physical environments in Mon R2: (i) the ionization front (IF) at (0",0"); two peaks in the molecular cloud (ii) MP1 at the offset (+15",-15") and (iii) MP2 at the farther offset (0",40"). In addition, we carried out extensive modeling to explain the chemical differences between the three observed regions. We detected more than thirty different species. We detected SO+ and C4H suggesting that UV radiation plays an important role in the molecular chemistry of this region. We detected the typical PDR molecules CN, HCN, HCO, C2H, and c-C3H2. While the IF and the MP1 have a chemistry similar to that found in high UV field and dense PDRs like the Orion Bar, the MP2 is more similar to lower UV/density PDRs like the Horsehead nebula. We also detected complex molecules that are not usually found in PDRs (CH3CN, H2CO, HC3N, CH3OH and CH3C2H). Sulfur compounds CS, HCS+, C2S, H2CS, SO and SO2 and the deuterated species DCN and C2D were also identified. [DCN]/[HCN]=0.03 and [C2D]/[C2H]=0.05, are among the highest in warm regions. Our results show that the high UV/dense PDRs present a different chemistry from that of the low UV case. Abundance ratios like [CO+]/[HCO+] or [HCO]/[HCO+] are good diagnostics to differentiate between them. In Mon R2 we have the two classes of PDRs, a high UV PDR towards the IF and the adjacent molecular bar and a low-UV PDR which extends towards the north-west following the border of the cloud., Comment: 31 pages

We present the first complete 55-671 um spectral scan of a low-mass Class 0 protostar (Serpens SMM1) taken with the PACS and SPIRE spectrometers on board Herschel. More than 145 lines have been detected, most of them rotationally excited lines of 12CO (full ladder from J=4-3 to 42-41), H2O, OH, 13CO, HCN and HCO+ . Bright [OI]63,145um and weaker [CII]158 and [CI]370,609um lines are also detected. Mid-IR spectra retrieved from the Spitzer archive are also first discussed here, they show clear detections of [NeII], [FeII], [SiII] and [SI] fine structure lines as well as weaker H2 S(1) and S(2) pure rotational lines. The observed line luminosity is dominated by CO (~54%), H2O (~22%), [OI] (~12%) and OH (~9%) emission. A non-LTE radiative transfer model allowed us to quantify the contribution of the 3 different temperature components suggested by the 12CO rotational ladder (Tk(hot)~800 K, Tk(warm)~375 K and Tk(cool)~150 K). Gas densities n(H2)~5x10^6 cm^-3 are needed to reproduce the observed far-IR lines arising from shocks in the inner protostellar envelope for which we derive upper limit abundances of x(CO)~10^-4, x(H2O)~0.2x10^-5 and x(OH)~10^-6. The lower energy submm 12CO and H2O lines show more extended emission that we associate with the cool entrained outflow gas. Fast dissociative J-shocks (v_s > 60 km s^-1) as well as lower velocity non-dissociative shocks (v_s

We present the first complete 55-671 um spectral scan of a low-mass Class 0 protostar (Serpens SMM1) taken with the PACS and SPIRE spectrometers on board Herschel. More than 145 lines have been detected, most of them rotationally excited lines of 12CO (full ladder from J=4-3 to 42-41), H2O, OH, 13CO, HCN and HCO+ . Bright [OI]63,145um and weaker [CII]158 and [CI]370,609um lines are also detected. Mid-IR spectra retrieved from the Spitzer archive are also first discussed here, they show clear detections of [NeII], [FeII], [SiII] and [SI] fine structure lines as well as weaker H2 S(1) and S(2) pure rotational lines. The observed line luminosity is dominated by CO (~54%), H2O (~22%), [OI] (~12%) and OH (~9%) emission. A non-LTE radiative transfer model allowed us to quantify the contribution of the 3 different temperature components suggested by the 12CO rotational ladder (Tk(hot)~800 K, Tk(warm)~375 K and Tk(cool)~150 K). Gas densities n(H2)~5x10^6 cm^-3 are needed to reproduce the observed far-IR lines arising from shocks in the inner protostellar envelope for which we derive upper limit abundances of x(CO)~10^-4, x(H2O)~0.2x10^-5 and x(OH)~10^-6. The lower energy submm 12CO and H2O lines show more extended emission that we associate with the cool entrained outflow gas. Fast dissociative J-shocks (v_s > 60 km s^-1) as well as lower velocity non-dissociative shocks (v_s

We present the first complete 55-671 um spectral scan of a low-mass Class 0 protostar (Serpens SMM1) taken with the PACS and SPIRE spectrometers on board Herschel. More than 145 lines have been detected, most of them rotationally excited lines of 12CO (full ladder from J=4-3 to 42-41), H2O, OH, 13CO, HCN and HCO+ . Bright [OI]63,145um and weaker [CII]158 and [CI]370,609um lines are also detected. Mid-IR spectra retrieved from the Spitzer archive are also first discussed here, they show clear detections of [NeII], [FeII], [SiII] and [SI] fine structure lines as well as weaker H2 S(1) and S(2) pure rotational lines. The observed line luminosity is dominated by CO (~54%), H2O (~22%), [OI] (~12%) and OH (~9%) emission. A non-LTE radiative transfer model allowed us to quantify the contribution of the 3 different temperature components suggested by the 12CO rotational ladder (Tk(hot)~800 K, Tk(warm)~375 K and Tk(cool)~150 K). Gas densities n(H2)~5x10^6 cm^-3 are needed to reproduce the observed far-IR lines arising from shocks in the inner protostellar envelope for which we derive upper limit abundances of x(CO)~10^-4, x(H2O)~0.2x10^-5 and x(OH)~10^-6. The lower energy submm 12CO and H2O lines show more extended emission that we associate with the cool entrained outflow gas. Fast dissociative J-shocks (v_s > 60 km s^-1) as well as lower velocity non-dissociative shocks (v_s

We present the first complete 55-671 um spectral scan of a low-mass Class 0 protostar (Serpens SMM1) taken with the PACS and SPIRE spectrometers on board Herschel. More than 145 lines have been detected, most of them rotationally excited lines of 12CO (full ladder from J=4-3 to 42-41), H2O, OH, 13CO, HCN and HCO+ . Bright [OI]63,145um and weaker [CII]158 and [CI]370,609um lines are also detected. Mid-IR spectra retrieved from the Spitzer archive are also first discussed here, they show clear detections of [NeII], [FeII], [SiII] and [SI] fine structure lines as well as weaker H2 S(1) and S(2) pure rotational lines. The observed line luminosity is dominated by CO (~54%), H2O (~22%), [OI] (~12%) and OH (~9%) emission. A non-LTE radiative transfer model allowed us to quantify the contribution of the 3 different temperature components suggested by the 12CO rotational ladder (Tk(hot)~800 K, Tk(warm)~375 K and Tk(cool)~150 K). Gas densities n(H2)~5x10^6 cm^-3 are needed to reproduce the observed far-IR lines arising from shocks in the inner protostellar envelope for which we derive upper limit abundances of x(CO)~10^-4, x(H2O)~0.2x10^-5 and x(OH)~10^-6. The lower energy submm 12CO and H2O lines show more extended emission that we associate with the cool entrained outflow gas. Fast dissociative J-shocks (v_s > 60 km s^-1) as well as lower velocity non-dissociative shocks (v_s

The aromatic infrared bands (AIBs) observed in the mid infrared spectrum are attributed to Polycyclic Aromatic Hydrocarbons (PAHs). We observe the NGC 7023-North West (NW) PDR in the mid-infrared (10 - 19.5 micron) using the Infrared Spectrometer (IRS), on board Spitzer. Clear variations are observed in the spectra, most notably the ratio of the 11.0 to 11.2 micron bands, the peak position of the 11.2 and 12.0 micron bands, and the degree of asymmetry of the 11.2 micron band. The observed variations appear to change as a function of position within the PDR. We aim to explain these variations by a change in the abundances of the emitting components of the PDR. A Blind Signal Separation (BSS) method, i.e. a Non-Negative Matrix Factorization algorithm is applied to separate the observed spectrum into components. Using the NASA Ames PAH IR Spectroscopic Database, these extracted signals are fit. The observed signals alone were also fit using the database and these components are compared to the BSS components. Three component signals were extracted from the observation using BSS. We attribute the three signals to ionized PAHs, neutral PAHs, and Very Small Grains (VSGs). The fit of the BSS extracted spectra with the PAH database further confirms the attribution to ionized and neutral PAHs and provides confidence in both methods for producing reliable results. The 11.0 micron feature is attributed to PAH cations while the 11.2 micron band is attributed to neutral PAHs. The VSG signal shows a characteristically asymmetric broad feature at 11.3 micron with an extended red wing. By combining the NASA Ames PAH IR Spectroscopic Database fit with the BSS method, the independent results of each method can be confirmed and some limitations of each method are overcome.

The aromatic infrared bands (AIBs) observed in the mid infrared spectrum are attributed to Polycyclic Aromatic Hydrocarbons (PAHs). We observe the NGC 7023-North West (NW) PDR in the mid-infrared (10 - 19.5 micron) using the Infrared Spectrometer (IRS), on board Spitzer. Clear variations are observed in the spectra, most notably the ratio of the 11.0 to 11.2 micron bands, the peak position of the 11.2 and 12.0 micron bands, and the degree of asymmetry of the 11.2 micron band. The observed variations appear to change as a function of position within the PDR. We aim to explain these variations by a change in the abundances of the emitting components of the PDR. A Blind Signal Separation (BSS) method, i.e. a Non-Negative Matrix Factorization algorithm is applied to separate the observed spectrum into components. Using the NASA Ames PAH IR Spectroscopic Database, these extracted signals are fit. The observed signals alone were also fit using the database and these components are compared to the BSS components. Three component signals were extracted from the observation using BSS. We attribute the three signals to ionized PAHs, neutral PAHs, and Very Small Grains (VSGs). The fit of the BSS extracted spectra with the PAH database further confirms the attribution to ionized and neutral PAHs and provides confidence in both methods for producing reliable results. The 11.0 micron feature is attributed to PAH cations while the 11.2 micron band is attributed to neutral PAHs. The VSG signal shows a characteristically asymmetric broad feature at 11.3 micron with an extended red wing. By combining the NASA Ames PAH IR Spectroscopic Database fit with the BSS method, the independent results of each method can be confirmed and some limitations of each method are overcome.

We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species., Comment: Accepted for publication in Astronomy and Astrophysics Letters (Herschel HIFI special issue), 5 pages, 5 figures