Large-scale Mapping Observations of DCN and DCO+ toward Orion KL

We present emission maps ( 1.′5×1.′5 scale, corresponding to 0.18 pc) of the DCN (J = 2 − 1) and DCO+ (J = 2 − 1) lines in the 2 mm band toward the Orion KL region obtained with the 2 mm receiver system named B4R installed on the Large Millimeter Telescope. The DCN emission shows a peak at the Orion KL hot core position, whereas no DCO+ emission has been detected there. The DCO+ emission shows enhancement at the west side of the hot core, which is well shielded from the UV radiation from OB massive stars in the Trapezium cluster. We have derived the abundance ratio of DCN/DCO+ at three representative positions where both species have been detected. The gas components with V LSR ≈ 7.5–8.7 km s−1 are associated with low abundance ratios of ∼4–6, whereas much higher abundance ratios (∼22–30) are derived for the gas components with V LSR ≈ 9.2–11.6 km s−1. We have compared the observed abundance ratio to our chemical models and found that the observed differences in the DCN/DCO+ abundance ratios are explained by different densities.

Deuterium fractionation is well known to start with an enhancement of the D/H ratios in several ions such as H + 3 , CH + 3 , and C 2 H + 2 by the following exothermic isotope-exchange reactions with HD; CH + 3 + HD ⇌ CH 2 D + + H 2 + ∆ E (390 K), (2) and C 2 H + 2 + HD ⇌ C 2 HD + + H 2 + ∆ E (550 K). (3) In cold conditions (T ≤ 10 K), the backward reactions of reactions (1) -( 3) are suppressed, and the deuterium is exceeded in these ions.These ions are precursors to other deuterium species, both molecular and ionic forms.For example, DCO + is considered to form from H 2 D + in low-temperature regions (T < 30 K; Ceccarelli et al. 2014) as follows: One main formation route of DCN is considered as the following reaction scheme (Millar et al. 1989;Turner 2001): followed by CHD + N → DCN + H.
Another pathway starts with the following reaction: followed by reaction (7).Thus, in cold conditions, DCO + and DCN are formed via H 2 D + and CH 2 D + , which are produced by reactions (1) and (2), respectively.
Several reactions can produce deuterium species in warm regions.For instance, the following reaction is responsible for the formation of DCO + in warm regions (> 30 K; Adams & Smith 1985): The deuterium fractionation mechanisms that proceed in cold and warm regions have been found to work in active star-forming regions.Parise et al. (2009) revealed the deuterium fractionation starting from CH 2 D + in warm regions by observations using the APEX and IRAM 30m telescopes toward two clumps in Orion Bar, which is one of the most famous photon-dominated regions (PDRs).They detected DCN, DCO + , and HDCO, and found significant deuterium fractionation for HCN and H 2 CO but a low fractionation in HCO + .They compared the observational results with their chemical simulations and concluded that the deuterium fractionation starting from CH 2 D + is important in this region.Gerner et al. (2015) investigated chemical evolution using several deuterium species toward 59 sources in high-mass star-forming regions using the Arizona Radio Observatory Submillimeter Telescope (SMT).They found that the DCN/HCN abundance ratio shows peaks at the hot molecular core stage, whereas the D/H ratios in HNC, HCO + , and N 2 H + decrease with time evolution.These results suggest that the deuterium fractionation of DCN can proceed efficiently even in warm or hot regions.Sakai et al. (2022) observed several deuterium species using the Atacama Large Millimeter/submillimeter Array (ALMA) toward the infrared dark cloud (IRDC) G 14.492-0.139and found that the emission of DCO + and DCN is signs of star formation activity in the IRDC.They also found that the observed DCO + /N 2 D + abundance ratios in this IRDC clump are lower than those in starless cores in low-mass star-forming regions, implying that the dense cores in this IRDC are warmer and denser than those in lowmass star-forming regions.Taniguchi et al. (2021) presented observational results of DCO + and DCN toward two low-mass protostars in the Ophiuchus star-forming region obtained with the ALMA.They have revealed that different deuterium fractionation mechanisms work around each protostar; deuterium fractionation proceeds in the warm regions (> 30 K) around Oph-emb9, whereas pathways that can proceed in cold regions (< 30 K) are efficient around Oph-emb5.Taniguchi et al. (2021) found that the UV radiation field may be responsible for the enhancement of these species around Oph-emb9.
The previous studies show that DCO + and DCN are likely key species to investigate present and past physical conditions in active star-forming regions because they are detected ubiquitously and behave differently in star-forming regions.Studies to investigate deuterium fractionation have mainly focused on small scales, i.e., cores and protoplanetary disks (e.g., Salinas et al. 2017;Öberg et al. 2021;Cataldi et al. 2021).On the other hand, studies of large-scale (i.e., cloud-scale) deuterium fractionation are scarce so far (e.g., N 2 H + /N 2 D + by Cosentino et al. 2023).
In this paper, we present the DCN and DCO + maps toward the Orion KL region (d = 418 ± 6 pc; Kim et al. 2008).Orion KL is located in the Orion Molecular Cloud-1 (OMC-1).Large-scale mapping observations toward the Orion region including OMC-1 were conducted with the Nobeyama 45m radio telescope as a part of a legacy project (Nakamura et al. 2019).OMC-1 is irradiated by the very strong far-ultraviolet (FUV: 6 eV < hν < 13.6 eV) radiation field but has the highest column densities of 13 CO and C 18 O among the Orion A region (Ishii et al. 2019).Recent studies with the James Webb Space Telescope (JWST) have revealed the 3D structure around this region.The Trapezium cluster irradiates the background OMC-1 creating the PDR Orion Bar (Peeters et al. 2023).Although the Trapezium cluster includes some massive OB stars, the effects from the most massive one, the θ1 Orionis C (θ 1 Ori C), seem to be dominant (Peeters et al. 2023).Orion KL is one of the best-studied hot cores rich with molecular lines from complex organic molecules and has been focused by many astrochemical studies (e.g., Schilke et al. 1997;Feng et al. 2015;Gong et al. 2015;Rizzo et al. 2017).
The data presented in this paper were obtained as a part of commissioning observations of the 2 mm band receiver system, which was named B4R 1 , installed on the Large Millimeter Telescope (LMT) or Gran Telescopio Milimétrico Alfonso Serrano (Hughes et al. 2020).This paper is organized as follows.In Section 2, we explain observational setups and data reduction procedures.Obtained maps and spectra are presented in Section 3. We investigate relationships between the observed DCN/DCO + abundance ratios and the velocity components in Section 4.1, and compare the observational results with the chemical simulations in Section 4.2.The main conclusions of this paper are summarized in Section 5.

OBSERVATIONS AND DATA REDUCTION
Observations presented in this paper were obtained during the commissioning of the B4R system from 20 to 29 November 2019.The data were obtained to make a beam map using the SiO maser line (J = 3 − 2, v = 1).The coordinate of the map center is (α J2000 , δ J2000 ) = (5 h 35 m 14. s 16, -5 • 22 ′ 21. ′′ 50).The details of the B4R system and results of its commissioning observations are presented in R. Kawabe et al. (2023, in prep.).The main beam efficiency and beam size at 140 GHz were 59 % and 10.6 ′′ , respectively (R. Kawabe et al., 2023).The frequency setup covers 128.9-131.4GHz and 142.6-145.1 GHz in LSB and USB, respectively.The DCO + (J = 2− 1; 144.077289GHz, E up = 10.4K) and DCN (J = 2 − 1; 144.8280015GHz, E up = 10.4K) lines were covered in the USB.The system noise temperatures (T sys ) were between 100 K and 150 K during the commissioning.
We conducted data reduction with the pipeline (b4rpipe)2 which was developed for the B4R observations.We used the data in the condition that the pointing accuracy was not degraded by wind during the observations to make final maps.The final maps were made in the fits format through the pipeline.The map size is 1.5 ′ × 1.5 ′ , corresponding to a 0.18-pc scale.A frequency resolution is 76.29 kHz, corresponding to ∼ 0.16 km s −1 in the velocity resolution at the observed frequency band.The beam size of the final maps is 11 ′′ , corresponding to 0.02 pc at the source distance (418 pc).
After making the fits files with the pipeline, we treated data and conducted further data reduction and analyses using the Common Astronomy Software Applications (CASA) package (CASA Team et al. 2022).We corrected the position offset.The 5th-order polynomial fitting calculated from line-free channels was applied for the baseline fitting.The detailed methods of the procedure to make final maps will be explained in T. Yonetsu et al. (in prep.).
When constructing moment maps in CASA, we applied the sigma-clipping method, because several velocity components are mixed in this region.We made masks for emissions above 3σ for the 3D cube fits using the immath task, and smoothed the masks with a single-beam scale (11 ′′ ) using the imsmooth task.The 3σ noise thresholds are 0.27 K and 0.25 K for the DCN and DCO + maps, respectively.We applied the generated masks for the original 3D cube data using the immath and obtained the 3D cube data.Using the new 3D cube data applied to the sigma-clipping method, we made moment maps presented in Section 3.1.

Maps of DCN and DCO +
Figure 1 shows integrated-intensity maps (moment 0 maps; top panels) and peak intensity maps (moment 8 maps; middle panels) of the DCN and DCO + lines, and a velocity map of DCN (moment 1 map; bottom panel).Based on the moment 0 maps, we identified five representative positions; DCN Peak, DCO + Peak1, DCO + Peak2, Tail, and Outer.The coordinate, H 2 column density, indicated as N (H 2 ), dust temperature, and visual extinction (A V ) at each position are summarized in Table 1.The values of N (H 2 ) and dust temperature are obtained from maps of Schuller et al. (2021).The H 2 column density map has an angular resolution of 8 ′′ .This map was made using SPIRE 160, 250, 350, and 500 µm data from Herschel Gould Belt Survey (HGBS; André et al. 2010), and ArTéMiS 350 and 450 µm data.The fits file of the dust temperature has an angular resolution of SPIRE 250 µm (18.2 ′′ ).This map was made by combining data of ArTéMiS 350 and 450 µm, SPIRE 350 and 500 µm, and SPIRE 160 and 250 µm data from HGBS.
We calculated the visual extinction (A v ) from the H 2 column density using the following formulae (Gratier et al. 2021): (10) Here, we assumed that all of the hydrogen nuclei are fixed in the molecular form (H 2 ).These five positions are indicated in the maps as black crosses.Figure 2 shows maps of N (H 2 ) and dust temperature obtained from Schuller et al. (2021), overlaid by moment 0 maps of DCN (black contours) and DCO + (gray contours).Moment 0 maps of the two deuterium species clearly show different spatial distributions.The DCN emission shows a peak associated with the hot core position and extends from northeast to southwest, similar to the N (H 2 ) map (left panel of Figure 2).On the other hand, the DCO + emission is mainly enhanced at the western parts from the hot core position.In addition, the DCO + emission is associated with Tail as well as DCN.Several recombination lines from some atoms (e.g., hydrogen and carbon), which trace H II regions, have been detected at the east side from the Orion KL hot core position (e.g., Goicoechea & Cuadrado 2021).Thus, the DCO + emission is enhanced in the shield regions against the UV radiation from the H II region.Enhancement in the temperature at the west side near DCO + Peak1 is likely heated by the Orion KL hot core source, rather than the UV radiation from the H II region.

Spectral Analyses
Black curves in Figure 3 show the observed spectra of the DCN (left) and DCO + (right) lines at the five representative positions.These spectra are extracted over the beam size (11 ′′ ), and the center coordinates are set at the values summarized in Table 1.The DCN (J = 2−1) line has been detected from all of the positions, whereas the DCO + (J = 2 − 1) line has been detected from the three positions except for DCN Peak and Outer, both of which are illuminated by the UV radiation from the massive star(s).
We analyzed spectra with the Markov Chain Monte Carlo (MCMC) method in the CASSIS software (Vastel et al. 2015).In this method, column density (N ), excitation temperature (T ex ), velocity component (V LSR ), line width (FWHM) are treated as free parameters.We assumed the thermodynamic equilibrium (LTE) condition.This assumption is reasonable here because lines of deuterium species are usually optically thin.We have only one transition line for each species, and we could not determine their column density and excitation temperature simultaneously.We then analyzed spectra changing excitation temperatures from 10 K to 50 K with the 5-K step.The lower limit of the temperature was determined because the deuterium fractionation prefers cold regions (Section 1), and its upper limit was determined from the dust temperature map (Figure 2 and Table 1).Regarding DCN Peak, we changed the excitation temperature from 40 K to 50 K, judging from the dust temperature map (right panel of Figure 2) and broad line widths.The red curves in Figure 3 show the best-fitting model with an assumed excitation temperature of 40 K.Even if we changed the assumed excitation temperatures from 10 K to 50 K, the best-fitting models show similar spectral features and there is no significant difference.We applied two-velocity-component fitting for the spectra at DCN Peak and DCO + Peak2.At DCO + Peak1, DCO + Peak2, and Tail, we analyzed the DCO + spectra at first and utilized the obtained velocity components to constrain the velocity ranges of the DCN line at each position for secure line identification.The hyperfine splitting due to nitrogen nuclei and deuterium is taken into account.
The results obtained by the MCMC method are summarized in Table 2.The velocity components are consistent between the DCN and DCO + lines at each position, suggesting that they trace similar regions.The DCN/DCO + abundance ratios at each position are summarized in the last column of Table 2. Table 3 summarizes fractional abundances denoted as X(molecules), defined by N (molecules)/N (H 2 ), at each position.The optical depths of the DCN line are less than 0.73, and those of the DCO + line are less than 0.29, in the assumed excitation temperature range and at all of the positions.Thus, these lines are optically thin in all of the cases.
The wing emission in the DCN spectra is likely caused by turbulence or outflows.We did not include these effects in the spectral fitting procedure.Another velocity component around V LSR ≈ 1 − 4 km s −1 is seen in the DCN spectra at DCO + Peak1 and Tail.We tried fitting these lines as another velocity component, but the lines could not be fitted simultaneously.Thus, these velocity components are contamination of another molecular line or emission from completely different regions.The following discussions and our conclusions are not affected by these components, because we will discuss mainly using the DCN/DCO + abundance ratios of the velocity components which are common in both species.Note-The H2 column density and temperature are obtained from fits files obtained from Schuller et al. (2021).The figures in parentheses refer to the standard deviation obtained from the fits files.

Note-
The figures in parentheses refer to the standard deviation.

Characteristics of the observed DCN/DCO + abundance ratios
We found that the DCN/DCO + abundance ratios are divided into two groups (Table 2); the ratio of ∼ 4 − 7 (low) and ∼ 22 − 30 (high).At DCO + Peak2, we found both components in the different velocity components, which means that we need to investigate the DCN/DCO + abundance ratios in the three dimensions.
Figure 4 shows a dependence of the DCN/DCO + abundance ratio on the velocity component.In Figure 4, we utilized the velocity components derived by fitting the DCO + line (Table 2).The lower abundance ratios are associated with the gas of V LSR ≈ 7.5 − 8.7 km s −1 , whereas the higher abundance ratios are accompanied with the velocity components of ≈ 9.2 − 11.6 km s −1 .These results imply that gases with different physical conditions are mixed.Pabst et al. (2019) found that adiabatic expansion of the hot gas has swept up the surrounding gas, and the velocity of the expanding shell has a velocity of 13 km s −1 .This expanding shell is formed by the most massive star in the Trapezium cluster (θ 1 Ori C).The massive star is located at (α J2000 , δ J2000 ) = (5 h 35 m 16. s 46, -5 • 23 ′ 22. ′′ 8).The velocity of the expanding shell is close to the gas components with the higher DCN/DCO + abundance ratios (V LSR ≈ 9.2 − 11.6 km s −1 ).On the other hand, the systemic velocity of the Orion Molecular Cloud (OMC) A is ∼ 8 km s −1 (Pabst et al. 2019), which is similar to the gases with the lower DCN/DCO + abundance ratios (V LSR ≈ 7.5 − 8.7 km −1 ).
In these different gas components, different deuterium fractionation pathways most likely proceed.Based on comparisons of the velocity components, the gas components with the higher velocity component (V LSR ≈ 9.2 − 11.6 km s −1 ) are affected by the Trapezium cluster more significantly, and fluxes of high energy particles (i.e., UV photons and cosmic rays) should be higher compared to the behind dense OMC-1 with the lower velocity component (V LSR ≈ 7.5 − 8.7 km s −1 ).The higher DCN/DCO + abundance ratios come from the gas affected by the expanding shell, while the lower abundance ratios are associated with the dense molecular clouds.We will discuss detailed deuterium fractionation of these two species using our latest chemical simulation in Section 4.2.2017).The DNAUTILUS.2.0 is capable of studying deuterium fractionation in both two-and three-phase modes, with or without ortho and para spin-states of major hydrogenated species, such as H 2 , H + 2 , and H + 3 , and their isotopologues.In this work, we employed the three-phase version of DNAUTILUS 2.0 without spin-state chemistry for computational efficiency and simplicity when running the parameter space summarized in Table 4.This version of DNAUTILUS 2.0 encompasses all the physicochemical processes and their corresponding equations derived from the state-of-the-art publicly available astrochemical code NAUTILUS, as detailed in Ruaud et al. (2016).
DNAUTILUS 2.0 includes 1606 gas species, 1472 grainsurface species, and 1472 grain-mantle species, connected by 83,715 gas-phase reactions, 10,967 reactions on grain surfaces, and 9,431 reactions in the grain mantles.This version was generated by applying our welltested deuteration routine to the gas, grain, and mantle reaction networks of NAUTILUS, where deuterons substitute for protons in the reactions, and branching ratios are calculated assuming complete scrambling.This methodology is similar to the one described in Albertsson et al. (2013).
The binding energies for neutral deuterated species are assumed to be similar to their non-deuterated counterparts, following previous studies (Taquet et al. 2013;Albertsson et al. 2013;Majumdar et al. 2017).This is unavoidable due to the limited data available on the binding energies of deuterated molecules on grain surfaces.The binding energy of DCN is set at 2050 K, similar to the standard value of HCN as considered in NAUTILUS (Ruaud et al. 2016).For initial abundances, the species are assumed to be initially in atomic form, similar to diffuse clouds, except for hydrogen and deuterium, which are initially in H 2 and HD forms, respectively.The initial elemental abundances are summarized in Table 4.

Physical parameters used in DNAUTILUS
The physical model we adopted for the simulation is similar to that of Sipilä et al. (2015).We chose to explore a broader parameter space, including temperature (T ), density (n H ), and the cosmic-ray (CR) ionization rate (ζ).The physical parameter spaces used in this paper are listed in Table 4.This choice is justified because isotopic ratios are observed in different regions of Orion KL with varying physical conditions.
The adapted density range covers the expanding shell (∼ 10 3 cm −3 ) and OMC-1 (10 4 − 10 5 cm −3 ) (Abel et al. 2019;Pabst et al. 2019).The eastern side of the Orion KL hot core is primarily exposed to CRs and UV flux from OB massive stars in the Trapezium cluster, while the western side is heavily shielded.Consequently, we can anticipate higher temperatures and secondary UV photon effects from CRs on the east side, which could impact the chemical outcomes.We ran models with different CR ionization rates from 10 −17 s −1 to 10 −15 s −1 to simulate ionization fronts in the hot core generated by the OB massive stars in the Trapezium cluster and the shielded regions.
We maintained a fixed visual extinction of 300 mag, which approximately corresponds to the average (330 mag) and median (296 mag) values obtained from the analysis of five positions (Table 1) presented in this paper.We kept the gas temperature and grain temperature constant.

Comparisons between the observations and models
We compared the modeled results with the temperatures of 10 -50 K to the observational results and confirmed that the models at 40 K reproduce the observational results best.The temperature of 40 K agrees with the dust temperature in this region (see Table 1).We thus show and discuss the modeled results at 40 K in the following parts.Figure 5 shows the modeled DCN/DCO + abundance ratios along with the observed values at 0.2 Myr, which corresponds to the age of the expanding bubble constrained by Pabst et al. (2019).We found that the CR ionization rates larger than 5×10 −17 s −1 cannot reproduce the observed abundance ratios and we exclude them in Figure 5. Major production and destruction reactions of DCN and DCO + for different densities and CR ionization rates at 40 K are summarized in Table 5.
Dissociative recombination reactions between molecular ions, DCNH + and HDCN + , and electrons (e − ) produce DCN in low-density regions (10 3 − 10 4 cm −3 ): In such conditions, DCO + is efficiently destroyed through dissociative recombination with e − : DCO + + e − → CO + D. (12) Consequently, the DCN/DCO + abundance ratio increases in the low-density regions where the electron density is high.In the dense regions (≳ 10 4 cm −3 ), freezing of DCN onto dust grains becomes important and decreases the gas-phase DCN abundances.This decreases the DCN/DCO + ratios in high density regions.In all these environments, the production of DCO + is primarily driven by the deuterium exchange reaction involving HCO + (reaction (9); HCO + + D → DCO + + H).We found that density is the key parameter to produce the observed differences in the DCN/DCO + abundance ratio.The CR ionization rate has little effect on the ratio in our assumed conditions.The low observed DCN/DCO + ratios (∼ 4 − 7) favor models with high densities of 10 4 − 10 5 cm −3 , which are consistent with the density in OMC-1 (Pabst et al. 2019).This situation is consistent with observational results, where gas components with V LSR ≈ 8 km s −1 , unaffected by the expanding shell, exhibit the low DCN/DCO + ratios.On the other hand, the high observed abundance ratios (DCN/DCO + ≈ 22 − 30) are well replicated by models with low densities of 10 3 − 10 4 cm −3 .Such lower densities agree with conditions of the expanding shell (Abel et al. 2019).This is consistent with the observational results that the high abundance ratios are found in V LSR ≈ 9.2 − 11.6 km s −1 , likely influenced by the expanding shell generated by the most massive star in the Trapezium cluster.

CONCLUSIONS
We have obtained the 1.5 ′ (0.18 pc)-square maps of the DCN and DCO + lines toward the Orion KL region with the B4R system installed on the Large Millimeter Telescope.The main conclusions of this paper are as follows.
1.The DCN emission shows a peak at the hot core position, whereas that of DCO + does not emit there.The DCO + emission is enhanced in the west side from the hot core which is well shielded from the UV photons and CRs from the OB stars in the Trapezium cluster.
2. We have derived the DCN/DCO + abundance ratios at the three representative positions where both species have been detected.The ratios are divided into two groups; DCN/DCO + ≈ 4 − 7 and ≈ 22 − 30.The former is associated with the gas components of V LSR ≈ 7.5 − 8.7 km s −1 , which is consistent with the systemic velocity of OMC-1.
On the other hand, the gas components with the higher abundance ratios have velocity components of 9.2-11.6km s −1 , which are likely affected by the expanding shell produced by the most massive star in the Trapezium cluster.
3. The two distinct abundance ratios can be reproduced by our models.Models featuring highdensity provide a better explanation for the low observed DCN/DCO + ratio.On the other hand, models with low-density reproduce the high observed abundance ratio.This condition is consistent with the presence of an expanding shell caused by the OB stars.This paper makes use of data taken by the Large Millimeter Telescope Alfonso Serrano (LMT) in Mexico.The LMT project is a joint effort of the Instituto Nacional de Astrófisica, Óptica, y Electrónica (INAOE) and the University of Massachusetts at Amherst (UMASS).We also appreciate the support of the technical staff and the support scientists of the LMT during the commissioning campaign of the B4R.Data analysis was in part carried out on the Multi-wavelength Data Analysis System operated by the Astronomy Data Center (ADC), National Astronomical Observatory of Japan.K.T. is supported by JSPS KAKENHI grant No. JP20K14523 and 21H01142.T.Y. is supported by JSPS KAKENHI grant No. JP22J22889.We thank the anonymous referee whose comments helped improve the paper.

Figure 1 .Figure 2 .
Figure 1.Integrated-intensity of DCN and DCO + (moment 0 maps; top panels), peak intensity maps of DCN and DCO + (moment 8 maps; middle panels), and velocity map of DCN (moment 1 map; bottom panel) toward the Orion KL region.The black contours in the DCN moment 1 map show the same one as shown in its moment 0 map.The black crosses indicate five representative positions where we analyze spectra.The names of each position are indicated in the DCN moment 1 map.

Figure 3 .
Figure3.Spectra of the DCN J = 2 − 1 line (left panels) and DCO + J = 2 − 1 line (right panels) at the five positions.Black and red curves are observed spectra and the best fitting results with an assumed excitation temperature of 40 K, respectively.

Figure 4 .
Figure 4. Comparison of the N (DCN)/N (DCO + ) ratio and the velocity component.The velocities are values derived by fitting the DCO + line.

4. 2 .
Comparison with the chemical modeling 4.2.1.Chemical parameters used in DNAUTILUS We compare the observed abundances of DCN and DCO + as well as the DCN/DCO + ratios around the Orion KL region with modeled values.To simulate the observed values around Orion KL, we have utilized our three-phase gas-grain chemical code called DNAUTILUS.2.0.This code calculates the timedependent abundances of chemical species, including deuterated ones, by solving rate equations.It is an updated version of the DNAUTILUS.1.0astrochemical code introduced in Majumdar et al. (

Figure 5 .
Figure5.Comparisons between the modeled DCN/DCO + ratios at T = 40 K with different CR ionization rates and the observed values.Shadow areas indicate the observed error ranges.The orange mesh area depicts the low-density region with a high DCN/DCO + ratio, while the blue mesh area depicts the high-density area with a low DCN/DCO + ratio.

Table 2 .
Results of the MCMC analysisPosition N (×10 12 cm −2 ) FWHM (km s −1 ) VLSR (km s −1 ) N (×10 12 cm −2 ) FWHM (km s −1 ) VLSR (km s −1 ) DCN/DCO + Note-The figures in parentheses refer to the standard deviation derived from the MCMC analysis, expressed in units of the last significant digits.The errors of the column densities include the absolute calibration error of 10%.

Table 3 .
Fractional abundances of DCN and DCO + at each position

Table 4 .
Physical parameters and initial abundances