Depletion and low gas temperature in the L183 prestellar core: the N2H+ - N2D+ tool

Context. The study of pre-stellar cores (PSCs) suffers from a lack of undepleted species to trace the gas physical properties in their very dense inner parts. Aims. We want to carry out detailed modelling of N2H+ and N2D+ cuts across the L183 main core to evaluate the depletion of these species and their usefulness as a probe of physical conditions in PSCs. Methods. We have developed a non-LTE (NLTE) Monte-Carlo code treating the 1D radiative transfer of both N2H+ and N2D+, making use of recently published collisional coefficients with He between individual hyperfine levels. The code includes line overlap between hyperfine transitions. An extensive set of core models is calculated and compared with observations. Special attention is paid to the issue of source coupling to the antenna beam. Results. The best fitting models indicate that i) gas in the core center is very cold (7$\pm$ 1 K) and thermalized with dust, ii) depletion of N2H+ does occur, starting at densities 5-7E5 cm−3 and reaching a factor of 6 (+13/−3) in abundance, iii) deuterium fractionation reaches ∼70% at the core center, and iv) the density profile is proportional to r^-1 out to ∼4000 AU, and to r^−2 beyond. Conclusions. Our NLTE code could be used to (re-)interpret recent and upcoming observations of N2H+ and N2D+ in many pre-stellar cores of interest, to obtain better temperature and abundance profiles

The Initial Conditions of Clustered Star Formation III. The Deuterium Fractionation of the Ophiuchus B2 Core

We present N2D+ 3-2 (IRAM) and H2D+ 1_11 - 1_10 and N2H+ 4-3 (JCMT) maps of the small cluster-forming Ophiuchus B2 core in the nearby Ophiuchus molecular cloud. In conjunction with previously published N2H+ 1-0 observations, the N2D+ data reveal the deuterium fractionation in the high density gas across Oph B2. The average deuterium fractionation R_D = N(N2D+)/N(N2H+) ~ 0.03 over Oph B2, with several small scale R_D peaks and a maximum R_D = 0.1. The mean R_D is consistent with previous results in isolated starless and protostellar cores. The column density distributions of both H2D+ and N2D+ show no correlation with total H2 column density. We find, however, an anticorrelation in deuterium fractionation with proximity to the embedded protostars in Oph B2 to distances >= 0.04 pc. Destruction mechanisms for deuterated molecules require gas temperatures greater than those previously determined through NH3 observations of Oph B2 to proceed. We present temperatures calculated for the dense core gas through the equating of non-thermal line widths for molecules (i.e., N2D+ and H2D+) expected to trace the same core regions, but the observed complex line structures in B2 preclude finding a reasonable result in many locations. This method may, however, work well in isolated cores with less complicated velocity structures. Finally, we use R_D and the H2D+ column density across Oph B2 to set a lower limit on the ionization fraction across the core, finding a mean x_e, lim >= few x 10^{-8}. Our results show that care must be taken when using deuterated species as a probe of the physical conditions of dense gas in star-forming regions.Comment: ApJ accepte

The puzzling detection of D_2CO in the molecular cloud L1689N

We present new observations of the D_2CO emission towards the small cloud L1689N in the ρ Ophiuchus complex. We surveyed five positions, three being a cut across a shock site and two probing the quiescent gas of the molecular cloud. We detected D_2CO emission in the first three positions. The measured [D_2CO] /[ H2CO] is about 3%, whereas it is ≤2% in the quiescent gas. We discuss the implications of these new observations, which suggest that the bulk of the D_2CO molecules is stored in grain mantles, and removed from the cold storage by the shock at the interface between the outflowing and quiescent gas. We review the predictions of the published models proposed to explain the observed high deuteration of formaldehyde. They fall in two basic schemes: gas phase and grain surface chemistry. None of the reviewed models is able to account for the observed [D_2CO] /[H_2CO] abundance ratio. A common characteristics shared by the models is apparently that all underestimate the atomic [D]/[H] ratio in the accreting gas

SOLIS XVII: Jet candidate unveiled in OMC-2 and its possible link to the enhanced cosmic-ray ionisation rate

The study of the early phases of star and planet formation is important to understand the physical and chemical history of stellar systems such as our own. In particular, protostars born in rich clusters are prototypes of the young Solar System. In the framework of the Seeds Of Life In Space (SOLIS) large observational project, the aim of the present work is to investigate the origin of the previously inferred high flux of energetic particles in the protocluster FIR4 of the Orion Molecular Cloud 2 (OMC-2), which appears asymmetric within the protocluster itself. Interferometric observations carried out with the IRAM NOEMA interferometer were used to map the silicon monoxide (SiO) emission around the FIR4 protocluster. Complementary archival data from the ALMA interferometer were also employed to help constrain excitation conditions. A physical-chemical model was implemented to characterise the particle acceleration along the protostellar jet candidate, along with a non-LTE analysis of the SiO emission along the jet. The emission morphology of the SiO rotational transitions hints for the first time at the presence of a collimated jet originating very close to the brightest protostar in the cluster, HOPS-108. The NOEMA observations unveiled a possible jet in the OMC-2 FIR4 protocluster propagating towards a previously measured enhanced cosmic-ray ionisation rate. This suggests that energetic particle acceleration by the jet shock close to the protostar might be at the origin of the enhanced cosmic-ray ionisation rate, as confirmed by modelling the protostellar jet.Comment: 8 pages, 7 figures, accepted for publications in A&

Deuterated water in the solar-type protostars NGC 1333 IRAS 4A and IRAS 4B

Aims. The aim of this paper is to study deuterated water in the solar-type protostars NGC1333 IRAS4A and IRAS4B, to compare their HDO abundance distribution with other star-forming regions, and to constrain their HDO/H2O ratios. Methods. Using the Herschel/HIFI instrument as well as ground-based telescopes, we observed several HDO lines covering a large excitation range (Eup/k=22-168 K) towards these protostars and an outflow position. Non-LTE radiative transfer codes were then used to determine the HDO abundance profiles in these sources. Results. The HDO fundamental line profiles show a very broad component, tracing the molecular outflows, in addition to a narrower emission component and a narrow absorbing component. In the protostellar envelope of NGC1333 IRAS4A, the HDO inner (T>100 K) and outer (T<100 K) abundances with respect to H2 are estimated at 7.5x10^{-9} and 1.2x10^{-11}, respectively, whereas, in NGC1333 IRAS4B, they are 1.0x10^{-8} and 1.2x10^{-10}, respectively. Similarly to the low-mass protostar IRAS16293-2422, an absorbing outer layer with an enhanced abundance of deuterated water is required to reproduce the absorbing components seen in the fundamental lines at 465 and 894 GHz in both sources. This water-rich layer is probably extended enough to encompass the two sources as well as parts of the outflows. In the outflows emanating from NGC1333 IRAS4A, the HDO column density is estimated at about (2-4)x10^{13} cm^{-2}, leading to an abundance of about (0.7-1.9)x10^{-9}. An HDO/H2O ratio between 7x10^{-4} and 9x10^{-2} is derived in the outflows. In the warm inner regions of these two sources, we estimate the HDO/H2O ratios at about 1x10^{-4}-4x10^{-3}. This ratio seems higher (a few %) in the cold envelope of IRAS4A, whose possible origin is discussed in relation to formation processes of HDO and H2O.Comment: 16 pages, 13 figure

Seeds of Life in Space (SOLIS): XI. First measurement of nitrogen fractionation in shocked clumps of the L1157 protostellar outflow

Context. The isotopic ratio of nitrogen presents a wide range of values in the Solar System: from 140 in meteorites and comets to 441 in the solar wind. In star-forming systems, we observe evena higher spread of ~150-1000. The origin of these differences is still unclear. Aims. Chemical reactions in the gas phase are one of the possible processes that could modify the 14N/15N ratio. We aim to investigate if and how the passage of a shock wave in the interstellar medium, which activates a rich chemistry, can affect the relative fraction of nitrogen isotopes. Theideal place for such a study is the chemically rich outflow powered by the L1157-mm protostar, where several shocked clumps are present. Methods. We present the first measurement of the 14N/15N ratio in the two shocked clumps, B1 and B0, of the protostellar outflow L1157. The measurement is derived from the interferometeric maps of the H13CN (1-0) and the HC15N (1-0) lines obtained with the NOrthern Extended Millimeter Array (NOEMA) interferometeras part of the Seeds of Life in Space (SOLIS) programme. Results. In B1, we find that the H13CN (1-0) and HC15N (1-0) emission traces the front of the clump, that is the apex of the shocked region, where the fast jet impacts the lower velocity medium with an averaged column density of N(H13CN) ~ 7 × 1012 cm-2 and N(HC15N) 2 × 1012 cm-2. In this region, the ratio H13CN (1-0)/HC15N (1-0) is almost uniform with an average value of ~5 ± 1. The same average value isalso measured in the smaller clump B0e. Assuming the standard 12C/13C = 68, we obtain 14N/15N = 340 ± 70. This ratio is similar to those usually found with the same species in prestellar cores and protostars. We analysed the prediction of a chemical shock model for several shock conditions and we found that the nitrogen and carbon fractionations do not vary much for the first period after the shock. The observed H13CN/HC15N can be reproduced by a non-dissociative, C-type shock with pre-shock density n(H) = 105 cm-3, shock velocity Vs between 20 and 40 km s-1, and cosmic-ray ionization rate of 3 × 10-16 s-1; this agrees with previous modelling of other chemical species in L1157-B1. Conclusions. Both observations and chemical models indicate that the rich chemistry activated by the shock propagation does not affect the nitrogen isotopic ratio, which remains similar to that measured in lower temperature gas in prestellar cores and protostellar envelopes

CO Depletion and Deuterium Fractionation in Prestellar Cores

We report the detection of D2CO in a sample of starless dense cores, in which we previously measured the degree of CO depletion. The deuterium fractionation is found extremely high, [D2CO]/[H2CO] ~ 1-10 %, similar to that reported in low-mass protostars. This provides convincing evidence that D2CO is formed in the cold pre-stellar cores, and later desorbed when the gas warms up in protostars. We find that the cores with the highest CO depletions have also the largest [D2CO]/[H2CO] ratios, supporting the theoretical prediction that deuteration increases with increasing CO depletion.Comment: 11 pages, 2 figures, accepted by ApJ Letter

Herschel/HIFI observations of interstellar OH+ and H2O+ towards W49N: a probe of diffuse clouds with a small molecular fraction

We report the detection of absorption by interstellar hydroxyl cations and water cations, along the sight-line to the bright continuum source W49N. We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 972 GHz N = 1 - 0 transition of OH+ and the 1115 GHz 1(11) - 0(00) transition of ortho-H2O+. The resultant spectra show absorption by ortho-H2O+, and strong absorption by OH+, in foreground material at velocities in the range 0 to 70 km/s with respect to the local standard of rest. The inferred OH+/H2O+ abundance ratio ranges from ~ 3 to ~ 15, implying that the observed OH+ arises in clouds of small molecular fraction, in the 2 - 8% range. This conclusion is confirmed by the distribution of OH+ and H2O+ in Doppler velocity space, which is similar to that of atomic hydrogen, as observed by means of 21 cm absorption measurements, and dissimilar from that typical of other molecular tracers. The observed OH+/H abundance ratio of a few E-8 suggests a cosmic ray ionization rate for atomic hydrogen of (0.6 - 2.4) E-16 s-1, in good agreement with estimates inferred previously for diffuse clouds in the Galactic disk from observations of interstellar H3+ and other species.Comment: Accepted for publication in A&A Letter

The puzzling detection of D_2CO in the molecular cloud L1689N

We present new observations of the D_2CO emission towards the small cloud L1689N in the ρ Ophiuchus complex. We surveyed five positions, three being a cut across a shock site and two probing the quiescent gas of the molecular cloud. We detected D_2CO emission in the first three positions. The measured [D_2CO] /[ H2CO] is about 3%, whereas it is ≤2% in the quiescent gas. We discuss the implications of these new observations, which suggest that the bulk of the D_2CO molecules is stored in grain mantles, and removed from the cold storage by the shock at the interface between the outflowing and quiescent gas. We review the predictions of the published models proposed to explain the observed high deuteration of formaldehyde. They fall in two basic schemes: gas phase and grain surface chemistry. None of the reviewed models is able to account for the observed [D_2CO] /[H_2CO] abundance ratio. A common characteristics shared by the models is apparently that all underestimate the atomic [D]/[H] ratio in the accreting gas

Reversal of infall in SgrB2(M) revealed by Herschel/HIFI observations of HCN lines at THz frequencies

To investigate the accretion and feedback processes in massive star formation, we analyze the shapes of emission lines from hot molecular cores, whose asymmetries trace infall and expansion motions. The high-mass star forming region SgrB2(M) was observed with Herschel/HIFI (HEXOS key project) in various lines of HCN and its isotopologues, complemented by APEX data. The observations are compared to spherically symmetric, centrally heated models with density power-law gradient and different velocity fields (infall or infall+expansion), using the radiative transfer code RATRAN. The HCN line profiles are asymmetric, with the emission peak shifting from blue to red with increasing J and decreasing line opacity (HCN to H$^{13}$CN). This is most evident in the HCN 12--11 line at 1062 GHz. These line shapes are reproduced by a model whose velocity field changes from infall in the outer part to expansion in the inner part. The qualitative reproduction of the HCN lines suggests that infall dominates in the colder, outer regions, but expansion dominates in the warmer, inner regions. We are thus witnessing the onset of feedback in massive star formation, starting to reverse the infall and finally disrupting the whole molecular cloud. To obtain our result, the THz lines uniquely covered by HIFI were critically important.Comment: A&A, HIFI special issue, accepte