Quantum study of reaction O(3 P) + H2 (v, j) → OH + H: OH formation in strongly UV-irradiated gas.

The reaction between atomic oxygen and molecular hydrogen is an important one in astrochemistry as it regulates the abundance of the hydroxyl radical and serves to open the chemistry of oxygen in diverse astronomical environments. However, the existence of a high activation barrier in the reaction with ground state oxygen atoms limits its efficiency in cold gas. In this study we calculate the dependence of the reaction rate coefficient on the rotational and vibrational state of H2 and evaluate the impact on the abundance of OH in interstellar regions strongly irradiated by far-UV photons, where H2 can be efficiently pumped to excited vibrational states. We use a recently calculated potential energy surface and carry out time-independent quantum mechanical scattering calculations to compute rate coefficients for the reaction O(3 P) + H2 (v, j) → OH + H, with H2 in vibrational states v = 0-7 and rotational states j = 0-10. We find that the reaction becomes significantly faster with increasing vibrational quantum number of H2, although even for high vibrational states of H2 (v = 4-5) for which the reaction is barrierless, the rate coefficient does not strictly attain the collision limit and still maintains a positive dependence with temperature. We implemented the calculated state-specific rate coefficients in the Meudon PDR code to model the Orion Bar PDR and evaluate the impact on the abundance of the OH radical. We find the fractional abundance of OH is enhanced by up to one order of magnitude in regions of the cloud corresponding to A V = 1.3-2.3, compared to the use of a thermal rate coefficient for O + H2, although the impact on the column density of OH is modest, of about 60%. The calculated rate coefficients will be useful to model and interpret JWST observations of OH in strongly UV-illuminated environments.

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.

Gas phase Elemental abundances in Molecular cloudS (GEMS) II. On the quest for the sulphur reservoir in molecular clouds: the H2S case.

CONTEXT: Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium and the identity of the main sulphur reservoir is still an open question. AIMS: Our goal is to investigate the H2S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. METHODS: Using millimeter observations of CS, SO, H2S, and their isotopologues, we determine the physical conditions and H2S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H2S abundance. RESULTS: Our modeling shows that chemical desorption is the main source of gas-phase H2S in dark cores. The measured H2S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when n H > 2 × 104. This change in the desorption rate is consistent with the formation of thick H2O and CO ice mantles on grain surfaces. The observed SO and H2S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of 5 - 10. Along the three cores, atomic S is predicted to be the main sulphur reservoir. CONCLUSIONS: The gaseous H2S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H2S. The behavior of the observed H2S abundance suggests a changing desorption efficiency, which would probe the snowline in these cold cores. Our model, however, highly overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, we can only conclude that our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.

Dynamics of cluster-forming hub-filament systems: The case of the high-mass star-forming complex Monoceros R2.

CONTEXT: High-mass stars and star clusters commonly form within hub-filament systems. Monoceros R2 (hereafter Mon R2), at a distance of 830 pc, harbors one of the closest such systems, making it an excellent target for case studies. AIMS: We investigate the morphology, stability and dynamical properties of the Mon R2 hub-filament system. METHODS: We employ observations of the 13CO and C18O 1→0 and 2→1 lines obtained with the IRAM-30m telescope. We also use H2 column density maps derived from Herschel dust emission observations. RESULTS: We identified the filamentary network in Mon R2 with the DisPerSE algorithm and characterized the individual filaments as either main (converging into the hub) or secondary (converging to a main filament) filaments. The main filaments have line masses of 30-100 M ⊙ pc-1 and show signs of fragmentation, while the secondary filaments have line masses of 12-60 M ⊙ pc-1 and show fragmentation only sporadically. In the context of Ostriker's hydrostatic filament model, the main filaments are thermally supercritical. If non-thermal motions are included, most of them are trans-critical. Most of the secondary filaments are roughly transcritical regardless of whether non-thermal motions are included or not. From the morphology and kinematics of the main filaments, we estimate a mass accretion rate of 10-4-10-3 M ⊙ yr-1 into the central hub. The secondary filaments accrete into the main filaments with a rate of 0.1-0.4×10-4 M ⊙ yr-1. The main filaments extend into the central hub. Their velocity gradients increase towards the hub, suggesting acceleration of the gas.We estimate that with the observed infall velocity, the mass-doubling time of the hub is ~ 2:5 Myr, ten times larger than the free-fall time, suggesting a dynamically old region. These timescales are comparable with the chemical age of the Hii region. Inside the hub, the main filaments show a ring- or a spiral-like morphology that exhibits rotation and infall motions. One possible explanation for the morphology is that gas is falling into the central cluster following a spiral-like pattern.

Abundances of sulphur molecules in the Horsehead nebula First NS+ detection in a photodissociation region.

CONTEXT: Sulphur is one of the most abundant elements in the Universe (S/H∼1.3×10 -5 ) and plays a crucial role in biological systems on Earth. The understanding of its chemistry is therefore of major importance. AIMS: Our goal is to complete the inventory of S-bearing molecules and their abundances in the prototypical photodissociation region (PDR) the Horsehead nebula to gain insight into sulphur chemistry in UV irradiated regions. Based on the WHISPER (Wide-band High-resolution Iram-30m Surveys at two positions with Emir Receivers) millimeter (mm) line survey, our goal is to provide an improved and more accurate description of sulphur species and their abundances towards the core and PDR positions in the Horsehead. METHODS: The Monte Carlo Markov Chain (MCMC) methodology and the molecular excitation and radiative transfer code RADEX were used to explore the parameter space and determine physical conditions and beam-averaged molecular abundances. RESULTS: A total of 13 S-bearing species (CS, SO, SO2, OCS, H2CS - both ortho and para - HDCS, C2S, HCS+, SO+, H2S, S2H, NS and NS+) have been detected in the two targeted positions. This is the first detection of SO+ in the Horsehead and the first detection of NS+ in any PDR. We find a differentiated chemical behaviour between C-S and O-S bearing species within the nebula. The C-S bearing species C2S and o-H2CS present fractional abundances a factor of > two higher in the core than in the PDR. In contrast, the O-S bearing molecules SO, SO2, and OCS present similar abundances towards both positions. A few molecules, SO+, NS, and NS+, are more abundant towards the PDR than towards the core, and could be considered as PDR tracers. CONCLUSIONS: This is the first complete study of S-bearing species towards a PDR. Our study shows that CS, SO, and H2S are the most abundant S-bearing molecules in the PDR with abundances of ∼ a few 10-9. We recall that SH, SH+, S, and S+ are not observable at the wavelengths covered by the WHISPER survey. At the spatial scale of our observations, the total abundance of S atoms locked in the detected species is < 10-8, only ∼0.1% of the cosmic sulphur abundance.

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.

Gas phase Elemental abundances in Molecular cloudS (GEMS): I. The prototypical dark cloud TMC 1.

GEMS is an IRAM 30m Large Program whose aim is determining the elemental depletions and the ionization fraction in a set of prototypical star-forming regions. This paper presents the first results from the prototypical dark cloud TMC 1. Extensive millimeter observations have been carried out with the IRAM 30m telescope (3 mm and 2 mm) and the 40m Yebes telescope (1.3 cm and 7 mm) to determine the fractional abundances of CO, HCO+, HCN, CS, SO, HCS+, and N2H+ in three cuts which intersect the dense filament at the well-known positions TMC 1-CP, TMC 1-NH3, and TMC 1-C, covering a visual extinction range from A V ~ 3 to ~20 mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of 1-5×103 cm-3; and ii) the dense phase, located at A V > 10 mag, with molecular hydrogen densities >104 cm-3. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C+/C/CO transition zone (A V ~ 3 mag) where C/H ~ 8×10-5 and C/O~0.8-1, until the beginning of the dense phase at A V ~ 10 mag. This is consistent with the grain temperatures being below the CO evaporation temperature in this region. In the case of sulfur, a strong depletion should occur before the translucent phase where we estimate a S/H ~ (0.4 - 2.2) ×10-6, an abundance ~7-40 times lower than the solar value. A second strong depletion must be present during the formation of the thick icy mantles to achieve the values of S/H measured in the dense cold cores (S/H ~8×10-8). Based on our chemical modeling, we constrain the value of ζ H2 to ~ (0.5 - 1.8) ×10-16 s-1 in the translucent cloud.

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.

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.

The Abundance of SiC2 in Carbon Star Envelopes: Evidence that SiC2 is a gas-phase precursor of SiC dust.

CONTEXT: Silicon carbide dust is ubiquitous in circumstellar envelopes around C-rich AGB stars. However, the main gas-phase precursors leading to the formation of SiC dust have not yet been identified. The most obvious candidates among the molecules containing an Si-C bond detected in C-rich AGB stars are SiC2, SiC, and Si2C. To date, the ring molecule SiC2 has been observed in a handful of evolved stars, while SiC and Si2C have only been detected in the C-star envelope IRC +10216. AIMS: We aim to study how widespread and abundant SiC2, SiC, and Si2C are in envelopes around C-rich AGB stars and whether or not these species play an active role as gas-phase precursors of silicon carbide dust in the ejecta of carbon stars. METHODS: We carried out sensitive observations with the IRAM 30m telescope of a sample of 25 C-rich AGB stars to search for emission lines of SiC2, SiC, and Si2C in the λ 2 mm band. We performed non-LTE excitation and radiative transfer calculations based on the LVG method to model the observed lines of SiC2 and to derive SiC2 fractional abundances in the observed envelopes. RESULTS: We detect SiC2 in most of the sources, SiC in about half of them, and do not detect Si2C in any source, at the exception of IRC +10216. Most of these detections are reported for the first time in this work. We find a positive correlation between the SiC and SiC2 line emission, which suggests that both species are chemically linked, the SiC radical probably being the photodissociation product of SiC2 in the external layer of the envelope. We find a clear trend in which the denser the envelope, the less abundant SiC2 is. The observed trend is interpreted as an evidence of efficient incorporation of SiC2 onto dust grains, a process which is favored at high densities owing to the higher rate at which collisions between particles take place. CONCLUSIONS: The observed behavior of a decline in the SiC2 abundance with increasing density strongly suggests that SiC2 is an important gas-phase precursor of SiC dust in envelopes around carbon stars.

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.

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.

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

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 H2 or hot gas (~500-1000 K) in photodissociation regions with high incident FUV radiation field. The excitation may also originate in dense gas (> 105 cm-3) followed by nonreactive collisions with H2, 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 H2, 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 H2. 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.

Trans-cis molecular photoswitching in interstellar Space.

As many organic molecules, formic acid (HCOOH) has two conformers (trans and cis). The energy barrier to internal conversion from trans to cis is much higher than the thermal energy available in molecular clouds. Thus, only the most stable conformer (trans) is expected to exist in detectable amounts. We report the first interstellar detection of cis-HCOOH. Its presence in ultraviolet (UV) irradiated gas exclusively (the Orion Bar photodissociation region), with a low trans-to-cis abundance ratio of 2.8 ± 1.0, supports a photoswitching mechanism: a given conformer absorbs a stellar photon that radiatively excites the molecule to electronic states above the interconversion barrier. Subsequent fluorescent decay leaves the molecule in a different conformer form. This mechanism, which we specifically study with ab initio quantum calculations, was not considered in Space before but likely induces structural changes of a variety of interstellar molecules submitted to UV radiation.

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.

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.