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.

IRC +10 216 in 3-D: morphology of a TP-AGB star envelope.

During their late pulsating phase, AGB stars expel most of their mass in the form of massive dusty envelopes, an event that largely controls the composition of interstellar matter. The envelopes, however, are distant and opaque to visible and NIR radiation: their structure remains poorly known and the mass-loss process poorly understood. Millimeter-wave interferometry, which combines the advantages of longer wavelength, high angular resolution and very high spectral resolution is the optimal investigative tool for this purpose. Mm waves pass through dust with almost no attenuation. Their spectrum is rich in molecular lines and hosts the fundamental lines of the ubiquitous CO molecule, allowing a tomographic reconstruction of the envelope structure. The circumstellar envelope IRC +10 216 and its central star, the C-rich TP-AGB star closest to the Sun, are the best objects for such an investigation. Two years ago, we reported the first detailed study of the CO(2-1) line emission in that envelope, made with the IRAM 30-m telescope. It revealed a series of dense gas shells, expanding at a uniform radial velocity. The limited resolution of the telescope (HPBW 11″) did not allow us to resolve the shell structure. We now report much higher angular resolution observations of CO(2-1), CO(1-0), CN(2-1) and C4H(24-23) made with the SMA, PdB and ALMA interferometers (with synthesized half-power beamwidths of 3″, 1″ and 0.3″, respectively). Although the envelope appears much more intricate at high resolution than with an 11″ beam, its prevailing structure remains a pattern of thin, nearly concentric shells. The average separation between the brightest CO shells is 16″ in the outer envelope, where it appears remarkably constant. Closer to the star (< 40″), the shell pattern is denser and less regular, showing intermediary arcs. Outside the small (r < 0.3″) dust formation zone, the gas appears to expand radially at a constant velocity, 14.5 km s-1, with small turbulent motions. Based on that property, we have reconstructed the 3-D structure of the outer envelope and have derived the gas temperature and density radial profiles in the inner (r < 25″) envelope. The shell-intershell density contrast is found to be typically 3. The over-dense shells have spherical or slightly oblate shapes and typically extend over a few steradians, implying isotropic mass loss. The regular spacing of shells in the outer envelope supports the model of a binary star system with a period of 700 years and a near face-on elliptical orbit. The companion fly-by triggers enhanced episodes of mass loss near periastron. The densification of the shell pattern observed in the central part of the envelope suggests a more complex scenario for the last few thousand years.

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.

ALMA observations of the young protostellar system Barnard 1b: signatures of an incipient hot corino in B1b-S.

The Barnard 1b core shows signatures of being at the earliest stages of low-mass star formation, with two extremely young and deeply embedded protostellar objects. Hence, this core is an ideal target to study the structure and chemistry of the first objects formed in the collapse of prestellar cores. We present ALMA Band 6 spectral line observations at ~0.6″ of angular resolution towards Barnard 1b. We have extracted the spectra towards both protostars, and used a Local Thermodynamic Equilibrium (LTE) model to reproduce the observed line profiles. B1b-S shows rich and complex spectra, with emission from high energy transitions of complex molecules, such as CH3OCOH and CH3CHO, including vibrational level transitions. We have tentatively detected for the first time in this source emission from NH2CN, NH2CHO, CH3CH2OH, CH2OHCHO, CH3CH2OCOH and both aGg' and gGg' conformers of (CH2OH)2. This is the first detection of ethyl formate (CH3CH2OCOH) towards a low-mass star forming region. On the other hand, the spectra of the FHSC candidate B1b-N are free of COMs emission. In order to fit the observed line profiles in B1b-S, we used a source model with two components: an inner hot and compact component (200 K, 0.35″) and an outer and colder one (60 K, 0.6″). The resulting COM abundances in B1b-S range from 10-13 for NH2CN and NH2CHO, up to 10-9 for CH3OCOH. Our ALMA Band 6 observations reveal the presence of a compact and hot component in B1b-S, with moderate abundances of complex organics. These results indicate that a hot corino is being formed in this very young Class 0 source.

Evidence for disks at an early stage in class 0 protostars?

AIMS: The formation epoch of protostellar disks is debated because of the competing roles of rotation, turbulence, and magnetic fields in the early stages of low-mass star formation. Magnetohydrodynamics simulations of collapsing cores predict that rotationally supported disks may form in strongly magnetized cores through ambipolar diffusion or misalignment between the rotation axis and the magnetic field orientation. Detailed studies of individual sources are needed to cross check the theoretical predictions. METHODS: We present 0.06 - 0.1 ″ resolution images at 350 GHz toward B1b-N and B1b-S, which are young class 0 protostars, possibly first hydrostatic cores. The images have been obtained with ALMA, and we compare these data with magnetohydrodynamics simulations of a collapsing turbulent and magnetized core. RESULTS: The submillimeter continuum emission is spatially resolved by ALMA. Compact structures with optically thick 350 GHz emission are detected toward both B1b-N and B1b-S, with 0.2 and 0.35″ radii (46 and 80 au at the Perseus distance of 230 pc), within a more extended envelope. The flux ratio between the compact structure and the envelope is lower in B1b-N than in B1b-S, in agreement with its earlier evolutionary status. The size and orientation of the compact structure are consistent with 0.2″ resolution 32 GHz observations obtained with the Very Large Array as a part of the VANDAM survey, suggesting that grains have grown through coagulation. The morphology, temperature, and densities of the compact structures are consistent with those of disks formed in numerical simulations of collapsing cores. Moreover, the properties of B1b-N are consistent with those of a very young protostar, possibly a first hydrostatic core. These observations provide support for the early formation of disks around low-mass protostars.

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.

Chemical segregation in the young protostars Barnard 1b-N and S Evidence of pseudo-disk rotation in Barnard 1b-S.

The extremely young Class 0 object B1b-S and the first hydrostatic core (FSHC) candidate, B1b-N, provide a unique opportunity to study the chemical changes produced in the elusive transition from the prestellar core to the protostellar phase. We present 40"×70" images of Barnard 1b in the 13CO 1→0, C18O 1→0, NH2D 11,1a→10,1s, and SO 32→21 lines obtained with the NOEMA interferometer. The observed chemical segregation allows us to unveil the physical structure of this young protostellar system down to scales of ∼500 au. The two protostellar objects are embedded in an elongated condensation, with a velocity gradient of ∼0.2-0.4 m s-1 au-1 in the east-west direction, reminiscent of an axial collapse. The NH2D data reveal cold and dense pseudo-disks (R∼500-1000 au) around each protostar. Moreover, we observe evidence of pseudo-disk rotation around B1b-S. We do not see any signature of the bipolar outflows associated with B1b-N and B1b-S, which were previously detected in H2CO and CH3OH, in any of the imaged species. The non-detection of SO constrains the SO/CH3OH abundance ratio in the high-velocity gas.

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

Ionization fraction and the enhanced sulfur chemistry in Barnard 1.

CONTEXT: Barnard B1b has revealed as one of the most interesting globules from the chemical and dynamical point of view. It presents a rich molecular chemistry characterized by large abundances of deuterated and complex molecules. Furthermore, it hosts an extremely young Class 0 object and one candidate to First Hydrostatic Core (FHSC) proving the youth of this star forming region. AIMS: Our aim is to determine the cosmic ray ionization rate, [Formula: see text], and the depletion factors in this extremely young star forming region. These parameteres determine the dynamical evolution of the core. METHODS: We carried out a spectral survey towards Barnard 1b as part of the IRAM Large program ASAI using the IRAM 30-m telescope at Pico Veleta (Spain). This provided a very complete inventory of neutral and ionic C-, N- and S- bearing species with, up to our knowledge, the first secure detections of the deuterated ions DCS+ and DOCO+. We use a state-of-the-art pseudo-time-dependent gas-phase chemical model that includes the ortho and para forms of [Formula: see text] and [Formula: see text] to determine the local value of the cosmic ray ionization rate and the depletion factors. RESULTS: Our model assumes n(H2)=105 cm-3 and T k =12 K, as derived from our previous works. The observational data are well fitted with ζH2 between 3×10-17 s-1 and 10-16 s-1, and the following elemental abundances: O/H=3 10-5, N/H=6.4-8 10-5, C/H=1.7 10-5 and S/H between 6.0 10-7 and 1.0 10-6. The large number of neutral/protonated species detected, allows us to derive the elemental abundances and cosmic ray ionization rate simultaneously. Elemental depletions are estimated to be ~10 for C and O, ~1 for N and ~25 for S. CONCLUSIONS: Barnard B1b presents similar depletions of C and O than those measured in pre-stellar cores. The depletion of sulfur is higher than that of C and O but not as extreme as in cold cores. In fact, it is similar to the values found in some bipolar outflows, hot cores and photon-dominated regions. Several scenarios are discussed to account for these peculiar abundances. We propose that it is the consequence of the initial conditions (important outflows and enhanced UV fields in the surroundings) and a rapid collapse (~0.1 Myr) that permits to maintain most S- and N-bearing species in gas phase to great optical depths. The interaction of the compact outflow associated with B1b-S with the surrounding material could enhance the abundances of S-bearing molecules, as well.

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.

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.

Nascent bipolar outflows associated with the first hydrostatic core candidates Barnard 1b-N and 1b-S.

In the theory of star formation, the first hydrostatic core (FHSC) phase is a critical step in which a condensed object emerges from a prestellar core. This step lasts about one thousand years, a very short time compared with the lifetime of prestellar cores, and therefore is hard to detect unambiguously. We present IRAM Plateau de Bure observations of the Barnard 1b dense molecular core, combining detections of H2CO and CH3OH spectral lines and dust continuum at 2.3" resolution (~ 500 AU). The two compact cores B1b-N and B1b-S are detected in the dust continuum at 2mm, with fluxes that agree with their spectral energy distribution. Molecular outflows associated with both cores are detected. They are inclined relative to the direction of the magnetic field, in agreement with predictions of collapse in turbulent and magnetized gas with a ratio of mass to magnetic flux somewhat higher than the critical value, µ ~ 2 - 7. The outflow associated with B1b-S presents sharp spatial structures, with ejection velocities of up to ~ 7 kms-1 from the mean velocity. Its dynamical age is estimated to be ~ 2000 yr. The B1b-N outflow is smaller and slower, with a short dynamical age of ~ 1000 yr. The B1b-N outflow mass, mass-loss rate, and mechanical luminosity agree well with theoretical predictions of FHSC. These observations confirm the early evolutionary stage of B1b-N and the slightly more evolved stage of B1b-S.