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.

We investigate inequalities in access to permanent positions in professional astronomy in France, focusing on the hiring stage. We use results from a national survey conducted on behalf of the French society of astronomy and astrophysics (SF2A) aimed at young astronomers holding a PhD obtained in France, and answered by over 300 researchers. We find that women are nearly two times less likely than men to be selected by the (national or local) committees attributing permanent positions ($p=0.06$). We also find that applicants who did their undergraduate studies in an elite school ("Grande \'Ecole"), where women are largely under-represented, rather than in a university, are nearly three times more likely to succeed in obtaining a position ($p=0.0026$). Our analysis suggests the existence of two biases in committees attributing permanent positions in astronomy in France: a gender bias, and a form of elitism. These biases against women in their professional life impacts their personal life as our survey shows that a larger fraction of them declare that having children can have a negative effect on their careers. They are half as many as men having children in the sample. National committees (such as the CNRS) have acknowledged this issue for several years now, hence one can hope that changes will be seen in the next decade.

This paper aims at providing a comprehensive framework to generate an astrophysical scene and to simulate realistic hyperspectral and multispectral data acquired by two JWST instruments, namely NIRCam Imager and NIRSpec IFU. We want to show that this simulation framework can be resorted to assess the benefits of fusing these images to recover an image of high spatial and spectral resolutions. To do so, we create a synthetic scene associated with a canonical infrared source, the Orion Bar. This scene combines pre-existing modelled spectra provided by the JWST Early Release Science Program 1288 and real high resolution spatial maps from the Hubble space and ALMA telescopes. We develop forward models including corresponding noises for the two JWST instruments based on their technical designs and physical features. JWST observations are then simulated by applying the forward models to the aforementioned synthetic scene. We test a dedicated fusion algorithm we developed on these simulated observations. We show the fusion process reconstructs the high spatio-spectral resolution scene with a good accuracy on most areas, and we identify some limitations of the method to be tackled in future works. The synthetic scene and observations presented in the paper are made publicly available and can be used for instance to evaluate instrument models (aboard the JWST or on the ground), pipelines, or more sophisticated algorithms dedicated to JWST data analysis. Besides, fusion methods such as the one presented in this paper are shown to be promising tools to fully exploit the unprecedented capabilities of the JWST.

Hyperspectral imaging has become a significant source of valuable data for astronomers over the past decades. Current instrumental and observing time constraints allow direct acquisition of multispectral images, with high spatial but low spectral resolution, and hyperspectral images, with low spatial but high spectral resolution. To enhance scientific interpretation of the data, we propose a data fusion method which combines the benefits of each image to recover a high spatio-spectral resolution datacube. The proposed inverse problem accounts for the specificities of astronomical instruments, such as spectrally variant blurs. We provide a fast implementation by solving the problem in the frequency domain and in a low-dimensional subspace to efficiently handle the convolution operators as well as the high dimensionality of the data. We conduct experiments on a realistic synthetic dataset of simulated observation of the upcoming James Webb Space Telescope, and we show that our fusion algorithm outperforms state-of-the-art methods commonly used in remote sensing for Earth observation.

The JWST will deliver large data sets of high-quality spectral data over the 0.6-28 $\mu$m range. It will combine sensitivity, spectral and spatial resolution. Specific tools are required to provide efficient scientific analysis of such large data sets. Our aim is to illustrate the potential of unsupervised learning methods to get insights into chemical variations in the populations that carry the aromatic infrared bands (AIBs), more specifically PAH species and carbonaceous very small grains (VSGs). We present a method based on linear fitting and blind signal separation (BSS) for extracting representative spectra for a spectral data set. The method is fast and robust, which ensures its applicability to JWST spectral cubes. We tested this method on a sample of ISO-SWS data, which resemble most closely the JWST spectra in terms of spectral resolution and coverage. Four representative spectra were extracted. Their main characteristics appear consistent with previous studies with populations dominated by cationic PAHs, neutral PAHs, evaporating VSGs, and large ionized PAHs, known as the PAH$^x$ population. In addition, the 3 $\mu$m range, which is considered here for the first time in a BSS method, reveals the presence of aliphatics connected to neutral PAHs. Each representative spectrum is found to carry second-order spectral signatures (e.g. small bands), which are connected with the underlying chemical diversity of populations. However, the precise attribution of theses signatures remains limited by the combined small size and heterogeneity of the sample of astronomical spectra available in this study. The upcoming JWST data will allow us to overcome this limitation. The large data sets of hyperspectral images provided by JWST analysed with the proposed method, which is fast and robust, will open promising perspectives for our understanding of the chemical evolution of the AIB carriers.

One serious challenge for planet formation is the rapid inward drift of pebble-sized dust particles in protoplanetary disks. Dust trapping at local maxima in the disk gas pressure has received much theoretical attention but still lacks observational support. The cold dust emission in the AB Aur disk forms an asymmetric ring at a radius of about 120 au, which is suggestive of dust trapping in a gas vortex. We present high spatial resolution (0".58x0".78 ~ 80x110 au) NOEMA observations of the 1.12 mm and 2.22 mm dust continuum emission from the AB Aur disk. Significant azimuthal variations of the flux ratio at both wavelengths indicate a size segregation of the large dust particles along the ring. Our continuum images also show that the intensity variations along the ring are smaller at 2.22 mm than at 1.12 mm, contrary to what dust trapping models with a gas vortex have predicted. Our two-fluid (gas+dust) hydrodynamical simulations demonstrate that this feature is well explained if the gas vortex has started to decay due to turbulent diffusion, and dust particles are thus losing the azimuthal trapping on different timescales depending on their size. The comparison between our observations and simulations allows us to constrain the size distribution and the total mass of solid particles in the ring, which we find to be of the order of 30 Earth masses, enough to form future rocky planets.

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.

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.

Protoplanetary disks undergo substantial mass-loss by photoevaporation, a mechanism which is crucial to their dynamical evolution. However, the processes regulating the gas energetics have not been well constrained by observations so far. We aim at studying the processes involved in disk photoevaporation when it is driven by far-UV photons. We present a unique Herschel survey and new ALMA observations of four externally-illuminated photoevaporating disks (a.k.a. proplyds). For the analysis of these data, we developed a 1D model of the photodissociation region (PDR) of a proplyd, based on the Meudon PDR code and computed the far infrared line emission. We successfully reproduce most of the observations and derive key physical parameters, i.e. densities at the disk surface of about $10^{6}$ cm$^{-3}$ and local gas temperatures of about 1000 K. Our modelling suggests that all studied disks are found in a transitional regime resulting from the interplay between several heating and cooling processes that we identify. These differ from those dominating in classical PDRs, i.e. grain photo-electric effect and cooling by [OI] and [CII] FIR lines. This energetic regime is associated to an equilibrium dynamical point of the photoevaporation flow: the mass-loss rate is self-regulated to set the envelope column density at a value that maintains the temperature at the disk surface around 1000 K. From our best-fit models, we estimate mass-loss rates - of the order of $10^{-7}$ $\mathrm{M}_\odot$/yr - that are in agreement with earlier spectroscopic observation of ionised gas tracers. This holds only if we assume an evaporation flow launched from the disk surface at sound speed (supercritical regime). We have identified the energetic regime regulating FUV-photoevaporation in proplyds. This regime could be implemented into models of the dynamical evolution of protoplanetary disks.

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.

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.

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.

The near- and mid-IR spectrum of many astronomical objects is dominated by emission bands due to UV-excited polycyclic aromatic hydrocarbons (PAH) and evaporating very small grains (eVSG). Previous studies with the ISO, Spitzer and AKARI space telescopes have shown that the spectral variations of these features are directly related to the local physical conditions that induce a photo-chemical evolution of the band carriers. Because of the limited sensitivity and spatial resolution, these studies have focused mainly on galactic star-forming regions. We discuss how the advent of JWST will allow to extend these studies to previously unresolved sources such as near-by galaxies, and how the analysis of the infrared signatures of PAHs and eVSGs can be used to determine their physical conditions and chemical composition.

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.

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.

Recent studies have confirmed the presence of buckminsterfullerene (C$_{60}$) in different interstellar and circumstellar environments. However, several aspects regarding C$_{60}$ in space are not well understood yet, such as the formation and excitation processes, and the connection between C$_{60}$ and other carbonaceous compounds in the interstellar medium, in particular polycyclic aromatic hydrocarbons (PAHs). In this paper we study several photodissociation regions (PDRs) where C$_{60}$ and PAHs are detected and the local physical conditions are reasonably well constrained, to provide observational insights into these questions. C$_{60}$ is found to emit in PDRs where the dust is cool ($T_d = 20-40$ K) and even in PDRs with cool stars. These results exclude the possibility for C$_{60}$ to be locked in grains at thermal equilibrium in these environments. We observe that PAH and C$_{60}$ emission are spatially uncorrelated and that C$_{60}$ is present in PDRs where the physical conditions (in terms of radiation field and hydrogen density) allow for full dehydrogenation of PAHs, with the exception of Ced 201. We also find trends indicative of an increase in C$_{60}$ abundance within individual PDRs, but these trends are not universal. These results support models where the dehydrogenation of carbonaceous species is the first step towards C$_{60}$ formation. However, this is not the only parameter involved and C$_{60}$ formation is likely affected by shocks and PDR age.

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.

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.

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.

The mid-infrared (IR; 5-15~$\mu$m) spectrum of a wide variety of astronomical objects exhibits a set of broad emission features at 6.2, 7.7, 8.6, 11.3 and 12.7 $\mu$m. About 30 years ago it was proposed that these signatures are due to emission from a family of UV heated nanometer-sized carbonaceous molecules known as polycyclic aromatic hydrocarbons (PAHs), causing them to be referred to as aromatic IR bands (AIBs). Today, the acceptance of the PAH model is far from settled, as the identification of a single PAH in space has not yet been successful and physically relevant theoretical models involving ``true'' PAH cross sections do not reproduce the AIBs in detail. In this paper, we use the NASA Ames PAH IR Spectroscopic Database, which contains over 500 quantum-computed spectra, in conjunction with a simple emission model, to show that the spectrum produced by any random mixture of at least 30 PAHs converges to the same 'kernel'-spectrum. This kernel-spectrum captures the essence of the PAH emission spectrum and is highly correlated with observations of AIBs, strongly supporting PAHs as their source. Also, the fact that a large number of molecules are required implies that spectroscopic signatures of the individual PAHs contributing to the AIBs spanning the visible, near-infrared, and far infrared spectral regions are weak, explaining why they have not yet been detected. An improved effort, joining laboratory, theoretical, and observational studies of the PAH emission process, will support the use of PAH features as a probe of physical and chemical conditions in the nearby and distant Universe.

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.

We report a Herschel detection of high-J rotational CO lines from a dense knot in the supernova remnant Cas A. Based on a combined analysis of these rotational lines, and previously observed ro-vibrational CO lines, we find the gas to be warm (two components at 400 and 2000 K) and dense (1e6-7 cm-3), with a CO column density of 5e17 cm-2. This, along with the broad line widths (400 kms-1), suggests that the CO emission originates in the post-shock region of the reverse shock. As the passage of the reverse shock dissociates any existing molecules, the CO has most likely reformed in the last few years, in the post-shock gas. The CO cooling time is comparable to the CO formation time, so possible heating sources (UV photons from the shock front, X-rays, electron conduction) to maintain the large column density of warm CO are discussed.

In the early 90s, C60+ was proposed as the carrier of two diffuse interstellar bands (DIBs) at 957.7 and 963.2 nm, but a firm identification still awaits gas-phase spectroscopic data. Neutral C60, on the other hand, was recently detected through its infrared emission bands in the interstellar medium and evolved stars. In this contribution, we present the detection of C60+ through its infrared vibrational bands in the NGC 7023 nebula, based on spectroscopic observations with the Spitzer space telescope, quantum chemistry calculation, and laboratory data from the literature. This detection supports the idea that C60+ could be a DIB carrier, and provides robust evidence that fullerenes exist in the gas-phase in the interstellar medium. Modeling efforts to design specific observations, combined with new gas-phase data, will be essential to confirm this proposal. A definitive attribution of the 957.7 and 963.2 nm DIBs to C60+ would represent a significant step forward in the field.

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.

We observationally investigate the relation between the photoelectric heating efficiency in PDRs and the charge of PAHs, which are considered to play a key role in photoelectric heating. Using PACS onboard Herschel, we observed six PDRs spanning a wide range of FUV radiation fields (G_0=100-10^5). To measure the photoelectric heating efficiency, we obtained the intensities of the main cooling lines, i.e., the [OI]63um, 145um, and [CII]158um, as well as the FIR continuum intensity. We used Spitzer/IRS spectroscopic mapping observations to investigate the MIR PAH features in the same regions. We decomposed the MIR PAH emission into that of neutral (PAH^0) and positively ionized (PAH^+) species to derive the fraction of the positively charged PAHs, and compare it to the photoelectric heating efficiency. The heating efficiency traced by ([OI]63um+[OI]145um+[CII]158um) / TIR, ranges between 0.1% and 0.9% in different sources, and the fraction of PAH^+ relative to (PAH^0 + PAH^+) spans from 0(+11)% to 87(+/-10)%. All positions with a high PAH^+ fraction show a low heating efficiency, and all positions with a high heating efficiency have a low PAH^+ fraction, supporting the scenario in which a positive grain charge results in a decreased heating efficiency. Theoretical estimates of the photoelectric heating efficiency show a stronger dependence on the charging parameter gamma=G_0 T^{1/2}/n_e than the observed efficiency reported in this study, and the discrepancy is significant at low gamma. The photoelectric heating efficiency on PAHs, traced by ([OI]63um+[OI]145um+[CII]158um) / (PAH+[OI]63um+[OI]145um+[CII]158um), shows a much better match between the observations and the theoretical estimates. The good agreement of the photoelectric heating efficiency on PAHs with a theoretical model indicates the dominant contribution of PAHs to the photoelectric heating. (abridged for arXiv)

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.

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.

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.

Buckminsterfullerene (C60) was recently confirmed to be the largest molecule identified in space. However, it remains unclear how, and where this molecule is formed. It is generally believed that C60 is formed from the build up of small carbonaceous compounds, in the hot and dense envelopes of evolved stars. Analyzing infrared observations, obtained by Spitzer and Herschel, we found that C60 is efficiently formed in the tenuous and cold environment of an interstellar cloud illuminated by strong ultraviolet (UV) radiation fields. This implies that another formation pathway, efficient at low densities, must exist. Based on recent laboratory and theoretical studies, we argue that Polycyclic Aromatic Hydrocarbons are converted into graphene, and subsequently C60, under UV irradiation from massive stars. This shows that alternative - top-down - routes are key to understanding the organic inventory in space.

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.

Massive stars influence their parental molecular cloud, and it has long been suspected that the development of hydrodynamical instabilities can compress or fragment the cloud. Identifying such instabilities has proved difficult. It has been suggested that elongated structures (such as the `pillars of creation') and other shapes arise because of instabilities, but alternative explanations are available. One key signature of an instability is a wave-like structure in the gas, which has hitherto not been seen. Here we report the presence of `waves' at the surface of the Orion molecular cloud near where massive stars are forming. The waves seem to be a Kelvin-Helmholtz instability that arises during the expansion of the nebula as gas heated and ionized by massive stars is blown over pre-existing molecular gas.

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.