Free Access
Issue
A&A
Volume 631, November 2019
Article Number A25
Number of page(s) 12
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201936012
Published online 15 October 2019

© ESO 2019

1 Introduction

How water evolves during star formation, from the molecular cloud core down to the protoplanetary disk, is a key question concerning the origin of the Solar System, the formation of exoplanets, and, ultimately, the possible emergence of life in plane- tary systems.

Water is observed during all stages of star formation: from molecular clouds, through the dense-core phase and the protoplanetary disk and finally in planetary systems such as our Solar System (van Dishoeck et al. 2014). In star-forming pre-stellar molecular clouds, water is predominantly present in the ice: the gas-phase abundance is constrained to be of the order Xgas (H2O) ~ 10−8–10−9 relative to H2 (Bergin & Snell 2002; Caselli et al. 2012) while absorption ice spectroscopy has revealed Xice (H2O) ~ 10−4–10−5 (Pontoppidan et al. 2004). It is an open question what processing, if any, water experiences during star formation. Classically, two scenarios are considered: inheritance or processing. In the inheritance scenario, the bulk of the water present in the planetary system after formation is inherited from the molecular cloud and is therefore representative of the chemistry in the molecular cloud before star formation begins (e.g., Visser et al. 2009; Drozdovskaya et al. 2016). In the alternative scenario, a substantial amount of water is destroyed and reformed in local processes within the envelope or protoplanetary disk during the formation process. In this case the final water chemistry is determined by local processes in the specific system (see, e.g., Cleeves et al. 2014, for an discussion of each scenario).

The deuterium fractionation of water is a useful proxy for the processing of water during star formation and can help distinguish between the two scenarios outlined above. The enrichment of deuterium is driven by several chemical processes, and their efficiency depend on physical conditions such as temperature, density, visual extinction, and ionization sources (Ceccarelli et al. 2014). At low temperatures the prominent pathway is the gas-phase exchange reaction H3+$_{3}^{+}$ + HD H2 D+ + H2 + ΔE, where ΔE ≈ 124 K1. This reaction is effectively one-way at low temperatures (T ≲50 K) due to the endothermicity of the backward reaction. This leads to an enrichment of H2 D+ which subsequently dissociatively recombines with free electrons to form atomic D, thus increasing the local atomic D/H ratio in the gas-phase and ultimately on dust grain surfaces where water and other molecules are formed through hydrogenation.

Measurements of water deuteration through the different stages of star formation can therefore be used to trace the chemical evolution of water from the molecular cloud down to the protoplanetary disk. Observations of water are challenging due to the high abundance of water in the Earths atmosphere which makes the atmosphere opaque to prominent water emission lines. One solution is to observe water from space, as was done with the Herschel Space Observatory (see, e.g., van Dishoeck et al. 2011). Observations from Herschel have greatly expanded our knowledge of the different origins of water emission toward protostars, i.e., cold and warm envelope gas, shocks and outflows (e.g., Kristensen et al. 2012; Coutens et al. 2013; Visser et al. 2013). Another approach is to target rarer water isotopologs such as H218$_{2}^{18}$O or the deuterated isotopologs HDO and D2O. These molecules have transitions which fall outside of the opaque water bands in the atmosphere and can be observed with ground based telescopes. First attempts to constrain the HDO/H2O ratio toward low-mass protostars utilized single-dish telescopes, thus observing a mixture of small- and large-scale emission and suffering from beam dilution. This led to a large discrepancy between measurements obtained with different telescopes and highlighted the importance of high spatial resolution to constrain the origin of the emission (e.g., Stark et al. 2004; Parise et al. 2005; Coutens et al. 2012).

Jørgensen & van Dishoeck (2010) reported the first interferometric determination of the HDO/H2O abundance ratio in the hot corino of NGC 1333 IRAS 4B. Since then water deuteration has been studied toward a number of low-mass Class 0 protostars using interferometers to resolve the water emission in the hot corino where T > 100 K and ice is entirely sublimated off the dust grains (Taquet et al. 2013a; Coutens et al. 2014; Persson et al. 2014). These measurements reveals varying degrees of deuterium fractionation on different spatial scales toward Class 0 protostars. On larger spatial scales, in cold gas, a high degree of deuterium fractionation has been detected with HDO/H2O and D2O/HDO ratios of the order ≿10−2 in the gas-phase (Coutens et al. 2014). Meanwhile, the water emission from the hot corino, i.e., small spatial scales, show lower HDO/H2O ratios in the range ~ 10−4–10−3. In comparison, the measured HDO/H2O ratios in the Solar System range from as low as the local ISM value2 ratio of ~ 4 × 10−5 to cometary values as high as [HDO/H2O] ~ 10−3; the Earth’s D/H ratio, as measured from Vienna Standard Mean Ocean Water, is D/H = 1.557 × 10−4 (de Laeter et al. 2003) which corresponds to HDO/H2O ~ 3 × 10−4.

A possible explanation for the observed variation of the HDO/H2O ratio on different spatial scales is that water and its deuterated version are not well mixed in the ice. The evolution of water deuteration during star formation and the effects of the layered ice structure have been the subject of recent modelling efforts (e.g., Cazaux et al. 2011; Taquet et al. 2013b; Furuya et al. 2016). Combining chemical models and the available observations of the HDO/H2O and D2O/HDO ratios, Furuya et al. (2016) proposed that water is primarily formed in the molecular cloud stage, before the dense core phase. This leads to a lower deuterium fractionation of water which initially freezes out onto the interstellar dust grains and constitutes the bulk of the water ice reservoir. Later, in the dense prestellar core phase, the deuterium fractionation is enhanced as the temperature drops and the visual extinction increases. In this phase highly deuterated ice is formed on top of the existing water ice on grain mantles. This scenario can explain the observed variation in the HDO/H2O ratio on different spatial scales. On larger scales, the observed gas-phase HDO/H2O ratio reflects the high deuterium fractionation of the outer layers of ice formed in the dense core phase as well as gas-phase synthesis of deuterated water through ion-neutral reactions in the cold outer regions of the envelope (Taquet et al. 2014; Furuya et al. 2016). In this region, the dust temperature is not sufficiently high to entirely sublimate the ice, and the inferred gas-phase abundance comes from a combination of the photodesorption of the outermost ice layers and gas-phase reactions with H2 D+ which can form HDO and D2O under these conditions. Meanwhile, on smaller scales the entire ice is sublimated off the grains and the ratio is lowered since the bulk of the water is formed with lower deuteration in the molecular cloud phase. Hence the ice abundances record the physical and chemical conditions of the different stages of star formation and the abundance of deuterated water could reflect the duration of the dense core phase. This in turn may depend on the local cloud conditions, since timescales of star formation is likely influenced heavily by external factors influencing the stability of the cloud, such as turbulence and radiation (e.g., Ward-Thompson et al. 2007).

Here we present observations of HDO and H218$_{2}^{18}$O toward three isolated Class 0 protostars, L483, BHR71–IRS1, and B335. We classify these sources as isolated as they are not associated with any known cloud complexes. This is in contrast to the previously targeted protostars which are identified as part of star-forming regions like NGC 1333 and Ophiucus and likely influenced by the dynamics within these stellar nurseries. The sources were selected to address if and how the deuterium fractionation of water varies as a function of the local cloud environment.

BHR71–IRS1 is a Bok globule located in the Southern Coalsack dark nebulae with a bolometric luminosity Lbol ≈ 15 L (Tobin et al. 2019). It is part of a wide binary system, with the companion located at 3200 au at a distance d ≈ 200 pc. Lynds 483, commonly referred to as L483, is an isolated dense core which harbors the infrared Class 0 source IRAS 18148–0440. Traditionally, L483 was associated with the Aquila Rift region at an inferred distance of 200 pc, however recent astrometry has revised the distance to the Aquila Rift up to d ≈ 436 ± 9 pc (Ortiz-León et al. 2018). Subsequent analysis based on stellar extinction and parallaxes from Gaia DR2 has shown that L483 is in fact located at a distance of 200–250 pc, i.e., not a part of the Aquila Rift complex (Jacobsen et al. 2019). We assume a distance of 200 pc in this paper. At this distance the estimated bolometric luminosity is 10–13 L (Tafalla et al. 2000; Shirley et al. 2000). B335 is a Bok globule located at d ≈ 100 pc (Olofsson & Olofsson 2009) and is the least luminous of our targets with a bolometric luminosity of Lbol ≈ 0.72 L (Evans et al. 2015). The central source is also known as IRAS 19347+0727 and identified as a low-mass Class 0 protostar. The central gas shows signs of infall and a rotational structure (e.g., Imai et al. 2019).

The paper is structured as follows. The observations and calibration procedures are presented in Sect. 2. Results are presented in Sect. 3. The HDO/H2O ratios are deduced and compared with previous measurements of water deuteration toward protostars in Sect. 4 along with a discussion of the implications. Finally, the results are summarized in Sect. 5 along with an outlook for the study of water deuteration and the impact of the local cloud environment.

Table 1

Observation log.

2 Observations

The low-mass embedded protostars L483, BHR71–IRS1, and B335 were observed with ALMA during Cycle 5 (PI: Jes K. Jørgensen, projectid: 2017.1.00693.S). For L483 the observations were centered on αJ2000 = 18:17:29.9, δJ2000 = –04:39:39.6, for BHR71–IRS1 αJ2000 = 12:01:36.5, δJ2000 = –65:08:49.3, and for B335 αJ2000 = 19:37:00.9, δJ2000 = 07:34:09.6. Information on the observation dates and calibration sources can be found in Table 1.

One spectral setup targeted the HDO 31,2 –22,1 and 21,1 –21,2 transitions at 225.8967 GHz (LSB) and 241.5616 GHz (USB) respectively in the ALMA band 6. Another spectral setup targeted the H218$_{2}^{18}$O 31,3 –22,0 transition at 203.4075 GHz in Band 5. Each spectral window contains 1920 channels with a width of 122 kHz (0.11 km s−1). The source velocities, estimated from the HDO 31,2 –22,1 transition are 4.5 km s−1, − 5.0 km s−1, 7.9 km s−1 for L483, BHR71–IRS1, and B335 respectively.

Each dataset was pipeline-calibrated using CASA 5.1 (McMullin et al. 2007). Phase self-calibration was performed using continuum channels for each dataset with CASA 5.4. For B335 a substantial improvement in the signal-to-noise ratio (S/N) was achieved through self-calibration, while only marginal gains were achieved for BHR71-IRS1 and no gains for L483. For the latter we opted to use the pipeline product, since the self-calibrated data offered no improvements. For the self-calibrated sources we performed continuum subtraction using CASA UVCONTSUB before inversion. The images were deconvolved using the TCLEAN algorithm with a robust parameter of − 0.5. For each source a continuum image was created at 202.7 GHz. The synthesized beam size range from 0.′′ 4 × 0.′′3 to 0.′′ 7 × 0.′′ 5 for the HDO spectral windows and 0.′′8 × 0.′′ 5 to 1.′′ 2 × 1.′′ 0 for the H218$_{2}^{18}$O spectral window (see Fig. 2).

3 Results

Figure 1 shows the continuum emission at 202.7 GHz above 3σrms toward the three sources. All sources are clearly detected and the continuum structure is resolved on linear scales of ~ 100 au. Toward L483 the continuum emission is extended along the north-west to south-east diagonal potentially tracing the cavity walls of an outflow directed perpendicular to this direction. Toward B335 the dust emission extends far out from the central source with a notable lack of emission in the east-west direction, likely driven by outflows along this direction. Since B335 is located closer than the other two sources more of the envelope emission is filtered out by the interferometer. Comparing the continuum toward B335 with Imai et al. (2016), who observed the source at similar angular resolution and continuum wavelength, we see good agreement: the 10σ contours (3 mJy beam−1) appear almost identical. Toward BHR71–IRS1 the continuum appears more circular with no evidence of the known outflows perturbing the dust in the plane of the sky.

The HDO emission lines are identified using data from the Jet Propulsion Laboratory (JPL, Pickett et al. 1998), referencing Messer et al. (1984), while the H218$_{2}^{18}$O transitional data originate from the Cologne Database for Molecular Spectroscopy (CDMS, Müller et al. 2001), with spectral details from de Lucia et al. (1972). All querying was done through the Splatalogue interface3.

The targeted water isotopologs are detected toward each of the sources: the HDO transitions are detected with high S/N > 10 while the H218$_{2}^{18}$O lines are slightly weaker detections (S/N ~ 5−10, Fig. 2). The targeted emission lines are presented in Fig. 3 along with fitted Gaussian profiles. The H218$_{2}^{18}$O data are rebinned by a factor of two for clarity, however this does not influence the abundances and HDO/H2O ratios derived in this paper.

Each spectrum is extracted toward the 202.7 GHz continuum peak of the source (see Fig. 1). For the spectral windows centered on the HDO transitions the image cube is convolved with a Gaussian kernel using the IMSMOOTH function in CASA. The size of the Gaussian kernel is determined such that the shape of the beam in the resulting image matches the larger beam of the H218$_{2}^{18}$O image cube at 203 GHz, allowing for direct comparison between the extracted spectra (see Fig. 2 for beam sizes before the convolution). Toward B335 and BHR71–IRS1 the emission lines are fitted using a single Gaussian profile. Toward L483 the emission lines show two components: a brighter redshifted component and weaker blueshifted component. In this case two Gaussian profiles were used to fit the lines. The double-peaked line profiles are consistent with previous high-resolution observations of hot corino species toward L483 presented in Oya et al. (2017) and Jacobsen et al. (2019) who independently found evidence of infalling rotating motion in the hot corino region toward the source.

Toward BHR71–IRS1 there is a slight hint of a broader outflow component in the wings of the emission lines, but this component is too weak to influence the fitted profile significantly. No extended outflow structure is seen when imaging the line-wings of any of the water transitions. The H218$_{2}^{18}$O line is partialy blended with dimethyl ether (CH3OCH3) toward BHR71-IRS1. In this case, six transitions of CH3OCH3 are fitted with Gaussian profiles using a fixed FWHM and system velocity and the fit is subsequently subtracted from the spectrum to remove the blending before the H218$_{2}^{18}$O line is fitted.

In Appendix C we present velocity field for the three transitions toward each source. The maps in Fig. C.1 confirm that no high-velocity outflow components contribute to the observed emission for the emission lines in question. Toward L483 the velocity field is consistent with the results of Jacobsen et al. (2019).

Figure 2 shows integrated emission maps for H218$_{2}^{18}$O and HDO toward the three sources. The detected emission originates from the central ~ 100 au for each source, consistent with the emission originating in the hot corino (T > 100 K). In both cases the central emission is unresolved and the slight variation in the extent of the emitting region can be explained by the differences in synthesised beam size. The peak positions of the continuum emission and water emission overlap down to scales of the synthesised beam size. Since we aim to determine the gas-phase abundances of HDO and H218$_2^{18}$O in the hot corino, it is essential that the observed emission lines originate from this region. We are confident this is the case for several reasons. The emission is compact, originating from scales of less than 100 au. This is confirmed by fitting two-dimensional Gaussian profiles to the emission maps. Here the FWHM extent on linear scales ranges from 50 au to 100 au. The fitted 2D profiles for the HDO 31,2 −22,1 transition toward the three sources are shown in the left column of Fig. 2. Additionally, the upper energy levels of the detected transitions lie around 100–200 K and the transitions are thus not easily excited in the cold outer envelope. Lastly, the line widths (FWHM ≲ 6 km s−1) of the emission lines are consistent with emission from the hot corino, with little or no evidence of outflow emission blending as described above and the emission is not extended along the known outflow directions, as in NGC 1333 IRAS 2A (Persson et al. 2012).

thumbnail Fig. 1

Continuum emission at 202.7 GHz toward the three sources. The map shows emission above 3σrms with white contours at 5σrms, 10σrms, and 30σrms. The maps are presented on same linear scale. The black cross marks the peak location and the grey arrows indicate the potential direction of outflows which may be perturbing the dust distribution. For BHR71–IRS1 the grey arrows indicates the outflow direction presented by Benedettini et al. (2017) and Yang et al. (2017).

thumbnail Fig. 2

Integrated emission for the targeted water transitions toward each of the sources. Left column: HDO 31,2 –22,1 transition at 225.9 GHz. The white shaded regions show the FWHM extent of a 2D Gaussian fitted to the data in the image plane and the white cross indicates the peak position of the fit. Middle column: HDO 21,1 –21,2 transition at 241.6 GHz. Right column: H218$_{2}^{18}$O 31,3 –22,0 transition at 203.4 GHz. Emission below 5σ is not included, where σ=σrms×Nchannels0.5×dv$\sigma = \sigma_{\textrm{rms}} \,{\times}\, N_{\textrm{channels}}^{0.5} \,{\times}\, dv$; dv is the channel width and N is the number of collapsed channels. The black cross marks the 202.7 GHz continuum peak position toward which the spectra are extracted.

thumbnail Fig. 3

Continuum-subtracted spectra for the three target transitions, extracted toward the continuum peak. Blue lines represent the Gaussian fits. For BHR71–IRS1, dimethyl ether lines (green) are fitted and subtracted from the spectrum before the H218$_{2}^{18}$O line (blue) is fitted. The H218$_{2}^{18}$O lines have been rebinned by a factor of two for clarity.

Table 2

Fitparameters for the targeted HDO and H218$_{2}^{18}$O transitions.

4 Analysis and discussion

4.1 Estimating the column densities of HDO and H218$_{{2}}^{{18}}$O and derivingthe HDO/H2O ratio

We estimated the column densities of HDO and H218$_{2}^{18}$O considering optically thin emission and local thermodynamic equilibrium (LTE) with an excitation temperature of 124 K. The choice of analysis was motivated by the aim to compare the present results with previous measurements presented in Persson et al. (2014) and Coutens et al. (2014) who determined the water deuteration in the hot corino toward a number of sources using the same method and excitation temperature. Choosing the same methodology, we can directly compare the water deuteration under similar assumptions. The results are summarized in Table 2. The estimated column densities range from (2− 7) × 1015 cm−2, comparable to the estimates of Persson et al. (2014) where the HDO column densities range from 1.5 × 1015 cm−2 for the faintest source, IRAS 4B, to 1.2 × 1016 cm−2 for IRAS 2A.

The estimated HDO/H2O ratios for L483, B335, and BHR71–IRS1 are (2.2 ± 0.4) × 10−3, (1.7 ± 0.3) × 10−3, and (1.8 ± 0.4) × 10−3, respectively. Uncertainties are derived from the statistical uncertainties of the fitted Gaussian profiles with an additional 10% uncertainty on the flux calibration. From the column density of H218$_{2}^{18}$O we infer the H2O water column densities by assuming the Galactic oxygen isotope ratio of 16O/18O = 560 (Wilson & Rood 1994). The HDO column density is determined as the weighted average of the two transitions. The HDO/H2O ratios toward the three sources show little scatter and are well within the uncertainties of one another which suggests similar chemical evolution of water in these systems.

The assumption of a fixed excitation temperature for the three sources only has a moderate effect on the estimated HDO/H2O ratios. Calculations presented by Persson et al. (2014) and Jørgensen & van Dishoeck (2010) have shown that varying the excitation temperature in the range 50–300 K has limited influence on the HDO column densities in the LTE approximation. We confirm this in Appendix B, where the HDO/H2O ratio is calculated for the measured line strengths and excitation temperatures in the range 30–300 K. Toward L483, an excitation temperature Tex ≤ 60 K is needed to bring the HDO/H2O ratio within the range of the four clustered sources reported in Persson et al. (2014) while B335 and BHR71–IRS1 require Tex ≤ 70 K. Such low excitation temperatures in the inner ~ 50–100 au are not expected. Since we observe two transitions of HDO for all three sources we can calculate the excitation temperature for each source under the assumption of LTE. The excitation temperatures are 127 ± 24 K, 174 ± 59 K, and 103 ± 21 K for L483, B335, and BHR71–IRS1; hence the choice of 124 K as the excitation temperature is consistent with the HDO line strengths toward the three sources. Meanwhile, the derived excitation temperature for HDO is inconsistent with the lower excitation temperatures needed to bring the HDO/H2O ratio down the values similar to the clustered sources, i.e., 60 K for L483 and 70 K for B335 and BHR71–IRS1. Computing the HDO/H2O ratio for the estimated excitation temperatures yields (2.2 ± 0.4) × 10−3 for L483, (2.1 ± 0.3) × 10−3 for B335, and (1.6 ± 0.3) × 10−3 for BHR71–IRS, all well within the error bars. Using the estimated excitation temperatures rather than the fixed value does not affect the conclusions presented here.

The difference between optically thin LTE calculations, similar to those presented here, and more advanced modelling for comparable observations of HDO and H218$_2^{18}$O have been studied in previous works to determine how well the former assumptions work in this regime. Persson et al. (2014) ran radiative transfer models for IRAS 16293–2422 and IRAS 2A and found that the water column densities were consistent with those derived from optically thin LTE calculations. Similarly, Coutens et al. (2014) ran non-LTE RADEX calculations and found column densities consistent with the LTE calculations. These results suggest that the water emission originating in the hot corino of these low-mass Class 0 protostars is well approximated by LTE calculations and that the emission is not optically thick. We note that if the densities in the inner envelopes of the isolated protostars are lower than the clustered counterparts then the emission from the water isotopologs could be sub-thermal and the LTE approximation invalid. We consider it unlikely to be the case; observational estimates of the overall envelope masses and bolometric luminosities for L483 and B335 are comparable to those of IRAS 2A, IRAS 4A, and IRAS 4B (see, e.g., Shirley et al. 2002; Jørgensen et al. 2007; Kristensen et al. 2012). With comparable envelope masses and luminosities it is unlikely that the inner density profile, i.e. on hot corino scales, would differ substantially between the isolated and clustered sources.

The advantage of more advanced models, as opposed to the LTE approach adopted here, is often diminished by the uncertainties of the physical source parameters such as the kinematics, the density profile, protostellar parameters etc. Attempting more advanced modelling may thus offer little improvement over the LTE approach and makes direct comparison between sources more challenging.

thumbnail Fig. 4

Comparison between the D/H ratio and HDO/H2O ratio for comets in the Solar System and hot corino observations toward Class 0 protostars. Values for IRAS 16293–2422, IRAS 2A, and IRAS 4B are from Persson et al. (2014), while IRAS 4A–NW has been adjusted from the value quoted in the paper as a mistake in the data analysis was discovered, which enhanced the HDO abundance by a factor of ~2. Errorbars show 1σ uncertainties. For the isolated sources the uncertainty is based on statistical errors from the fitted Gaussian profiles and a flux calibration uncertainty of 10%. On the right axis the corresponding D/H ratio is shown. The references for the Oort Cloud Comets (OCC) and Jupiter Family Comets (JFC) can be found in Appendix A. The colored regions show the standard deviation for each class of objects. Note that the HDO/H2O ratio for the protostars are derived from observations of HDO and H218$_{2}^{18}$O while some cometary values are derived from other proxies for the D/H ratio.

4.2 Water deuteration: comparison with existing observations

Figure 4 shows the calculated HDO/H2O ratios for the three isolated protostars along with existing data for a number of low-mass embedded protostars in clustered star-forming regions and cometary values from the Solar System. The protostellar values in Fig. 4 are all derived from interferometric observations with high spatial resolution to probe the hot corino emission where the ice is sublimated entirely off the dust grains. The three sources observed in this work have a similar degree of deuterium fractionation with HDO/H2O ratios in the range (1.7−2.2) × 10−3. Meanwhile, the four previously observed protostars lie in the range (5.5−9.2) × 10−4. This suggests that the sources presented here, which are all isolated protostars, have a distinct chemical history compared to the clustered protostars previously targeted. Comets for which the D/H ratio have been determined generally show lower values than the isolated protostars as seen in Fig. 4. Meanwhile, the clustered protostars show reasonable agreement with the Oort Cloud Comets, which have led to the suggestion that comets form from gas that is chemically similar to the gas observed in hot corinos, with little or no processing after this stage (Persson et al. 2014). Assuming this is the case, and that the D/H ratio of comets have not changed significantly after their formation, this would suggest that the Solar System was formed in a clustered region of star formation. Such a scenario is also supported by evidence such as the abundance of short-lived radionuclides, high-eccentricity orbits of small Solar System bodies, and low occurrence rate of isolated protostars (Adams 2010).

The apparent differentiation between clustered and isolated protostars can be understood in the framework for water formation and deuterium fractionation proposed by Furuya et al. (2016). They propose that water is primarily formed in the molecular cloud with limited deuteration. Later on, in the dense-core phase, deuteration is enhanced due to the low temperature and high shielding leading to the freeze-out of CO and a higher D/H ratio in the gas-phase. In this scenario, the duration of the dense-core phase determines the amount of deuteration, with a longer dense-core phase resulting in an enhanced deuterium enrichment in the ice. Conversely, a prolonged molecular cloud phase could decrease the deuterium fractionation and the HDO/H2O ratio is thusrelated to the ratio of the life-time of the two stages.

The star-formation process is recognized to be heterogenous in nature, with some stars born in dense clusters while others are formed in relatively isolated regions of molecular clouds (e.g., Ward-Thompson et al. 2007). Stars born in dense clusters are likely to collapse on shorter timescales; the mean free-fall time scales with the inverse square root of the mean density ρ0, tff=(3π/32Gρ0)1/2$t_{\mathrm{ff}} = (3\pi/32G\rho_{0})^{1/2}$ (e.g., Padoan et al. 2014; Krumholz 2014). Furthermore, external pressure from nearby massive stars in the local cloud region can trigger and potentially accelerate the collapse process in dense star forming regions. Accordingly, isolated protostars could experience alonger dense-core phase, enhancing the deuterium fractionation and potentially also enhancing the abundance of complex organic molecules formed from CO ice. The isolated protostars presented here all exhibit a higher deuterium fractionation than comparable counterparts in more dense, star-forming regions, in agreement with the theoretical expectation as outlined above. Assuming that the water formation processes are understood, this enhanced deuterium fractionation can imply two things: either the timescale of the dense-core phase is longer for the isolated cores or the molecular cloud phase, where most H2O is formed, isshorter for clustered star formation. The latter option does not appear likely since both theoretical and observational data imply that isolated cores do not collapse more easily than their clustered counterparts (Ward-Thompson et al. 2007). Given that the isolated protostars presented here exhibit higher deuterium fractionation compared to clustered protostars, this implies that the ratio between the duration of the dense-core phase and the molecular cloud phase is higher for these sources.

Another possible explanation for the differences between the HDO/H2O ratios toward clustered and isolated protostars could be variations in the temperatures between the different regions. For example, higher initial gas temperatures in the clustered regions, e.g., due to irradiation by neighboring young stars, would also reduce the efficiency of deuteration enrichment in these regions. Further, observational constraints and numerical studies are needed to address the relative importance of the collapse time-scales and gas temperatures.

The observations presented here support the hypothesis that the local cloud environment influences the early physical evolution and ultimately the chemistry of young stellar systems. A consequence of this is that the water deuteration can be an important proxy for both the chemical and the physical history of protostars.

Regarding the question of inheritance or local processing of water during star formation, the results presented here favor the inheritance scenario, at least at the earliest protostellar phase, since all three isolated protostars show similar HDO/H2O ratios. The lack of pronounced variation between the sources indicate a similar physical and chemical evolution for the systems, with little impact from local variations in, e.g., protostellar luminosity or accretion bursts. This conclusion is compatible with previous studies which modeled the water evolution from the collapse of an isolated pre-stellar cores to circumstellar disk (Visser et al. 2009; Cleeves et al. 2014; Drozdovskaya et al. 2016; Furuya et al. 2017).

5 Summary and outlook

In this paper we present the first ALMA Band 5 observations of the H218$_{2}^{18}$O 31,3 –22,0 transition toward three isolated low-mass Class 0 protostars. Combined with observations of the HDO 31,2 –22,1 and 21,1 –21,2 transitions we have determined the HDO/H2O ratio for the sources and compared with previous determinations of the water deuteration toward low-mass Class 0 protostars.

  • 1.

    The targeted water emission is detected in the hot corino toward the targeted sources with a high S/N on angular scales of 0.′′ 3–1.′′1 corresponding to linear scales of ~50–150 au. The column densities of HDO and H218$_2^{18}$O have been determined assuming optically thin emission under local thermodynamic equilibrium with an excitation temperature of 124 K.

  • 2.

    From the estimated column densities the derived HDO/H2O ratios for L483, B335, and BHR71–IRS1 are (2.2 ± 0.4) × 10−3, (1.7 ± 0.3) × 10−3, and (1.8 ± 0.4) × 10−3, respectively.

  • 3.

    The three isolated protostars show a factor of ~2 higher deuterium fractionation than the previously targeted sources, which are associated with larger cloud complexes. This observed differentiation in water deuteration is significant and can be explained by variations in the collapse timescales or initial gas temperatures depending on the cloud environment. In clustered regions external perturbations from nearby stars may either accelerate the collapse process through turbulence or heat the gas, leading to a lower deuterium fractionation. If this is the case then the degree of deuterium fractionation correlates with the local cloud environment, providing a new proxy for the early evolutionary history of young stars.

  • 4.

    The similarity of the HDO/H2O ratio toward the three isolated protostars could indicate that little processing of the water has occurred from cold cloud to hot core and suggests that the conditions in the dense core phase, before the onset of the collapse, determines the deuterium fractionation at later stages.

These observations present the first measurements of water deuteration targeting specific cloud environments and have nearly doubled the number of protostars for which the hot corino water deuteration has been measured. The observations indicate that isolated protostars have a distinct chemical history and further exploration of the relationship between cloud environment and deuterium fractionation could strengthen our understanding of the physical and chemical evolution during star formation. A natural progression is to determine the water deuteration toward more clustered or high-mass protostars, which should show lower fractionation levels according to the chemical evolution outlined above. Another option is to target doubly deuterated water, D2O. At this stage the D2O column density has only been determined toward the hot corino region of one source, NGC 1333 IRAS 2A (Coutens et al. 2014). Expanding the number of sources for which the D2O/HDO ratio is measured would test if the trend shown in Fig. 4 is real. Should the trend not be present for the doubly deuterated water isotopolog then we are missing important details in the current chemical models.

A third option to constrain the importance of the local cloud environment for the chemical evolution is to determine the deuteration of other molecules, such as methanol, in the hot corino. If the relationship between local cloud environment and molecular deuteration proposed here is correct, then the deuteration should be enhanced for all molecules.

In parallel with observational efforts, numerical modeling is needed to improve our understanding of the effects of the local cloud environment on the chemistry of young stellar systems. So far, little work has been done to model the chemical evolution in the context of a dynamic molecular cloud environment, i.e., collapse models which include the influence of the surrounding cloud environment with the local differences in temperature, density, UV-field, and turbulence. Such modeling would strengthen our understanding of the link between chemistry, particularly the deuterium fractionation, and the local cloud environment.

Acknowledgements

We thank the anonymous referee for valuable comments, which improved the manuscript. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2017.1.00693.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The group of J.K.J. acknowledges support from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 646908) through ERC Consolidator Grant “S4F”. Research at the Centre for Star and Planet Formation is funded by the Danish National Research Foundation. Astrochemistry in Leiden is supported by the Netherlands Research School for Astronomy (NOVA). A.C. postdoctoral grant is funded by the ERC Starting Grant 3DICE (grant agreement 336474).

Appendix A Table of D/H and HDO/H2O ratios

Values reported in Fig. 4 are presented in Table A.1 along with references and information of the tracers used to determine the water deuterium fractionation.

Table A.1

Measured HDO/H2O and D/H ratios for comets and protostars.

Appendix B Computing HDO/H2O as a function of excitation temperature in the LTE approximation

In Fig. B.1 we show the HDO/H2O ratio as a function of excitation temperature in the range 30–300 K under the assumption of optically thin LTE emission. All values are calculated using the line strengths from the fitted line profiles presented in Fig. 3. Evidently the exact choice of excitation temperature in the regime 100–200 K has limited effect on the corresponding HDO/H2O ratio. The green shaded region indicates the upper limit for the clustered protostars in Persson et al. (2014), HDO/H2O = 1.18 × 10−3 for IRAS 16293–2422. For L483 an excitation temperature of ~60 K is needed to lower the HDO/H2O to 1.18 × 10−3 while the threshold for B335 and BHR71–IRS1 lies around 70 K.

thumbnail Fig. B.1

HDO/H2O ratio for the three sources assuming excitation temperatures in the range 30–300 K and optically thin LTE emission. The green shaded area show the upper limit for the HDO/H2O ratio for the clustered protostars presented in Persson et al. (2014). The blue circle marks the value at the assumed excitation temperature of 124 K value.

Appendix C Velocity fields for HDO and H218$_{2}^{18}$O

Moment 1 velocity fields computed in CASA to verify that no high-velocity components are present for any of the targeted water transitions.

thumbnail Fig. C.1

Moment 1 velocity fields toward each of the sources for the three water transitions. The maps show no signs of outflow emission contributing to the emission for the targeted transitions. Emission below 5σrms is excluded in the integration.

References

  1. Adams, F. C. 2010, ARA&A, 48, 47 [NASA ADS] [CrossRef] [Google Scholar]
  2. Altwegg, K., Balsiger, H., Bar-Nun, A., et al. 2015, Science, 347, 1261952 [Google Scholar]
  3. Benedettini, M., Gusdorf, A., Nisini, B., et al. 2017, A&A, 598, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Bergin, E. A.,& Snell, R. L. 2002, ApJ, 581, L105 [NASA ADS] [CrossRef] [Google Scholar]
  5. Biver, N., Bockelée-Morvan, D., Crovisier, J., et al. 2006, A&A, 449, 1255 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Biver, N., Moreno, R., Bockelée-Morvan, D., et al. 2016, A&A, 589, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Bockelée-Morvan, D., Gautier, D., Lis, D. C., et al. 1998, Icarus, 133, 147 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bockelée-Morvan, D., Biver, N., Swinyard, B., et al. 2012, A&A, 544, L15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Brown, R. H., Lauretta, D. S., Schmidt, B., & Moores, J. 2012, Planet. Space Sci., 60, 166 [NASA ADS] [CrossRef] [Google Scholar]
  10. Caselli, P., Keto, E., Bergin, E. A., et al. 2012, ApJ, 759, L37 [Google Scholar]
  11. Cazaux, S., Caselli, P., & Spaans, M. 2011, ApJ, 741, L34 [NASA ADS] [CrossRef] [Google Scholar]
  12. Ceccarelli, C., Caselli, P., Bockelée-Morvan, D., et al. 2014, in Protostars and Planets VI, eds. H. Beuther, R. S. Klessen, C. P. Dullemond, & T. Henning (Tucson: University of Arizona Press), 859 [Google Scholar]
  13. Cleeves, L. I., Bergin, E. A., Alexander, C. M. O. D., et al. 2014, Science, 345, 1590 [NASA ADS] [CrossRef] [Google Scholar]
  14. Coutens, A., Vastel, C., Caux, E., et al. 2012, A&A, 539, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Coutens, A., Vastel, C., Cabrit, S., et al. 2013, A&A, 560, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Coutens, A., Jørgensen, J. K., Persson, M. V., et al. 2014, ApJ, 792, L5 [NASA ADS] [CrossRef] [Google Scholar]
  17. de Laeter, J. R., Böhlke, J. K., Bièvre, P. D., et al. 2003, Pure Appl. Chem., 75, 683 [CrossRef] [Google Scholar]
  18. de Lucia, F. C., Helminger, P., Cook, R. L., & Gordy, W. 1972, Phys. Rev. A, 6, 1324 [NASA ADS] [CrossRef] [Google Scholar]
  19. Drozdovskaya,M. N., Walsh, C., van Dishoeck, E. F., et al. 2016, MNRAS, 462, 977 [NASA ADS] [CrossRef] [Google Scholar]
  20. Evans, II, N. J., Di Francesco, J., Lee, J.-E., et al. 2015, ApJ, 814, 22 [CrossRef] [Google Scholar]
  21. Furuya, K., van Dishoeck, E. F., & Aikawa, Y. 2016, A&A, 586, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Furuya, K., Drozdovskaya, M. N., Visser, R., et al. 2017, A&A, 599, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Gibb, E. L., Bonev, B. P., Villanueva, G., et al. 2012, ApJ, 750, 102 [NASA ADS] [CrossRef] [Google Scholar]
  24. Hartogh, P., Lis, D. C., Bockelée-Morvan, D., et al. 2011, Nature, 478, 218 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  25. Hutsemékers, D., Manfroid, J., Jehin, E., Zucconi, J. M., & Arpigny, C. 2008, A&A, 490, L31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Imai, M., Sakai, N., Oya, Y., et al. 2016, ApJ, 830, L37 [NASA ADS] [CrossRef] [Google Scholar]
  27. Imai, M., Oya, Y., Sakai, N., et al. 2019, ApJ, 873, L21 [NASA ADS] [CrossRef] [Google Scholar]
  28. Jacobsen, S. K., Jørgensen, J. K., Di Francesco, J., et al. 2019, A&A, 629, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Jørgensen, J. K., & van Dishoeck, E. F. 2010, ApJ, 725, L172 [NASA ADS] [CrossRef] [Google Scholar]
  30. Jørgensen, J. K., Bourke, T. L., Myers, P. C., et al. 2007, ApJ, 659, 479 [NASA ADS] [CrossRef] [Google Scholar]
  31. Kristensen, L. E., van Dishoeck, E. F., Bergin, E. A., et al. 2012, A&A, 542, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Krumholz, M. R. 2014, Phys. Rep., 539, 49 [NASA ADS] [CrossRef] [Google Scholar]
  33. Lis, D. C., Biver, N., Bockelée-Morvan, D., et al. 2013, ApJ, 774, L3 [NASA ADS] [CrossRef] [Google Scholar]
  34. Lis, D. C., Bockelée-Morvan, D., Güsten, R., et al. 2019, A&A, 625, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, in Astronomical Data Analysis Software and Systems XVI, eds. R. A. Shaw, F. Hill, & D. J. Bell, ASP Conf. Ser., 376, 127 [NASA ADS] [Google Scholar]
  36. Meier, R., Owen, T. C., Matthews, H. E., et al. 1998, Science, 279, 842 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  37. Messer, J., Lucia, F. C. D., & Helminger, P. 1984, J. Mol. Spectr., 105, 139 [NASA ADS] [CrossRef] [Google Scholar]
  38. Müller, H. S. P., Thorwirth, S., Roth, D. A., & Winnewisser, G. 2001, A&A, 370, L49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Olofsson, S., & Olofsson, G. 2009, A&A, 498, 455 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Ortiz-León, G. N., Loinard, L., Dzib, S. A., et al. 2018, ApJ, 869, L33 [NASA ADS] [CrossRef] [Google Scholar]
  41. Oya, Y., Sakai, N., Watanabe, Y., et al. 2017, ApJ, 837, 174 [NASA ADS] [CrossRef] [Google Scholar]
  42. Padoan, P., Federrath, C., Chabrier, G., et al. 2014, in Protostars and Planets VI, eds. H. Beuther, R. S. Klessen, C. P. Dullemond, & T. Henning (Tucson: University of Arizona Press), 77 [Google Scholar]
  43. Parise, B., Caux, E., Castets, A., et al. 2005, A&A, 431, 547 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Persson, M. V., Jørgensen, J. K., & van Dishoeck, E. F. 2012, A&A, 541, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Persson, M. V., Jørgensen, J. K., & van Dishoeck, E. F., & Harsono, D. 2014, A&A, 563, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Pickett, H. M., Poynter, R. L., Cohen, E. A., et al. 1998, J. Quant. Spectr. Rad. Transf., 60, 883 [NASA ADS] [CrossRef] [Google Scholar]
  47. Pontoppidan, K. M., van Dishoeck, E. F., & Dartois, E. 2004, A&A, 426, 925 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Prodanović, T., Steigman, G., & Fields, B. D. 2010, MNRAS, 406, 1108 [NASA ADS] [Google Scholar]
  49. Shirley, Y. L., Evans, N. J., II, Rawlings, J. M. C., & Gregersen, E. M. 2000, ApJS, 131, 249 [NASA ADS] [CrossRef] [Google Scholar]
  50. Shirley, Y. L., Evans, N. J., II, & Rawlings, J. M. C. 2002, ApJ, 575, 337 [NASA ADS] [CrossRef] [Google Scholar]
  51. Stark, R., Sandell, G., Beck, S. C., et al. 2004, ApJ, 608, 341 [NASA ADS] [CrossRef] [Google Scholar]
  52. Tafalla, M., Myers, P. C., Mardones, D., & Bachiller, R. 2000, A&A, 359, 967 [NASA ADS] [Google Scholar]
  53. Taquet, V., López-Sepulcre, A., Ceccarelli, C., et al. 2013a, ApJ, 768, L29 [NASA ADS] [CrossRef] [Google Scholar]
  54. Taquet, V., Peters, P. S., Kahane, C., et al. 2013b, A&A, 550, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Taquet, V., Charnley, S. B., & Sipilä, O. 2014, ApJ, 791, 1 [NASA ADS] [CrossRef] [Google Scholar]
  56. Tobin, J. J., Bourke, T. L., Mader, S., et al. 2019, ApJ, 870, 81 [NASA ADS] [CrossRef] [Google Scholar]
  57. van Dishoeck, E. F., Kristensen, L. E., Benz, A. O., et al. 2011, PASP, 123, 138 [NASA ADS] [CrossRef] [Google Scholar]
  58. van Dishoeck, E. F., Bergin, E. A., Lis, D. C., & Lunine, J. I. 2014, in Protostars and Planets VI, eds. H. Beuther, R. S. Klessen, C. P. Dullemond, & T. Henning (Tucson: University of Arizona Press), 835 [Google Scholar]
  59. Villanueva, G. L., Mumma, M. J., Bonev, B. P., et al. 2009, ApJ, 690, L5 [NASA ADS] [CrossRef] [Google Scholar]
  60. Visser, R., van Dishoeck, E. F., Doty, S. D., & Dullemond, C. P. 2009, A&A, 495, 881 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Visser, R., Jørgensen, J. K., Kristensen, L. E., van Dishoeck, E. F., & Bergin, E. A. 2013, ApJ, 769, 19 [NASA ADS] [CrossRef] [Google Scholar]
  62. Ward-Thompson, D., André, P., Crutcher, R., et al. 2007, in Protostars and Planets V, eds. B. Reipurth, D. Jewitt, & K. Keil (Tucson: University of Arizona Press), 33 [Google Scholar]
  63. Wilson, T. L.,& Rood, R. 1994, ARA&A, 32, 191 [NASA ADS] [CrossRef] [Google Scholar]
  64. Yang, Y.-L., Evans, II, N. J., Green, J. D., Dunham, M. M., & Jørgensen, J. K. 2017, ApJ, 835, 259 [NASA ADS] [CrossRef] [Google Scholar]

1

The exact value of ΔE depends on the spin state of the involved reactants.

2

[HDO/H2O] = 2 × [D/H] and [D/H]ISM = 2 × 10−5 (Prodanović et al. 2010).

All Tables

Table 1

Observation log.

Table 2

Fitparameters for the targeted HDO and H218$_{2}^{18}$O transitions.

Table A.1

Measured HDO/H2O and D/H ratios for comets and protostars.

All Figures

thumbnail Fig. 1

Continuum emission at 202.7 GHz toward the three sources. The map shows emission above 3σrms with white contours at 5σrms, 10σrms, and 30σrms. The maps are presented on same linear scale. The black cross marks the peak location and the grey arrows indicate the potential direction of outflows which may be perturbing the dust distribution. For BHR71–IRS1 the grey arrows indicates the outflow direction presented by Benedettini et al. (2017) and Yang et al. (2017).

In the text
thumbnail Fig. 2

Integrated emission for the targeted water transitions toward each of the sources. Left column: HDO 31,2 –22,1 transition at 225.9 GHz. The white shaded regions show the FWHM extent of a 2D Gaussian fitted to the data in the image plane and the white cross indicates the peak position of the fit. Middle column: HDO 21,1 –21,2 transition at 241.6 GHz. Right column: H218$_{2}^{18}$O 31,3 –22,0 transition at 203.4 GHz. Emission below 5σ is not included, where σ=σrms×Nchannels0.5×dv$\sigma = \sigma_{\textrm{rms}} \,{\times}\, N_{\textrm{channels}}^{0.5} \,{\times}\, dv$; dv is the channel width and N is the number of collapsed channels. The black cross marks the 202.7 GHz continuum peak position toward which the spectra are extracted.

In the text
thumbnail Fig. 3

Continuum-subtracted spectra for the three target transitions, extracted toward the continuum peak. Blue lines represent the Gaussian fits. For BHR71–IRS1, dimethyl ether lines (green) are fitted and subtracted from the spectrum before the H218$_{2}^{18}$O line (blue) is fitted. The H218$_{2}^{18}$O lines have been rebinned by a factor of two for clarity.

In the text
thumbnail Fig. 4

Comparison between the D/H ratio and HDO/H2O ratio for comets in the Solar System and hot corino observations toward Class 0 protostars. Values for IRAS 16293–2422, IRAS 2A, and IRAS 4B are from Persson et al. (2014), while IRAS 4A–NW has been adjusted from the value quoted in the paper as a mistake in the data analysis was discovered, which enhanced the HDO abundance by a factor of ~2. Errorbars show 1σ uncertainties. For the isolated sources the uncertainty is based on statistical errors from the fitted Gaussian profiles and a flux calibration uncertainty of 10%. On the right axis the corresponding D/H ratio is shown. The references for the Oort Cloud Comets (OCC) and Jupiter Family Comets (JFC) can be found in Appendix A. The colored regions show the standard deviation for each class of objects. Note that the HDO/H2O ratio for the protostars are derived from observations of HDO and H218$_{2}^{18}$O while some cometary values are derived from other proxies for the D/H ratio.

In the text
thumbnail Fig. B.1

HDO/H2O ratio for the three sources assuming excitation temperatures in the range 30–300 K and optically thin LTE emission. The green shaded area show the upper limit for the HDO/H2O ratio for the clustered protostars presented in Persson et al. (2014). The blue circle marks the value at the assumed excitation temperature of 124 K value.

In the text
thumbnail Fig. C.1

Moment 1 velocity fields toward each of the sources for the three water transitions. The maps show no signs of outflow emission contributing to the emission for the targeted transitions. Emission below 5σrms is excluded in the integration.

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.