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.

Time-dependent molecular emission in IRC+10216.

CONTEXT: The variability in IRC+10216, the envelope of the asymptotic giant branch (AGB) star CW Leo, has attracted increasing attention in recent years. Studying the details of this variability in the molecular emission required a systematic observation program. AIMS: We aim to reveal and characterize the periodical variability of the rotational lines from several molecules and radicals in IRC+10216, and to compare it with previously reported IR variability. METHODS: We carried out systematic monitoring within the ~80 to 116 GHz frequency range with the IRAM 30m telescope. RESULTS: We report on the periodical variability in IRC+10216 of several rotational lines from the following molecules and radicals: HC3N, HC5N, CCH, C4H, C5H, and CN. The analysis of the variable molecular lines provides periods that are consistent with previously reported IR variability, and interesting phase lags are revealed that point toward radiative transfer and pumping, rather than chemical effects. CONCLUSIONS: This study indicates that observations of several lines of a given molecule have to be performed simultaneously or at least at the same phase in order to avoid erroneous interpretation of the data. In particular, merging ALMA data from different epochs may prove to be difficult, as shown by the example of the variability we studied here. Moreover, radiative transfer codes have to incorporate the effect of population variability in the rotational levels in CW Leo.

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.

High-velocity hot CO emission close to Sgr A*: Herschel/HIFI☠, ☠☠submillimeter spectral survey toward Sgr A.

The properties of molecular gas, the fuel that forms stars, inside the cavity of the circumnuclear disk (CND) are not well constrained. We present results of a velocity-resolved submillimeter scan (~480 to 1250 GHz) and [C ii] 158 µm line observations carried out with Herschel/HIFI toward Sgr A*; these results are complemented by a ~2'×2' 12CO (J=3-2) map taken with the IRAM 30 m telescope at ~7″ resolution. We report the presence of high positive-velocity emission (up to about +300 km s-1) detected in the wings of 12CO J=5-4 to 10-9 lines. This wing component is also seen in H2O (11,0-10,1), a tracer of hot molecular gas; in [C ii]158 µm, an unambiguous tracer of UV radiation; but not in [C i] 492, 806 GHz. This first measurement of the high-velocity 12CO rotational ladder toward Sgr A* adds more evidence that hot molecular gas exists inside the cavity of the CND, relatively close to the supermassive black hole (< 1 pc). Observed by ALMA, this velocity range appears as a collection of 12CO (J=3-2) cloudlets lying in a very harsh environment that is pervaded by intense UV radiation fields, shocks, and affected by strong gravitational shears. We constrain the physical conditions of the high positive-velocity CO gas component by comparing with non-LTE excitation and radiative transfer models. We infer T k≃400 K to 2000 K for n H≃(0.2-1.0)·105 cm-3. These results point toward the important role of stellar UV radiation, but we show that radiative heating alone cannot explain the excitation of this ~10-60 M ⊙ component of hot molecular gas inside the central cavity. Instead, strongly irradiated shocks are promising candidates.

Using radio astronomical receivers for molecular spectroscopic characterization in astrochemical laboratory simulations: A proof of concept.

We present a proof of concept on the coupling of radio astronomical receivers and spectrometers with chemical reactors and the performances of the resulting setup for spectroscopy and chemical simulations in laboratory astrophysics. Several experiments including cold plasma generation and UV photochemistry were performed in a 40 cm long gas cell placed in the beam path of the Aries 40 m radio telescope receivers operating in the 41-49 GHz frequency range interfaced with fast Fourier transform spectrometers providing 2 GHz bandwidth and 38 kHz resolution. The impedance matching of the cell windows has been studied using different materials. The choice of the material and its thickness was critical to obtain a sensitivity identical to that of standard radio astronomical observations. Spectroscopic signals arising from very low partial pressures of CH3OH, CH3CH2OH, HCOOH, OCS, CS, SO2 (<10-3 mbar) were detected in a few seconds. Fast data acquisition was achieved allowing for kinetic measurements in fragmentation experiments using electron impact or UV irradiation. Time evolution of chemical reactions involving OCS, O2 and CS2 was also observed demonstrating that reactive species, such as CS, can be maintained with high abundance in the gas phase during these experiments.

[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.

Herschel/HIFI observations of the circumstellar ammonia lines in IRC+10216.

CONTEXT: A discrepancy exists between the abundance of ammonia (NH3) derived previously for the circumstellar envelope (CSE) of IRC+10216 from far-IR submillimeter rotational lines and that inferred from radio inversion or mid-infrared (MIR) absorption transitions. AIMS: To address the discrepancy described above, new high-resolution far-infrared (FIR) observations of both ortho- and para-NH3 transitions toward IRC+10216 were obtained with Herschel, with the goal of determining the ammonia abundance and constraining the distribution of NH3 in the envelope of IRC+10216. METHODS: We used the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel to observe all rotational transitions up to the J = 3 level (three ortho- and six para-NH3 lines). We conducted non-LTE multilevel radiative transfer modelling, including the effects of near-infrared (NIR) radiative pumping through vibrational transitions. The computed emission line profiles are compared with the new HIFI data, the radio inversion transitions, and the MIR absorption lines in the ν2 band taken from the literature. RESULTS: We found that NIR pumping is of key importance for understanding the excitation of rotational levels of NH3. The derived NH3 abundances relative to molecular hydrogen were (2.8 ± 0.5) × 10-8 for ortho-NH3 and [Formula: see text] for para-NH3, consistent with an ortho/para ratio of 1. These values are in a rough agreement with abundances derived from the inversion transitions, as well as with the total abundance of NH3 inferred from the MIR absorption lines. To explain the observed rotational transitions, ammonia must be formed near to the central star at a radius close to the end of the wind acceleration region, but no larger than about 20 stellar radii (1σ confidence level).

VELOCITY-RESOLVED [C ii] EMISSION AND [C ii]/FIR MAPPING ALONG ORION WITH HERSCHEL.

We present the first ~7.5'×11.5' velocity-resolved (~0.2 km s-1) map of the [C ii] 158 µm line toward the Orion molecular cloud 1 (OMC 1) taken with the Herschel/HIFI instrument. In combination with far-infrared (FIR) photometric images and velocity-resolved maps of the H41α hydrogen recombination and CO J=2-1 lines, this data set provides an unprecedented view of the intricate small-scale kinematics of the ionized/PDR/molecular gas interfaces and of the radiative feedback from massive stars. The main contribution to the [C ii] luminosity (~85 %) is from the extended, FUV-illuminated face of the cloud (G0>500, nH>5×103 cm-3) and from dense PDRs (G≳104, nH≳105 cm-3) at the interface between OMC 1 and the H ii region surrounding the Trapezium cluster. Around ~15 % of the [C ii] emission arises from a different gas component without CO counterpart. The [C ii] excitation, PDR gas turbulence, line opacity (from [13C ii]) and role of the geometry of the illuminating stars with respect to the cloud are investigated. We construct maps of the L[C ii]/LFIR and LFIR/MGas ratios and show that L[C ii]/LFIR decreases from the extended cloud component (~10-2-10-3) to the more opaque star-forming cores (~10-3-10-4). The lowest values are reminiscent of the "[C ii] deficit" seen in local ultra-luminous IR galaxies hosting vigorous star formation. Spatial correlation analysis shows that the decreasing L[C ii]/LFIR ratio correlates better with the column density of dust through the molecular cloud than with LFIR/MGas. We conclude that the [C ii] emitting column relative to the total dust column along each line of sight is responsible for the observed L[C ii]/LFIR variations through the cloud.

Discovery of Time Variation of the Intensity of Molecular Lines in IRC+10216 in The Submillimeter and Far Infrared Domains.

We report on the discovery of strong intensity variations in the high rotational lines of abundant molecular species towards the archetypical circumstellar envelope of IRC+10216. The observations have been carried out with the HIFI instrument on board Herschel and with the IRAM 30-m telescope. They cover several observing periods spreading over 3 years. The line intensity variations for molecules produced in the external layers of the envelope most probably result from time variations in the infrared pumping rates. We analyze the main implications this discovery has on the interpretation of molecular line emission in the envelopes of Mira-type stars. Radiative transfer calculations have to take into account both the time variability of infrared pumping and the possible variation of the dust and gas temperatures with stellar phase in order to reproduce the observation of molecular lines at different epochs. The effect of gas temperature variations with stellar phase could be particularly important for lines produced in the innermost regions of the envelope. Each layer of the circumstellar envelope sees the stellar light radiation with a different lag time (phase). Our results show that this effect must be included in the models. The sub-mm and FIR lines of AGB stars cannot anymore be considered as safe intensity calibrators.