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

Multi-line Observations, Models, and Data Needed to Understand the Nature of UV-irradiated Interstellar Matter

Far-ultraviolet photons from OB-type massive stars regulate the heating, ionization, and chemistry of much of the neutral interstellar gas in star-forming galaxies. The interaction of FUV radiation and interstellar matter takes place in environments broadly known as photodissociation regions (PDRs). PDR line diagnostics are the smoking gun of the radiative feedback from massive stars. Improving our understanding of stellar feedback in the ISM requires quantifying the energy budget, gas dynamics, and chemical composition of PDR environments. This goal demands astronomical instrumentation able to deliver multi-line spectroscopic images of the ISM (of the Milky Way and nearby galaxies). It also requires interdisciplinary collaborations to obtain the rate coefficients and cross sections of the many microphysical processes that occur in the ISM and that are included in models such as the Meudon PDR code.

Bringing high spatial resolution to the far-infrared: A giant leap for astrophysics

The far-infrared (FIR) regime is one of the wavelength ranges where no astronomical data with sub-arcsecond spatial resolution exist. None of the medium-term satellite projects like SPICA, Millimetron, or the Origins Space Telescope will resolve this malady. For many research areas, however, information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, from highly excited carbon monoxide (CO), light hydrides, and especially from water lines would open the door for transformative science. A main theme will be to trace the role of water in proto-planetary discs, to observationally advance our understanding of the planet formation process and, intimately related to that, the pathways to habitable planets and the emergence of life. Furthermore, key observations will zoom into the physics and chemistry of the star-formation process in our own Galaxy, as well as in external galaxies. The FIR provides unique tools to investigate in particular the energetics of heating, cooling, and shocks. The velocity-resolved data in these tracers will reveal the detailed dynamics engrained in these processes in a spatially resolved fashion, and will deliver the perfect synergy with ground-based molecular line data for the colder dense gas.

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.

Direct estimation of electron density in the Orion Bar PDR from mm-wave carbon recombination lines

CONTEXT

A significant fraction of the molecular gas in star-forming regions is irradiated by stellar UV photons. In these environments, the electron density (n e) plays a critical role in the gas dynamics, chemistry, and collisional excitation of certain molecules.

AIMS

We determine n e in the prototypical strongly irradiated photodissociation region (PDR), the Orion Bar, from the detection of new millimeter-wave carbon recombination lines (mmCRLs) and existing far-IR [13Cii] hyperfine line observations.

METHODS

We detect 12 mmCRLs (including α, β, and γ transitions) observed with the IRAM 30m telescope, at ~ 25″ angular resolution, toward the H/H2 dissociation front (DF) of the Bar. We also present a mmCRL emission cut across the PDR.

RESULTS

These lines trace the C+/C/CO gas transition layer. As the much lower frequency carbon radio recombination lines, mmCRLs arise from neutral PDR gas and not from ionized gas in the adjacent Hii region. This is readily seen from their narrow line profiles (Δv = 2.6 ± 0.4 km s-1) and line peak velocities (ν LSR = +10.7 ± 0.2 km s-1). Optically thin [13Cii] hyperfine lines and molecular lines - emitted close to the DF by trace species such as reactive ions CO+ and HOC+ - show the same line profiles. We use non-LTE excitation models of [13Cii] and mmCRLs and derive n e = 60 - 100 cm-3 and T e = 500 - 600 K toward the DF.

CONCLUSIONS

The inferred electron densities are high, up to an order of magnitude higher than previously thought. They provide a lower limit to the gas thermal pressure at the PDR edge without using molecular tracers. We obtain P th ≥ (2 - 4)·108 cm-3 K assuming that the electron abundance is equal to or lower than the gas-phase elemental abundance of carbon. Such elevated thermal pressures leave little room for magnetic pressure support and agree with a scenario in which the PDR photoevaporates.

Molecular tracers of radiative feedback in Orion (OMC-1) Widespread CH+ (J = 1-0), CO (10-9), HCN (6-5), and HCO+ (6-5) emission

Young massive stars regulate the physical conditions, ionization, and fate of their natal molecular cloud and surroundings. It is important to find tracers that help quantifying the stellar feedback processes that take place at different spatial scales. We present ~85 arcmin2 (~1.3 pc2) velocity-resolved maps of several submillimeter molecular lines, taken with Herschel/HIFI, toward the closest high-mass star-forming region, the Orion molecular cloud 1 core (OMC-1). The observed rotational lines include probes of warm and dense molecular gas that are difficult, if not impossible, to detect from ground-based telescopes: CH+ (J = 1-0), CO (J = 10-9), HCO+ (J = 6-5) and HCN (J = 6-5), and CH (N, J =1, 3/2-1, 1/2). These lines trace an extended but thin layer (A V ≃3-6 mag or ~1016 cm) of molecular gas at high thermal pressure, P th = n H · T k ≈ 107 - 109 cm-3 K, associated with the far ultraviolet (FUV) irradiated surface of OMC-1. The intense FUV radiation field, emerging from massive stars in the Trapezium cluster, heats, compresses and photoevaporates the cloud edge. It also triggers the formation of specific reactive molecules such as CH+. We find that the CH+ (J = 1-0) emission spatially correlates with the flux of FUV photons impinging the cloud: G 0 from ~103 to ~105. This correlation is supported by constant-pressure photodissociation region (PDR) models in the parameter space P th/G 0 ≈ [5 · 103 - 8 · 104] cm-3 K where many observed PDRs seem to lie. The CH+ (J = 1-0) emission spatially correlates with the extended infrared emission from vibrationally excited H2 (v ≥ 1), and with that of [C ii] 158 μm and CO J = 10-9, all emerging from FUV-irradiated gas. These correlations link the presence of CH+ to the availability of C+ ions and of FUV-pumped H2 (v ≥ 1) molecules. We conclude that the parsec-scale CH+ emission and narrow-line (Δv ≃ 3 km s-1) mid-J CO emission arises from extended PDR gas and not from fast shocks. PDR line tracers are the smoking gun of the stellar feedback from young massive stars. The PDR cloud surface component in OMC-1, with a mass density of 120-240 M ⊙ pc-2, represents ~5% to ~10% of the total gas mass, however, it dominates the emitted line luminosity; the average CO J = 10-9 surface luminosity in the mapped region being ~35 times brighter than that of CO J = 2-1. These results provide insights into the source of submillimeter CH+ and mid-J CO emission from distant star-forming galaxies.

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 stars^1,2). Stellar winds, supernova explosions and ionization by ultraviolet photons control the lifetimes of molecular clouds^3-7. Theoretical studies predict that momentum injection by radiation should dominate that by stellar winds⁸, 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 θ¹ 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.

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.

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

Complex organic molecules in strongly UV-irradiated gas

We investigate the presence of complex organic molecules (COMs) in strongly UV-irradiated interstellar molecular gas. We have carried out a complete millimetre (mm) line survey using the IRAM 30 m telescope towards the edge of the Orion Bar photodissociation region (PDR), close to the H2 dissociation front, a position irradiated by a very intense far-UV (FUV) radiation field. These observations have been complemented with 8.5″ resolution maps of the H2CO JKa,Kc = 51,5 → 41,4 and C18O J = 3 → 2 emission at 0.9 mm. Despite being a harsh environment, we detect more than 250 lines from COMs and related precursors: H2CO, CH3OH, HCO, H2CCO, CH3CHO, H2CS, HCOOH, CH3CN, CH2NH, HNCO, [Formula: see text] and HC3N (in decreasing order of abundance). For each species, the large number of detected lines allowed us to accurately constrain their rotational temperatures (Trot) and column densities (N). Owing to subthermal excitation and intricate spectroscopy of some COMs (symmetric- and asymmetric-top molecules such as CH3CN and H2CO, respectively), a correct determination of N and Trot requires building rotational population diagrams of their rotational ladders separately. The inferred column densities are in the 1011 - 1013cm-2 range. We also provide accurate upper limit abundances for chemically related molecules that might have been expected, but are not conclusively detected at the edge of the PDR (HDCO, CH3O, CH3NC, CH3CCH, CH3OCH3, HCOOCH3, CH3CH2OH, CH3CH2CN, and CH2CHCN). A non-thermodynamic equilibrium excitation analysis for molecules with known collisional rate coefficients suggests that some COMs arise from different PDR layers but we cannot resolve them spatially. In particular, H2CO and CH3CN survive in the extended gas directly exposed to the strong FUV flux (Tk = 150 - 250 K and Td ≳ 60 K), whereas CH3OH only arises from denser and cooler gas clumps in the more shielded PDR interior (Tk = 40 - 50 K). The non-detection of HDCO towards the PDR edge is consistent with the minor role of pure gas-phase deuteration at very high temperatures. We find a HCO/H2CO/CH3OH ≃ 1/5/3 abundance ratio. These ratios are different from those inferred in hot cores and shocks. Taking into account the elevated gas and dust temperatures at the edge of the Bar (mostly mantle-free grains), we suggest the following scenarios for the formation of COMs: (i) hot gas-phase reactions not included in current models; (ii) warm grain-surface chemistry; or (iii) the PDR dynamics is such that COMs or precursors formed in cold icy grains deeper inside the molecular cloud desorb and advect into the PDR.

Spatial distribution of FIR rotationally excited CH+ and OH emission lines in the Orion Bar PDR

CONTEXT

The methylidyne cation (CH⁺) and hydroxyl (OH) are key molecules in the warm interstellar chemistry, but their formation and excitation mechanisms are not well understood. Their abundance and excitation are predicted to be enhanced by the presence of vibrationally excited H₂ or hot gas (~500-1000 K) in photodissociation regions with high incident FUV radiation field. The excitation may also originate in dense gas (> 10⁵ cm^-3) followed by nonreactive collisions with H₂, H, and electrons. Previous observations of the Orion Bar suggest that the rotationally excited CH⁺ and OH correlate with the excited CO, a tracer of dense and warm gas, and formation pumping contributes to CH⁺ excitation.

AIMS

Our goal is to examine the spatial distribution of the rotationally excited CH⁺ and OH emission lines in the Orion Bar in order to establish their physical origin and main formation and excitation mechanisms.

METHODS

We present spatially sampled maps of the CH⁺ J=3-2 transition at 119.8 µm and the OH Λ-doublet at 84 µm in the Orion Bar over an area of 110″×110″ with Herschel (PACS). We compare the spatial distribution of these molecules with those of their chemical precursors, C⁺, O and H₂, and tracers of warm and dense gas (high-J CO). We assess the spatial variation of CH⁺ J=2-1 velocity-resolved line profile at 1669 GHz with Herschel HIFI spectrometer observations.

RESULTS

The OH and especially CH⁺ lines correlate well with the high-J CO emission and delineate the warm and dense molecular region at the edge of the Bar. While notably similar, the differences in the CH⁺ and OH morphologies indicate that CH⁺ formation and excitation are strongly related to the observed vibrationally excited H₂. This, together with the observed broad CH⁺ line widths, indicates that formation pumping contributes to the excitation of this reactive molecular ion. Interestingly, the peak of the rotationally excited OH 84 µm emission coincides with a bright young object, proplyd 244-440, which shows that OH can be an excellent tracer of UV-irradiated dense gas.

CONCLUSIONS

The spatial distribution of CH⁺ and OH revealed in our maps is consistent with previous modeling studies. Both formation pumping and nonreactive collisions in a UV-irradiated dense gas are important CH⁺ J=3-2 excitation processes. The excitation of the OH Λ-doublet at 84 µm is mainly sensitive to the temperature and density.

Analysis of the Herschel/HEXOS Spectral Survey Towards Orion South: A massive protostellar envelope with strong external irradiation

We present results from a comprehensive submillimeter spectral survey toward the source Orion South, based on data obtained with the HIFI instrument aboard the Herschel Space Observatory, covering the frequency range 480 to 1900 GHz. We detect 685 spectral lines with S/N > 3σ, originating from 52 different molecular and atomic species. We model each of the detected species assuming conditions of Local Thermodynamic Equilibrium. This analysis provides an estimate of the physical conditions of Orion South (column density, temperature, source size, & V LSR ). We find evidence for three different cloud components: a cool (T ex ~ 20 - 40 K), spatially extended (> 60″), and quiescent (ΔVFWHM ~ 4 km s -1) component; a warmer (T ex ~ 80 - 100 K), less spatially extended (~ 30″), and dynamic (ΔVFWHM ~ 8 km s -1) component, which is likely affected by embedded outflows; and a kinematically distinct region (T ex > 100 K; V LSR ~ 8 km s -1), dominated by emission from species which trace ultraviolet irradiation, likely at the surface of the cloud. We find little evidence for the existence of a chemically distinct "hot core" component, likely due to the small filling factor of the hot core or hot cores within the Herschel beam. We find that the chemical composition of the gas in the cooler, quiescent component of Orion South more closely resembles that of the quiescent ridge in Orion-KL. The gas in the warmer, dynamic component, however, more closely resembles that of the Compact Ridge and Plateau regions of Orion-KL, suggesting that higher temperatures and shocks also have an influence on the overall chemistry of Orion South.

Compression and ablation of the photo-irradiated molecular cloud the Orion Bar

The Orion Bar is the archetypal edge-on molecular cloud surface illuminated by strong ultraviolet radiation from nearby massive stars. Owing to the close distance to Orion (about 1,350 light-year), the effects of stellar feedback on the parental cloud can be studied in detail. Visible-light observations of the Bar1 show that the transition between the hot ionised gas and the warm neutral atomic gas (the ionisation front) is spatially well separated from the transition from atomic to molecular gas (the dissociation front): about 15 arcseconds or 6,200 astronomical units. (One astronomical unit is the Earth-Sun distance.) Static equilibrium models2,3 used to interpret previous far-infrared and radio observations of the neutral gas in the Bar4,5,6 (typically at 10-20 arcsecond resolution) predict an inhomogeneous cloud structure consisting of dense clumps embedded in a lower density extended gas component. Here we report 1 arcsecond resolution millimetre-wave images that allow us to resolve the molecular cloud surface and constrain the gas density and temperature structures at small spatial scales. In contrast to stationary model predictions7,8,9, there is no appreciable offset between the peak of the H₂ vibrational emission (delineating the H/H₂ transition) and the edge of the observed CO and HCO⁺ emission. This implies that the H/H₂ and C⁺/C/CO transition zones are very close. These observations reveal a fragmented ridge of high-density substructures, photo-ablative gas flows and instabilities at the molecular cloud surface. They suggest that the cloud edge has been compressed by a high-pressure wave that currently moves into the molecular cloud. The images demonstrate that dynamical and nonequilibrium effects are important. Thus, they should be included in any realistic description of irradiated interstellar matter.