Issue |
A&A
Volume 516, June-July 2010
|
|
---|---|---|
Article Number | A105 | |
Number of page(s) | 8 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/200913806 | |
Published online | 22 July 2010 |
The puzzling behavior of HNCO isomers in molecular clouds![[*]](/icons/foot_motif.png)
N. Marcelino1 - S. Brünken2 - J. Cernicharo1 - D. Quan3 - E. Roueff4 - E. Herbst5 - P. Thaddeus6
1 - Centro de Astrobiología (CSIC-INTA), Laboratorio de Astrofísica Molecular,
Ctra de Ajalvir Km 4, 28850, Torrejón de Ardoz, Madrid, Spain
2 -
I. Physikalisches Institut, Universität zu Köln, Zülpicher Str. 77,
50937 Köln, Germany
3 -
Chemical Physics Program, The Ohio State University, Columbus, OH 43210, USA
4 -
Observatoire de Paris-Meudon, LUTH UMR 8102, 5 Place Jules Janssen,
F-92195 Meudon Cedex, France
5 -
Departments of Physics, Astronomy, and Chemistry, The Ohio State University,
Columbus, OH 43210, USA
6 - Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA
Received 4 December 2009 / Accepted 26 March 2010
Abstract
Context. Isocyanic acid (HNCO) has been observed in
different physical environments in the interstellar medium (ISM) and in
external galaxies. HNCO has several metastable isomers with a ground
electronic singlet state: HOCN, HCNO, and HONC. The recent detection of
fulminic acid (HCNO) in prestellar and protostellar cores and cyanic
acid (HOCN) in warm molecular sources (e.g. hot cores) in the Galactic
center proves that these species could also be common constituents of
the ISM.
Aims. To shed some light on the possible formation pathways of
these species, we searched for HCNO in the sources where HOCN has been
previously detected and vice versa. We have also included the low-mass
protostar IRAS 16293-2422, where HNCO is found to be prominent.
Methods. Using the new EMIR receivers at the IRAM 30-m
telescope, we performed deep searches for three rotational transitions
of HOCN and four of HCNO.
Results. We report the detection of HOCN in four sources - three
dense cores and the lukewarm corino L1527 - where HCNO has been
previously observed. HOCN is tentatively detected toward the
protostellar binary IRAS 16293-2422. However, HCNO has been detected
neither in this source nor in the sources of the Galactic center where
HOCN has been previously reported. The derived abundance ratios
HCNO/HOCN are close to unity in quiescent clouds, while they are less
than 0.01-0.1 in warm clouds. We attempt to explain these results by
using both gas-phase and gas-grain chemical models.
Key words: astrochemistry - line: identification - ISM: abundances - ISM: clouds - ISM: molecules
1 Introduction
Isomeric forms of molecular species are of interest in astrochemistry because they can help constrain the chemical processes leading to their formation and depletion, whether they occur in the gas phase or on the surfaces of dust particles. Furthermore, for some isomers, the stablest one is not necessarily the most abundant in molecular clouds, as is the case for HCN and HNC in cold sources. This pair of isomers, for which gas phase chemical models predict similar abundances at low temperature (see Herbst et al. 2000, and references therein), show however a peculiar behavior depending on the kinetic temperature of the source. Whereas the abundance of the metastable isomer HNC is similar to, and even higher than, that of HCN in dark cloud cores, in warm molecular clouds the abundance ratio HNC/HCN is much less than unity (Hirota et al. 1998; Irvine & Schloerb 1984; Churchwell et al. 1984; Schilke et al. 1992). Indeed, the HCN abundance does not change considerably between cold and warm clouds, but the HNC abundance decreases. This is puzzling since both isomers are mainly formed from the dissociative recombination of HCNH+. According to Watson & Walmsley (1982), these observational results require either a production mechanism that produces more HNC than HCN and a destruction process that primarily destroys HNC when the temperature increases, or a more effective isomerization process in clouds with higher temperatures. Regarding the last point, Talbi et al. (1996) showed that the HNC/HCN isomerization reactions with atomic hydrogen indeed become more effective at higher temperatures owing to activation energy barriers in both directions. Thus, the relative abundance of different isomers in molecular clouds is probably not thermodynamic but kinetic in origin.
HNCO has several metastable isomers with a singlet ground electronic state: HOCN, HCNO, and HONC.
These isomers lie at increasingly higher energy than HNCO according to quantum calculations
(24.7, 70.7, and 84.1 kcal mol-1, respectively; see Schuurman et al. 2004).
Since there is not an obvious common precursor for them, their relative abundances in
molecular clouds might provide some clues concerning the respective contributions of
gas-phase and grain processes to the chemistry.
Two isomeric forms of HNCO have recently been detected in space, HCNO (Marcelino et al. 2009) and HOCN (Brünken et al. 2010).
Fulminic acid (HCNO) is observed in three dark cloud cores and toward the low-mass star-forming
region L1527, with an abundance ratio of HNCO/HCNO between 20 and 60. However, HCNO
has not been detected toward the direction of the cyanopolyyne peak of TMC-1 or in the Orion
Hot Core region.
Cyanic acid (HOCN) is detected in several positions in the giant molecular cloud SgrB2
(Brünken et al. 2010,2009a), after its characterization in the laboratory in the cm- and mm-wave bands
and a tentative detection in SgrB2(OH) based on previous spectral line surveys (Brünken et al. 2009a).
The reported HNCO/HOCN abundance ratio in the Galactic center (GC) clouds is between 40-400.
Brünken et al. (2009a) also obtained a tentative detection of HOCN in TMC-1 through the fundamental
transition at 20.9 GHz. The derived column density is
cm-2 and the
HNCO/HOCN abundance ratio is
100.
The chemistry of the CHNO isomers, where the term refers to all four species, seems to be rather complex. While in Marcelino et al. (2009), the observed abundances of HCNO and HNCO were reproduced well with a steady-state, gas-phase model, Brünken et al. (2010,2009a) proposed models with both gas-phase and grain-surface chemistry to explain the isomeric abundances of HOCN and HNCO in the molecular clouds of the GC. In this paper we report observations of HOCN and HCNO toward different sources in order to shed some light on the possible formation pathways of these species in both the gas phase and on dust surfaces. We looked for HCNO in several positions in SgrB2 associated with hot or warm environments, where HOCN had previously been detected (see Brünken et al. 2010), while HOCN was searched for in the cold to lukewarm dense cores observed in Marcelino et al. (2009). We have also included in our search the well-studied low-mass protostar IRAS 16293-2422, where HNCO has been found to be prominent (Bisschop et al. 2008), and a rich organic chemistry has been widely observed as a result of mantle grain desorption and, possibly, warm gas chemistry (see, e.g., Garrod et al. 2008; Kuan et al. 2004; Garrod & Herbst 2006; Cazaux et al. 2003, and references therein). The observed positions are shown in Table 1.
Table 1: Observed sources.
2 Observations
The observations were performed with the IRAM 30-m telescope (Granada, Spain) in April and May 2009.
We used the new EMIR receivers at 3, 2, and 1 mm in combinations E090/E150 and E090/E230
(3 and 2 mm, 3 and 1 mm, respectively). All receivers were tuned in single-side band mode
with image rejections of 13 dB, and with both polarizations at the same frequency.
System temperatures were between 70-130 K at 3 and 2 mm (with the highest value reflecting
worse weather conditions) and 180 K at 1 mm,
except for the HCNO (
505-404) transition at 114 GHz, where the temperature was higher
owing to the poorer atmospheric transmission (
K).
Intensity calibration was performed using two absorbers at different temperatures.
The atmospheric opacity at 225 GHz, obtained from the measurement of the sky emissivity
and the use of the ATM code (Cernicharo 1985), was typically
0.1 at 1 mm
(water vapor column 1-3 mm) during night time, and between 0.1-0.2 in the afternoon
(water vapor column
5 mm).
![]() |
Figure 1: Line profiles for the observed transitions of HOCN (black solid lines) and simulation from LVG calculations (red solid lines). Upper limits of the intensity scale for the central position of L1527 ( second row of panels) have been reduced by the factor 0.5. For IRAS 16293-2422 ( lowest panels), vertical dashed lines show the observed velocity components (see text). |
Open with DEXTER |
Table 2: HOCN observed line parameters.
Observing modes and spectral resolutions were selected depending on the sources.
For the cold cores, we observed in frequency switching mode, since the lines are
sufficiently narrow, while we used the wobbler switching mode for IRAS 16293-2422 and SgrB2.
The off positions were checked to be free of HCNO and HOCN emission.
The autocorrelator VESPA was used as a backend, with 40 kHz (0.08-0.14 km s-1),
80 kHz (
0.11-0.28 km s-1), and 320 kHz (
0.45-1 km s-1) of
resolution for the dense cores, IRAS 16293-2422, and SgrB2, respectively.
Pointing and focus were checked regularly (every 1.5 and 3 h, respectively)
on strong and nearby sources.
At the observed frequencies, the beamwidth of the antenna is in the range
29''-17''and the main beam efficiency is
0.78-0.64.
All spectra have been calibrated in main beam temperature scale.
The observed line profiles for HOCN are shown in Fig. 1 and
the derived line parameters, obtained from Gaussian fits using the GILDAS
package,
are shown in Table 2.
A total of 50 h were devoted to completion of the whole project. Three and four molecular
transitions arising from HOCN and HCNO, respectively, were observed in different sources.
Owing to the expected low intensity (few tens of mK) of the HOCN lines in dark clouds, we performed
deep integrations on them.
The lowest energy observed transition - the
404-303 (
MHz,
K) line -
was searched for in all dense cores shown in Table 1 (including a second position for L1527)
for 2.5 to 5 h of on-source observing time. In those cores where this transition was detected - B1-b,
L1527 (0,0), L1544, and L183 - we also searched for the
505-404 (
MHz,
K)
line for 2 to 4 h per position. We observed the
707-606 (
MHz,
K) line
simultaneously with the former or the latter transition, resulting in a total observing time per source
of 5-8 h.
Toward IRAS 16293-2422, we used 2 h of observing time each for the HOCN
404-303 and
505-404 lines, and 4 h for the
707-606 transition.
Four HCNO lines - namely the 4-3 (
MHz,
K), 5-4 (
MHz,
K), 6-5 (
MHz,
K), and 9-8 (
MHz,
K) transitions - were searched for in several positions toward SgrB2
(see Table 1), using
1 h per position and transition.
To compare our intensities with the previous HOCN observations (Brünken et al. 2010), we observed
the
404-303 and
707-606 lines in all SgrB2 positions. For that purpose, we used 40 min per
source and transition. HCNO (4-3) and (9-8) were observed in IRAS 16293-2422, but they were not detected
after 1 h of on-source observing time.
3 Results
Cyanic acid (HOCN) has been detected towards three cold dense cores and the center of the low-mass
``lukewarm'' corino L1527 (see Fig. 1).
Three rotational transitions were detected toward B1-b, which shows the strongest intensities,
while in the cold cores L183 and L1544 two and one transitions were detected, respectively.
The HOCN (
404-303) transition is marginally detected in the central position of L1527, while
no lines were detected toward L1527 (20'', -20'') with an rms of 9 mK, where one transition of
HCNO has been previously observed (Marcelino et al. 2009).
Although in some cases we only have one line to support the identification of HOCN in dark clouds,
we are fully confident of the assignment as the density of lines in these sources is very low
(10 lines/GHz), and the agreement in frequency between laboratory and astronomical observations
is better than 20 kHz; i.e., near laboratory accuracy.
Despite the non-detection of HCNO towards the cyanopolyyne peak of TMC-1, we have searched for HOCN
to check the behavior of both isomers, but the
404-303 transition was not detected with an
rms of 6 mK.
The observed line intensities are similar for HOCN and HCNO towards the dense cores observed in
Marcelino et al. (2009), except for L1527, where the intensity of the
404-303 transition of HOCN is
half that of the HCNO J = 4-3 one (see Fig. 1 in Marcelino et al. 2009).
Since the dipole moments of cyanic and fulminic acid are similar (HOCN:
D,
D;
HCNO:
D), i.e., the a-type line strengths of the observed transitions in the bent cyanic
acid are very close to the observed line strengths of the linear fulminic acid, similar line intensities
should reflect similar abundances (see below).
Three rotational transitions of cyanic acid (HOCN) were detected toward the low-mass, star-forming region
IRAS 16293-2422. On the other hand, HCNO was not detected, with an observed rms of 10 mK for the 4-3 line.
Although each HOCN individual spectrum has a lower S/N ratio than in the other sources, the three observed
transitions show the same double-peaked line profile and consistent
(lowest panels in
Fig. 1), supporting a positive detection.
IRAS 16293-2422 is known to be a proto-binary system whose components, A and B, separated by 5''
(Wootten 1989), show a complex velocity pattern (see, e.g., Kuan et al. 2004). HNCO has been observed in
IRAS 16293-2422 with emission at a systemic velocity of 2.5 km s-1 arising from both components
but being more intense in component A, which also shows a strong velocity feature at
km s-1 (see Bisschop et al. 2008).
It is therefore likely that most of the HOCN emission also comes from the A component.
The double-peaked profile has also been interpreted as arising from a disk-like structure
(see Stark et al. 2004; Huang et al. 2005, and references therein).
Table 3: Derived column densities for HOCN and HCNO, and abundance ratios.
Four different transitions arising from fulminic acid (HCNO) were searched for in the SgrB2 molecular complex. The observed positions - SgrB2M, SgrB2M (20'', 100''), SgrB2N, and SgrB2S - were selected on the basis of previous detections of HOCN (see Brünken et al. 2010). However, no HCNO emission was found at any position, with rms values between 5-14 mK. Although a line was observed near 91.751 GHz, it has been recently identified as a transition arising from HSCN (see Halfen et al. 2009; Brünken et al. 2009b).
The three detected transitions towards B1-b and IRAS 16293-2422 arising from HOCN allow us to derive
rotational temperatures of
K and
K for B1-b and IRAS 16293-2422,
respectively; in the latter, a source size of 10'' has been assumed (van Dishoeck et al. 1995).
With these temperatures, we obtained the column densities given in Table 3.
For the other cold cores, those with HOCN detections and non-detections, we adopted the temperatures
obtained for HNCO (see Marcelino et al. 2009) in order to compute column densities and upper limits to the
column density (for the offset position of L1527 and TMC-1).
The obtained upper limits to the column density and to the HNCO/HOCN abundance ratio in TMC-1 are
similar to those observed by Brünken et al. (2009a).
For IRAS 16293-2422 and the GC molecular clouds we calculated upper limits to the column density of HCNO.
In the former, we used the HOCN rotational temperature derived above, while for the latter we assumed
the parameters (source sizes and temperatures) used by Brünken et al. (2010) for HOCN.
Column densities, column density ratios, and fractional abundances for the assorted isomers are shown in
Table 3 for cold cores and GC sources based on this work and prior studies. The fractional
abundances for the GC sources are particularly uncertain mainly because of uncertain H2 column densities and
densities. For Sgr B2N and Sgr B2M, we distinguished between the actual hot core and the extended gas
surrounding the hot core.
Our salient results are that the HOCN and HCNO column densities are similar in all the cold dense cores,
while the warmer molecular clouds in SgrB2 show a very low abundance ratio for HCNO/HOCN.
To better constrain the observed column densities, we also
performed LVG calculations. Since collisional rates for HOCN are not known,
we adopted those of HNCO, as we did for HCNO (see Marcelino et al. 2009).
We assumed a kinetic temperature of 10 K for all dark clouds. Although
the derived rotational temperature for B1-b is lower, the observed intensities
are better reproduced using
K.
For both positions in L1527, we adopted a kinetic temperature of 15 K
(see Marcelino et al. 2009).
In the cores where only one transition was detected and also to obtain
the upper limits for non detections, we used the volume density obtained in Marcelino et al. (2009).
For IRAS 16293-2422, we used the temperature derived above from the rotational diagrams
to calculate both the HOCN column density and the upper limit to that of HCNO.
The derived column densities and ratios, together with the volume densities,
are also shown in Table 3.
The expected intensities obtained with the LVG aproximation for non-detected
lines in L1544, L183, and L1527 (0, 0) are consistent with the observed rms noise
in the spectra (see Fig. 1).
The uncertainties in the LVG calculations are a factor 2 in density and
2 K in
kinetic temperature. The corresponding errors in column density are dominated by
the noise in the observations.
4 Discussion
The results shown in Table 3 indicate a change in the
HCNO/HOCN abundance ratio between cold clouds (quiescent or showing some
indications of low-mass star formation activity) and sources in the GC,
which are assumed to be warmer. In the latter, HCNO is not detected
and the HCNO/HOCN abundance ratio is <
,
while in
cold clouds the ratio is close to unity (0.7-1.7). Nevertheless, the
derived HNCO/HOCN abundance ratio is between 20-320 in all clouds.
It therefore seems that the nature of the cloud does not influence the relative abundance
of HOCN strongly with respect to the most abundant isomer HNCO. The situation is
completely different for fulminic acid, HCNO. The abundance of HNCO with respect
to fulminic acid was found by Marcelino et al. (2009) to range between 20 and 90 in
cold clouds and to be much larger in warm clouds. This result is confirmed by
the upper limits we obtain for this isomer towards the clouds in the GC.
We used two models to explain these results. The first is an extended pseudo-time-dependent approach based on the gas-phase network used in the steady-state results reported in Marcelino et al. (2009), while the second is a time-dependent, gas-grain approach (Hassel et al. 2008). In the pseudo-time-dependent approach, the physical conditions are homogeneous and time-independent while the chemical process is occurring.
4.1 Gas-phase model
Little attention has been paid to gas-phase processes involving the elements H, N, C, and O in the context of interstellar chemistry. Iglesias (1977) had a pioneer paper where chemical processes were proposed for NCO and HNCO. The major synthetic reaction for NCO is the reaction between CN and O2, which is sometimes referred to as a prototypical radical-radical reaction (Smith 1995) and has been studied from 13 K to 3720 K. The NCO radical is a first step towards isocyanic acid and several of its isomers, which form by protonation, hydrogenation, and dissociative recombination. In our preceding paper (Marcelino et al. 2009), we considered exothermic neutral-neutral, ion-molecule, and dissociative recombination reactions to study the chemistry of HNCO, HOCN, and HCNO. Here, we reinvestigated the chemical processes and introduced the additional isomer isofulminic acid (HONC), which is the most energetic stable isomer of isocyanic acid, at an energy calculated to be 84.1 kcal/mol above that of HNCO (Schuurman et al. 2004), and whose rotational spectrum has recently been obtained in the laboratory (Mladenovic et al. 2009). To preserve the consistency of the chemical network, we considered new ions arising from the protonation of NCO and CNO, which may be progenitors of the four different CHNO isomers: i.e., HNCO+, HOCN+, HCNO+, and HONC+. The protonated CHNO ions have been theoretically studied by Ijjaali et al. (2001), and the present chemical network includes their 5 most stable structures S1 to S5 (H2NCO+, HNCOH+, HCNOH+, H2OCN+, and H2CNO+). These authors also determined transition states linking the stable protonated ions, among which two transition states, TS1-2 and TS3-5, connect stable isomers with transition energies of 95.2 and 172.5 kcal/mol, respectively. These values are obviously too high to play a role even in the warm environments of Sg B2.
Because little experimental information on the ion-molecule chemistry is available, we used the most
conservative hypotheses by including only known exothermic channels and assuming that the structures of the
reactants and products do not change dramatically in a chemical reaction. For example, in dissociative
recombination reactions, we only allow those product channels in which the remaining heavy species has
the same skeletal structure of the reactant ion upon separation of an H or H2; viz.,
![]() |
(1) | ||
![]() |
(2) | ||
![]() |
(3) | ||
![]() |
(4) | ||
![]() |
(5) | ||
![]() |
(6) | ||
![]() |
(7) | ||
![]() |
(8) | ||
![]() |
(9) |
When several product channels arise, we used the same value of the branching ratios in the absence of any theoretical or experimental argument.
Uniquely among the CHNO isomers, fulminic acid can possibly be formed by an efficient neutral-neutral reaction (Marcelino et al. 2009):
![]() |
(10) |
The reaction between CH2 and NO has been studied both theoretically and experimentally by a number of investigators, without a unique result concerning the product branching fractions. Some detailed studies show that HCNO is a dominant product (Roggenbuck & Temps 1998; Eschenko et al. 2002; Fikri et al. 2001), while others do not (Zhang et al. 2004). In our standard model, based on the previous work of Marcelino et al. (2009), we consider HCNO to be a dominant product, while we do not in an alternative model. In general, the calculated abundances for HCNO in the alternative model are reduced by a wide range of factors from just a few to two orders of magnitude.
We display in Tables 4 (standard model) and 5 (alternative model)
the fractional abundances relative to molecular hydrogen of the four CHNO isomers, as well as three isomeric
ratios for a variety of different densities at two kinetic temperatures (10 K and 50 K). These temperatures
represent, respectively, those of a cold dark cloud and the envelope of a hot core. Results are listed both
for so-called early time, a time of 0.1 Myr, at which the HNCO abundance peaks, and at steady-state,
which are reached around 10 Myr.
We also ran the models at 200 K (hot core case), but the results do not differ much from those
obtained at 50 K at the densities used, which do not range as high as our estimated central hot core densities
in Table 3. We used the early-time results in comparison with observation, and the former are
higher than the steady-state values by up to a factor of
5-10.
Table 4: Fractional abundances and abundance ratios from the standard gas-phase modela.
Table 5: Fractional abundances and abundance ratios from the alternative gas-phase modela.
In these models, the cosmic ionization rate is taken as
s-1, and elemental abundances
corresponding to the ``low metal'' case have been used. Some of the results for the lowest densities
are unusual in that they are characteristic of the high ionization phase (HIP) (Le Bourlot et al. 1993), with
vanishingly low values for saturated molecules, so the displayed ratios are not significant. The other
results are characteristic of the better known, more standard low ionization phase.
The fractional abundances calculated for HNCO and HOCN for the cold sources tend to be too high by an
order of magnitude or so, while for HCNO this is true only for the standard model. The results of the
alternative model, on the other hand, underestimate the observed HCNO fractional abundances in the cold
sources. In the GC sources, the H2 column densities and densities are more uncertain, but it appears
that the 50 K theoretical results at a density of
cm-3 for HNCO and HOCN reproduce the
hot core envelopes well, but the higher-density results are significantly lower than the observed values for
the hot cores themselves.
If we focus on the abundance ratios among the HNCO isomers, we note that the observed HNCO/HOCN ratio in the
cold clouds (50-100) is fit closely at the appropriate densities by the early-time theoretical results,
whereas the HCNO/HOCN ratio (
1) is fit in the standard case but is underproduced by the alternative
case. We conclude that the standard model is needed and that HCNO must be produced in the gas-phase model by
the reaction between CH2 and NO.
In the GC clouds and IRAS 16293-2422, the HNCO/HOCN ratio is also reproduced reasonably well. On the other hand, the low HCNO/HOCN ratios are not reproduced by the standard model, which tends to be 1-2 orders of magnitude higher. The alternative model, on the other hand, does somewhat better. This conflict between the standard and alternative models, depending upon whether the clouds are cold or warm, suggests that more detailed gas-grain models should be considered, which for warm sources include the actual warm-up stage.
4.2 Gas-grain model
To simulate the different regions where the CHNO isomers are detected, we utilized the Ohio State
gas-grain network (Garrod et al. 2008; Garrod & Herbst 2006) with a variety of physical conditions. In this code, both
gas-phase and grain-surface chemistry are included, and solutions are found as functions of time.
Even more so than in the pure gas-phase modeling discussed above, many of the rate coefficients
involving these isomers are estimated. In addition to a model for a cold core at 10 K and an
H2 density of 104 cm-3, we computed three different warm-up models, which undergo
a cold phase (10 K) that lasts for 105 yr followed by a warm-up that lasts
yr
where the temperature increases quadratically with time (Viti et al. 2004) to an asymptotic value depending
on the physical object being simulated. The asymptotic temperatures for a lukewarm corino (e.g. L1527),
a warm envelope surrounding a hot core, and the hot core itself are 30 K, 50 K, and 200 K, respectively
(Hassel et al. 2008; Garrod & Herbst 2006).
The H2 densities in the three models are 106, 105, and 106 cm-3, respectively.
As in the gas-phase calculations, two models have been run, depending upon the products of the CH2 + NO
reaction. We once again use the terms ``standard'' and ``alternative'' for the two models.
In general, the results of the gas-grain model calculations do not totally solve the theoretical problem of accounting for the comparable abundances of HCNO and HOCN in cold regions and the much lower abundance of HCNO with respect to HOCN in warmer ones. Since the gas-grain results will be the subject of a future paper (Quan, Osamura, & Herbst, in preparation), we concentrate here on one type of warm-up model, the warm envelope surrounding a hot core.
Figure 2 shows calculated gas-phase fractional abundances for the four CHNO isomers
in the gas phase as functions of time for the standard and alternative warm-up envelope models.
At times up to 105 yr, the calculated abundances are similar to those of the cold core model.
Here, the HCNO isomer is predicted to be significantly more abundant than HOCN, in disagreement with
observations. During the subsequent warm-up period, the abundances of HNCO and HOCN increase markedly,
whereas those of HCNO and HONC do not. Still, at the start of the warm-up era, there is little HOCN,
so that despite its subsequent increase, it achieves a higher abundance than HCNO only for a brief period.
After the asymptotic temperature of 50 K is reached, all abundances decrease sharply or gradually, and the
abundance of HCNO remains higher than that of HOCN, even with the alternative model. During the brief
period near the time of peak abundances, which occurs near the end of the warm-up period, the fractional
abundance of HNCO equals (alternative model) or exceeds (standard model) 10-8, while those of HOCN and
HCNO lie in the vicinity of 10-11 and 10-11-10-12, respectively. The peak calculated
abundance of HNCO exceeds its observed value, while that of HOCN reasonably agrees with observation
(see Fig. 2). The predicted minimum HCNO/HOCN ratio is 0.5 for the standard case
and
0.1 for the alternative one; the latter is in closer agreement with observations.
The observed higher upper limit for the HCNO/HOCN ratio in IRAS 16293-2422, which is similar in
temperature to the environs of a hot core, is reproduced slightly before or slightly after the
peak abundances in our envelope model, after which the HCNO abundance becomes larger than for HOCN.
In the standard model, an abundance ratio of unity occurs at 2 and
yr, while
in the alternative model, such an abundance ratio occurs at slightly later times.
![]() |
Figure 2: Fractional abundances (X) for the gaseous CHNO isomers with respect to molecular hydrogen computed with the gas-grain model for the case of a warm envelope surrounding a hot core. The observed values are in the vicinity of those observed for the extended environments considered here. |
Open with DEXTER |
5 Conclusions
In this paper we have presented observations of the HNCO metastable isomers, HOCN and HCNO, in several sources including quiescent prestellar cores, low-mass protostellar objects, and hot cores and their environs, in particular the GC sources. While HOCN was observed in all sources (except TMC-1) with a similar HNCO/HOCN abundance ratio, HCNO was not detected towards the low-mass protostar IRAS 16293-2422 or the warm GC clouds. Because HCNO is higher in energy than HOCN, this behavior of the CHNO isomers is similar to that of HCN and HNC, where the latter seems to be depleted when the kinetic temperature of the cloud increases.
In an attempt to explain the fractional and relative abundances of the CHNO isomers, we used both a gas-phase, pseudo-time-dependent model and a gas-grain model that specifically considers a warm-up phase in between cold and warm eras. The results of both models are mixed and, for HCNO, depend to some extent on whether or not the gas-phase reaction between CH2 and NO produces HCNO.
In summary, much of the chemistry is still poorly understood, and it is not clear in particular why HCNO is so underabundant with respect to HOCN in the warm sources but not in the cold ones. It is probable that we are missing an important destruction mechanism for HCNO. Although we considered several destruction channels (in particular reactions with atomic C, O, and H), we are not able to reproduce the low observed abundances at high temperatures. One possibility would be a mechanism that converts HCNO into HOCN, in analogy to the conversion of HNC into HCN in warm sources.
Even so, the chemistry of the CHNO isomers that occurs on the dust grains may play a key role in the relative abundances of HNCO/HCNO/HOCN. Laboratory experiments on the surface formation of these isomers via hydrogenation of NCO and CNO could add quantitative and qualitative information to our models. One basis of both our gas-phase and gas-grain chemistry is the radical NCO. Detection or non-detection of this species in space would also help to constrain our chemical models.
AcknowledgementsWe thank the referee for pointing out additional information on the reaction between CH2 and NO and for other useful comments. We thank M. Gerin for useful discussions concerning the chemistry. This work has been supported by Spanish Ministerio de Ciencia e Innovación through grants AYA2006-14876 and ESP2007-65812-CO2-01, and by DGU of the Madrid community government under IV-PRICIT project S-0505/ESP-0237 (ASTROCAM). E.H. acknowledges the support of the National Science Foundation for his astrochemistry program through grant AST-0702876 and for his program on chemical kinetics through the Center for the Chemistry of the Universe. E.H. also acknowledges the support of NASA for studies in the evolution of pre-planetary matter.
References
- Bisschop, S. E., Jørgensen, J. K., Bourke, T. L., Bottinelli, S., & van Dishoeck, E. F. 2008, A&A, 488, 959 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Belloche, A., Menten, K. M., Comito, C., et al. 1990, A&A, 482, 179 [Google Scholar]
- Bottinelli, S., Ceccarelli, C., Neri, R., et al. 2004, ApJ, 617, 69 [Google Scholar]
- Brünken, S., Gottlieb, C. A., McCarthy, M. C., & Thaddeus, P. 2009a, ApJ, 697, 880 [NASA ADS] [CrossRef] [Google Scholar]
- Brünken, S., Yu, Z., Gottlieb, C. A., McCarthy, M. C., & Thaddeus, P. 2009b, ApJ, 706, 1588 [NASA ADS] [CrossRef] [Google Scholar]
- Brünken, S., Belloche, A., Martín, S., Verheyen, L., & Menten, K. M. 2010, A&A, 516, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cazaux, S., Tielens, A. G. G. M., Ceccarelli, C., et al. 2003, ApJ, 593, 51 [Google Scholar]
- Cernicharo, J. 1985, IRAM report No. 52 (Granada: IRAM) [Google Scholar]
- Cernicharo J., & Guélin M. 1987, A&A, 176, 299 [NASA ADS] [Google Scholar]
- Churchwell, E., Nash, A. G., & Walmsley, C. M. 1984, ApJ, 287, 681 [NASA ADS] [CrossRef] [Google Scholar]
- Eschenko, G., Koecher, T., Kerst, C., & Temps, F. 2002, Chem. Phys. Lett., 356, 181 [NASA ADS] [CrossRef] [Google Scholar]
- Fikri, M., Meyer, S., Roggenbuck, J., & Temps, F. 2001, Faraday Disc., 119, 223 [Google Scholar]
- Garrod, R. T., & Herbst, E. 2006, A&A, 457, 927 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Garrod, R. T., Widicus Weaver, S. L, & Herbst, E. 2008, ApJ, 682, 283 [NASA ADS] [CrossRef] [Google Scholar]
- Halfen, D. T., Ziurys, L. M., Brünken, S., et al. 2009, ApJ, 702, L124 [NASA ADS] [CrossRef] [Google Scholar]
- Hassel, G. E., Herbst, E., & Garrod, R. T. 2008, ApJ, 681, 1385 [NASA ADS] [CrossRef] [Google Scholar]
- Herbst, E., Terzieva, R., & Talbi, D. 2000, MNRAS, 311, 869 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hirano, N., Kamazaki, T., Mikami, H., Ohashi, N., & Umemoto, T. 1999, in Star Formation 1999, ed. T. Nakamoto, NRO, 181 [Google Scholar]
- Hirota, T., Yamamoto, S., Mikami, H., & Ohishi, M. 1998, ApJ, 503, 717 [NASA ADS] [CrossRef] [Google Scholar]
- Huang, H.-C., Kuan, Y.-J., Charnley, S. B., et al. 2005, Adv. Space Res., 36, 146 [NASA ADS] [CrossRef] [Google Scholar]
- Iglesias E. 1977, ApJ, 218, 697 [NASA ADS] [CrossRef] [Google Scholar]
- Ijjaali, F., Alcami, M., Mó, O., & Yáñez, M. 2001, Mol. Phys., 99, 1129 [NASA ADS] [CrossRef] [Google Scholar]
- Irvine, W. M., & Schloerb, F. P. 1984, ApJ, 282, 516 [Google Scholar]
- Jørgensen, J. K., Schöier, F. L., & van Dishoeck, E. F. 2002, A&A, 389, 908 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kuan, Y.-J., Huang, H.-C., Charnley, S. B., et al. 2004, ApJ, 616, 27 [Google Scholar]
- Le Bourlot, J., Pineau des Forets, G., Roueff, E., & Schilke, P. 1993, ApJ, 416, L87 [NASA ADS] [CrossRef] [Google Scholar]
- Lis, D. C., & Goldsmith, P. F. 1990, ApJ, 356, 195 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, S.-Y., & Snyder, L. E. 1999, ApJ, 523, 683 [NASA ADS] [CrossRef] [Google Scholar]
- Marcelino, N., Cernicharo, J., Tercero, B., & Roueff, E. 2009, ApJ, 690, 27 [Google Scholar]
- Mladenovic M., Lewerenz M., McCarthy M. C., & Thaddeus P. 2009, J. Chem. Phys., 131, 174308 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Reid, M. J. 1993, ARA&A, 31, 345 [NASA ADS] [CrossRef] [Google Scholar]
- Roggenbuck, J., & Temps, F. 1998, Chem. Phys. Lett., 285, 422 [NASA ADS] [CrossRef] [Google Scholar]
- Schilke, P., Walmsley, C. M., Pineau des Forêts, G., et al. 1992, A&A, 256, 595 [NASA ADS] [Google Scholar]
- Schuurman M. S., Muir, S. R., Allen W. D., & Shaefer III H. F. 2004, J. Chem. Phys., 120, 11586 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Smith I. W. M. 1995, Int. J. Mass. Spectr. Ion Proc., 149, 231 [NASA ADS] [CrossRef] [Google Scholar]
- Stark, R., Sandell, G., Beck, S. C., et al. 2004, ApJ, 608, 341 [NASA ADS] [CrossRef] [Google Scholar]
- Swade, D. A. 1989, ApJ, 345, 828 [NASA ADS] [CrossRef] [Google Scholar]
- Talbi, D., Ellinger, Y., & Herbst, E. 1996, A&A, 314, 688 [NASA ADS] [Google Scholar]
- van Dishoeck, E. F., Blake, G. A., Jansen, D. J., & Groesbeck, T. D. 1995, ApJ, 447, 760 [NASA ADS] [CrossRef] [Google Scholar]
- Viti, S., Collings, M. P., Dever, J. W., McCoustra, M. R. S., & Williams, D. A. 2004, MNRAS, 354, 1141 [NASA ADS] [CrossRef] [Google Scholar]
- Ward-Thompson, D., Motte, F., & André, P. 1999, MNRAS, 305, 143 [NASA ADS] [CrossRef] [Google Scholar]
- Watson, W. D., & Walmsley, C. M. 1982, in Regions of Recent Star Formation, ed. R. S. Roger, & E. Dewdney (Dordrecht: Reidel), 357 [Google Scholar]
- Wootten, A. 1989, ApJ, 337, 858 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, W., Du, B., & Feng, C. 2004, J. Mol. Struc. (Theochem), 679, 121 [CrossRef] [Google Scholar]
Footnotes
- ... clouds
- This work was based on observations carried out with the IRAM 30-m telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain).
- ...
package
- http://www.iram.fr/IRAMFR/GILDAS
All Tables
Table 1: Observed sources.
Table 2: HOCN observed line parameters.
Table 3: Derived column densities for HOCN and HCNO, and abundance ratios.
Table 4: Fractional abundances and abundance ratios from the standard gas-phase modela.
Table 5: Fractional abundances and abundance ratios from the alternative gas-phase modela.
All Figures
![]() |
Figure 1: Line profiles for the observed transitions of HOCN (black solid lines) and simulation from LVG calculations (red solid lines). Upper limits of the intensity scale for the central position of L1527 ( second row of panels) have been reduced by the factor 0.5. For IRAS 16293-2422 ( lowest panels), vertical dashed lines show the observed velocity components (see text). |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Fractional abundances (X) for the gaseous CHNO isomers with respect to molecular hydrogen computed with the gas-grain model for the case of a warm envelope surrounding a hot core. The observed values are in the vicinity of those observed for the extended environments considered here. |
Open with DEXTER | |
In the text |
Copyright ESO 2010
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.