CO Depletion Factor Research Articles

The ALMA Survey of 70 μm Dark High-mass Clumps in Early Stages (ASHES). VI. The Core-scale CO Depletion

Studying the physical and chemical properties of cold and dense molecular clouds is crucial for the understanding of how stars form. Under the typical conditions of infrared dark clouds, CO is removed from the gas phase and trapped onto the surface of dust grains by the so-called depletion process. This suggests that the CO-depletion factor (f D ) can be a useful chemical indicator for identifying cold and dense regions (i.e., prestellar cores). We have used the 1.3 mm continuum and C18O (2–1) data observed at the resolution of ∼5000 au in the ALMA Survey of 70 μm Dark High-mass Clumps in Early Stages (ASHES) to construct averaged maps of f D in 12 clumps to characterize the earliest stages of the high-mass star formation process. The average f D determined for 277 of the 294 ASHES cores follows an unexpected increase from the prestellar to the protostellar stage. If we exclude the temperature effect due to the slight variations in the NH3 kinetic temperature among different cores, we explain this result as a dependence primarily on the average gas density, which increases in cores where protostellar conditions prevail. This shows that f D determined in high-mass star-forming regions at the core scale is insufficient to distinguish among prestellar and protostellar conditions for the individual cores and should be complemented by information provided by additional tracers. However, we confirm that the clump-averaged f D values correlate with the luminosity-to-mass ratio of each source, which is known to trace the evolution of the star formation process.

Astrochemical modelling of infrared dark clouds

Context. Infrared dark clouds (IRDCs) are cold, dense regions of the interstellar medium (ISM) that are likely to represent the initial conditions for massive star and star cluster formation. It is thus important to study the physical and chemical conditions of IRDCs to provide constraints and inputs for theoretical models of these processes. Aims. We aim to determine the astrochemical conditions, especially the cosmic ray ionisation rate (CRIR) and chemical age, in different regions of the massive IRDC G28.37+00.07 by comparing observed abundances of multiple molecules and molecular ions with the predictions of astrochemical models. Methods. We have computed a series of single-zone, time-dependent, astrochemical models with a gas-grain network that systematically explores the parameter space of the density, temperature, CRIR, and visual extinction. We have also investigated the effects of choices of CO ice binding energy and temperatures achieved in the transient heating of grains when struck by cosmic rays. We selected ten positions across the IRDC that are known to have a variety of star formation activity. We utilised mid-infrared extinction maps and sub-millimetre (sub-mm) emission maps to measure the mass surface densities of these regions needed for abundance and volume density estimates. The sub-mm emission maps were also used to measure temperatures. We then used Instituto de Radioas-tromía Milimétrica (IRAM) 30 m observations of various tracers, especially C18O(1-0), H13CO+(1-0), HC18O+(1-0), and N2H+(1-0), to estimate column densities and thus abundances. Finally, we investigated the range of astrochemical conditions that are consistent with the observed abundances. Results. The typical physical conditions of the IRDC regions are nH ~ 3 × 104 to 105 cm−3 and T ≃ 10 to 15 K. Strong emission of H13CO+(1-0) and N2H+(1-0) is detected towards all the positions and these species are used to define relatively narrow velocity ranges of the IRDC regions, which are used for estimates of CO abundances, via C18O(1-0). We would like to note that CO depletion factors are estimated to be in the range fD ~ 3 to 10. Using estimates of the abundances of CO, HCO+, and N2H+, we find consistency with astrochemical models that have relatively low CRIRs of ζ ~ 10−18 to ~10−17 s−1, with no evidence for systematic variation with the level of star formation activity. Astrochemical ages, which are defined with a reference to an initial condition of all H in H2, all C in CO, and all other species in atomic form, are found to be <1 Myr. We also explore the effects of using other detected species, that is HCN, HNC, HNCO, CH3OH, and H2CO, to constrain the models. These generally lead to implied conditions with higher levels of CRIRs and older chemical ages. Considering the observed fD versus nH relation of the ten positions, which we find to have relatively little scatter, we discuss potential ways in which the astrochemical models can match such a relation as a quasi-equilibrium limit valid at ages of at least a few free-fall times, that is ≳0.3 Myr, including the effect of CO envelope contamination, small variations in temperature history near 15 K, CO-ice binding energy uncertainties, and CR-induced desorption. We find general consistency with the data of ~0.5 Myr-old models that have ζ ~ 2-5 × 10−18 s−1 and CO abundances set by a balance of freeze-out with CR-induced desorption. Conclusions. We have constrained the astrochemical conditions in ten regions in a massive IRDC, finding evidence for relatively low values of CRIR compared to diffuse ISM levels. We have not seen clear evidence for variation in the CRIR with the level of star formation activity. We favour models that involve relatively low CRIRs (≲10−17 s−1) and relatively old chemical ages (≳0.3 Myr, i.e. ≳3tff). We discuss potential sources of systematic uncertainties in these results and the overall implications for IRDC evolutionary history and astrochemical models.

Deuterium chemistry and D/H ratios in Class 0/I proto-brown dwarfs

ABSTRACT We have conducted the first extensive observational survey of several deuterated species in 16 Class 0/I proto-brown dwarfs (proto-BDs) and 4 Class Flat/Class II brown dwarfs. Observations were obtained with the IRAM 30m telescope in the DCO+ (3–2), DCN (3–2), DNC (3–2), and N2D+ (3–2) lines. The DCO+/H13CO+, DCN/H13CN, and DNC/HN13C ratios are comparatively higher and show a narrower range than the DCO+/HCO+, DCN/HCN, and DNC/HNC ratios, respectively. The mean D/H ratios for the proto-BDs derived from these molecules are in the range of ∼0.02–3. Both low-temperature gas-phase ion-molecule deuteron transfer and grain surface reactions are required to explain the enhanced deuterium fractionation. The very dense and cold ($n_{H_{2}} \ge 10^{6}$ cm−3, T ≤10 K) interior of the proto-BDs provide the suitable conditions for efficient deuterium fractionation in these cores. There is no correlation between the D/H ratios and the CO depletion factor, with the exception of the DCN/HCN ratios that show a strong anti-correlation possibly due to the difference in the peak emitting regions of the DCN and HCN molecules. Over a wide range in the bolometric luminosities spanning ∼0.002–40 L⊙, we find a trend of higher DCO+/HCO+ and DCN/HCN ratios, nearly constant DNC/HNC and DNC/HN13C ratios, and lower N2D+/N2H+ ratios in the proto-BDs compared to protostars. Only one Class II brown dwarf shows emission in the DCO+ (3–2) line. No correlation is seen between the D/H ratios and the evolutionary stage.

Deuterium chemodynamics of massive pre-stellar cores

ABSTRACT High levels of deuterium fractionation of N2H+ (i.e. $\mathrm{D^{{\mathrm{N_2H^+}}}_{\text{frac}}}$≳ 0.1) are often observed in pre-stellar cores (PSCs) and detection of N2D+ is a promising method to identify elusive massive PSCs. However, the physical and chemical conditions required to reach such high levels of deuteration are still uncertain, as is the diagnostic utility of N2H+ and N2D+ observations of PSCs. We perform 3D magnetohydrodynamics simulations of a massive, turbulent, magnetized PSC, coupled with a sophisticated deuteration astrochemical network. Although the core has some magnetic/turbulent support, it collapses under gravity in about one free-fall time, which marks the end of the simulations. Our fiducial model achieves relatively low $\mathrm{D^{{\mathrm{N_2H^+}}}_{\text{frac}}}$∼0.002 during this time. We then investigate effects of initial ortho-para ratio of H2 ($\mathrm{OPR^{H_2}}$), temperature, cosmic ray (CR) ionization rate, CO and N-species depletion factors, and prior PSC chemical evolution. We find that high CR ionization rates and high depletion factors allow the simulated $\mathrm{D^{{\mathrm{N_2H^+}}}_{\text{frac}}}$ and absolute abundances to match observational values within one free-fall time. For $\mathrm{OPR^{H_2}}$, while a lower initial value helps the growth of $\mathrm{D^{{\mathrm{N_2H^+}}}_{\text{frac}}}$, the spatial structure of deuteration is too widespread compared to observed systems. For an example model with elevated CR ionization rates and significant heavy element depletion, we then study the kinematic and dynamic properties of the core as traced by its N2D+ emission. The core, undergoing quite rapid collapse, exhibits disturbed kinematics in its average velocity map. Still, because of magnetic support, the core often appears kinematically subvirial based on its N2D+ velocity dispersion.

Survey of ortho-H2D+ in high-mass star-forming regions

Context. Deuteration has been suggested to be a reliable chemical clock of star-forming regions due to its strong dependence on density and temperature changes during cloud contraction. In particular, the H3+ isotopologues (e.g. ortho-H2D+) seem to act as good proxies of the evolutionary stages of the star formation process. While this has been widely explored in low-mass star-forming regions, in the high-mass counterparts only a few studies have been pursued, and the reliability of deuteration as a chemical clock remains inconclusive. Aims. We present a large sample of o-H2D+ observations in high-mass star-forming regions and discuss possible empirical correlations with relevant physical quantities to assess its role as a chronometer of star-forming regions through different evolutionary stages. Methods. APEX observations of the ground-state transition of o-H2D+ were analysed in a large sample of high-mass clumps selected from the ATLASGAL survey at different evolutionary stages. Column densities and beam-averaged abundances of o-H2D+ with respect to H2, X(o-H2D+), were obtained by modelling the spectra under the assumption of local thermodynamic equilibrium. Results. We detect 16 sources in o-H2D+ and find clear correlations between X(o-H2D+) and the clump bolometric luminosity and the dust temperature, while only a mild correlation is found with the CO-depletion factor. In addition, we see a clear correlation with the luminosity-to-mass ratio, which is known to trace the evolution of the star formation process. This would indicate that the deuterated forms of H3+ are more abundant in the very early stages of the star formation process and that deuteration is influenced by the time evolution of the clumps. In this respect, our findings would suggest that the X(o-H2D+) abundance is mainly affected by the thermal changes rather than density changes in the gas. We have employed these findings together with observations of H13CO+, DCO+, and C17O to provide an estimate of the cosmic-ray ionisation rate in a sub-sample of eight clumps based on recent analytical work. Conclusions. Our study presents the largest sample of o-H2D+ in star-forming regions to date. The results confirm that the deuteration process is strongly affected by temperature and suggests that o-H2D+ can be considered a reliable chemical clock during the star formation processes, as proved by its strong temporal dependence.

The chemical structure of the very young starless core L1521E

Context. L1521E is a dense starless core in Taurus that was found to have relatively low molecular depletion by earlier studies, thus suggesting a recent formation. Aims. We aim to characterize the chemical structure of L1521E and compare it to the more evolved L1544 pre-stellar core. Methods. We have obtained ~2.5 × 2.5 arcminute maps toward L1521E using the IRAM-30 m telescope in transitions of various species, including C17O, CH3OH, c-C3H2, CN, SO, H2CS, and CH3CCH. We derived abundances for the observed species and compared them to those obtained toward L1544. We estimated CO depletion factors using the C17O IRAM-30 m map, an N(H2) map derived from Herschel/SPIRE data and a 1.2 mm dust continuum emission map obtained with the IRAM-30 m telescope. Results. Similarly to L1544, c-C3H2 and CH3OH peak at different positions. Most species peak toward the c-C3H2 peak including C2S, C3S, HCS+, HC3N, H2CS, CH3CCH, and C34S. C17O and SO peak close to both the c-C3H2 and the CH3OH peaks. CN and N2H+ peak close to the Herschel dust peak. We found evidence of CO depletion toward L1521E. The lower limit of the CO depletion factor derived toward the Herschel dust peak is 4.3±1.6, which is about a factor of three lower than toward L1544. We derived abundances for several species toward the dust peaks of L1521E and L1544. The abundances of most sulfur-bearing molecules such as C2S, HCS+, C34S, C33S, and HCS+ are higher toward L1521E than toward L1544 by factors of ~2–20, compared to the abundance of A-CH3OH. The abundance of methanol is very similar toward the two cores. Conclusions. The fact that the abundances of sulfur-bearing species toward L1521E are higher than toward L1544 suggests that significant sulfur depletion takes place during the dynamical evolution of dense cores, from the starless to pre-stellar stage. The CO depletion factor measured toward L1521E suggests that CO is more depleted than previously found. Similar CH3OH abundances between L1521E and L1544 hint that methanol is forming at specific physical conditions in the Taurus Molecular Cloud Complex, characterized by densities of a few ×104 cm−3 and N(H2) ≳ 1022 cm−2, when CO starts to catastrophically freeze-out, while water can still be significantly photodissociated, so that the surfaces of dust grains become rich in solid CO and CH3OH, as already found toward L1544. Methanol can thus provide selective crucial information about the transition region between dense cores and the surrounding parent cloud.

A Study of the c-C3HD/c-C3H2 Ratio in Low-mass Star-forming Regions*

Abstract We use the deuteration of c-C3H2 to probe the physical parameters of starless and protostellar cores, related to their evolutionary states, and compare it to the N2H+-deuteration in order to study possible differences between the deuteration of C- and N-bearing species. We observed the main species c-C3H2, the singly and doubly deuterated species c-C3HD and c-C3D2, as well as the isotopologue c-H13CC2H toward 10 starless cores and five protostars in the Taurus and Perseus complexes. We examined the correlation between the N(c-C3HD)/N(c-C3H2) ratio and the dust temperature along with the H2 column density and the CO depletion factor. The resulting N(c-C3HD)/N(c-C3H2) ratio is, within error bars, consistent with 10% in all starless cores with detected c-C3HD. This also accounts for the protostars except for the source HH211, where we measure a high deuteration level of 23%. The deuteration of N2H+ follows the same trend but is considerably higher in the dynamically evolved core L1544. We find no significant correlation between the deuteration of c-C3H2 and the CO depletion factor among the starless and protostellar cores. Toward the latter the coolest objects show the largest deuterium fraction in c-C3H2. We show that the deuteration of c-C3H2 can trace the early phases of star formation and is comparable to that of N2H+. However, the largest c-C3H2 deuteration level is found toward protostellar cores, suggesting that while c-C3H2 is mainly frozen onto dust grains in the central regions of starless cores, active deuteration is taking place on ice.

Star Formation Conditions in a Planck Galactic Cold Clump, G108.84–00.81

Abstract We present the results from a series of ground-based radio observations toward a Planck Galactic Cold Clump (PGCC), PGCC G108.84–00.81, which is located in one curved filamentary cloud in the vicinity of an extended H ii region Sh2-152 and SNR G109.1-1.0. PGCC G108.84–00.81 is mainly composed of two clumps, “G108–N” and “G108–S”. In the 850 μm dust continuum emission map, G108–N is shown as one component while G108–S is fragmented into four components. There is no infrared source associated with G108–N, while there are two infrared sources (IRS 1 and IRS 2) associated with G108–S. The total mass of G108–N is larger than the Jeans mass, suggesting that G108–N is gravitationally unstable and a potential place for a future star formation. The clump properties of G108–N and G108–S such as the gas temperature and the column density, are not distinctly different. However, G108–S is slightly more evolved than G108–N, if considering the CO depletion factor, molecular abundances, and association with infrared sources. G108–S seems to be affected by the compression from Sh2-152, while G108–N is relatively protected from the external effect.

Characterising the physical and chemical properties of a young Class 0 protostellar core embedded in the Orion B9 filament

The present study aims to characterise the physical and chemical properties of the protostellar core Orion B9-SMM3. The APEX telescope was used to perform a follow-up molecular line survey of SMM3. The following species were identified from the frequency range 218.2-222.2 GHz: $^{13}$CO, C$^{18}$O, SO, para-H$_2$CO, and E$_1$-type CH$_3$OH. The on-the-fly mapping observations at 215.1-219.1 GHz revealed that SMM3 is associated with a dense gas core as traced by DCO$^+$ and p-H$_2$CO. Altogether three different p-H$_2$CO transitions were detected with clearly broadened linewidths (8.2-11 km s$^{-1}$ in FWHM). The derived p-H$_2$CO rotational temperature, $64\pm15$ K, indicates the presence of warm gas. We also detected a narrow p-H$_2$CO line (FWHM=0.42 km s$^{-1}$) at the systemic velocity. The p-H$_2$CO abundance for the broad component appears to be enhanced by two orders of magnitude with respect to the narrow line value ($\sim3\times10^{-9}$ versus $\sim2\times10^{-11}$). The detected methanol line shows a linewidth similar to those of the broad p-H$_2$CO lines, which indicates their coexistence. The CO isotopologue data suggest that the CO depletion factor decreases from $\sim27\pm2$ towards the core centre to a value of $\sim8\pm1$ towards the core edge. In the latter position, the N$_2$D$^+$/N$_2$H$^+$ ratio is revised down to $0.14\pm0.06$. The origin of the subfragments inside the SMM3 core we found previously can be understood in terms of the Jeans instability if non-thermal motions are taken into account. The estimated fragmentation timescale, and the derived chemical abundances suggest that SMM3 is a few times $10^5$ yr old, in good agreement with its Class 0 classification inferred from the spectral energy distribution analysis. The broad p-H$_2$CO and CH$_3$OH lines, and the associated warm gas provide the first clear evidence of a molecular outflow driven by SMM3.

Deuterium fractionation in the Ophiuchus molecular cloud

Aims. We measure the deuterium fraction, RD, and the CO-depletion factor, fd, toward a number of starless and protostellar cores in the L1688 region of the Ophiuchus molecular cloud complex and search for variations based upon environmental differences across L1688. The kinematic properties of the dense gas traced by the N2H+ and N2D+ (1-0) lines are also discussed. Methods. RD has been measured via observations of the J=1-0 transition of N2H+ and N2D+ toward 33 dense cores in different regions of L1688. fd estimates have been done using C17O(1-0) and 850 micron dust continuum emission from the SCUBA survey. All line observations were carried out with the IRAM 30 meter antenna. Results. The dense cores show large (2-40%) deuterium fractions, with significant variations between the sub-regions of L1688. The CO-depletion factor also varies from one region to another (1-7). Two different correlations are found between deuterium fraction and CO-depletion factor: cores in regions A, B2 and I show increasing RD with increasing fd, similar to previous studies of deuterium fraction in pre-stellar cores; cores in regions B1, B1B2, C, E, F and H show a steeper RD-fd correlation, with large deuterium fractions occurring in fairly quiescent gas with relatively low CO freeze-out factors. These are probably recently formed, centrally concentrated starless cores which have not yet started the contraction phase toward protostellar formation. We also find that the deuterium fraction is affected by the amount of turbulence, dust temperature and distance from heating sources in all regions of L1688, although no clear trend is found.

Gas kinematics and excitation in the filamentary IRDC G035.39-00.33

Some theories of dense molecular cloud formation involve dynamical environments driven by converging atomic flows or collisions between preexisting molecular clouds. The determination of the dynamics and physical conditions of the gas in clouds at the early stages of their evolution is essential to establish the dynamical imprints of such collisions, and to infer the processes involved in their formation. We present multi-transition 13CO and C18O maps toward the IRDC G035.39-00.33, believed to be at the earliest stages of evolution. The 13CO and C18O gas is distributed in three filaments (Filaments 1, 2 and 3), where the most massive cores are preferentially found at the intersecting regions between them. The filaments have a similar kinematic structure with smooth velocity gradients of ~0.4-0.8 km s-1 pc-1. Several scenarios are proposed to explain these gradients, including cloud rotation, gas accretion along the filaments, global gravitational collapse, and unresolved sub-filament structures. These results are complemented by HCO+, HNC, H13CO+ and HN13C single-pointing data to search for gas infall signatures. The 13CO and C18O gas motions are supersonic across G035.39-00.33, with the emission showing broader linewidths toward the edges of the IRDC. This could be due to energy dissipation at the densest regions in the cloud. The average H2 densities are ~5000-7000 cm-3, with Filaments 2 and 3 being denser and more massive than Filament 1. The C18O data unveils three regions with high CO depletion factors (f_D~5-12), similar to those found in massive starless cores.

GASEOUS CO ABUNDANCE—AN EVOLUTIONARY TRACER FOR MOLECULAR CLOUDS

Planck cold clumps are among the most promising objects to investigate the initial conditions of the evolution of molecular clouds. In this work, by combing the dust emission data from the survey of Planck satellite with the molecular data of $^{12}$CO/$^{13}$CO (1-0) lines from observations with the Purple Mountain Observatory (PMO) 13.7 m telescope, we investigate the CO abundance, CO depletion and CO-to-H$_{2}$ conversion factor of 674 clumps in the early cold cores (ECC) sample. The median and mean values of the CO abundance are 0.89$\times10^{-4}$ and 1.28$\times10^{-4}$, respectively. The mean and median of CO depletion factor are 1.7 and 0.9, respectively. The median value of $X_{CO-to-H_{2}}$ for the whole sample is $2.8\times10^{20}$ cm$^{-2}$K$^{-1}$km$^{-1}$ s. The CO abundance, CO depletion factor and CO-to-H$_{2}$ conversion factor are strongly (anti-)correlated to other physical parameters (e.g. dust temperature, dust emissivity spectral index, column density, volume density and luminosity-to-mass ratio). To conclude, the gaseous CO abundance can be used as an evolutionary tracer for molecular clouds.

Dynamics, CO depletion, and deuterium fractionation of the dense condensations within the fragmented prestellar core Orion B9–SMM 6

We present APEX observations of C17O(2-1), N2H+(3-2), and N2D+(3-2) towards the subfragments inside the prestellar core SMM 6 in Orion B9. We combined these spectral line data with our previous SABOCA 350-{\mu}m dust continuum map of the source. The subfragments are characterised by subsonic internal non-thermal motions ({\sigma}NT~0.5cs), and most of them appear to be gravitationally bound. The dispersion of the N2H+ velocity centroids among the condensations is very low (0.02 km/s). The CO depletion factors we derive, fD=0.8+/-0.4 - 3.6+/-1.5, do not suggest any significant CO freeze-out but this may be due to the canonical CO abundance we adopt. The fractional abundances of N2H+ and N2D+ with respect to H2 are found to be ~0.9-2.3x10^-9 and ~4.9-9.9x10^-10, respectively. The deuterium fractionation of N2H+, or the N2D+/N2H+ column density ratio, lies in the range 0.30+/-0.07 - 0.43+/-0.09. The detected substructure inside SMM 6 is likely the result of cylindrical Jeans-type gravitational fragmentation. We estimate the timescale for this fragmentation to be ~1.8x10^5 yr. The condensations are unlikely to be able to interact with one another and coalesce before local gravitational collapse ensues. Moreover, significant mass growth of the condensations via competitive-like accretion from the parent core seems unfeasible. The high level of molecular deuteration in the condensations suggests that gas-phase CO should be strongly depleted. It also points towards an advanced stage of chemical evolution. The subfragments of SMM 6 might therefore be near the onset of gravitational collapse, but whether they can form protostellar or substellar objects (brown dwarfs) depends on the local star formation efficiency and remains to be clarified.

Deuterium Fractionation and Ionization Degree in Massive Protostellar/cluster Cores

AbstractWe have conducted a survey of deuterium fractionation of N2H+, RD(N2H+) ≡ N(N2D+)/N(N2H+), with the Arizona Radio Observatory (ARO) Submillimeter Telescope (SMT) to assess the use of RD(N2H+) as an evolutionary tracer among massive protostellar/cluster cores in early stages. Our sample includes 32 dense cores in various evolutionary stages, from high-mass starless cores (HMSCs), high-mass protostellar objects (HMPOs), to ultra-compact (UC) HII regions, in infrared dark clouds (IRDCs) and high infrared extinction clouds. The results show a decreasing trend in deuterium fractionation with evolutionary stage traced by gas temperature and line width (Fig. 1). A moderate increasing trend of deuterium fractionation with the CO depletion factor is also found among cores in IRDCs and HMSCs. These suggest a general chemical behavior of deuterated species in low- and high-mass protostellar candidates. Upper limits to the ionization degree are also estimated to be in the range of 4 × 10−8 − 5 × 10−6.

High CO depletion in southern infrared dark clouds

Infrared-dark high-mass clumps are among the most promising objects to study the initial conditions of the formation process of high-mass stars and rich stellar clusters. In this work, we have observed the (3-2) rotational transition of C18O with the APEX telescope, and the (1,1) and (2,2) inversion transitions of NH3 with the Australia Telescope Compact Array in 21 infrared-dark clouds already mapped in the 1.2 mm continuum, with the aim of measuring basic chemical and physical parameters such as the CO depletion factor (fD), the gas kinetic temperature and the gas mass. In particular, the C18O (3-2) line allows us to derive fD in gas at densities higher than that traced by the (1-0) and (2-1) lines, typically used in previous works. We have detected NH3 and C18O in all targets. The clumps possess mass, H2 column and surface densities consistent with being potentially the birthplace of high-mass stars. We have measured fD in between 5 and 78, with a mean value of 32 and a median of 29. These values are, to our knowledge, larger than the typical CO depletion factors measured towards infrared-dark clouds and high-mass dense cores, and are comparable to or larger than the values measured in low-mass pre-stellar cores close to the onset of the gravitational collapse. This result suggests that the earliest phases of the high-mass star and stellar cluster formation process are characterised by fD larger than in low-mass pre-stellar cores. Thirteen out of 21 clumps are undetected in the 24 {\mu}m Spitzer images, and have slightly lower kinetic temperatures, masses and H2 column densities with respect to the eight Spitzer-bright sources. This could indicate that the Spitzer-dark clumps are either less evolved or are going to form less massive objects.

A (sub)millimetre study of dense cores in Orion B9

We aim to further constrain the properties and evolutionary stages of dense cores in Orion B9. The central part of Orion B9 was mapped at 350 micron with APEX/SABOCA. A sample of nine cores in the region were observed in C17O(2-1), H13CO+(4-3) (towards 3 sources), DCO+(4-3), N2H+(3-2), and N2D+(3-2) with APEX/SHFI. These data are used in conjunction with our previous APEX/LABOCA 870-micron dust continuum data. Many of the LABOCA cores show evidence of substructure in the higher-resolution SABOCA image. In particular, we report on the discovery of multiple very low-mass condensations in the prestellar core SMM 6. Based on the 350-to-870 micron flux density ratios, we determine dust temperatures of ~7.9-10.8 K, and dust emissivity indices of ~0.5-1.8. The CO depletion factors are in the range ~1.6-10.8. The degree of deuteration in N2H+ is ~0.04-0.99, where the highest value (seen towards the prestellar core SMM 1) is, to our knowledge, the most extreme level of N2H+ deuteration reported so far. The level of HCO+ deuteration is about 1-2%. We also detected D2CO towards two sources. The detection of subcondensations within SMM 6 shows that core fragmentation can already take place during the prestellar phase. The origin of this substructure is likely caused by thermal Jeans fragmentation of the elongated parent core. A low depletion factor and the presence of gas-phase D2CO in SMM 1 suggest that the core chemistry is affected by the nearby outflow. The very high N2H+ deuteration in SMM 1 is likely to be remnant of the earlier CO-depleted phase.

DEUTERIUM FRACTIONATION AS AN EVOLUTIONARY PROBE IN MASSIVE PROTOSTELLAR/CLUSTER CORES

Clouds of high infrared extinction are promising sites of massive star/cluster formation. A large number of cloud cores discovered in recent years allows investigation of possible evolutionary sequence among cores in early phases. We have conducted a survey of deuterium fractionation toward 15 dense cores in various evolutionary stages, from high-mass starless cores to ultracompact Hii regions, in the massive star-forming clouds of high extinction, G34.43+0.24, IRAS 18151-1208, and IRAS 18223-1243, with the Submillimeter Telescope (SMT). Spectra of N2H+ (3 - 2), N2D+ (3 - 2), and C18O (2 - 1) were observed to derive the deuterium fractionation of N2H+, Dfrac \equiv N(N2D+)/N(N2H+), as well as the CO depletion factor for every selected core. Our results show a decreasing trend in Dfrac with both gas temperature and linewidth. Since colder and quiescent gas is likely to be associated with less evolved cores, larger Dfrac appears to correlate with early phases of core evolution. Such decreasing trend resembles the behavior of Dfrac in the low-mass protostellar cores and is consistent with several earlier studies in high-mass protostellar cores. We also find a moderate increasing trend of Dfrac with the CO depletion factor, suggesting that sublimation of ice mantles alters the competition in the chemical reactions and reduces Dfrac. Our findings suggest a general chemical behavior of deuterated species in both low- and high-mass proto-stellar candidates at early stages. In addition, upper limits to the ionization degree are estimated to be within 2 \times 10^-7 and 5 \times 10^-6. The four quiescent cores have marginal field-neutral coupling and perhaps favor turbulent cooling flows.

Survey of ortho-H$\mathsf{_{2}}$D$\mathsf{^{+}}$ (1$\mathsf{_{1,0}}$–1$\mathsf{_{1,1}}$) in dense cloud cores

Aims. We present a survey of the ortho-H2D + (11,0−11,1) line toward a sample of 10 starless cores and 6 protostellar cores, carried out at the Caltech Submillimeter Observatory. The high diagnostic power of this line is revealed for the study of the chemistry, and the evolutionary and dynamical status of low-mass dense cores. Methods. The derived ortho-H2D + column densities (N(ortho-H2D + )) are compared with predictions from simple chemical models of centrally concentrated cloud cores. Results. The line is detected in 7 starless cores and in 4 protostellar cores. N(ortho-H2D + ) ranges between 2 and 40 × 10 12 cm −2 in starless cores and between 2 and 9 × 10 12 cm −2 in protostellar cores. The brightest lines are detected toward the densest and most centrally concentrated starless cores, where the CO depletion factor and the deuterium fractionation are also largest. The large scatter observed in plots of N(ortho-H2D + ) vs. the observed deuterium fractionation and vs. the CO depletion factor is likely to be due to variations in the ortho-to-para (o/p) ratio of H2D + from >0.5 for Tkin < 10 K gas in pre-stellar cores to � 0.03 (consistent with Tkin � 15 K for protostellar cores). The two Ophiuchus cores in our sample also require a relatively low o/p ratio (� 0.3). Other parameters, including the cosmic-ray ionization rate, the CO depletion factor (or, more in general, the depletion factor of neutral species), the volume density, the fraction of dust grains and PAHs also largely affect the ortho-H2D + abundance. In particular, gas temperatures above 15 K, low CO depletion factors and large abundance of negatively charged small dust grains or PAHs drastically reduce the deuterium fractionations to values inconsistent with those observed toward pre-stellar and protostellar cores. The most deuterated and H2D + -rich objects (L 429, L 1544, L 694-2 and L 183) are reproduced by chemical models of centrally concentrated (central densties � 10 6 cm −3 ) cores with chemical ages between 10 4 and 10 6 yr. Upper limits of the para-H3O + (1 − −2 + ) and para-D2H + (11,0−10,1) lines are also given. The upper limit to the para-H3O + fractional abundance is � 10 −8 and we find an upper limit to the para-D2H + /orthoH2D + column density ratio equal to 1, consistent with chemical model predictions of high density (2 × 10 6 cm −3 ) and low temperature (Tkin < 10 K) clouds. Conclusions. Our results point out the need for better determinations of temperature and density profiles in dense cores as well as for observations of para-H2D + .

TMC‐1C: An Accreting Starless Core

We have mapped the starless core TMC-1C in a variety of molecular lines with the IRAM 30m telescope. High density tracers show clear signs of self-absorption and sub-sonic infall asymmetries are present in N2H+ (1-0) and DCO+ (2-1) lines. The inward velocity profile in N2H+ (1-0) is extended over a region of about 7,000 AU in radius around the dust continuum peak, which is the most extended ``infalling'' region observed in a starless core with this tracer. The kinetic temperature (~12 K) measured from C17O and C18O suggests that their emission comes from a shell outside the colder interior traced by the mm continuum dust. The C18O (2-1) excitation temperature drops from 12 K to ~10 K away from the center. This is consistent with a volume density drop of the gas traced by the C18O lines, from ~4x10^4 cm^-3 towards the dust peak to ~6x10^3 cm^-3 at a projected distance from the dust peak of 80 (or 11,000 AU). The column density implied by the gas and dust show similar N2H+ and CO depletion factors (f_D < 6). This can be explained with a simple scenario in which: (i) the TMC-1C core is embedded in a relatively dense environment (H2 ~10^4 cm^-3), where CO is mostly in the gas phase and the N2H+ abundance had time to reach equilibrium values; (ii) the surrounding material (rich in CO and N2H+) is accreting onto the dense core nucleus; (iii) TMC-1C is older than 3x10^5 yr, to account for the observed abundance of N2H+ across the core (~10^-10 w.r.t. H2); and (iv) the core nucleus is either much younger (~10^4 yr) or ``undepleted'' material from the surrounding envelope has fallen towards it in the past 10,000 yr.

Searching for massive pre-stellar cores through observations of N$_\mathsf{2}$H+and N$_\mathsf{2}$D+

Aims. We have measured the deuterium fractionation and the CO depletion factor (ratio between expected and observed CO abundance) in a sample of high-mass protostellar candidates, in order to understand whether the earliest evolutionary stages of high-mass stars have chemical characteristics similar to those of low-mass ones. It has been found that low-mass starless cores on the verge of star formation have large values both of the column density ratio $N({\rm N_2D^+})/N({\rm N_2H^+})$ and of the CO depletion factor. Methods. With the IRAM-30 m telescope and the JCMT we have observed two rotational lines of N 2 H + and N 2 D + , the (2–1) line of C 17 O and DCO + , and the sub-millimeter continuum towards a sample of 10 high-mass protostellar candidates. Results. We have detected N 2 D + emission in 7 of the 10 sources of our sample, and found an average value $N({\rm N_2D^+})/N({\rm N_2H^+})\sim 0.015$. This value is ~3 orders of magnitude larger than the interstellar D/H ratio, indicating the presence of cold and dense gas, in which the physical-chemical conditions are similar to those observed in low-mass pre-stellar cores. The integrated CO depletion factors show that in the majority of the sources the expected CO abundances are larger than the observed values, with a median ratio of 3.2. Conclusions. In principle, the cold gas that generates the N 2 D + emission can be the remnant of the massive molecular core in which the high-mass (proto-)star was born, not yet heated up by the central object. If so, our results indicate that the chemical properties of the clouds in which high-mass stars are born are similar to their low-mass counterparts. Alternatively, this cold gas could be located in one (or more) starless core (cores) near the protostellar object. Due to the poor angular resolution of our data, we cannot distinguish between the two scenarios.