A&A 456, 1037-1043 (2006)
DOI: 10.1051/0004-6361:20040294

HCN and HNC mapping of the protostellar core Chamaeleon-MMS1[*]

P. P. Tennekes1,2 - J. Harju1 - M. Juvela1 - L. V. Tóth3,4


1 - Observatory, PO Box 14, 00014 University of Helsinki, Finland
2 - Julius Institute, Utrecht University, The Netherlands
3 - Department of Astronomy of the Loránd Eötvös University, Pázmány Péter sétány 1, 1117 Budapest, Hungary
4 - Konkoly Observatory, PO Box 67, 1525 Budapest, Hungary

Received 18 February 2004 / Accepted 19 June 2006

Abstract
Aims. The purpose of this study is to investigate the distributions of the isomeric molecules HCN and HNC and estimate their abundance ratio in the protostellar core Cha-MMS1 located in Chamaeleon I.
Methods. The core was mapped in the J=1-0 rotational lines of \ensuremath {{\rm HCN}}, \ensuremath {{\rm HNC}}, and \ensuremath {{\rm HN^{13}C}}. The column densities of \ensuremath{{\rm H^{13}CN}}, \ensuremath {{\rm HN^{13}C}}, \ensuremath{{\rm H^{15}NC}} and \ensuremath {{\rm NH_3}} were estimated towards the centre of the core.
Results. The core is well delineated in all three maps. The kinetic temperature in the core, derived from the \ensuremath {{\rm NH_3}} (1, 1) and (2, 2) inversion lines, is $12.1 \pm 0.1$ K. The \ensuremath {{\rm HN^{13}C}}/ \ensuremath{{\rm H^{13}CN}} column density ratio is between 3 and 4, i.e. similar to values found in several other cold cores. The \ensuremath {{\rm HN^{13}C}}/ \ensuremath{{\rm H^{15}NC}} column density ratio is $\sim $7. In case no 15N fractionation occurs in \ensuremath {{\rm HNC}} (as suggested by recent modelling results), the \ensuremath {{\rm HNC}}/ \ensuremath {{\rm HN^{13}C}} abundance ratio is in the range 30-40, which indicates a high degree of 13C fractionation in \ensuremath {{\rm HNC}}. Assuming no differential 13C fractionation the \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} abundances are estimated to be ${\sim} 7 \times 10^{-10}$ and ${\sim} 2\times 10^{-9}$, respectively, the former being nearly two orders of magnitude smaller than that of \ensuremath {{\rm NH_3}}. Using also previously determined column densities in Cha-MMS1, we can put the most commonly observed nitrogenous molecules in the following order according to their fractional abundances: $\chi(\ensuremath{{\rm NH_3}} ) > \chi(\ensuremath{{\rm HC_3N}} ) > \chi(\ensure...
...{\rm HNC}} ) > \chi(\ensuremath{{\rm HCN}} ) > \chi(\ensuremath{{\rm N_2H^+}} )$.
Conclusions. The relationships between molecular abundances suggest that Cha-MMS1 represents an evolved chemical stage, experiencing at present the "late-time'' cyanopolyyne peak. The possibility that the relatively high HNC/HCN ratio derived here is only valid for the 13C isotopic substitutes cannot be excluded on the basis of the present and other available data.

Key words: ISM: individual objects: Chamaeleon-MMS1 - ISM: abundances - ISM: molecules

1 Introduction

Hydrogen cyanide, \ensuremath {{\rm HCN}}, and its metastable isomer hydrogen isocyanide, \ensuremath {{\rm HNC}}, are commonly used tracers of dense gas in molecular clouds. The \ensuremath {{\rm HNC}} column density has been observed to be similar or larger than that of \ensuremath {{\rm HCN}} in cold regions, whereas in warm GMC cores the [HNC]/[HCN] ratio is typically much smaller than unity (e.g., Goldsmith et al. 1981; Churchwell et al. 1984; Irvine & Schloerb 1984; Schilke et al. 1992; Hirota et al. 1998).

Gas-phase chemistry models, including both ion-molecule and neutral-neutral reactions, predict that the \ensuremath {{\rm HNC}}/ \ensuremath {{\rm HCN}} abundance ratio should be close to unity in cold gas (see e.g. Herbst et al. 2000, where also previous work is reviewed). The two isomers are thought to form primarily via the dissociative recombination of HCNH+, which in turn is produced by a reaction between \ensuremath {{\rm NH_3}} and C+. The recombination reaction regulates the abundance ratio despite the initial production mechanisms of the two isomers, provided that they are efficiently protonated via reactions with ions, such as H3+and HCO+ (Churchwell et al. 1984; Brown et al. 1989). Recent theoretical studies corroborate the equal branching ratios for \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} in dissociative recombination (Hickman et al. 2005; Ishii et al. 2006).

\ensuremath {{\rm HNC}}/ \ensuremath {{\rm HCN}} column density ratios clearly larger than unity have been reported: Churchwell et al. 1984 find that [ \ensuremath {{\rm HNC}}]/[ \ensuremath {{\rm HCN}}] may be as high as $\sim $10 in some dark clouds cores, and Hirota et al. (1998) find a ratio of $\sim $5 in L1498. The first result quoted is uncertain because it depends on the assumed degree of 13C fractionation in HNC, whereas the second value refers to \ensuremath {{\rm HN^{13}C}} and \ensuremath{{\rm H^{13}CN}} measurements. Nevertheless, both results suggest that the \ensuremath {{\rm HCN}}/ \ensuremath {{\rm HNC}} chemistry is not fully understood. Further determinations of the relative abundances of these molecules towards cores with well-defined chemical and physical characteristics seem therefore warranted.

In this paper we report on observations of \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} and some of their isotopologues towards the dense core Cha-MMS1, which is situated near the reflection nebula Ced 110 in the central part of the Chamaeleon I dark cloud. We also derive the gas kinetic temperature in the core by using the (1, 1) and (2, 2) inversion lines of \ensuremath {{\rm NH_3}}. Cha-MMS1 was discovered in the 1.3 mm dust continuum by Reipurth et al. (1996), and was studied in several molecules by Kontinen et al. (2000). The object is embodied in one of the most massive "clumps'' in Chamaeleon I identified in the large scale \ensuremath{{\rm C^{18}O}} survey of Haikala et al. 2005 (clump No. 3, ${\sim} 12~\ensuremath{{M}_\odot} $). Reipurth et al. (1996) suggested that Cha-MMS1 contains a Class 0 protostar. A FIR source (Ced 110 IRS10) was detected near the centre of the core by Lehtinen et al. (2001). However, no centimetre continuum nor near-infrared sources have been found in its neighbourhood (Lehtinen et al. 2003, and references therein), and the core therefore represents a very early stage of star formation. At the same time, its chemical composition has probably reached an advanced stage (Kontinen et al. 2000).

In Sect. 2 of this paper we describe the observations. In Sect. 3 we present the direct observational results in the form of spectra and maps, and derive the \ensuremath{{\rm H^{13}CN}}, \ensuremath {{\rm HN^{13}C}}, and \ensuremath{{\rm H^{15}NC}} column densities. By combining our results with previous observations we derive in Sect. 4 the fractional abundances of several nitrogen containing molecules in Cha-MMS1, and discuss briefly the chemical state of the core. In Sect. 5 we summarize our results.

2 Observations

The J=1-0 transitions of \ensuremath {{\rm HCN}}, \ensuremath{{\rm H^{13}CN}}, \ensuremath {{\rm HNC}}, \ensuremath {{\rm HN^{13}C}} and \ensuremath{{\rm H^{15}NC}} at about 90 GHz were observed with the Swedish-ESO-Submillimeter-Telescope, SEST, located on La Silla in Chile. The observations took place in December 1990. A 3 mm dual polarization Schottky receiver was used in the frequency switching mode. The backend was a 2000 channel acousto optical spectrometer with an 86 MHz bandwidth. The velocity resolution with this configuration is about $0.15~\ensuremath{{\rm km~s^{-1}}} $ at 90 GHz. Typical values for the single sideband system temperatures ranged from 350 K to 450 K. The half-power beam width of the SEST is about $55\hbox{$^{\prime\prime}$ }$ at the frequencies observed. Calibration was achieved by the chopper wheel method. Pointing was checked every 2-3 h towards a nearby SiO maser source. We estimate the pointing accuracy to have been better than  $5\hbox{$^{\prime\prime}$ }$ during the observations. The focusing was done using a strong SiO maser.

The (1, 1) and (2, 2) inversion lines of \ensuremath {{\rm NH_3}} at about 23 GHz were observed with the Parkes 64-m radio telescope of the Australian Telescope National Facility (ATNF) in April 2000. The FWHM of the telescope is about  $80\hbox{$^{\prime\prime}$ }$ at this frequency. The two lines were observed simultaneously in two orthogonal polarizations in the frequency switching mode, using the 8192 channel autocorrelator. The calibration was checked by observing standard extragalactic calibration sources at different elevations, and by comparing the ammonia line intensities towards some strong galactic sources also visible from the Effelsberg 100-m telescope, taking the different beam sizes into account.

The electric dipole moments of the molecules and the line frequencies used in this paper are listed in Table 1. Except for the frequencies of the \ensuremath {{\rm HNC}} and \ensuremath {{\rm HN^{13}C}} hyperfine transitions adopted from Frerking et al. (1979), the parameters are obtained from the Jet Propulsion Laboratory molecular spectroscopy data base (http://spec.jpl.nasa.gov), where also the original references are given. The hyperfine components of \ensuremath {{\rm NH_3}} are not listed. The frequencies given represent the centre frequencies of 18 and 21 individual components of the (1, 1) and (2, 2) lines, respectively, distributed over a range of about 4 MHz.

Table 1: The permanent dipole moments of the molecules and the rest frequencies of the observed transitions.


  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{0294fig1.eps}
\end{figure} Figure 1: Integrated \ensuremath {{\rm HCN}}(J=1-0) line intensity map. The integration is performed over all three hyperfine structure components (F=1-1, F=2-1 and F=0-1). The map is centered on the position of Cha-MMS1a (denoted by a cross), $\alpha_{2000} = 11 \ensuremath{^{\rm h}}06 \ensuremath{^{\rm m}} 31\hbox{$.\!\!^{\rm s}$ }7$, $\delta _{2000} = -77 \hbox {$^\circ $ }23\hbox {$^\prime $ }32\hbox {$^{\prime \prime }$ }$, and the offsets are in arcseconds. The contour levels indicated with ten shades of grey and their values can be read in the bottom right corner. The intensity unit is K km s-1. The positions used for the column density estimates ( $-20\hbox {$^{\prime \prime }$ },20\hbox {$^{\prime \prime }$ }$) and $T_{\rm ex}$estimates (see Table 2) are marked with a black dot and open circles, respectively. The grid spacing used in the map is  $30\hbox {$^{\prime \prime }$ }$.
Open with DEXTER

3 Maps and column densities

3.1 Maps

The integrated \ensuremath {{\rm HCN}}(J=1-0), \ensuremath {{\rm HNC}}(J=1-0) and \ensuremath {{\rm HN^{13}C}}(J=1-0)intensity maps of Cha-MMS1 are presented as contour plots in Figs. 1-3, respectively. The values of the contour levels are given in the bottom right of each figure. The position of the 1.3 mm dust emission maximum, Cha-MMS1a (Reipurth et al. 1996), is chosen as the map centre. The coordinates of Cha-MMS1a are $\alpha_{2000} = 11 \ensuremath{^{\rm h}}06 \ensuremath{^{\rm m}} 31\hbox{$.\!\!^{\rm s}$ }7$, $\delta _{2000} = -77 \hbox {$^\circ $ }23\hbox {$^\prime $ }32\hbox {$^{\prime \prime }$ }$. The grid spacing used in the maps was $30\hbox {$^{\prime \prime }$ }$. Also indicated in these figures are the location of the \ensuremath {{\rm HN^{13}C}} maximum ( $\Delta\alpha,\Delta\delta = -20\hbox{$^{\prime\prime}$ },20\hbox{$^{\prime\prime}$ }$, black dot), and the positions used for the \ensuremath {{\rm HCN}}(J=1-0) excitation temperature estimates in the LTE method (open circles, see below). Long integrations towards the \ensuremath {{\rm HN^{13}C}} maximum were made besides \ensuremath {{\rm HNC}} and \ensuremath {{\rm HCN}}  also in the J=1-0 lines of \ensuremath{{\rm H^{13}CN}} and \ensuremath{{\rm H^{15}NC}}. The five spectra obtained towards this position are shown in Fig. 5.

The three maps have differences, but they all show a relatively compact distribution peaking near Cha-MMS1a. The \ensuremath {{\rm HN^{13}C}} emission with the lowest opacity reflects probably best the distribution dense gas.


  \begin{figure}
\par\includegraphics[width=7.6cm,clip]{0294fig2.eps}
\end{figure} Figure 2: Integrated \ensuremath {{\rm HNC}}(J=1-0) line intensity map. The coordinates and markers as in Fig. 1.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{0294fig3.eps}
\end{figure} Figure 3: Integrated \ensuremath {{\rm HN^{13}C}}(J=1-0) line intensity map. The coordinates and markers as in Fig. 1.
Open with DEXTER

3.2 LTE column densities

The column densities of the molecules have been determined with two different methods. First, we have derived the optical thicknesses of the lines and the total column densities of the molecules by assuming local thermodynamic equilibrium, LTE. The LTE method involves idealistic assumptions about the homogenity of the cloud. In the next subsection we present modelling of the observed spectra with a Monte-Carlo program using a realistic density distribution.

In the following, we assume that the rotational levels of all five isotopologues are characterized by the same  $T_{\rm ex}$, and that the core is homogeneous, i.e. the same  $T_{\rm ex}$ is valid for different locations. The former assumption can be justified by the similar dipole moments and rotational constants of the molecules. The second assumption infers that the excitation conditions remain constant along the line of sight and also between close-lying positions in the core.

The excitation temperature is derived by fitting first the optical thickness, $\tau$, to the observed hyperfine ratios of a "clean'' \ensuremath {{\rm HCN}} spectrum, i.e. a spectrum with a good S/N and hyperfine intensity ratios in accordance with the LTE assumption. The $T_{\rm ex}$ is solved from the antenna equation using observed peak main beam brightness temperature, $T_{\rm MB}$ and assuming uniform beam filling. The fit is possible towards the four positions listed in Table 2. In other positions with strong signal, like the \ensuremath {{\rm HN^{13}C}} maximum $(-20\hbox {$^{\prime \prime }$ },+20\hbox {$^{\prime \prime }$ })$, or the \ensuremath {{\rm HCN}} maximum at $(-20\hbox{$^{\prime\prime}$ },-20\hbox{$^{\prime\prime}$ })$ the \ensuremath {{\rm HCN}} spectra show "anomalous'' hyperfine ratios. The smallest error is obtained towards the offset ( $40\hbox {$^{\prime \prime }$ },20\hbox {$^{\prime \prime }$ }$) near the edge of the core with $T_{\rm ex}\approx 5$ K. This \ensuremath {{\rm HCN}} spectrum is shown in Fig. 4. The other positions listed in Table 2 give similar values of  $T_{\rm ex}$ but the errors are larger. Also \ensuremath {{\rm HCN}} column density estimates using the values for $\tau$ and $\Delta v$ obtained from hyperfine fits are given in Table 2.

Table 2: Estimates of the \ensuremath {{\rm HCN}}(J=1-0) excitation temperatures and the \ensuremath {{\rm HCN}} column densities towards selected positions using hyperfine structure fits. The coordinates are offsets from the position of Cha-MMS1a.


  \begin{figure}
\par\includegraphics[angle=270,width=7.1cm,clip]{0294fig4.eps}
\end{figure} Figure 4: The \ensuremath {{\rm HCN}}(J=1-0) spectrum towards offset ( $40\hbox {$^{\prime \prime }$ },20\hbox {$^{\prime \prime }$ }$) used to determine the $T_{\rm ex}$. The dotted line represents a fit to the hyperfine structure.
Open with DEXTER

The \ensuremath {{\rm HCN}}, \ensuremath{{\rm H^{13}CN}}, \ensuremath {{\rm HN^{13}C}} and \ensuremath{{\rm H^{15}NC}} column densities towards the \ensuremath {{\rm HN^{13}C}} maximum, $(-20\hbox {$^{\prime \prime }$ },+20\hbox {$^{\prime \prime }$ })$, have been derived from the integrated J=1-0 line intensities and the  $T_{\rm ex}$ estimate derived from the \ensuremath {{\rm HCN}} spectrum shown in Fig. 4. For \ensuremath {{\rm HCN}} the weakest of the well separated hyperfine components, F=0-1, with the least optical thickness has been used in this calculation by taking its relative strength into account. For \ensuremath {{\rm HNC}} we did not attemp to do a similar estimate because the hyperfine components overlap and are very likely optically thick. Overlapping hyperfine components hamper also the \ensuremath {{\rm HN^{13}C}} column density estimate, but their optical thicknesses should be lower than those of \ensuremath {{\rm HNC}} by an order of magnitude or more. \ensuremath{{\rm H^{15}NC}} has no hyperfine structure due to zero nuclear spin of 15N. The results of these calculations are presented in Table 3.


  \begin{figure}
\par\includegraphics[width=6.5cm,clip]{0294fig5.eps}
\end{figure} Figure 5: Observations used for the column density estimates. The position corresponds to the \ensuremath {{\rm HN^{13}C}} maximum, which lies at the offset $(-20\hbox {$^{\prime \prime }$ },+20\hbox {$^{\prime \prime }$ })$ from the dust continuum peak Cha-MMS1a (marked with a black dot in Figs. 1-3).
Open with DEXTER

Table 3: Column densities towards the \ensuremath {{\rm HN^{13}C}} maximum derived using the LTE method. The assumed value for $T_{\rm ex}$ is 5.2+1.0-0.6 K.   a The intensity value is for the weakest hyperfine component F=0-1 which comprises one ninth of the total intensity of the J=1-0 line.

The following column density ratios can be derived from Table 3:

$\displaystyle \frac{N[\ensuremath{{\rm HN^{13}C}} ]}{N[\ensuremath{{\rm H^{13}CN}} ]}$ = $\displaystyle 3.3~ \pm0.4$  
$\displaystyle \frac{N[\ensuremath{{\rm HN^{13}C}} ]}{N[\ensuremath{{\rm H^{15}NC}} ]}$ = $\displaystyle 7.4~ \pm 1.8$  
$\displaystyle \frac{N[\ensuremath{{\rm HCN}} ]}{N[\ensuremath{{\rm H^{13}CN}} ]}$ = $\displaystyle 17 ~ \pm 2 .$  

The \ensuremath {{\rm HCN}}/ \ensuremath{{\rm H^{13}CN}} column density ratio, 17, is much lower than the isotopic 12C/13C ratios usually observed in nearby molecular clouds or local diffuse ISM ($\sim $60, e.g. Lucas & Liszt 1998; the terrestrial value is 89), and is likely to be affected by absorption or scattering in the less dense envelope of the core. Radiative transfer effects are even more obvious in the case of the \ensuremath {{\rm HNC}}/ \ensuremath {{\rm HN^{13}C}} intensity ratio, 1.7. We do not use therefore the values derived from the \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} spectra towards the core centre. The \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} spectra are discussed at the end of Sect. 4.

The ammonia column density, $N(\ensuremath{{\rm NH_3}} )$, and the gas kinetic temperature, $T_{\rm kin}$, were derived using hyperfine fits to the \ensuremath {{\rm NH_3}} inversion lines arising from the rotational states (J,K) = (1,1) and (2,2). The spectra are shown in Fig. 6. The standard method of analysis described e.g. in Ho et al. (1979) gave the following values: $N(\ensuremath{{\rm NH_3}} ) = (1.4 \pm
0.2) \times 10^{15}$ \ensuremath{{\rm cm^{-2}}} and $T_{\rm kin} = 12.1 \pm 0.1$ K. Here it is assumed that the relative populations of all metastable levels of both ortho- and para- \ensuremath {{\rm NH_3}} are determined by thermal equilibrium at the (2,2)/(1,1) rotational temperature.


  \begin{figure}
\par\includegraphics[width=7cm,clip]{0294fig6.ps}
\end{figure} Figure 6: \ensuremath {{\rm NH_3}} (J,K) = (1,1) and (2,2) inversion spectra towards Cha-MMS1a.
Open with DEXTER

3.3 Monte-Carlo modelling

\ensuremath{{\rm H^{15}NC}}, \ensuremath {{\rm HN^{13}C}}, \ensuremath {{\rm HCN}} and \ensuremath{{\rm H^{13}CN}} spectra were calculated with the Monte Carlo radiative transfer program of Juvela (1997). A constant kinetic temperature and a spherically symmetric density distribution, $n\sim r^{-1.5}$, were assumed. The selected density power law conforms with observational results towards star-forming cores (e.g., Tatematsu et al. 2004). In order to avoid divergence in the cloud centre, we assume a constant density within a distance corresponding to 6% of the outer radius of the model cloud. Compared with the beam size, the constant density region is small and the results are not very sensitive to the selected value of its radius. The spectra were convolved with a Gaussian that corresponds to the size of the beam used in the observations.The simulated spectra were fitted to observations by adjusting four parameters: the outer radius of the core, the scaling of gas density, and the fractional abundances of \ensuremath{{\rm H^{13}CN}} and \ensuremath {{\rm HN^{13}C}}. The \ensuremath {{\rm HCN}} fit was optimized for the weakest hyperfine component F=0-1. The abundances of the other species were determined by assumed fixed ratios [12C]/[13C] = 20, as suggested by the LTE column density determination, and [14N]/[15N] = 280 (terrestrial). The resulting fractional \ensuremath{{\rm H^{13}CN}} and \ensuremath {{\rm HN^{13}C}} abundances do not, however, depend on these assumptions about isotopic ratios. The calculations were done with $T_{\rm kin}$ = 12 K, which was derived from the ammonia observations, and $T_{\rm kin}$ = 20 K.The lower temperature resulted in better fit although the difference in the $\chi^2$ values was no more than 25%. The results for the 12 K are shown in Fig. 7. The outer radius of the model is 60 $^{\prime\prime}$. The volume density increases from $2.3\times 10^4$ cm-3 at the outer edge to $1.4\times 10^6$ cm-3 in the centre, and the molecular hydrogen column density towards the centre is $N({\rm {H}}_2)=1.9\times 10^{22}$ cm-2.

The fractional \ensuremath {{\rm HN^{13}C}} and \ensuremath{{\rm H^{13}CN}} abundances resulting from the fit are $2.1\times 10^{-10}$ and $5.5\times 10^{-11}$, respectively. The \ensuremath {{\rm HN^{13}C}}/ \ensuremath{{\rm H^{13}CN}} abundance ratio obtained from the simulations is thus 3.8, i.e. close to the value obtained from the LTE method.

The largest discrepancy between the observed and modelled profiles is seen in \ensuremath {{\rm HCN}}. It can be explained if the optical depth is larger than in this particular model. For example, a foreground layer of cold gas would readily decrease the intensity of the \ensuremath {{\rm HCN}} lines and would bring the ratios between the hyperfine components closer to the observed values. The hyperfine ratios are, however, sensitive to the actual structure of the source (Gonzáles-Alfonso & Cernicharo 1993).

The optical thickness of the \ensuremath {{\rm HN^{13}C}} main component is $\sim $0.7. This moderate value suggests that the method used in the previous subsection results in a slight underestimatation of column densities. The \ensuremath{{\rm H^{13}CN}} and \ensuremath {{\rm HN^{13}C}} column densities (fractional abundances multiplied by $N({\rm {H}}_2)$) are indeed about two times larger than those presented in the previous section: $N(\ensuremath{{\rm H^{13}CN}} ) = 1.1 \times 10^{12}
~\ensuremath{{\rm cm^{-2}}} $, $N(\ensuremath{{\rm HN^{13}C}} ) = 4.0 \times 10^{12}~ \ensuremath{{\rm cm^{-2}}} $. These results suggest, however, that the LTE method does not lead to crude underestimates of the column densities of rarer isotopologues of HCN and HNC.


  \begin{figure}
\par\includegraphics[width=6.5cm,clip]{0294fig7.eps}
\end{figure} Figure 7: Simulated spectra fitted to the observations towards the \ensuremath {{\rm HN^{13}C}} maximum. The assumed kinetic temperature, $T_{\rm kin}=12$ K, is based on ammonia observations. The following isotopic ratios were used in the model calculations: [12C]/[13C] = 20, and [14N]/[15N] = 280.
Open with DEXTER

4 Discussion

4.1 Fractional abundances

The present observations and the LTE method give reasonable estimates for the column densities of \ensuremath{{\rm H^{13}CN}}, \ensuremath {{\rm HN^{13}C}}, \ensuremath{{\rm H^{15}NC}} and \ensuremath {{\rm NH_3}} near the centre of Cha-MMS1. Unfortunately the column densities of the common isotopologues of \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} cannot be directly determined because of the large optical thickness of their lines. \ensuremath{{\rm H^{15}NC}} is, however, likely to be useful for this purpose. According to the modelling results of Terzieva & Herbst (2000) hardly any 15N fractionation occurs in dense interstellar clouds. The isotopic 14N/15N ratio lies probably between 240 and 280, i.e. the values derived for the local ISM (Lucas & Liszt 1998), and the terrestrial atmosphere. Using this range and the \ensuremath{{\rm H^{15}NC}} column density listed in Table 3 we obtain $N(\ensuremath{{\rm HNC}} )
\sim 6 \times 10^{13} ~ \ensuremath{{\rm cm^{-2}}} $. Assuming that no differential 13C fractionation occurs, we get from the \ensuremath {{\rm HN^{13}C}}/ \ensuremath{{\rm H^{13}CN}} ratio derived above that $N(\ensuremath{{\rm HCN}} ) \sim 2 \times 10^{13}$ \ensuremath{{\rm cm^{-2}}}.

The dust continuum emission provides the most reliable estimate for the \ensuremath{{\rm H_2}} column density. Reipurth et al. (1996) measured towards Cha-MMS1 a total 1.3 mm flux density of 950 mJy from a nearly circular region with a diameter of about  $50\hbox{$^{\prime\prime}$ }$ (bordered by the lowest white contour in Fig. 8). As this size is similar to the SEST beamsize at 3 mm, the resulting average intensity and average  $N(\ensuremath{{\rm H_2}} )$ can be compared with the column densities derived here. The dust temperature as derived from ISO observations by Lehtinen et al. (2001) is 20 K (Cha-MMS1 = Ced 110 IRS10). Due to crowding of sources in this region and the limited angular resolution of ISOPHOT, this value may be affected by emission from other sources in its vicinity (see Fig. 2 of Lehtinen et al. 2001). Therefore, value 20 K is probably an upper limit. A lower limit is provided by the gas kinetic temperature $\sim $12 K derived above. By using the two temperatures and Eq. (1) of Motte et al. 1998 with $\kappa_{\rm 1.3 mm} = 0.01$ g  \ensuremath{{\rm cm^{-2}}} appropriate for circumstellar envelopes around Class 0 and Class I protostars, we get $N(\ensuremath{{\rm H_2}} ) = 1.9 {-} 4.0 \times 10^{22}~ {\ensuremath{{\rm cm^{-2}}} }$. The corresponding mass range of the dust core obtained from the total flux given by Reipurth et al. (1996) is $0.4{-}0.8~ \ensuremath{{M}_\odot} $.

The median of the reasonable H2 column density range, $N(\ensuremath{{\rm H_2}} ) =
3.0 \times 10^{22}~ {\ensuremath{{\rm cm^{-2}}} }$, yields the following values for the fractional \ensuremath {{\rm HCN}}, \ensuremath {{\rm HNC}} and \ensuremath {{\rm NH_3}} abundances: $\chi(\ensuremath{{\rm HCN}} ) \sim 7 \times 10^{-10}$, $\chi(\ensuremath{{\rm HNC}} ) \sim 2 \times 10^{-9}$, and $\chi(\ensuremath{{\rm NH_3}} ) \sim 5 \times 10^{-8}$. The column densities of two other nitrogen bearing molecules towards Cha-MMS1 were determined by Kontinen et al. (2000): $N(\ensuremath{{\rm HC_3N}} ) = (4.5 \pm 0.8) \times 10^{14}$  \ensuremath{{\rm cm^{-2}}}, and $N(\ensuremath{{\rm N_2H^+}} ) = (1.4
\pm 0.1) \times 10^{13}$  \ensuremath{{\rm cm^{-2}}}, which imply the fractional abundances $\chi(\ensuremath{{\rm HC_3N}} ) \sim 2 \times 10^{-8}$ and $\chi(\ensuremath{{\rm N_2H^+}} ) \sim 5 \times 10^{-10}$.

We find that both \ensuremath {{\rm HCN}} and \ensuremath{{\rm N_2H^+}} are about 100 times less abundant than \ensuremath {{\rm NH_3}}, whereas \ensuremath {{\rm HC_3N}} is almost equally abundant as \ensuremath {{\rm NH_3}}. The derived \ensuremath{{\rm N_2H^+}}/ \ensuremath {{\rm NH_3}} column density ratio is in accordance with previous observational results (e.g., Tafalla et al. 2002; Hotzel et al. 2004). The abundance ratios are discussed in Sect. 4.3. The \ensuremath{{\rm C^{18}O}} column density in this direction is $(1.5 \pm 0.1) \times 10^{15}$ \ensuremath{{\rm cm^{-2}}} (Haikala et al. 2005). The fractional \ensuremath{{\rm C^{18}O}} abundance is therefore $5 \times 10^{-8}$ or smaller. This upper limit falls below the "standard'' value ${\sim} 2\times 10 ^{-7}$ of Frerking et al. (1982), and suggests CO depletion.

4.2 Comparison with previous maps


  \begin{figure}
\par\includegraphics[width=7.2cm,clip]{0294fig8.ps}
\end{figure} Figure 8: The \ensuremath {{\rm HN^{13}C}} map (dark contours) superposed on the \ensuremath {{\rm HC_3N}} map from Kontinen et al. (2000) (greyscale) and the 1.3 mm dust continuum map of Reipurth et al. (1996) (white contours).
Open with DEXTER

The superposition of our \ensuremath {{\rm HN^{13}C}} map, the \ensuremath {{\rm HC_3N}} map of Kontinen et al. (2000), and the 1.3 mm dust continuum map of Reipurth et al. (1996) is shown in Fig. 8. The two molecular line maps have similar angular resolution (about  $60\hbox{$^{\prime\prime}$ }$) whereas the resolution of the continuum map is  $22\hbox{$^{\prime\prime}$ }$. Both molecules show rather compact, symmetric distributions around the dust continuum source, \ensuremath {{\rm HC_3N}} agreeing slightly better with dust than \ensuremath {{\rm HN^{13}C}}. Nevertheless, the dense core can be clearly localized with the aid of \ensuremath {{\rm HN^{13}C}} as well as \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}}. This is in contrast with SO, which shows two peaks on both sides of the core (Kontinen et al. 2000), and with C18O which has a flat distribution in this region (Haikala et al. 2005), probably due to depletion. This result agrees with the results of Hirota et al. (2003) and Nikolic et al. (2003), which suggest that \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} are less affected by the accretion onto dust grains than some other commonly observed molecules.

4.3 The chemical state of Cha-MMS1 and the HNC/HCN abundance ratio

CO depletion and a large \ensuremath {{\rm NH_3}} abundance point towards an advanced stage of chemical evolution. The \ensuremath {{\rm NH_3}}/ \ensuremath{{\rm N_2H^+}} abundance ratio ${\sim}
100$ agrees with model predictions for dense cores, but does not put severe constraints on the evolutionary stage (Aikawa et al. 2005). According to the observational results of Hotzel et al. (2004) this value should be more appropriate for a star-forming than for a prestellar core. However, due to slightly different beamsizes the accuracy of the present determination of the \ensuremath {{\rm NH_3}}/ \ensuremath{{\rm N_2H^+}} ratio is not sufficient for distinguishing between the two cases. Kontinen et al. (2000) suggested that the large \ensuremath {{\rm HC_3N}} abundance observed in Cha-MMS1 indicates either a so called "late-time'' cyanopolyyne peak associated with a large degree of depletion according to the model of Ruffle et al. (1997), or carbon chemistry revival by the influence of a newly born, so far undetected star. A large ammonia abundance can be understood in either case. Ammonia production benefits from the diminishing of CO and can resist depletion longer (e.g. Nejad & Wagenblast 1999). On the other hand, ammonia can possibly form on grain surfaces and be released into the gas phase in shocks (Nejad et al. 1990; Aikawa et al. 2003). As there is no clear evidence for star-cloud interaction in Cha-MMS1, the alternative that the core - or the portion traced by the present observations - represents matured pre-stellar chemistry is more likely.

The behaviour of \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} under these circumstances is not clear. The main gas-phase production pathway via the HCNH+, which is effective especially at later stages, requires \ensuremath {{\rm NH_3}} and C+ions. The latter can be produced also in deep interiors of molecular clouds via the cosmic ray ionization of C atoms, or in the reaction between CO and the He+ ion, the latter mechanism becoming less important with CO depletion (Nejad & Wagenblast 1999). To our knowledge the production of \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} on interstellar grains and their desorption has not been studied.

The models of Aikawa et al. (2005) for collapsing prestellar cores include the gas-phase production of \ensuremath {{\rm HCN}}, \ensuremath {{\rm HC_3N}}, \ensuremath {{\rm NH_3}}, and \ensuremath{{\rm N_2H^+}}. Although the order of fractional abundances, $\chi(\ensuremath{{\rm NH_3}} ) > \chi(\ensuremath{{\rm HCN}} )
> \chi(\ensuremath{{\rm N_2H^+}} )$, and the \ensuremath {{\rm NH_3}}/ \ensuremath{{\rm N_2H^+}} abundance ratio derived from the present observations are well reproduced by this model, especially for the gravity-to-pressure ratio $\alpha=4$ impying a rapid collapse and less marked depletion, the predicted \ensuremath {{\rm HCN}} column densities relative to \ensuremath {{\rm NH_3}} and \ensuremath{{\rm N_2H^+}} are an order of magnitude larger than those observed here. A likely explanation is that \ensuremath {{\rm HCN}} is not as centrally peaked than \ensuremath {{\rm NH_3}} and \ensuremath{{\rm N_2H^+}} (see Fig. 3 of Aikawa et al. 2005). Therefore a large fraction of \ensuremath {{\rm HCN}} resides in on the outskirts of the core, where the density is not sufficient for the collisional excitation of its rotational levels, and the molecule is not seen in emission (see also the next section). Another discrepancy is that the model predicts much lower \ensuremath {{\rm HC_3N}} column densities relative to ammonia than observed here. The agreement is slightly better for the model with overwhelmingly large gravitational potential ($\alpha=4$) than for the slowly evolving model near hydrostatic equilibrium ( $\alpha=1.1$) which results in very high degrees of depletion. Unfortunately the spectral resolution and S/N ratio of the present data do not allow to study the core kinematics.

As compared to the "principal'' N-bearing molecules \ensuremath {{\rm NH_3}} and \ensuremath {{\rm HC_3N}}, \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} have clearly smaller abundances, and probably only marginal importance to the nitrogen chemistry in the core. Nevertheless, the fact that \ensuremath {{\rm HNC}}/ \ensuremath {{\rm HCN}} abundance ratio is once again observed to be larger than unity remains a problem. Talbi & Herbst (1998) have shown that even though the reaction between \ensuremath {{\rm NH_3}} and C+ may produce a significant amount of the metastable ion ${\rm H_2NC^+}$, which should produce only \ensuremath {{\rm HNC}} in electron recombination, the energy released in the former reaction leads to the efficient transformation into the linear isomer ${\rm
HCNH}^+$ before the recombination. Likewise, according to Herbst et al. (2000), the vibrational energies of \ensuremath {{\rm HNC}} and \ensuremath {{\rm HCN}} formed in the neutral reactions ${\rm C} + {\rm NH_2} \rightarrow \ensuremath{{\rm HNC}} + {\rm H}$, and ${\rm N} + {\rm CH_2} \rightarrow \ensuremath{{\rm HCN}} + {\rm H}$ are sufficient to overcome the isomerization barrier, and roughly equal amounts of both isomers are produced.

The theoretical results presented in Talbi & Herbst (1998) and Herbst et al. (2000), and more recently by Hickman et al. (2005) and Ishii et al. (2006), show quite convincingly that the \ensuremath {{\rm HNC}}/ \ensuremath {{\rm HCN}} abundance ratio resulting from the chemical reactions known to be effective in dark clouds cannot differ significantly from unity. There are two possible explanations for the discrepancy between the observations and theoretical predictions. Either so far unknown processes favour \ensuremath {{\rm HNC}} in the cost of \ensuremath {{\rm HCN}}, or then the derived \ensuremath {{\rm HN^{13}C}}/ \ensuremath{{\rm H^{13}CN}} column density ratio is affected by different degrees of 13C fractionation in \ensuremath {{\rm HNC}} and \ensuremath {{\rm HCN}}. In the latter case 13C would about three times more enhanced in \ensuremath {{\rm HNC}} than in \ensuremath {{\rm HCN}} in Cha-MMS1.

The \ensuremath {{\rm HN^{13}C}}/ \ensuremath{{\rm H^{15}NC}} column density ratio 7.4 together with the above quoted range of 14N/15N ratios imply that the \ensuremath {{\rm HNC}}/ \ensuremath {{\rm HN^{13}C}} ratio in Cha-MMS1 is in the range 30-40. As the typical 12C/13C isotopic ratio observed in nearby clouds is 60 (Lucas & Liszt 1998) the result suggests considerable 13C enhancement in \ensuremath {{\rm HNC}}. An \ensuremath{{\rm HC^{15}N}} column density determination would readily show, according to the results of Terzieva & Herbst (2000), whether also \ensuremath {{\rm HCN}} has a similar degree of 13C fractionation.

In the only survey so far including all these isotopologues, Irvine & Schloerb (1984) found very similar \ensuremath {{\rm HN^{13}C}}/ \ensuremath{{\rm H^{15}NC}} and \ensuremath{{\rm H^{13}CN}}/ \ensuremath{{\rm HC^{15}N}} intensity ratios ($\sim $1.5) towards TMC-1. This observation suggests similar 13C fractionation in \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}}. On the other hand, this result might not be universally valid. Ikeda et al. 2002 determined substantially different \ensuremath{{\rm H^{13}CN}}/ \ensuremath{{\rm HC^{15}N}} column density ratios in the starless dark cloud cores L1521E, TMC-1, and L1498, listed in order of increasing \ensuremath{{\rm H^{13}CN}}/ \ensuremath{{\rm HC^{15}N}} ratio (2.6, 6.7, and 11). The \ensuremath {{\rm HN^{13}C}}/ \ensuremath{{\rm H^{13}CN}} column density ratios in these sources increase in the same order (0.54, 1.7, and 4.5; Hirota et al. 1998). The values obtained towards TMC-1 are consistent with those of Irvine & Schloerb (1984), but the results towards L1521E and L1498 suggest (if the possibility of 15N fractionation is neglected) that the 13C fractionation in \ensuremath {{\rm HCN}} changes from source to source, and perhaps the variation in the case of \ensuremath {{\rm HNC}} is even larger. L1521E is suggested to represent youthful chemistry with little molecular depletion (Hirota et al. 2002; Tafalla & Santiago 2004) whereas the position observed in L1498 is known to be heavily depleted (Willacy et al. 1998). Thus, it appears that the freeze-out of molecules favour the 13C substituted isotopologues of \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}}. It seems warranted to further examine the possibility of their differential 13C fractionation, although no theoretical arguments are in favour of different gas-phase production and destruction rates of \ensuremath {{\rm HN^{13}C}} and \ensuremath{{\rm H^{13}CN}}, or their different freezing rates. As no empirical data on this matter is currently available, the best way to proceed is probably to perform similar observations as those done by Irvine & Schloerb (1984) towards a sample of cores with different degree of depletion.

4.4 Scattering envelope

The low \ensuremath {{\rm HCN}}/ \ensuremath{{\rm H^{13}CN}} and \ensuremath {{\rm HNC}}/ \ensuremath {{\rm HN^{13}C}} intensity ratios and the line shapes in the \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} spectra observed towards the centre of the core require some attention. The intensity ratios suggest very large optical thicknesses for the \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} lines, and consequently, deep self-absorption features and "anomalous'' hyperfine component ratios are to be expected. A clear absorption feature can be seen in the simulated \ensuremath {{\rm HCN}} spectrum shown in Fig. 7, although a very low 12C/13C isotopic ratio is used in order to avoid inverted hyperfine ratios, which would disagree with the observations. The observed \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} spectra are, however, free from absorption dips, and they look more like weakened "normal'' lines.

This weakening can probably be attributed to less dense gas surrounding the core. The \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} molecules in the envelope, where the density lies below the critical density of their J=1-0lines, function as scatterers of line emission arising from the core, i.e. they re-emit absorbed photons (into random directions) before colliding with other particles. The effect is less marked for 13C isotopologues, and the resulting \ensuremath {{\rm HCN}}/ \ensuremath{{\rm H^{13}CN}} and \ensuremath {{\rm HNC}}/ \ensuremath {{\rm HN^{13}C}} intensity ratios decrease, without \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} loosing their characteristic line patterns. This model was used by Cernicharo et al. (1984) to explain \ensuremath {{\rm HCN}} hyperfine anomalies, particularly the intesity of the intrinsically weakest component F=0-1 observed towards dark clouds. The \ensuremath {{\rm HCN}} spectrum shown in Fig. 5 conforms with this explanation: the F=0-1component on the left is stronger than F=1-1 component on the right, although the optically thin LTE ratio is 1:3.

The scattering of \ensuremath {{\rm HCN}} emission from the core has probably lead to underestimation of the excitation temperature and column densities derived using the LTE method in Sect. 3.2. For example, an increase of 1 K in  $T_{\rm ex}$, would mean 30% larger column densities, and 3 K would correspond to a 100% increase in the column density. Quantitative estimates of the scattering effects are, however, beyond the scope of this paper, and here we merely note this additional source of error.

5 Summary

A mapping in the J=1-0 rotational lines of \ensuremath {{\rm HCN}}, \ensuremath {{\rm HNC}} and \ensuremath {{\rm HN^{13}C}} shows a compact structure in the direction of the protostellar core Cha-MMS1. The \ensuremath {{\rm HN^{13}C}} map agrees best with the previous 1.3 mm dust continuum map of Reipurth et al. (1996), and the \ensuremath {{\rm HC_3N}} map of Kontinen et al. (2000).

The column density determinations of \ensuremath{{\rm H^{13}CN}}, \ensuremath {{\rm HN^{13}C}}, \ensuremath{{\rm H^{15}NC}} and \ensuremath {{\rm NH_3}}, and the results of the quoted two studies allow us to estimate that the fractional \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} abundances are $7 \times 10^{-10}$ and $2 \times 10^{-9}$, respectively. These are slightly larger than that of \ensuremath{{\rm N_2H^+}}, but clearly smaller than the \ensuremath {{\rm NH_3}} and \ensuremath {{\rm HC_3N}} abundances ( $3
\times 10^{-8}$ and $1 \times 10^{-8}$, respectively). In accordance with recent results in other dense cores we find that \ensuremath {{\rm NH_3}} is about 100 times more abundant than \ensuremath{{\rm N_2H^+}}. The kinetic temperature derived from the \ensuremath {{\rm NH_3}} lines is about 12 K. The derived fractional abundances conform with the previous suggestion of Kontinen et al. (2000), that Cha-MMS1 represents an advanced chemical state, although it has a large \ensuremath {{\rm HC_3N}} abundance, which usually is thought to indicate early times. This situation is likely to have developed under the circumstances of CO depletion (Ruffle et al. 1997).

The following column density ratios are found: $\ensuremath{{\rm HN^{13}C}} /\ensuremath{{\rm H^{13}CN}} ~
\sim 3{-}4$ and $\ensuremath{{\rm HN^{13}C}} /\ensuremath{{\rm H^{15}NC}} ~ \sim 7$. Using the latter ratio, and the assumption based on the modelling results of Terzieva & Herbst (2000), we estimate that the[ \ensuremath {{\rm HNC}}]/[ \ensuremath {{\rm HN^{13}C}}] ratio is between 30 and 40. As there is no \ensuremath{{\rm HC^{15}N}} column density determination available towards this source, the possibility that differential 13C fractionation causes the observed relatively high [ \ensuremath {{\rm HN^{13}C}}]/[ \ensuremath{{\rm H^{13}CN}}] ratio cannot be entirely ruled out.

The \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} lines are not much brighter than those of the 13C substituted species, but still show no absorption dips. These lines are probably weakened by scattering by \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} molecules residing in less dense gas around the core. As suggested by Cernicharo et al. (1984), the same process can cause the so called \ensuremath {{\rm HCN}} hyperfine anomalies in dark clouds, and is probably also responsible for the "too'' bright F=0-1 component of HCN(J=1-0) towards the centre of Cha-MMS1. The evidence for scattering and the relatively low \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} column densities derived towards the core suggest that the \ensuremath {{\rm HCN}} and \ensuremath {{\rm HNC}} distributions are less centrally peaked than those of \ensuremath {{\rm NH_3}} and \ensuremath{{\rm N_2H^+}}, which would be in agreement with the model of Aikawa et al. (2005) for a rapidly collapsing core.

Acknowledgements
We thank Wolf Geppert for useful discussions and the anonymous referee for valuable comments on the manuscript. This project was supported by the Academy of Finland, grant Nos. 173727 and 174854. P.P.T. thanks the ERASMUS programme of the European Union, Olga Koningfonds, Utrecht, and Kapteyn Fonds, Groningen.

References

 

Copyright ESO 2006