Expanding bubbles in Orion A: [C II] observations of M42, M43, and NGC 1977.

CONTEXT: 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. Stellar feedback is one of the most crucial cosmological parameters that determine the properties and evolution of the interstellar medium in galaxies. AIMS: We aim to understand the role that feedback by stellar winds and radiation play in the evolution of the interstellar medium. Velocity-resolved observations of the [C II] 158µm fine-structure line allow us to study the kinematics of UV-illuminated gas. Here, we present a square-degree-sized map of [C II] emission from the Orion Nebula complex at a spatial resolution of 16″ and high spectral resolution of 0.2kms-1, covering the entire Orion Nebula (M42) plus M43 and the nebulae NGC 1973, 1975, and 1977 to the north. We compare the stellar characteristics of these three regions with the kinematics of the expanding bubbles surrounding them. METHODS: We use [C II] 158µm line observations over an area of 1.2deg2 in the Orion Nebula complex obtained by the upGREAT instrument onboard SOFIA. RESULTS: The bubble blown by the O7V star θ 1 Ori C in the Orion Nebula expands rapidly, at 13kms-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 Myr. M43 with the B0.5V star NU Ori also exhibits an expanding bubble structure, with an expansion velocity of 6kms-1. Comparison with analytical models for the pressure-driven expansion of H II regions gives an age estimate of 0.02 Myr. The bubble surrounding NGC 1973, 1975, and 1977 with the central B1V star 42 Orionis expands at 1.5kms-1, likely due to the over-pressurized ionized gas as in the case of M43. We derive an age of 0.4 Myr for this structure. CONCLUSIONS: We conclude that the bubble of the Orion Nebula is driven by the mechanical energy input by the strong stellar wind from θ 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.

Disruption of the Orion molecular core 1 by wind from the massive star θ1 Orionis C.

Massive stars inject mechanical and radiative energy into the surrounding environment, which stirs it up, heats the gas, produces cloud and intercloud phases in the interstellar medium, and disrupts molecular clouds (the birth sites of new stars1,2). Stellar winds, supernova explosions and ionization by ultraviolet photons control the lifetimes of molecular clouds3-7. Theoretical studies predict that momentum injection by radiation should dominate that by stellar winds8, but this has been difficult to assess observationally. Velocity-resolved large-scale images in the fine-structure line of ionized carbon ([C II]) provide an observational diagnostic for the radiative energy input and the dynamics of the interstellar medium around massive stars. Here we report observations of a one-square-degree region (about 7 parsecs in diameter) of Orion molecular core 1-the region nearest to Earth that exhibits massive-star formation-at a resolution of 16 arcseconds (0.03 parsecs) in the [C II] line at 1.9 terahertz (158 micrometres). The results reveal that the stellar wind originating from the massive star θ1 Orionis C has swept up the surrounding material to create a 'bubble' roughly four parsecs in diameter with a 2,600-solar-mass shell, which is expanding at 13 kilometres per second. This finding demonstrates that the mechanical energy from the stellar wind is converted very efficiently into kinetic energy of the shell and causes more disruption of the Orion molecular core 1 than do photo-ionization and evaporation or future supernova explosions.

Learning mid-IR emission spectra of polycyclic aromatic hydrocarbon populations from observations.

CONTEXT: The James Webb Space Telescope (JWST) will deliver an unprecedented quantity of high-quality spectral data over the 0.6-28 µm range. It will combine sensitivity, spectral resolution, and spatial resolution. Specific tools are required to provide efficient scientific analysis of such large data sets. AIMS: 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 polycyclic aromatic hydrocarbon (PAH) species and carbonaceous very small grains (VSGs). METHODS: We present a method based on linear fitting and blind signal separation 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. RESULTS: 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 µm range, which is considered here for the first time in a blind signal separation (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. CONCLUSIONS: 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.

Structure of photodissociation fronts in star-forming regions revealed by observations of high-J CO emission lines with Herschel.

CONTEXT: In bright photodissociation regions (PDRs) associated to massive star formation, the presence of dense "clumps" that are immersed in a less dense interclump medium is often proposed to explain the difficulty of models to account for the observed gas emission in high-excitation lines. AIMS: We aim at presenting a comprehensive view of the modeling of the CO rotational ladder in PDRs, including the high-J lines that trace warm molecular gas at PDR interfaces. METHODS: We observed the 12CO and 13CO ladders in two prototypical PDRs, the Orion Bar and NGC 7023 NW using the instruments onboard Herschel. We also considered line emission from key species in the gas cooling of PDRs (C+, O, H2) and other tracers of PDR edges such as OH and CH+. All the intensities are collected from Herschel observations, the literature and the Spitzer archive and are analyzed using the Meudon PDR code. RESULTS: A grid of models was run to explore the parameter space of only two parameters: thermal gas pressure and a global scaling factor that corrects for approximations in the assumed geometry. We conclude that the emission in the high-J CO lines, which were observed up to J up =23 in the Orion Bar (J up =19 in NGC 7023), can only originate from small structures of typical thickness of a few 10-3 pc and at high thermal pressures (Pth ~ 108 K cm-3). CONCLUSIONS: Compiling data from the literature, we found that the gas thermal pressure increases with the intensity of the UV radiation field given by G0, following a trend in line with recent simulations of the photoevaporation of illuminated edges of molecular clouds. This relation can help rationalising the analysis of high-J CO emission in massive star formation and provides an observational constraint for models that study stellar feedback on molecular clouds.

Herschel survey and modelling of externally-illuminated photoevaporating protoplanetary disks.

CONTEXT: 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. AIMS: We aim at studying the processes involved in disk photoevaporation when it is driven by far-UV photons (i.e. 6 < E < 13.6 eV). METHODS: 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 we computed the far infrared line emission. RESULTS: With this model, we successfully reproduce most of the observations and derive key physical parameters, i.e. densities at the disk surface of about 106 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 specific energetic regime is associated to an equilibrium dynamical point of the photoevaporation flow: the mass-loss rate is self-regulated to keep the envelope column density at a value that maintains the temperature at the disk surface around 1000 K. From the physical parameters derived from our best-fit models, we estimate mass-loss rates - of the order of 10-7 M⊙/yr - that are in agreement with earlier spectroscopic observation of ionised gas tracers. This holds only if we assume photoevaporation in the supercritical regime where the evaporation flow is launched from the disk surface at sound speed. CONCLUSIONS: 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.

[Cii] emission from L1630 in the Orion B molecular cloud.

CONTEXT: L1630 in the Orion B molecular cloud, which includes the iconic Horsehead Nebula, illuminated by the star system σ Ori, is an example of a photodissociation region (PDR). In PDRs, stellar radiation impinges on the surface of dense material, often a molecular cloud, thereby inducing a complex network of chemical reactions and physical processes. AIMS: Observations toward L1630 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. Our goal is to relate the [Cii] fine-structure line emission to the physical conditions predominant in L1630 and compare it to studies of OMC1. METHODS: The [Cii] 158 µm line emission of L1630 around the Horsehead Nebula, an area of 12' × 17', was observed using the upgraded German Receiver for Astronomy at Terahertz Frequencies (upGREAT) onboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). RESULTS: Of the [Cii] emission from the mapped area 95%, 13 L⊙, originates from the molecular cloud; the adjacent Hii region contributes only 5%, that is, 1 L⊙. 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 nH ∼ 3 · 103 cm-3, with surface layers, including the Horsehead Nebula, having a density of up to nH ∼ 4 · 104 cm-3. The temperature of the surface gas is T ∼ 100 K. The average [Cii] cooling efficiency within the molecular cloud is 1.3 · 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. CONCLUSIONS: In L1630 only a small fraction of the gas mass is traced by [Cii]. Most of the [Cii] emission in the mapped area stems from PDR surfaces. The layered edge-on structure of the molecular cloud and limitations in spatial resolution put constraints on our ability to relate different tracers to each other and to the physical conditions. From our study, we conclude that the relation between [Cii] emission and physical conditions is likely to be more complicated than often assumed. The theoretical heating efficiency is higher than the one we calculate from the observed [Cii] emission in the L1630 molecular cloud.

Detection of Buckminsterfullerene emission in the diffuse interstellar medium.

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 Å to electronic transitions of the buckminsterfullerene cation (i.e. [Formula: see text]) 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 µm emission bands commonly attributed to vibrational bands of neutral C60. According to classical models that compute the charge state of large molecules in space, C60 is expected to be mostly neutral in the diffuse ISM. This is in agreement with the abundances of diffuse C60 we derive here from observations. We also find that C60 is less abundant in the diffuse ISM than in star-forming regions, supporting the theory that C60 can be formed in these regions.

The first CO+ image: I. Probing the HI/H2 layer around the ultracompact HII region Mon R2.

The CO+ reactive ion is thought to be a tracer of the boundary between a HII region and the hot molecular gas. In this study, we present the spatial distribution of the CO+ rotational emission toward the Mon R2 star-forming region. The CO+ emission presents a clumpy ring-like morphology, arising from a narrow dense layer around the HII region. We compare the CO+ distribution with other species present in photon-dominated regions (PDR), such as [CII] 158 µm, H2 S(3) rotational line at 9.3 µm, polycyclic aromatic hydrocarbons (PAHs) and HCO+. We find that the CO+ emission is spatially coincident with the PAHs and [CII] emission. This confirms that the CO+ emission arises from a narrow dense layer of the HI/H2 interface. We have determined the CO+ fractional abundance, relative to C+ toward three positions. The abundances range from 0.1 to 1.9 ×10-10 and are in good agreement with previous chemical model, which predicts that the production of CO+ in PDRs only occurs in dense regions with high UV fields. The CO+ linewidth is larger than those found in molecular gas tracers, and their central velocity are blue-shifted with respect to the molecular gas velocity. We interpret this as a hint that the CO+ is probing photo-evaporating clump surfaces.

High spatial resolution imaging of SO and H2CO in AB Auriga: The first SO image in a transitional disk.

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 NOEMA (NOrthern Extended Millimeter Array) interferometer to image the 1mm continuum dust emission and the 13CO J=2 →1, C18OJ=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.

Polycyclic aromatic hydrocarbons and molecular hydrogen in oxygen-rich planetary nebulae: the case of NGC 6720.

Evolved stars are primary sources for the formation of polycyclic aromatic hydrocarbons (PAHs) and dust grains. Their circumstellar chemistry is usually designated as either oxygen-rich or carbon-rich, although dual-dust chemistry objects, whose infrared spectra reveal both silicate- and carbon-dust features, are also known. The exact origin and nature of this dual-dust chemistry is not yet understood. Spitzer-IRS mid-infrared spectroscopic imaging of the nearby, oxygen-rich planetary nebula NGC 6720 reveals the presence of the 11.3 µm aromatic (PAH) emission band. It is attributed to emission from neutral PAHs, since no band is observed in the 7-8 µm range. The spatial distribution of PAHs is found to closely follow that of the warm clumpy molecular hydrogen emission. Emission from both neutral PAHs and warm H2 is likely to arise from photo-dissociation regions associated with dense knots that are located within the main ring. The presence of PAHs together with the previously derived high abundance of free carbon (relative to CO) suggest that the local conditions in an oxygen-rich environment can also become conducive to in-situ formation of large carbonaceous molecules, such as PAHs, via a bottom-up chemical pathway. In this scenario, the same stellar source can enrich the interstellar medium with both oxygen-rich dust and large carbonaceous molecules.

Top-down formation of fullerenes in the interstellar medium.

Fullerenes have been recently detected in various circumstellar and interstellar environments, raising the question of their formation pathway. It has been proposed that they can form at the low densities found in the interstellar medium by the photo-chemical processing of large polycyclic aromatic hydrocarbons (PAHs). Following our previous work on the evolution of PAHs in the NGC 7023 reflection nebula, we evaluate, using photochemical modeling, the possibility that the PAH C66H20 (i.e. circumovalene) can lead to the formation of C60 upon irradiation by ultraviolet photons. The chemical pathway involves full dehydrogenation of C66H20, folding into a floppy closed cage and shrinking of the cage by loss of C2 units until it reaches the symmetric C60 molecule. At 10" from the illuminating star and with realistic molecular parameters, the model predicts that 100% of C66H20 is converted into C60 in ~ 105 years, a timescale comparable to the age of the nebula. Shrinking appears to be the kinetically limiting step of the whole process. Hence, PAHs larger than C66H20 are unlikely to contribute significantly to the formation of C60, while PAHs containing between 60 and 66 C atoms should contribute to the formation of C60 with shorter timescales, and PAHs containing less than 60 C atoms will be destroyed. Assuming a classical size distribution for the PAH precursors, our model predicts absolute abundances of C60 are up to several 10-4 of the elemental carbon, i.e. less than a percent of the typical interstellar PAH abundance, which is consistent with observational studies. According to our model, once formed, C60 can survive much longer (> 107 years for radiation fields below G0 = 104) than other fullerenes because of the remarkable stability of the C60 molecule at high internal energies. Hence, a natural consequence is that C60 is more abundant than other fullerenes in highly irradiated environments.