No Arabic abstract
We have measured the deuterium fractionation (through the column density ratio N(N2D+)/N(N2H+)) and the CO depletion factor (ratio between expected and observed CO abundance) in a sample of 10 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. The observations were carried out with the IRAM-30m telescope and the JCMT. We have detected N2D+ emission in 7 of the 10 sources of our sample, and found an average value N(N2D+)/N(N2H+)~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. Also, 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. In principle, the cold gas that gives origin to the N2D+ 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 can be located into one (or more) starless core (cores) near the protostellar object. Due to the poor angular resolution of our data, we cannot decide which is the correct scenario.
We have undertaken a survey of N2H+ and N2D+ towards 31 low-mass starless cores using the IRAM 30m telescope. Our main objective has been to determine the abundance ratio of N2D+ and N2H+ towards the nuclei of these cores and thus to obtain estimates of the degree of deuterium enrichment, a symptom of advanced chemical evolution according to current models. We find that the N(N2D+)/N(N2H+) ratio is larger in more centrally concentrated cores with larger peak H2 and N2H+ column density than the sample mean. The deuterium enrichment in starless cores is presently ascribed to depletion of CO in the high density (> 3*10^4 cm-3) core nucleus. To substantiate this picture, we compare our results with observations in dust emission at 1.2 mm and in two transitions of C18O. We find a good correlation between deuterium fractionation and N(C18O)/N(H2) for the nuclei of 14 starless cores. We, thus, identified a set of properties that characterize the most evolved, or pre-stellar, starless cores. These are: higher N2H+ and N2D+ column densities, higher N(N2D+)/N(N2H+), more pronounced CO depletion, broader N2H+ lines with infall asymmetry, higher central H2 column densities and a more compact density profile than in the average core. We conclude that this combination of properties gives a reliable indication of the evolutionary state of the core. Seven cores in our sample (L1521F, OphD, L429, L694, L183, L1544 and TMC2) show the majority of these features and thus are believed to be closer to forming a protostar than are the other members of our sample. Finally, we note that the subsample of Taurus cores behaves more homogeneously than the total sample, an indication that the external environment could play an important role in the core evolution.
Chemical models predict that the deuterated fraction (the column density ratio between a molecule containing D and its counterpart containing H) of N2H+, Dfrac(N2H+), is high in massive pre-protostellar cores and rapidly drops of an order of magnitude after the protostar birth, while that of HNC, Dfrac(HNC), remains constant for much longer. We tested these predictions by deriving Dfrac(HNC) in 22 high-mass star forming cores divided in three different evolutionary stages, from high-mass starless core candidates (HMSCs, 8) to high-mass protostellar objects (HMPOs, 7) to Ultracompact HII regions (UCHIIs, 7). For all of them, Dfrac (N2H+) was already determined through IRAM-30m Telescope observations, which confirmed the theoretical rapid decrease of Dfrac(N2H+) after protostar birth (Fontani et al. 2011). Therefore our comparative study is not affected by biases introduced by the source selection. We have found average Dfrac(HNC) of 0.012, 0.009 and 0.008 in HMSCs, HMPOs and UCHIIs, respectively, with no statistically significant differences among the three evolutionary groups. These findings confirm the predictions of the chemical models, and indicate that large values of Dfrac(N2H+) are more suitable than large values of Dfrac(HNC) to identify cores on the verge of forming high-mass stars, likewise what found in the low-mass regime.
Context : Dynamical studies of prestellar cores search for small velocity differences between different tracers. The highest radiation frequency precision is therefore required for each of these species. Aims : We want to adjust the frequency of the first three rotational transitions of N2H+ and N2D+ and extrapolate to the next three transitions. Methods : N2H+ and N2D+ are compared to NH3 the frequency of which is more accurately known and which has the advantage to be spatially coexistent with N2H+ and N2D+ in dark cloud cores. With lines among the narrowests, and N2H+ and NH3 emitting region among the largests, L183 is a good candidate to compare these species. Results : A correction of ~10 kHz for the N2H+ (J:1-0) transition has been found (~0.03 km/s) and similar corrections, from a few m/s up to ~0.05 km/s are reported for the other transitions (N2H+ J:3-2 and N2D+ J:1-0, J:2-1, and J:3-2) compared to previous astronomical determinations. Einstein spontaneous decay coefficients (Aul) are included.
High levels of deuterium fractionation of $rm N_2H^+$ (i.e., $rm D_{frac}^{N_2H^+} gtrsim 0.1$) are often observed in pre-stellar cores (PSCs) and detection of $rm N_2D^+$ 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 $rm N_2H^+$ and $rm N_2D^+$ observations of PSCs. We perform 3D magnetohydrodynamics simulations of a massive, turbulent, magnetised PSC, coupled with a sophisticated deuteration astrochemical network. Although the core has some magnetic/turbulent support, it collapses under gravity in about one freefall time, which marks the end of the simulations. Our fiducial model achieves relatively low $rm D_{frac}^{N_2H^+} sim 0.002$ during this time. We then investigate effects of initial ortho-para ratio of $rm H_2$ ($rm 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 $rm D_{frac}^{N_2H^+}$ and absolute abundances to match observational values within one freefall time. For $rm OPR^{H_2}$, while a lower initial value helps the growth of $rm D_{frac}^{N_2H^+}$, 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 $rm N_2D^+$ 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 sub-virial based on its $rm N_2D^+$ velocity dispersion.
High levels of deuterium fraction in N$_2$H$^+$ are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values ($D_{rm frac}^{{rm N}_2{rm H}^+} equiv {rm N}_2{rm D}^+/{rm N}_2{rm H}^+ gtrsim 0.1$) is longer than the free-fall time, possibly ten times longer. Here, we explore the deuteration of turbulent, magnetized cores with 3D magnetohydrodynamics simulations. We use an approximate chemical model to follow the growth in abundances of N$_2$H$^+$ and N$_2$D$^+$. We then examine the dynamics of the core using each tracer for comparison to observations. We find that the velocity dispersion of the core as traced by N$_2$D$^+$ appears slightly sub-virial compared to predictions of the Turbulent Core Model of McKee & Tan, except at late times just before the onset of protostar formation. By varying the initial mass surface density, the magnetic energy, the chemical age, and the ortho-to-para ratio of H$_2$, we also determine the physical and temporal properties required for high deuteration. We find that low initial ortho-to-para ratios ($lesssim 0.01$) and/or multiple free-fall times ($gtrsim 3$) of prior chemical evolution are necessary to reach the observed values of deuterium fraction in pre-stellar cores.