EDP Sciences
Free Access
Issue
A&A
Volume 621, January 2019
Article Number L7
Number of page(s) 7
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/201834602
Published online 09 January 2019

© ESO 2019

1. Introduction

Depletion of heavy elements proceeds progressively in dense and cold material (nH2 > few × 104 cm−3, T ≲ 25 K), rendering heavy molecular species unobservable, and thus not viable to trace gas in the earliest phases of the star-formation process. It has been demonstrated that deuterated species become more abundant under these conditions, and that there is a correlation between the degree of deuteration and the depletion of CO (Caselli et al. 1999, 2002; Bacmann et al. 2003). Large values for the deuterium fraction are also found towards protostellar sources (e.g. Ceccarelli et al. 1998; Parise et al. 2002).

Pivotal to the deuteration process is the exothermic reaction chain (e.g. Walmsley et al. 2004; Flower et al. 2004). H2D+, D2H+ , and react with other species, producing the vastly enhanced deuterium fractions in molecules at the core of the ion-neutral chemistry. N2 in particular can react with H2D+ to form N2D+ via

(1)

and more efficiently with D2H+ and . The molecule N2D+ is frequently used as a reliable tracer for gas affected by heavy CO depletion (e.g. Caselli et al. 2002; Barnes et al. 2016). Because H2D+ and N2D+ are chemically linked, and because they are abundant only under specific circumstances, they could be used to determine the evolutionary stage of a source in early phases (e.g. Caselli et al. 2008; Emprechtinger et al. 2009; Fontani et al. 2015), something of particular interest in the high-mass regime.

Pillai et al. (2012) have obtained maps for part of the DR21 complex of the o-H2D+(11, 0 − 11, 1) and N2D+(3–2) transitions, observed with the James-Clerk-Maxwell Telescope (JCMT) and the Submillimeter Array (SMA), respectively. They find very extended o-H2D+ emission as do Vastel et al. (2006) in the low-mass source L1544. They also find that this species mainly traces gas that is not seen in dust continuum emission or in the interferometric N2D+ data. H2D+ may therefore be sensitive to gas that would elude detection in the most commonly used tracers, and can represent an even earlier stage in the process of star formation.

These results raise two important questions that we try to address in this work: Firstly, is the lack of correlation between o-H2D+(11, 0 − 11, 1) and N2D+(3–2) found by Pillai et al. (2012) real, or is it an effect of interferometric filtering? Secondly, how do the abundances of these molecules evolve with time, and why? In this Letter we first describe our observations, present the results, and finally discuss the importance of our findings.

2. Source selection and observations

G351.77–0.51 (hereafter G351) is the closest and most massive filament (D = 1 kpc, M ∼ 2000 M, Leurini et al. 2019) identified in the 870 μm APEX Telescope Large Survey of the Galaxy (ATLASGAL) (Schuller et al. 2009). A three colour image of the region is shown in Fig. 1, emphasising the large amount of dense and cold gas in its massive ridge.

thumbnail Fig. 1.

Three-colour image (red: ATLASGAL 870 μm, green: Hi-GAL 250 μm, blue: Hi-GAL 70 μm) of the G351 complex. Contours of the ATLASGAL images are also indicated in grey (0.15, 0.5, 1.5, 5 Jy beam−1). Clumps 2, 5, and 7 are indicated in red. In the top-left corner we list the clumps from the least- to the most evolved, and the species we observed in each source.

Open with DEXTER

G351 offers several advantages for our purposes: 1) it is close-by, 2) it hosts massive clumps that exceed the Kauffmann & Pillai (2010) threshold for high-mass star formation, and 3) the clumps are in different evolutionary stages. Therefore, sources along the spine of the G351 filament are extremely well suited to studying the evolution of chemistry in high-mass clumps, not only because they are nearby and all at the same distance, but also because they have similar properties and share the same initial chemical conditions.

We select three clumps in this complex (Clumps 2, 5, and 7, following Leurini et al. 2011) in different evolutionary stages, based on their IR properties: Clump 7 is the least evolved, being still quiescent at 70 μm; Clumps 5 and 2 are both bright at this wavelength, with the former being weaker than the latter (Fig. 1). These sources are the ones with the largest column densities and masses found in G351 in each evolutionary stage. Their mass and peak column densities are also within a factor of two, using a temperature of 25 K for Clump 2, and 10 K for Clumps 5 and 7 (cf. Table 4 in Leurini et al. 2011).

The dust emission peak for each of the three clumps has been targeted by APEX 12m submillimeter telescope observations (project M-099.F-9508A), using the PI230 and FLASH+ receivers, to cover the N2D+(3–2) and o-H2D+(11, 0 − 11, 1) lines at 231.3216 and 372.42134 GHz, respectively. Observations were performed between 2017 May 31 and 2017 September 19, and completed on 2018 July 1–2. The rms noise on the main-beam brightness temperature scale, TMB, is ∼20 mK at 231 GHz and 30 mK at 372 GHz, for a spectral resolution of 0.45 Km s−1. We converted the antenna temperature to TMB, using ηMB = 0.60 for o-H2D+(11, 0 − 11, 1) and ηMB = 0.69 for N2D+(3–2).

3. Results

In our APEX observations Clumps 5 and 7 were detected in o-H2D+(11, 0 − 11, 1), while Clumps 2 and 5 were detected in N2D+(3–2). Figure 2 shows the spectra for all sources observed, highlighting their different features in line emission.

thumbnail Fig. 2.

Spectra of the three clumps for o-H2D+(11, 0 − 11, 1) (left column) and N2D+(3–2) (right column). The best fit is indicated in red, while upper limits are drawn in blue.

Open with DEXTER

To give a physical explanation of this behaviour, we first characterise the properties of the clumps.

3.1. Dust spectrum

We perform aperture photometry on the clumps to extract dust continuum fluxes. We use images from ATLASGAL (870 μm Schuller et al. 2009), HiGAL (350 μm, 250 μm, 160 μm and 70 μm Molinari et al. 2010), MIPSGAL (24 μm Carey et al. 2009) and MSX (Egan et al. 2003). All images were smoothed to 28″, the resolution of the APEX 230 GHz observations, and the fluxes were extracted in the central beam, after removing the median background estimated in an annulus with 2 × θbeam <  Rannulus <  5 × θbeam, for each clump. The apertures and annuli are centred on the positions of the clumps, as determined in Leurini et al. (2011).

The dust temperature and column density are estimated via a greybody fit at λ ≥ 70 μm. The bolometric luminosity is computed as the sum of the integral of the best-fit greybody curve and the integral at λ ≤ 70 μm, computed with the trapezoidal rule in log-log space. For the greybody we use κ870 μm = 1.85 cm2g−1 and β = 1.75, and consider the 70 μm flux as an upper limit.

To convert the dust column density to that of molecular hydrogen, we use a gas-to-dust ratio γ = 120, as calculated from Eq. (2) in Giannetti et al. (2017a), and RGC = 7.4 kpc (Leurini et al. 2019). The impression that Clump 7 is the least evolved of the three sources considered, and that Clump 2 is the most evolved one is confirmed by Td and the L/M ratio (see Table 1), two efficient indicators of evolution (e.g. Saraceno et al. 1996; Molinari et al. 2008; König et al. 2017; Urquhart et al. 2018).

Table 1.

Properties of the clumps derived from dust continuum emission.

The average volume density along the line of sight listed in Table 1 is estimated as nH2 = NH2/ℓ, where ℓ is the size of the clump obtained by Leurini et al. (2011). We note however that the filament width is nearly constant (∼0.2 pc; Leurini et al. 2019), even at the location of the clumps, so the difference in mean volume density may be overestimated.

3.2. Molecular column densities and abundances

To calculate the column densities, a spectral-line fit is performed, under the assumption of local thermodynamic equilibrium (LTE) (cf. Vastel et al. 2012), with MCWeeds (Giannetti et al. 2017b). This allows the uncertainty on this parameter to be obtained as well (Table 2). We use the partition function for o-H2D+ (in the relevant temperature range, 9.375 K: 10.3375, 18.750 K: 12.5068, 37.500 K: 15.5054) from CDMS (Müller et al. 2001). Because the observed clumps are sufficiently dense to attain thermal coupling between gas and dust (cf. Table 1), and because we lack an estimate of the excitation temperature for o-H2D+ (11, 0 − 11, 1) and N2D+(3–2), we use Td, assuming that Td = Tex. The fit is performed using Monte Carlo Markov chains. We used an adaptive Metropolis-Hastings sampler (Haario et al. 2001), with 100 000 total samples, a burn-in period and a delay for the adaptive sampling of 10 000 and 5000 iterations, respectively, and a thinning factor of 20 (see Giannetti et al. 2017b, for more details on these parameters). Convergence and independence of the samples are ensured with the Raftery-Lewis (Raftery & Lewis 1995), Geweke (Geweke 1992a) and Gelman-Rubin (Geweke 1992b) tests.

Table 2.

Results of the line-fitting procedure with MCWeeds.

We include in the budget a Gaussian calibration uncertainty with σ = 5%. Multiple tests have been performed on the priors to make sure that the choice of the latter is not crucial for the results.

Beam-averaged abundances can be estimated from the peak column density of H2 derived from the dust continuum emission. In the following, the same value of NH2 is used for both o-H2D+ and N2D+, despite the difference in the angular resolution of the spectral-line observations (18–28″, respectively). Smoothing the continuum maps to the resolution of the ATLASGAL data (excluding the HiGAL 350 μm image) leads to an increase in the peak H2 column density of ≲15%, well within the uncertainties; this has the effect of lowering X(H2D+) in Clump 5, thus increasing the difference with Clump 7 (see Fig. 3, left).

thumbnail Fig. 3.

Left panel: calculated o-H2D+ (blue) and N2D+ (orange) abundances. Right panel: relative abundance of the two species. The clumps are ordered by evolutionary stage, from the least- to the most evolved. The 95% credible interval is indicated for abundances and ratios; we point out that for non-detections, only the upper limit is shown.

Open with DEXTER

If X(o−H2D+) does not vary significantly over scales of 18000-28000 AU (18–28″ at 1 kpc), one can also estimate the relative abundance of the two species. In Fig. 3 one can see that o-H2D+ becomes rarer with evolution, while N2D+ shows the opposite behaviour. Considering the abundance variation and its uncertainty, o-H2D+ is much more sensitive to the clump evolution than N2D+. The relative abundance (panel b) of the two species offers another way of looking at the same finding, but with the advantage of being independent of the H2 column density, then removing an additional source of uncertainty. From the figure, we see that the o-H2D+/N2D+ ratio progressively decreases in time and o-H2D+ is more abundant than N2D+ by a factor ≳10 in the youngest Clump 5 and 7.

4. Discussion

In the previous section we show that o-H2D+ and N2D+ have opposite trends in abundance as a function of evolution. In the following we discuss a few possible explanations for the abundance differences between the three clumps.

4.1. Deuteration of and availability of gas-phase N2

H2D+ rapidly forms in cold and dense gas, where the reaction

(2)

only proceeds in the forward direction, unless there is a substantial fraction of ortho-H2 (e.g. Gerlich et al. 2002), the only case in which ΔE can be negative. H2D+ then drives the deuteration of the neutral molecules that have a deuteron affinity larger than molecular hydrogen (Millar et al. 1989). The o-H2D+ abundance derived in Clumps 5 and 7 is of the order of ∼3 × 10−11, similar to the values reported by Pillai et al. (2012) for their o-H2D+ peaks in the DR 21 region, and by Vastel et al. (2006) for the outskirts of L1544. For the temperatures of our sources, 13 K and 15.5 K for Clumps 7 and 5, respectively, the ortho-to-para ratio (OPR) predicted by Flower et al. (2004) for H2D+ is in the range 0.05–0.1. Recently, we also performed 3D magnetohydrodynamic simulations of collapsing clumps and filaments, which included chemistry under the assumption of total depletion (Körtgen et al. 2017, 2018). In these simulations we also follow the OPR for H2D+, finding that it rapidly decreases to values below unity, in agreement with Flower et al. (2004), on time scales of the order of 105 yr. This OPR can be used to estimate the total H2D+ abundance, which would be in the range 10−10 − 10−9. Caselli et al. (2003) show how in L1544, H2D+ is nearly as abundant as electrons (10−9 vs. 2 × 10−9), indicating that this species is a major ion, which only happens when the depletion is extremely high. Our results, when compared to these findings, also suggest that regions exist in the clump where heavy elements are vastly depleted, which is also supported by our ongoing study of CO depletion for the entire filament (Sabatini et al., in prep.).

Clump 7, however, not only has a high abundance of o-H2D+, but also shows low values of X(N2D+). One possible explanation for this is the time lag that is needed to form N2D+ from H2D+. As an example, Fig. 1 in Sipilä et al. (2015) shows that in the first 105 yr, according to their model, o-H2D+ is much more abundant than N2D+, and the difference progressively decreases. This provides, qualitatively, an upper limit to the age of the clump of 105 years, comparable to the current free-fall time (tff ∼ 1.2 × 105 yr). The clump is indeed sub-virial (α ∼ 0.41), and shows a very weak stellar activity in its centre. Clump 7 is also the only one to show a blue-skewed profile in HCO+(3–2), which indicates that the clump is collapsing (see Fig. A.1). We note that if the H2 ortho-to-para conversion occurs on the surface of dust grains, the chemical timescale could be even shorter (Bovino et al. 2017).

In Clump 5, on the other hand, X(N2D+) is two to three times larger than in Clump 7. In this clump, being more evolved, N2D+ has had time to form, while the stellar activity is not yet sufficient to warm up the gas to temperatures high enough to significantly alter the chemistry on large scales.

Interestingly, Clump 2 still has a large fraction of N2D+ with respect to H2, but at the same time X(H2D+) is at least an order of magnitude lower than X(N2D+). Because N2D+ forms via deuteron exchange from H2D+, one would expect to see a correlation in the abundance of the two molecules, after they have had time to form. Furthermore, Clump 2 hosts luminous young stellar objects (YSOs), possibly close to the zero-age main sequence (ZAMS), as indicated by its L/M ratio (Giannetti et al. 2017b), and is therefore a relatively evolved source.

What is keeping the level of N2D+ high, even higher than in Clump 5, and at the same time destroying o-H2D+? The beam- and line-of-sight-averaged depletion of CO does not vary significantly from Clump 7 to Clump 2, and a relatively high depletion is one of the key ingredients for the efficient formation of N2D+ because carbon monoxide acts as a destroyer of both H2D+ and N2D+ (e.g. Walmsley et al. 2004; Emprechtinger et al. 2009). The other fundamental ingredient for the formation of N2D+ is the presence of N2 in the gas phase, which can react with H2D+, D2H+, and . Close to the YSOs, where the temperature is high and CO is in the gas phase, H2D+ and N2D+ are both actively destroyed and Eq. (2) is also fast in the backward direction. However, at larger distances from luminous YSOs, the material is still dense and cold, and CO remains frozen out onto dust grains, contrary to N2, either because the latter has not been significantly depleted (due to its longer formation timescale compared to CO and the fact that atomic nitrogen could have a lower sticking coefficient than N2Flower et al. 2006), or because it evaporates faster. In fact, according to the most recent estimates (Wakelam et al. 2017) N2 and CO binding energies differ by a few hundred Kelvin (1100 K vs. 1300 K), enough to cause a significant difference in the evaporation timescales. The latter, evaluated by employing the standard , are ∼5 × 108 yr, and ∼2 × 104 yr, for CO and N2, respectively. In the equation above, ν0 = 1012 s−1 is the typical harmonic frequency (Hasegawa et al. 1992), ED is the binding energy, kB is the Boltzmann’s constant, and Td is the dust temperature (see Table 1). More importantly, with time, a large amount of H2D+ is converted to D2H+ and , decreasing its abundance. N2D+, on the other hand, is more efficiently formed by the reaction of D2H+ and with N2, compared to H2D+. The combination of these processes could succeed in reproducing the observed reduction in abundance of o-H2D+, and in maintaining a high column density of N2D+. In even more evolved sources, X(N2D+) decreases, as expected (Fontani et al. 2015).

If this is the case, we propose that the combined abundances of o-H2D+ and N2D+ represent a more efficient evolutionary indicator for the first stages of the high-mass star-formation process, compared to their individual values (e.g. Caselli et al. 2008; Fontani et al. 2015), in the same way that they are both needed to trace the full reservoir of gas (Pillai et al. 2012).

4.2. Alternative explanations

One alternative possibility to explain the emission in Clump 7 could be that the region in the clump where heavy elements (including N2) are completely depleted is extended, strongly affecting the abundance measurement. Contamination in the beam and along the line-of-sight would explain the weak emission of high-density tracers towards this source (see also Leurini et al. 2011, for N2H+ J = 1 → 0).

The critical density of the o-H2D+(11, 0 − 11, 1) is ∼105 cm−3, roughly one order of magnitude lower than that of N2D+(3–2), and similar to the mean values of n(H2) reported in Table 1. Although this is not always a direct indication that emission is coming from denser gas (see e.g. Kauffmann et al. 2017), it is possible that in Clump 7 the region dense enough to excite N2D+(3–2) is smaller than in the other two clumps.

In both cases, the abundance of o-H2D+ relative to N2D+ keeps its potential as an evolutionary indicator, on the one hand because if a region of complete depletion exists at the centre of the clump, the warm-up by the YSOs will progressively erode it, and on the other hand because density increases as a result of the clump collapse, before feedback becomes important at the clump scale and starts to dissipate it.

For Clumps 2 and 5, the N2D+(3–2) emission is likely to arise in cold gas surrounding the star-forming cores, as observed by Fontani et al. (2009) for example. If o-H2D+ is abundant in these compact structures, its emission may be completely saturated and diluted in the APEX beam for Clump 2; o-H2D+(11, 0 − 11, 1) must also be weak in the lower-density gas enveloping these cores. This is in contrast with the observations of extended emission in DR 21, and those of L1544 (Vastel et al. 2006), although the envelope of high-mass star forming regions could be warmer due to the stronger external illumination from close-by high-mass stars; this increases the ortho-to-para H2 ratio, which suppresses D-fractionation (e.g. Kong et al. 2015). Additional high-resolution observations are needed to investigate these possibilities.

5. Summary and conclusions

We observed the o-H2D+(11, 0 − 11, 1) and N2D+(3–2) lines in three clumps along the spine of G351, the most massive filament within 1 kpc (Leurini et al. 2011) in ATLASGAL. These observations provide for the first time the possibility to investigate the variation with evolution of the abundances of o-H2D+ and N2D+ in high-mass sources, from a clump that is still quiescent at 70 μm, to one which hosts luminous YSOs, likely close to the ZAMS. The selected clumps not only belong to the same complex, ensuring the best conditions for a comparison of the chemistry, but also have a comparable peak column density of H2.

The abundance of N2D+ progressively increases with evolution, while an opposite trend is found for o-H2D+. We propose that the chemical evolution of the clumps causes this behaviour, and that the relative abundance of these species is a good indicator of time evolution.

First, after the clump reaches a high-enough density to allow for an efficient CO depletion, H2D+ starts to form, but to attain a considerable amount of N2D+ would require a longer time (e.g. Sipilä et al. 2015). Clump 7 still has X(N2D+)≲9 × 10−13, which also suggests an age of no more than 105 yr, and a fast collapse time. As time passes, more N2D+ forms via reaction (1) in cold and dense gas. High degrees of deuteration are commonly observed in relatively young and massive clumps (e.g. Caselli et al. 2002; Fontani et al. 2011). This is the situation of Clump 5, where both H2D+ and N2D+ are abundant, despite the presence of YSOs (Fig. 1), which are not yet luminous enough to warm up a significant portion of the gas in the clump. In Clump 2, X(N2D+) increases by a factor of ∼2 and X(H2D+) decreases by a factor ≳10 compared to Clump 5. Because this source is the most evolved, N2 has had more time to form. Also, the line-of-sight and beam-averaged dust temperature in Clump 2 is ∼20 K, not high enough to release CO into the gas phase on large scales (see Sabatini et al., in prep.), but sufficient to reduce the evaporation timescale of N2 to ∼2 × 104 yr. Therefore, N2 is abundant in the gas phase throughout the clump, whether or not it was depleted from the gas phase earlier. H2D+ may be progressively transformed into D2H+ and , reducing the abundance of H2D+, and, at the same time, boosting the production of N2D+. This behaviour could be confirmed with observations of p-D2H+.

A weak N2D+(3–2) line in Clump 7 may alternatively be caused by excitation effects, because it has a critical density an order of magnitude larger than the o-H2D+ line (see also Pillai et al. 2012) and larger than the mean n(H2) of the source. A similar effect could result from a large part of the gas being affected by complete depletion of heavy elements. In both cases the potential of the relative abundance of o-H2D+ and N2D+ in tracing evolution remains unaltered, because density increases with time, and the region of complete depletion is progressively eroded by the warm-up caused by the YSOs.

Without high-resolution observations, however, we cannot exclude that H2D+ is also abundant in compact N2D+-emitting regions for sources like Clump 2. Here, o-H2D+(11, 0 − 11, 1) could be completely saturated, and dilution in the APEX beam could cause the non-detection of the line.


1

Assuming a homogeneous clump (see MacLaren et al. 1988), the radius from Table 1, and the line width from o-H2D+(11, 0 − 11, 1).

Acknowledgments

SB is financially supported by CONICYT Fondecyt Iniciación (project code 11170268), CONICYT programa de Astronomia Fondo Quimal 2017 QUIMAL170001, and BASAL Centro de Astrofisica y Tecnologias Afines (CATA) AFB-17002. DRGS is financially supported by CONICYT Fondecyt regular (project code 1161247) and the international collaboration project PII20150171. AG thanks the Astronomy Department, University of Concepción for having supported his visit during July 2018 (via Redes Internacionales project number REDI170093). B.K. acknowledges funding from the German Science Foundation (DFG) within the Priority Programm “The Physics of the ISM” (SPP 1573) via the grant BA 3706/3-2. This work was partly supported by the Collaborative Research Council 956, sub-project A6, funded by the Deutsche Forschungsgemeinschaft (DFG). This paper is based on data acquired with the Atacama Pathfinder EXperiment (APEX). APEX is a collaboration between the Max Planck Institute for Radioastronomy, the European Southern Observatory, and the Onsala Space Observatory. This research made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration 2013, http://www.astropy.org), of NASA’s Astrophysics Data System, and of Matplotlib (Hunter 2007). MCWeeds makes use of the PyMC package (Patil et al. 2010).

References

Appendix A: Additional spectra

thumbnail Fig. A.1.

Spectra of the three clumps for HCO+ J = 3 → 2 (black) and H13CO+ J = 3 → 2, multiplied by a factor of three (red).

Open with DEXTER

All Tables

Table 1.

Properties of the clumps derived from dust continuum emission.

Table 2.

Results of the line-fitting procedure with MCWeeds.

All Figures

thumbnail Fig. 1.

Three-colour image (red: ATLASGAL 870 μm, green: Hi-GAL 250 μm, blue: Hi-GAL 70 μm) of the G351 complex. Contours of the ATLASGAL images are also indicated in grey (0.15, 0.5, 1.5, 5 Jy beam−1). Clumps 2, 5, and 7 are indicated in red. In the top-left corner we list the clumps from the least- to the most evolved, and the species we observed in each source.

Open with DEXTER
In the text
thumbnail Fig. 2.

Spectra of the three clumps for o-H2D+(11, 0 − 11, 1) (left column) and N2D+(3–2) (right column). The best fit is indicated in red, while upper limits are drawn in blue.

Open with DEXTER
In the text
thumbnail Fig. 3.

Left panel: calculated o-H2D+ (blue) and N2D+ (orange) abundances. Right panel: relative abundance of the two species. The clumps are ordered by evolutionary stage, from the least- to the most evolved. The 95% credible interval is indicated for abundances and ratios; we point out that for non-detections, only the upper limit is shown.

Open with DEXTER
In the text
thumbnail Fig. A.1.

Spectra of the three clumps for HCO+ J = 3 → 2 (black) and H13CO+ J = 3 → 2, multiplied by a factor of three (red).

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.