Issue |
A&A
Volume 516, June-July 2010
|
|
---|---|---|
Article Number | A98 | |
Number of page(s) | 8 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/201013997 | |
Published online | 21 July 2010 |
HNCO enhancement by shocks in the L1157 molecular outflow
N. J. Rodríguez-Fernández1 - M. Tafalla2 - F. Gueth1 - R. Bachiller2
1 - IRAM, 300 rue de la Piscine, 38406 St. Martin d'Heres, France
2 - Observatorio Astronomico Nacional, IGN,
calle Alfonso XII 3, 28014 Madrid, Spain
Received 5 January 2010 / Accepted 23 March 2010
Abstract
Context. The isocyanic acid (HNCO) has an extended
distribution in the centers of the Milky Way and the spiral galaxy
IC342. Based on both the morphology of the emission and the HNCO
abundance with respect to H2, several authors hypothesized that HNCO could be a good tracer of interstellar shocks.
Aims. We test this hypothesis by observing a well-known Galactic source whose chemistry is dominated by shocks.
Methods. We observed several transitions of HNCO towards
L1157-mm and two positions (B1 and B2) in the blue lobe of the
molecular outflow.
Results. The HNCO line profiles exhibit the same characteristics as other well-known shock tracers such as CH3OH, H2CO, SO, or SO2. The three molecules HNCO, SO2,
and OCS are the only ones detected so far whose emission is more
intense in B2 than in B1, making these species valuable probes of
chemical differences along the outflow. The HNCO abundance with respect
to H2 is (0.4-
in B1 and (0.3-
in B2. These abundances are the highest ever measured, and imply an
increment with respect to L1157-mm of up to a factor of 83,
demonstrating that this molecule is a good shock tracer.
Conclusions. we demonstrate that shocks can produce the HNCO
abundance measured in galactic nuclei and even higher values. We
propose that the gas phase abundance of HNCO is produced by both grain
mantle erosion by the shock waves and by neutral-neutral reactions in
gas phase involving CN and O2. The observed anticorrelation
between CN and HNCO fluxes supports this scenario. The observed
similarities between the HNCO emission and sulfur-bearing molecules may
be caused by formation pathways that also involve O2.
Key words: ISM: individual objects: L1157 - ISM: jets and outflows - ISM: molecules - stars: formation
1 Introduction
Interstellar isocyanic acid (HNCO) was first detected towards Sgr B2 (Kuan & Snyder 1996; Snyder & Buhl 1972; Churchwell et al. 1986). Since the first detection in this source, the molecule has been observed in other hot cores around massive (MacDonald et al. 1996; Blake et al. 1987) and low mass protostars (Bisschop et al. 2008; van Dishoeck et al. 1995). It has also been detected in translucent clouds (Turner et al. 1999) and in the dense regions of Galactic molecular clouds (Jackson et al. 1984; Zinchenko et al. 2000), including those at the Galactic center (Dahmen et al. 1997; Lindqvist et al. 1995; Rizzo et al. 2000; Martín et al. 2008; Hüttemeister 1993; Minh et al. 2005). HNCO has also been detected in some extragalactic sources (Martín et al. 2009; Meier & Turner 2005; Nguyen-Q-Rieu et al. 1991). The isotopologue HCNO (fulminic acid) was detected by Marcelino et al. (2009).
The emission in HNCO lines with
K-1 = 0
is clearly extended in the two Galactic nuclei where
it has been mapped: IC342 (Meier & Turner 2005) and the Milky Way (Dahmen et al. 1997; Lindqvist et al. 1995).
In the Galactic center, the spatial distribution of the
HNCO emission differs from that of most other tracers.
In particular, there is a HNCO peak without any associated CO peak at Galactic longitude
(Dahmen et al. 1997).
This very special distribution suggests that HNCO may be an important
tracer of some physical processes that are not well revealed by other
molecules.
It is possible that the molecule is tracing shocks, which are
understood to occur at this Galactic longitude due to the gas dynamics
in the barred potential
of the Milky Way (Rodriguez-Fernandez & Combes 2008).
At the Galactic center, the highest gas phase abundances of
SiO (a well-known shock tracer) are measured in this region
(Rodriguez-Fernandez et al. 2006; Hüttemeister et al. 1998).
Because of the special morphology of the HNCO emission in Sgr B2,
Minh & Irvine (2006) suggested that HNCO is enhanced by shocks.
Regarding IC342, the HNCO emission resembles that of CH3OH. In particular,
HNCO is detected not only in the nuclear ring but also in
the inner spiral arms or dustlanes.
Since the dustlanes are understood to be
the locus of strong shocks in barred galaxies
(see for instance Rodriguez-Fernandez & Combes 2008, and references therein),
Meier & Turner (2005) suggested that HNCO, as CH3OH, could trace large-scale
shocks.
Nevertheless, the hypothesis that HNCO is a good shock tracer on the scale of galaxies still needs to be probed since the HNCO emission has never been studied in well-known Galactic templates of interstellar shocks. To understand more clearly the excitation and the chemistry of this promising molecule, we observed the protostar L1157 and its associated molecular outflow. This outflow presents the morphological signature of shocks (Gueth et al. 1998; Codella et al. 2009) and is frequently used to benchmark numerical models of shocks (Gusdorf et al. 2008a,b). The L1157 outflow is the clearest example of ``chemically active outflow'' (Bachiller et al. 2001) and a template of shock chemistry, since many species have significantly higher abundances than the protostar (Bachiller & Perez Gutierrez 1997). Therefore, it is a source of choice to characterize the emission of a given molecule in a shocked environment (see for instance Arce et al. 2008; Benedettini et al. 2007; Bachiller et al. 2001; Bachiller & Perez Gutierrez 1997; Codella et al. 2009).
![]() |
Figure 1: Integrated CO(2-1) emission from the L1157 outflow (Bachiller et al. 2001). We show the position of L1157-mm (star), B1 (empty square), and B2 (filled square). B1 and B2 are the positions of the shocks as traced, for instance, by SiO(3-2) (see Fig. 6 of Bachiller et al. 2001). The CO emission peaks behind the true positions of the shocks. The different circles represent the half maximum contours of the IRAM 30m beam at the frequencies of the lines discussed in this paper (see Table 1). |
Open with DEXTER |
2 Observations
We observed three lines of the
K-1 = 0 ladder of HNCO
towards the two main shocks in the southern lobe of the L1157 molecular outflow,
the B1 and B2 positions of Bachiller & Perez Gutierrez (1997).
In addition, as reference we observed the continuum source L1157-mm
(the protostar).
The equatorial coordinates of L1157-mm are RA =
,
Dec =
,
(J2000).
The offsets of B1 and B2 with respect to L1157-mm are
(
20'', -60'') and (
35'', -95''), respectively.
The observations were performed with the IRAM 30m telescope in Pico Veleta in July 2007. The line quantum numbers, frequencies, and the telescope parameters are listed in Table 1. As backends, we used the 100 KHz filterbank and the VESPA autocorrelator with a channel resolution of 20-40 MHz, which allows us to study in detail the line wings. The observations were carried out in position-switching mode with the OFF position located 24' NE of L1157-mm. Typical system temperatures and rms noise of the spectra are given in Table 1.
Table 2 indicates
the integrated flux of the different lines computed from spectra in
units (forward and beam efficiencies are given in Table 1).
The uncertainties in the calibration of the IRAM 30m telescope are within
10% at the frequency of the lines discussed here (Mauersberger et al. 1989).
Table 1: Observational parameters.
Table 2: Measured fluxes.
![]() |
Figure 2: Observed HNCO spectra. From bottom to top, the different panels show the HNCO 404-303, 505-404, and 606-505 spectra towards L1157 mm, B1, and B2. |
Open with DEXTER |
![]() |
Figure 3: Comparison of the HNCO(505-404) line profiles (in red) with the spectra of other species published by Bachiller & Perez-Gutierrez (1997), in black lines. The HNCO peak temperature has been scaled to the peak temperature of the other molecule to compare profiles. |
Open with DEXTER |
3 Results
3.1 Line profiles: differences between the observed positions
Figure 2 shows the spectra. The 404-303 and 505-404 lines were detected at all the positions. In addition the 606-505 line was also detected at L 1157 mm. The lines are narrow towards L1157-mm (linewidths of 0.7-0.9 km s-1) but the linewidth increases by more than a factor of 3 towards the B1 and B2 positions. In particular, they have prominent blue-wings towards B1.
Figure 3 compares of the HNCO (505-404) line profile measured towards the three observed positions and the spectra of other molecules observed by Bachiller & Perez Gutierrez (1997). The line profiles are in overall good agreement with the other molecules. In particular, the HNCO profiles are very similar to those of CH3OH, H2CO, SO, and SO2, which are all shock tracers whose abundance has increased by more than a factor of 10 in the shocked gas with respect to the quiescent gas in L1157-mm. The HNCO line wings are narrower than SiO lines, in agreement with the findings of Zinchenko et al. (2000) in Galactic dense cores, but the SiO line profiles are particularly broad in the L1157 molecular outflow, being broader than those of all the other molecules.
As discussed by Bachiller & Perez Gutierrez (1997), most of the molecules exhibit lines of similar intensity in B2 and B1. The exceptions are CN, which is very intense towards B1 but not detected towards B2, and the sulfur-bearing molecules SO2 and OCS, whose lines are far more intense in B2 than in B1. We note that HNCO, whose lines are a factor of 2 more intense in B2 than in B1, is indeed the only non-sulfur-bearing molecule showing this behavior.
3.2 Enhanced HNCO abundances in the shocked gas
The critical density of the HNCO transitions discussed here is high (

3.2.1 LTE analysis
We computed HNCO column densities assuming optically thin emission and LTE excitation. We also computed the excitation temperatures of the 505and 404 populations (T





3.2.2 RADEX analysis
In addition, we studied the HNCO excitation using RADEX (van der Tak et al. 2007), which is a non-LTE excitation and radiative transfer code that decouples the statistical equilibrium and the radiative transfer equations using the escape probability method. The calculations assume lines of a given width and rectangular shape. Therefore, the variation in the optical depth over the profile is not taken into account. This is not a problem if the opacity of the lines is not very high. We used the HNCO-H2 collisional coefficients listed in the LAMDA database (Schöier et al. 2005). The computations use a background temperature of 2.73 K. We modeled the line fluxes (or integrated intensities) as a function of the hydrogen volume density (



![]() |
Figure 4:
RADEX results for the 505-404 to 404-303
flux ratio
(horizontal curves) and the 404-303 line flux
(vertical curves) computed with a
brightness temperature for a source size of 10''
as a function of the hydrogen density, the HNCO column density and the kinetic
temperature (30 K in the upper panel and 250 K in the lower panel). A typical
uncertainty of 10 |
Open with DEXTER |
![]() |
Figure 5:
Same than Fig. 4 but using
|
Open with DEXTER |
We based our analysis on the two lowest lines, which are
detected towards the three observed positions.
Figure 4
shows the model predictions for the 505-404 to 404-303 line ratio and
the 404-303 brightness temperature (for a source size of 10'')
as a function of
and
.
The same diagrams for a source size equal to the
telescope main beam (using
)
are shown in Fig. 5.
The kinetic temperature cannot be constrained since
equally good solutions can be found by
just changing the hydrogen density by 0.5 dex.
For a given source size and kinetic temperature, the H2 density is similar for B1 and B2.
The density varies from 104-104.4 cm-3 for a source size of 10'' and TK=250 K, to
104.6-105.2 cm-3 for sources filling the main beam and TK=30 K.
In L1157-mm,
is lower than that in B1 and B2 by 0.2-0.3 dex.
The total HNCO columns densities are 0.2-1.3, 0.6-2.5, and 1.3-5.0
(in units of 1013 cm-2) for mm, B1, and B2, respectively (Table 3).
The lower limits are similar to those derived using LTE,
while the upper limits for B1 and B2 are somewhat higher because opacity
effects begin to play a role.
Table 3 also shows HNCO abundances computed with the RADEX
estimations of the HNCO column density and the total H2 column density
derived from CO by Bachiller & Perez Gutierrez (1997).
The HNCO abundances are
(0.3-
for L1157-mm, (4.3-
for L1157-B1,
and (25-
for L1157-B2 (Table 3).
Therefore, there is a clear increase in the HNCO abundance in the shocked regions,
by a factor of 6-14 in B1 and by a factor of 34-83 in B2.
This is the first direct evidence that HNCO is indeed a good tracer of shocked gas.
In the range of physical conditions and HNCO column densities derived from the 505-404 and 404-303 lines, the model predictions for the 606-505 to 505-404 ratio is 0.6-0.7. In contrast, the measured ratio is 0.3-0.4 for L1157-mm and lower than 0.3 for B1 and B2. As already mentioned above when discussing the LTE analysis, a multi-transition survey would be needed to fully understand the excitation of the molecule. In addition, HNCO maps with high angular resolution will be needed to determine the positions and sizes of the different HNCO clumps, and in particular to ascertain whether is any HNCO emission outside the IRAM-30 m beam at 2 mm.
3.2.3 Comparison with other Galactic sources
Table 4 shows HNCO abundances with respect to H2 measured for other Galactic sources such
as translucent clouds, photon-dominated regions (PDRs), dense cores, and hot cores.
The highest HNCO abundance is found in hot cores and dense cores
up to
.
The HNCO abundance in L1157-mm is similar to that measured in
these dense and hot cores with the lowest HNCO abundance
and in the only ``hot corino'' where the molecule has been previously
detected, IRAS 16293 (van Dishoeck et al. 1995).
The HNCO abundance in the L1157-B1 shock is
comparable, but can even be higher than that in hot cores and dense cores.
For the L1157-B2 shock, the HNCO abundance is at least 4-11 times higher than that in hot cores and dense cores. The HNCO abundance in the L1157-B2 shock are indeed the highest HNCO abundances
ever measured.
Table 3: Results.
Table 4: HNCO abundance with respect to H2 in different astronomical environments.
4 Discussion
4.1 HNCO formation in molecular clouds and hot cores
The possible formation pathways of HNCO in shocks have never been modeled. In contrast, HNCO has been included in some models of dark clouds and hot core chemistry. In the models by Iglesias (1977), HNCO is produced by the ion-neutral reaction





In the first phase, a dark cloud of
cm-3 undergoes
an isothermal collapse.
The collapse phase is halted once
cm-3, which occurs
after approximately
years (free-fall time).
The typical temperature at this phase is 10-20 K and surface chemistry is very important then.
Surface chemistry models (Tideswell et al. 2010; Garrod et al. 2008) produce
HNCO as a secondary radical formed via the reaction
.
On the grain surface,
HNCO can reach a maximum abundance
of 10-5 (Tideswell et al. 2010) before being
destroyed by new reactions with primary radicals (H, CH3, HCO, NH, ...)
giving complex species as HNCHO, HNCOCHO, or CH3CONH (Garrod et al. 2008).
The cloud-collapse phase is followed by a warm-up phase in which
the temperature increases to 200 K.
The timescale of the warm-up phase can be as short as
yr
in hot cores but can reach 106 yr in their low mass
equivalent (hot corinos),
giving a similar although slightly different chemistry (Garrod et al. 2008).
During the warm-up phase, molecules are evaporated from the grain mantles. With respect to HNCO, the models by Garrod et al. (2008) and Tideswell et al. (2010) indicate that its abundance increases monotonously with time in the gas phase, which implies that HNCO is not directly ejected from the grains. Instead, HNCO is formed by the destruction in the gas phase of the complex molecules formed from HNCO on the grain surface (Tideswell et al. 2010; Garrod et al. 2008). This mechanism can explain the HNCO abundances measured in hot cores of 10-9-10-8. As already mentioned, the contribution of the Turner et al. (1999) and Iglesias (1977) formation pathways is negligible at the moderate temperature of a hot core.
The abundance of HNCO in L1157-mm is similar to that measured in the hot corino of the only other low mass protostar where HNCO has been observed (IRAS16293-2422, see Table 4). Therefore, the HNCO abundance measured in L1157-mm can also be explained in the context of the hot-cores models discussed above. This would imply that L1157-mm can be considered as a hot corino, which is indeed in agreement with the intense emission of water detected by Herschel towards mm (Nisini et al., in prep). However, more observations will be needed to clearly establish the hot corino character of L1157-mm, in particular observations of complex organic molecules.
4.2 HNCO formation in shocks
What is the formation pathway of HNCO in the L1157 outflow shocks? The most likely scenario is a combination of grain surface and gas phase chemistry. In the context of shocks, molecules formed on the grain surfaces can be ejected to the gas phase due to grain sputtering instead of thermal evaporation as in hot cores. Therefore, the chemistry will be sensitive to the grain mantle composition at the time of the shock arrival. No models have yet been developed to study the HNCO abundance in dark clouds including grain surface chemistry. Dark cloud models by Tideswell et al. (2010) only consider gas phase reactions. Therefore, the exact grain mantle composition before the arrival of the shock in L1157-B1 and L1157-B2 is not known. Nevertheless, one can compare with the cloud-collapse phase of hot core models (the main difference is that the density during the collapse phase reaches higher values than in dark clouds). As discussed in the previous section, during the collapse phase the HNCO abundance on the grain surfaces reaches a peak of 10-5, grain sputtering at that time immediately producing a very high HNCO abundance in gas phase. Even the very high HNCO abundance in L1157-B2 (
On the other hand, gas phase chemistry in shocks differs considerably from
that in hot cores because the temperature in the shocked gas can be much
higher than that in a hot core.
As already pointed out by Zinchenko et al. (2000),
in shocks the reaction
can be efficient in spite of the activation barrier.
In addition, the reaction
will be assisted
by the high O2 abundances that are predicted in post-shock gas
(Gusdorf et al. 2008a,b; Bergin et al. 1998).
Using Tideswell et al. (2010) results, it is possible to verify whether
gas phase chemistry alone could explain the HNCO abundance measured in
the L1157 shocks.
In one of their dark cloud models (DC1), they
studied the HNCO formation exclusively via the gas phase neutral-neutral reactions
of Turner et al. (1999) as a function of the rate coefficient for the reaction
.
Their Fig. 1 shows that a rate coefficient k higher than
10-14 cm3 s-1 is needed to explain the observed abundances in L1157
of 10-8-10-7.
However, to produce a high HNCO abundance (a few 10-8) in a short time
(several 103 yr) the rate coefficient must be higher than 10-12 cm3 s-1.
Taking into account the formula
given by Tideswell et al. (2010),
k=10-12 cm3 s-1implies a temperature higher than 1600 K.
In contrast, the minimum temperature would be 428 K if the
activation barrier is only 1160 K (Turner et al. 1999).
Gusdorf et al. (2008b) computed shock models to compare the modelled
SiO and H2 emission to observations of L1157.
The closest agreement between models and observations is found
for a preshock density of 104 cm-3 and a shock speed of 20 km s-1.
In this shock, the gas temperature is higher than 430-1600 K
for a few thousand years, which is consistent with both
the dynamical age of the B1 and B2 shocks (2000-3000 yr, Gueth et al. 1998)
and the chemical time necessary
to produce a HNCO abundance of a few 10-8only with gas phase reactions (Fig. 1 of Tideswell et al. 2010).
Therefore, gas phase chemistry alone could explain the HNCO abundance in L1157-B1.
However, it is unclear whether a HNCO abundance as high as that
measured in L1157-B2 (
10-7) can be attained
in a few thousand years
only with gas phase reactions.
Therefore, we reckon that the most likely explanation of the high HNCO abundances in L1157 is dust grain mantles processing by the shock waves followed by neutral-neutral reactions in gas phase.
4.3 Chemical differences between B1 and B2 and the role of the O2 chemistry
We have discussed the line profiles and the line intensities of different species in Sect. 3.1. The HNCO profiles are very similar to those of SO and SO2. In addition, the lines of most of the molecules have similar intensities in B1 and B2 with the exception of CN, which is more intense in B1, and HNCO, SO2, and OCS, whose lines are more intense emission in B2 than in B1. These similarities between the HNCO and the sulfured molecules is somewhat puzzling.
The true reason for the chemical differences in B1 and B2 is not clear. The present gas density in B1 and B2 is similar (see Figs. 4 and 5 and Nisini et al. 2007; Bachiller & Perez Gutierrez 1997). In contrast, based on the line profiles, the shock velocity may be higher in B1 than in B2. Taking into account the extreme SiO line wings, the difference could reach 10 km s-1, which according to the Gusdorf et al. (2008a) models could be significant.
Bachiller & Perez Gutierrez (1997) suggested that the chemical
differences between B1 and B2 could alternatively
be due to different shock ages.
In contrast to SO2 and OCS,
the H2S abundance is indeed lower in B2 than in B1.
This is logical since H2S is a parent molecule for other
sulfured species.
H2S is formed on the grain
surfaces and released to the gas
phase by the effect of the shock waves.
Other sulfur-bearing molecules such as
SO and SO2 are produced in the gas phase very quickly (few 103 yr)
via reactions with H, OH, and O2 (Wakelam et al. 2005; Pineau des Forêts et al. 1993).
One possible explanation of the differences between B1 and B2 is that the
B2 shock could be older than that in B1
and that H2S has been converted into SO and SO2.
Even if there are uncertainties in the sulfur chemistry
(see for instance Codella et al. 2005, and references therein), sulfur-bearing
molecules such as H2S, SO, SO2, H2CS, and OCS have been invoked to be
potential valuable
probes of chemical evolution (Wakelam et al. 2005; Herpin et al. 2009; Codella et al. 2005).
In particular, Wakelam et al. (2005) found that the SO2/SO ratio
increases with the shock age.
The higher SO2/SO ratio in B2 than in
B1 is also indicative of an older shock in B2.
Bachiller & Perez Gutierrez (1997) also
proposed that the different CN abundance in B1 and B2
can be accounted for in
this scenario of an older shock in B2,
since after an enhancement of CN in the shock, reactions such as
CN+O
CO+N could be very efficient in destroying CN.
Our data globally agree with this scenario, nevertheless
the clear anticorrelation of the CN and HNCO intensities and the
similar abundances of CN and HNCO in L1157-B1 suggest that
all the CN can indeed be transformed into HNCO once the shock
has increased the temperature and the neutral-neutral reactions
form NCO from CN and O2 and HNCO from NCO (see Sect. 4.2).
Therefore, the chemical differences between B1 and B2 do support
the scenario proposed in Sect. 4.2 to explain the
HNCO abundances in the L1157 shocks.
In addition, the scheme proposed in Sect. 4.2
also provides insight into the possible
link between HNCO and the sulfured molecules, in particular
SO and SO2.
In shocks, these molecules can be produced by the reactions
and
(Wakelam et al. 2005; Pineau des Forêts et al. 1993; Charnley 1997).
Therefore, the common link with HNCO would be formation pathways
involving O2.
Another important shock tracer that
could be linked to the O2 chemistry is SiO.
Nevertheless, the situation regarding this molecule is more complex.
First, there are models that can explain the SiO emission in shocks
without a significant contribution of Si oxidation in gas phase,
either by sputtering in gas-grain collisions
if there is already SiO in the grain mantles (Gusdorf et al. 2008b)
or by dust vaporization in grain-grain collisions
(Guillet et al. 2009).
Second, if one assumes that SiO is formed in gas phase, there is
a threshold shock velocity of 25 km s-1 to
eject Si from the grain cores by sputtering (Gusdorf et al. 2008a).
In contrast, other molecules such as the organics are inside the grain mantles
and there is no such high velocity threshold for them to be
ejected into the gas phase.
Therefore, if SiO comes directly from the grains (Guillet et al. 2009; Gusdorf et al. 2008b)
or if it is formed in the gas phase (Gusdorf et al. 2008a), because of the shock velocity
threshold, we do not expect to find similarities (and they are not observed) between the SiO emission and that of SO, SO2, and HNCO.
We note that the impressive SiO blue wing in B1 (Bachiller & Perez Gutierrez 1997) is consistent with SiO formation having occurred, in comparison to other molecules, in the higher velocity gas. In this context, chemical differences between B1 and B2 may also be due to different shock velocities. If the shock velocity in B2 were actually lower than in B1 and the SiO formation less efficient, more O2 would be available in B2 to form HNCO, SO, and SO2.
4.4 On the origin of the HNCO emission in galactic nuclei
The HNCO abundance in the molecular clouds in the center of the Milky Way
and in the nuclei of starburst galaxies is similar or a bit higher (at most
by a factor of 2) than in Galactic hot cores (Table 4).
HNCO becomes another piece in the well known puzzle of
the chemistry of the Galactic center molecular clouds.
This puzzle can be summarized as follows.
The abundance of SiO and complex organic molecules in the Galacic center
is as high as in hot cores
(Martín-Pintado et al. 1997; Rodriguez-Fernandez et al. 2006; Requena-Torres et al. 2006).
In contrast, (i) the emission in the Galactic center
is extended over the central 300 pc and does not resemble
a collection of discrete sources with the size of a hot core (0.1 pc);
(ii) the gas density (103-104 cm-3)
is much lower than in hot cores; and (iii) the dust temperatures (<30 K) are also much lower than in hot cores.
The Galactic center clouds present a ``hot core chemistry without hot cores''
(Requena-Torres et al. 2006).
The origin of this chemistry is not known although it is understood to be
triggered by
some type of mechanical processes such as shock waves
(Rodríguez-Fernández et al. 2004; Martín-Pintado et al. 1997,2001).
The origin of the shocks may be related to the complex dynamics in
the inner regions of the Galaxy
(Rodriguez-Fernandez et al. 2006; Rodriguez-Fernandez & Combes 2008; Hüttemeister et al. 1998).
Based on the spatial distribution of the HNCO and the comparison with other species, such as CH3OH, Meier & Turner (2005) also suggested that the HNCO emission in IC 342 traces shocks. Martín et al. (2009,2008) also proposed that shocks may explain the high HNCO abundances measured in galactic nuclei.
Discussing the precise origin of shocks in Galactic nuclei is beyond the scope of this paper. Nevertheless, our observations support the scenario of HNCO tracing shocks in galactic nuclei since our L1157 results probe a medium of moderate H2 density, where the HNCO abundance is indeed high because of the grain processing and gas heating by shock waves.
5 Summary
We have observed three lines of HNCO towards the protostar L1157 and its
associated molecular outflow.
Our aim has been to characterize the HNCO properties in shocked gas.
HNCO was well detected in the shocked gas,
where the abundance was found to increase by a factor of 6-34
with respect to the abundance in the protostar L1157-mm.
The abundance in B1 and B2 is (0.4-
and (0.3-
,
respectively.
The abundance in B2 is the
highest ever measured, considerably higher than those in hot cores
(a few 10-9) and galactic nuclei (10-9-10-8).
Our results demonstrate that the HNCO abundances measured in galactic nuclei
can easily be attained in shocked gas, providing a solid
basis to previous suggestions that the extended HNCO in galactic
nuclei could trace large-scale shocks (Minh & Irvine 2006; Meier & Turner 2005).
The dominant formation pathway of HNCO in hot cores is grain mantle evaporation of complex molecules formed from HNCO and subsequent dissociation to produce again HNCO. In addition, there is contribution from gas phase reactions, but this is minor because of the high activation barriers of some reactions (Tideswell et al. 2010). In the L1157 shocks, we propose a relatively similar formation pathway, although with different relative importances of the dust surface and gas phase chemistry. First there should be grain mantle erosion by the shock waves that will increase the HNCO abundance in the gas phase. Second, in contrast to hot cores, the neutral-neutral reactions in the gas phase can be the dominant formation pathway due to the higher gas temperatures in the shocked gas with respect to the typical temperature of hot cores. These reactions start with the formation of OCN from CN and O2, which is expected to be very abundant in the gas phase of a magneto-hydrodynamic shock of relatively low velocity (20-30 km s-1) such as the L1157 shocks (Gusdorf et al. 2008a). The observed anticorrelation of CN and HNCO fluxes also support this scenario.
HNCO line profiles are very similar to those of CH3OH and H2CO lines, which is not surprising since all those molecules have, at least partially, a common origin linked to the grain surface chemistry and subsequent desorption. In addition, HNCO line profiles are also similar to some sulfur-bearing molecules such as SO and SO2. Another common property of HNCO and the sulfur-bearing molecules is that they are the only species that exhibit a more intense emission in B2 than B1. We propose that the common link between HNCO and the sulfured molecules is that in both cases there are formation pathways with reactions involving O2. HNCO and the sulfured molecules are reliable probes of chemical differentiation in shocks. In particular, in the case of L1157 the higher abundances in B2 than in B1 could be due to a lower shock velocity and/or an older shock in B2 than in B1.
AcknowledgementsN.J.R.-F. acknowledges fruitful discussions on chemistry and molecular outflows with S. Cabrit, A. Fuente and C. Codella. This work has been partially funded by the grant ANR-09-BLAN-0231-01 from the French Agence Nationale de la Recherche as part of the SCHISM project. We thank an anonymous referee for a very detailed report that helped to correct and clarify a number of issues in the original version of the manuscript.
References
- Arce, H. G., Santiago-García, J., Jørgensen, J. K., Tafalla, M., & Bachiller, R. 2008, ApJ, 681, L21 [NASA ADS] [CrossRef] [Google Scholar]
- Bachiller, R., & Perez Gutierrez, M. 1997, ApJ, 487, L93 [NASA ADS] [CrossRef] [Google Scholar]
- Bachiller, R., Pérez Gutiérrez, M., Kumar, M. S. N., & Tafalla, M. 2001, A&A, 372, 899 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Benedettini, M., Viti, S., Codella, C., et al. 2007, MNRAS, 381, 1127 [NASA ADS] [CrossRef] [Google Scholar]
- Bergin, E. A., Neufeld, D. A., & Melnick, G. J. 1998, ApJ, 499, 777 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- 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]
- Blake, G. A., Sutton, E. C., Masson, C. R., & Phillips, T. G. 1987, ApJ, 315, 621 [NASA ADS] [CrossRef] [Google Scholar]
- Bockelée-Morvan, D., Lis, D. C., Wink, J. E., et al. 2000, A&A, 353, 1101 [NASA ADS] [Google Scholar]
- Caselli, P., Hasegawa, T. I., & Herbst, E. 1993, ApJ, 408, 548 [NASA ADS] [CrossRef] [Google Scholar]
- Charnley, S. B. 1997, ApJ, 481, 396 [NASA ADS] [CrossRef] [Google Scholar]
- Churchwell, E., Wood, D., Myers, P. C., & Myers, R. V. 1986, ApJ, 305, 405 [NASA ADS] [CrossRef] [Google Scholar]
- Codella, C., Bachiller, R., Benedettini, M., et al. 2005, MNRAS, 361, 244 [NASA ADS] [CrossRef] [Google Scholar]
- Codella, C., Benedettini, M., Beltrán, M. T., et al. 2009, A&A, 507, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dahmen, G., Huettemeister, S., Wilson, T. L., et al. 1997, A&AS, 126, 197 [Google Scholar]
- Demyk, K., Dartois, E., D'Hendecourt, L., et al. 1998, A&A, 339, 553 [NASA ADS] [Google Scholar]
- Garrod, R. T., Weaver, S. L. W., & Herbst, E. 2008, ApJ, 682, 283 [NASA ADS] [CrossRef] [Google Scholar]
- Gueth, F., Guilloteau, S., & Bachiller, R. 1998, A&A, 333, 287 [NASA ADS] [Google Scholar]
- Guillet, V., Jones, A. P., & Pineau Des Forêts, G. 2009, A&A, 497, 145 [Google Scholar]
- Gusdorf, A., Cabrit, S., Flower, D. R., & Pineau Des Forêts, G. 2008a, A&A, 482, 809 [Google Scholar]
- Gusdorf, A., Pineau Des Forêts, G., Cabrit, S., & Flower, D. R. 2008b, A&A, 490, 695 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hatchell, J., Thompson, M. A., Millar, T. J., & MacDonald, G. H. 1998, A&A, 338, 713 [NASA ADS] [Google Scholar]
- Helmich, F. P., & van Dishoeck, E. F. 1997, A&AS, 124, 205 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Herpin, F., Marseille, M., Wakelam, V., Bontemps, S., & Lis, D. C. 2009, A&A, 504, 853 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hüttemeister, S. 1993, PhD Thesis, Dissertation, Bonn University [Google Scholar]
- Hüttemeister, S., Dahmen, G., Mauersberger, R., et al. 1998, A&A, 334, 646 [NASA ADS] [Google Scholar]
- Iglesias, E. 1977, ApJ, 218, 697 [NASA ADS] [CrossRef] [Google Scholar]
- Jackson, J. M., Armstrong, J. T., & Barrett, A. H. 1984, ApJ, 280, 608 [NASA ADS] [CrossRef] [Google Scholar]
- Jansen, D. J., Spaans, M., Hogerheijde, M. R., & van Dishoeck, E. F. 1995, A&A, 303, 541 [NASA ADS] [Google Scholar]
- Kuan, Y., & Snyder, L. E. 1996, ApJ, 470, 981 [NASA ADS] [CrossRef] [Google Scholar]
- Larsson, B., Liseau, R., Pagani, L., et al. 2007, A&A, 466, 999 [Google Scholar]
- Lindqvist, M., Sandqvist, A., Winnberg, A., Johansson, L. E. B., & Nyman, L. 1995, A&AS, 113, 257 [NASA ADS] [Google Scholar]
- MacDonald, G. H., Gibb, A. G., Habing, R. J., & Millar, T. J. 1996, A&AS, 119, 333 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Marcelino, N., Cernicharo, J., Tercero, B., & Roueff, E. 2009, ApJ, 690, L27 [NASA ADS] [CrossRef] [Google Scholar]
- Martín-Pintado, J., de Vicente, P., Fuente, A., & Planesas, P. 1997, ApJ, 482, L45 [NASA ADS] [CrossRef] [Google Scholar]
- Martín-Pintado, J., Rizzo, J. R., de Vicente, P., Rodríguez-Fernández, N. J., & Fuente, A. 2001, ApJ, 548, L65 [NASA ADS] [CrossRef] [Google Scholar]
- Martín, S., Requena-Torres, M. A., Martín-Pintado, J., & Mauersberger, R. 2008, ApJ, 678, 245 [NASA ADS] [CrossRef] [Google Scholar]
- Martín, S., Martín-Pintado, J., & Mauersberger, R. 2009, ApJ, 694, 610 [NASA ADS] [CrossRef] [Google Scholar]
- Mauersberger, R., Guelin, M., Martin-Pintado, J., et al. 1989, A&AS, 79, 217 [NASA ADS] [Google Scholar]
- Meier, D. S., & Turner, J. L. 2005, ApJ, 618, 259 [NASA ADS] [CrossRef] [Google Scholar]
- Minh, Y. C., & Irvine, W. M. 2006, New Astron., 11, 594 [NASA ADS] [CrossRef] [Google Scholar]
- Minh, Y. C., Haikala, L., Hjalmarson, A., & Irvine, W. M. 1998, ApJ, 498, 261 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Minh, Y. C., Kim, S., Pak, S., et al. 2005, New Astron., 10, 425 [NASA ADS] [CrossRef] [Google Scholar]
- Nguyen-Q-Rieu, Henkel, C., Jackson, J. M., & Mauersberger, R.1991, A&A, 241, L33 [NASA ADS] [Google Scholar]
- Nisini, B., Codella, C., Giannini, T., et al. 2007, A&A, 462, 163 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pineau des Forêts, G., Roueff, E., Schilke, P., & Flower, D. R. 1993, MNRAS, 262, 915 [NASA ADS] [CrossRef] [Google Scholar]
- Requena-Torres, M. A., Martín-Pintado, J., Rodríguez-Franco, A., et al. 2006, A&A, 455, 971 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rizzo, D., Hüttemeister, S., & Dahmen, G. 2000, in Molecules in Space and in the Laboratory, ed. I. Porceddu & S. Aiello, 141 [Google Scholar]
- Rodriguez-Fernandez, N. J., & Combes, F. 2008, A&A, 489, 115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rodríguez-Fernández, N. J., Martín-Pintado, J., Fuente, A., & Wilson, T. L. 2004, A&A, 427, 217 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rodriguez-Fernandez, N. J., Combes, F., Martin-Pintado, J., Wilson, T. L., & Apponi, A. 2006, A&A, 455, 963 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schöier, F. L., van der Tak, F. F. S., van Dishoeck, E. F., & Black, J. H. 2005, A&A, 432, 369 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Snyder, L. E., & Buhl, D. 1972, ApJ, 177, 619 [NASA ADS] [CrossRef] [Google Scholar]
- Tideswell, D. M., Fuller, G. A., Millar, T. J., & Markwick, A. J. 2010, A&A, 510, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Turner, B. E., Terzieva, R., & Herbst, E. 1999, ApJ, 518, 699 [NASA ADS] [CrossRef] [Google Scholar]
- van Broekhuizen, F. A., Keane, J. V., & Schutte, W. A. 2004, A&A, 415, 425 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van der Tak, F. F. S., Black, J. H., Schöier, F. L., Jansen, D. J., & van Dishoeck, E. F. 2007, A&A, 468, 627 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van Dishoeck, E. F., Blake, G. A., Jansen, D. J., & Groesbeck, T. D. 1995, ApJ, 447, 760 [NASA ADS] [CrossRef] [Google Scholar]
- Wakelam, V., Ceccarelli, C., Castets, A., et al. 2005, A&A, 437, 149 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zinchenko, I., Henkel, C., & Mao, R. Q. 2000, A&A, 361, 1079 [NASA ADS] [Google Scholar]
Footnotes
- ... 0
- HNCO is a planar, nearly linear, slightly asymmetric prolate rotor. The notation for a given level is JK-1K1. The K-1 = 0 are usually excited thermally, while the excited K ladders ( K-1 > 0) are probably excited by FIR radiation (Churchwell et al. 1986).
- ...
.
- HNCO formation on the grain surfaces could also occur by means of reactions of OCN- with NH4+ or H3O+. The inverse reactions (HNCO with NH3 or H2O) are invoked to explain the presence of OCN- ices in the dust grains, (Demyk et al. 1998; van Broekhuizen et al. 2004). To our knowledge, these reactions have never been included in chemical models.
- ...
). - This does not mean that the current O2 abundance in the L1157 shocks is high, but it may have been high just after the passage of a shock-wave as predicted by shock models. Since the chemistry in shocked regions proceeds rapidly, O2can be rapidly converted into other species. Therefore, this is not inconsistent with the O2 abundance in the interstellar medium being found to be low (<10-7 Larsson et al. 2007).
All Tables
Table 1: Observational parameters.
Table 2: Measured fluxes.
Table 3: Results.
Table 4: HNCO abundance with respect to H2 in different astronomical environments.
All Figures
![]() |
Figure 1: Integrated CO(2-1) emission from the L1157 outflow (Bachiller et al. 2001). We show the position of L1157-mm (star), B1 (empty square), and B2 (filled square). B1 and B2 are the positions of the shocks as traced, for instance, by SiO(3-2) (see Fig. 6 of Bachiller et al. 2001). The CO emission peaks behind the true positions of the shocks. The different circles represent the half maximum contours of the IRAM 30m beam at the frequencies of the lines discussed in this paper (see Table 1). |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Observed HNCO spectra. From bottom to top, the different panels show the HNCO 404-303, 505-404, and 606-505 spectra towards L1157 mm, B1, and B2. |
Open with DEXTER | |
In the text |
![]() |
Figure 3: Comparison of the HNCO(505-404) line profiles (in red) with the spectra of other species published by Bachiller & Perez-Gutierrez (1997), in black lines. The HNCO peak temperature has been scaled to the peak temperature of the other molecule to compare profiles. |
Open with DEXTER | |
In the text |
![]() |
Figure 4:
RADEX results for the 505-404 to 404-303
flux ratio
(horizontal curves) and the 404-303 line flux
(vertical curves) computed with a
brightness temperature for a source size of 10''
as a function of the hydrogen density, the HNCO column density and the kinetic
temperature (30 K in the upper panel and 250 K in the lower panel). A typical
uncertainty of 10 |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Same than Fig. 4 but using
|
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.