Issue |
A&A
Volume 507, Number 2, November IV 2009
|
|
---|---|---|
Page(s) | L25 - L28 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/200913340 | |
Published online | 27 October 2009 |
A&A 507, L25-L28 (2009)
LETTER TO THE EDITOR
Methyl cyanide as tracer of bow shocks in L1157-B1
C. Codella1 - M. Benedettini2 - M. T. Beltrán1 - F. Gueth3 - S. Viti4 - R. Bachiller5 - M. Tafalla5 - S. Cabrit6 - A. Fuente5 - B. Lefloch7
1 - INAF - Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy
2 -
INAF, Istituto di Fisica dello Spazio Interplanetario, via Fosso del Cavaliere 100, 00133 Roma, Italy
3 -
IRAM, 300 rue de la Piscine, 38406 Saint-Martin-d'Hères, France
4 -
Department of Physics and Astronomy, University College London, Gower Street, London WC1E6BT, UK
5 -
Observatorio Astrónomico Nacional (IGN), Alfonso XII, 28014 Madrid, Spain
6 -
LERMA, UMR 8112 du CNRS, Observatoire de Paris, 61 Av. de l'Observatoire, 75014 Paris, France
7 -
Laboratoire d'Astrophysique de l'Observatoire de Grenoble, BP 53, 38041 Grenoble Cedex, France
Received 22 September 2009 / Accepted 20 October 2009
Abstract
Context. L1157-mm is a low-mass protostar driving an outflow
which is considered the proto-type of the so-called chemically active
outflows. In particular, the blue-shifted lobe B1 stands out for
its rich mm-wave spectrum.
Aims. Our aim is to infer the physical conditions of the molecular gas within L1157-B1.
Methods. We carried out CH3CN(8K-7K) observations at 2 mm with the IRAM Plateau de Bure Interferometer to map L1157-B1.
Results. The CH3CN image shows a clumpy structure
superimposed to the classical B1 arch-like shape, displaying a
unique continuous structure tracing the propagation of a large bow
shock. The CH3CN spatial distribution well agrees with that of CH3OH, one of the most abundant species evaporated from dust grain mantles. Thanks to the present CH3CN
observations at unprecedented spatial resolution, we determine for the
first time the temperature distribution of the B1 bow shock. We
obtain values between 57 and 132 K, with the highest
temperature at the shock apex, thus confirming the association of CH3CN with gas affected by the passage of a shock wave. The
column densities are
8-40
1012 cm-2, depending on the position, leading to CH3CN/CH3OH abundance ratios
0.2-1.3
10-3, i.e. values on average smaller than those found in hot-corinos. This suggests that a significant amount of CH3CN is formed in the gas phase: in the shocked L1157-B1 regions the CH3CN abundance may not have reached its maximum yet.
Key words: ISM: individual objects: L1157 - ISM: molecules - stars: formation
1 Introduction
During the earliest evolutionary stages, young stellar objects (YSOs)
produce high-velocity jets which drive shocks travelling through the
surrounding high-density medium. Jets compress, and with time, the
ambient medium creates slower molecular outflows formed by swept-up
material. Shocks trigger several processes such as endothermic chemical
reactions and ice grain mantle sublimation, which lead to a dramatic
enhancement of the abundance of several species (see e.g. van
Dishoeck & Blake 1998,
and references therein). This has been demonstrated by multiline
observations at millimeter wavelengths towards a number of outflows
(e.g. Garay et al. 1998; Bachiller & Pérez Gutiérrez 1999; Jörgensen et al. 2007).
In particular, the L1157 outflow, powered by IRAS 20386+6751,
may be regarded as the typical laboratory to observe the effects of
shocks on the gas chemistry (Bachiller et al. 2001). IRAS 20386+6751 is a low-luminosity (11 )
class 0 protostar, located at 440 pc driving a precessing
powerful molecular outflow that shows several bow shocks seen in CO
(Gueth et al. 1996) as well as in NIR H2 images (Davis & Eislöffel 1995).
Interferometric PdBI and VLA observations (Tafalla & Bachiller 1995; Gueth et al. 1998; Benedettini et al. 2007,
hereafter Paper I) have mapped the brightest blue-shifted
bow-shock, called B1, in (i) molecules thought to be released by
dust mantles
such as CH3OH and NH3; (ii) simple species,
mainly sulphuretted such as CS and SO, expected to form, although less
efficiently, also in the ambient gas; and (iii) SiO, which is the
typical tracer of high-speed shocks, formed after the sputtering of
silicon from dust (Gusdorf et al. 2008a,b).
These observations revealed a rich and clumpy structure, the clumps
being located at the wall of the cavity with an arch-shape.
Recently, Arce et al. (2008) have reported the detection of complex organic molecules such as HCOOCH3, HCOOH, C2H5OH, and CH3CN towards L1157-B1. Such detections, made with the 30-m IRAM single-dish telescope, suggest that these complex molecules, usually considered as suitable hot-core tracers, can also increase in abundance due to the impact of outflows on the surrounding medium and not only due to a protostar's radiation. This finding is particularly important for CH3CN (methyl cyanide), whose emission has been extensively used to investigate the surroundings of high-mass (the so called hot-cores; e.g. Churchwell et al. 1992; Olmi et al. 1993, 1996a,b) and lower mass (hot-corinos; e.g. Cazaux et al. 2003; Bottinelli et al. 2004a,b; Bisschop et al. 2008) YSOs. In particular, CH3CN emission has been used to study the kinematics of hot-cores (Beltrán et al. 2004, 2005), revealing rotating massive toroids expected to host obscured accretion disks.
In this letter, we report for the first time the high-resolution image of CH3CN towards the L1157-B1, showing that this emission is issued from warm gas (60-130 K) tracing the bow structure created by the passage of shocks.
![]() |
Figure 1:
Chemical differentiation in L1157-B1: the maps are centred at the following coordinates:
|
Open with DEXTER |
2 Observations
The CH3CN observations of L1157-B1 were obtained in
August 2008 and March 2009 with the IRAM Plateau de Bure
Interferometer (PdBI) in France. Two configurations were used: the
6-element C configuration (baselines from 24 m up to
176 m), and the 5-element D configuration (from 24 m
to 97 m). The CH3CN(8K-7K) spectral pattern at 147 174.594 MHz (K = 0) was observed with a 320 MHz bandwidth and a resolution of 2.5 MHz (5.1 km s-1).
Amplitude and phase were calibrated by observing 1849+670 and 2021+614,
whereas the flux density scale was derived by observing 3C 84 and
3C 454.3, with an uncertainty of
20%. Images were produced using natural weighting, and were restored with a clean beam of
(PA = 10
).
3 CH3CN as tracer of shocked gas
The contours of Fig. 1 show the map of the CH3CN(8K-7K) emission integrated over the K = 0-2 components. The CH3CN structure is well in agreement with that of other molecular species observed by Tafalla & Bachiller (1995), Gueth et al. (1998), and Paper I with the VLA and PdB interferometers, revealing a clumpy structure superimposed to the classical B1 arch-like shape seen e.g. in CO. Beside the already detected B1a, B1b, B1c clumps, the CH3CN emission shows further emission peaks in B1, here labelled B1e and B1f. In addition, we have detected (i) the north-eastern B0e clump and (ii) a new clump, labelled B0f, towards the north-western portion of the map. In summary, the present CH3CN map clearly confirms B1 together with the northern B0 clumps as a unique continuous structure tracing the propagation of the large bow shock which created the large cavity called C2 (see e.g. Fig. 7 of Paper I) which was detected in CO by Gueth et al. (1996). The walls of the cavity are outlined well both by the B0e and B0f clumps and by the weaker CH3CN eastern emission between the clumps. In this context, what we have so far called clumps should not necessarily be considered isolated roundish bodies, but rather emission peaks belonging to a more elongated structure.
Figure 1 compares the CH3CN emission with those of CH3OH, SiO and CS (Paper I; Gueth et al. 1998), mapped with similar spatial resolution. These previous observations have revealed a dichotomy suggesting an inhomogeneity in the physical conditions and/or chemical composition along B1: a group of molecules (CS, SiO, and NH3, not shown here) peaking at the east-side and another group of species (CH3OH, SO, and OCS, not shown here) peaking at the west-side of the cavity. The brightest peak of CH3CN is B1b, thus well coinciding with the western CH3OH peak. Bright CH3CN emission is also found in the elongated structure composed by the B1a, B1e, and B1f peaks. On the other hand, CH3CN is anticorrelated with both the apex of the bow structure as well as the so-called SiO and CS finger which, according to Gueth et al. (1998), could trace a magnetic precursor of the shock.
![]() |
Figure 2:
Left panel: CH3CN J = 8K-7K spectra (in |
Open with DEXTER |
In conclusion, the present findings confirm the complexity of a region interested by the passage of a bow shock: for B1, in particular, CH3CN seems to introduce a third group of molecules peaking in the northern part of the B1 structure, quite apart from the effects seen in the west and east dichotomy and the CS/SiO finger. This finding suggests temporal effects, with CH3CN mostly confined to the back part of the shock. The possible implication on the formation mechanism of CH3CN is reported in Sect. 5.
4 The temperature structure in L1157-B1
The temperature estimates of L1157-B1 so far obtained are based on single-dish observations of CO, NH3, and CH3CN (Umemoto et al. 1999; Hirano & Taniguchi 2001; Arce et al. 2008).
The values are in the 50-170 K range, and they have to be
considered as measurements averaged on half power beam widths (HPBWs) of 8
-14
(CO), 72
(NH3), and 10
(CH3CN). Also the multi-lines survey of a classical shock tracer like SiO (Nisini et al. 2007), which shows temperatures between 100 and 300 K, has been performed with HPBWs
9
-27
.
An estimate based on a better spatial resolution (5-6
)
has been obtained with NH3 emission imaged at VLA by Tafalla & Bachiller (1995), but a unique value (80 K) for the whole L1157-B1 structure has been reported.
Figure 2 (left panels) shows the CH3CN(8K-7K) spectra integrated over the area of the L1157-B1 clumps. Several K-components (up to K = 4)
are clearly detected, which have been fitted by fixing their separation
in frequency in the laboratory values and forcing their line width to
be identical. For a symmetric-top molecule such as CH3CN, the population of one K-ladder relative to another is established by collisions and is characterised by a rotation temperature (
)
which is approximately equal to the kinetic temperature (e.g. Loren & Mundy 1984). We derived estimates for the rotational temperature,
,
and the total CH3CN column density,
,
by means of the standard rotation diagram method.
Note that the error bars refer to the line-fit uncertainty and do not
take into account the absolute flux scale error (20%), which is
affecting all lines simultaneously. This means that the temperatures,
identified by the slope of the fit, are not affected by the error on
absolute flux. Figure 2
(right panels) shows the rotation diagrams: for the first time,
it has been possible to estimate a temperature for several
L1157-B1 clumps, thus allowing to derive a temperature distribution
covering the whole B1 structure. The rotational temperatures are
in the 57-132 K range. The higher
has
been measured at the B1c position, near the apex of the bow
structure and the CS/SiO finger, which should lie along the jet
direction and represent the magnetic precursor of the shock (Gueth
et al. 1998).
In addition, the temperature estimates are lower at the edge of the
shock (67 K for B1a and 73 K for B1b). Also the B0e and
B1f clumps, where only 4 K-components have been detected (up to K = 3), show a relatively low temperature: 57 K and 55 K, respectively; note that the upper limits to the K = 4
integrated intensity are consistent with a temperature less than
60 K, confirming the result. On the other hand,
is higher in B1e (92 K), which is located on the jet axis, and that is where the highest blue-shifted velocities (
10 km s-1
)
have been observed (Gueth et al. 1998;
Paper I). These findings seem to confirm the presence of a
bow-shock with a temperature gradient decreasing with distance from the
apex, supporting the association of CH3CN with shocked gas.
Interestingly, a fit of the B1c rotation diagram without the K = 4 component leads to
= 72 K (instead of 132 K). By excluding the K =
4 line also the B1b temperature decreases from 73 K to
55 K and that of B1e from 92 K to 81 K; on the opposite
the B1a rotation diagram keeps the same solution (67 K) with
or without the K = 4 emission. These results clearly
illustrate the importance of probing temperatures in regions associated
with high temperatures, such as shocked regions, by using high
excitation transitions. These results suggest two gas components with
different temperatures and could be related to the SiO results
obtained by Nisini et al.
(2007), which show evidence for
velocity components having different excitation conditions, with the
denser and/or warmer gas associated with higher speeds: further high
spectral resolution observations would be instructive.
5 On the origin of CH3CN in bow-shocks
The origin of CH3CN molecules in star forming regions is still a matter of debate (see e.g. Bottinelli et al. 2007; and Bisschop et al. 2008, for detailed reviews). One possibility is that they form on the surface of grains from reactions of CH3 with CN and are then injected in the gas phase, together with species such as H2O, H2CO, and CH3OH, once the
icy grain mantles evaporate. Otherwise, CH3CN
may belong to the so-called second generation species, i.e. the
molecules formed in the warm gas-phase chemically enriched after the
mantle evaporation; in this case CH3CN would form from HCN with CH3+.
Both scenarios need an increase of temperature, due to stellar
radiation in case of hot-core and hot-corinos or the passage of a shock
in the present case, confirmed by the high
measured
here. As chemical models for hot-cores and corinos have shown, the
warm-up phase plays a major role in both scenarios (Viti et al. 2004; Garrod & Herbst 2006).
Under the assumption of the more realistic gradual warm-up of
temperatures instead of an instantaneous jump from 10 to several
hundred kelvins, then (i) the changes in the chemical evolution of
each species due to differential desorption are important and several
species, including CH3OH and CH3CN show a strong time
dependence (Viti et al. 2004) and (ii) gas-phase and grain-surface reactions are strongly coupled during the warm-up phase (Garrod & Herbst 2006). For example, CN and CH3 may become more mobile on the grains enhancing the production of CH3CN.
On the other hand, if instead the evaporation is caused by the passage
of the shock, it will be effectively instantaneous due to
sputtering; hence only surface reactions during the freeze-out process
and gas phase reactions would be relevant for the production of CH3CN.
Table 1: The L1157-B1 clumps: CH3CN rotational temperatures, and CH3CN and CH3OH column densities.
The
column densities, reported in Table 1, are in the 8-40
1012 cm-2 range. These values are larger than those given by Arce et al. (2008), 1
1012 cm-2: this is not surprising, since their estimate is based on a single-dish spectrum producing a CH3CN brightness and temperature averaged with a HPBW
10
,
affected by beam dilution due to the clumpy structure of the emission.
We derived the abundance ratio of CH3CN with respect to CH3OH, which is a well known abundant molecule injected in the gas phase from warmed dust mantles. By using the
values calculated in Paper I, we obtain [CH3CN]/[CH3OH] abundance ratios in the 0.2-1.3
10-3 range, depending on the clump. Even with higher
values with respect to Arce et al. (2008), the CH3CN-to-CH3OH ratios are still on average smaller with respect to those found in hot-corinos (0.3-1
10-1; Bottinelli et al. 2007; Bisschop et al. 2008), as well as those found in massive hot-cores (
10-2-10-1; Bisschop et al. 2007). This difference suggests that asignificant amount of CH3CN may form in the gas phase and that in the shocked L1157-B1 regions the CH3CN abundance
has not yet reached its maximum, given that the age of the warmgas in
L1157-B1 has been estimated to be less than 2000 yr (Arce
et al. 2008). This possibility is supported by the relative CH3CN weakness of the apex region of the shock compared to the rear of the bow. The [CH3CN]/[CH3OH] ratio,
however, is found to vary by several orders of magnitude among
hot-cores/corinos, depending on their physical properties as well as
their evolutionary stage.
6 Conclusions
The present high spatial resolution map shows that CH3CN
can be efficiently produced in regions far from the protostar and
affected by the passage of a shock wave. Thanks to the analysis of the
CH3CN(8K-7K) spectral
pattern, we obtained for the first time estimates of the temperature at
six different positions covering the whole B1 structure. This is
essential for the correct analysis of the forthcoming spectral survey
of L1157-B1 that will be carried out by the ESA Herschel satellite in
Guaranteed Time. The survey will cover the 500-2000 GHz spectral
window to obtain an accurate census of the molecular content in star
forming regions. The knowledge of the physical properties of the
L1157-B1 clumps, used as proto-type of shocked regions,
is needed to interpret the Herschel spectra, which will be
obtained with a 10
-40
spatial resolution.
We thank R. Cesaroni for helpful discussion.
References
- Arce, H. G., Santiago-García, J., Jörgensen, J. K., Tafalla, M., & Bachiller, R. 2008, ApJ, 681, L21 [CrossRef] [NASA ADS]
- Bachiller, R., & Peréz Gutiérrez, M. 1999, ApJ, 487, L93 [CrossRef] [NASA ADS]
- Bachiller, R., Peréz Gutiérrez, M., Kumar, M. S. N., & Tafalla, M. 2001, A&A, 372, 899 [EDP Sciences] [CrossRef] [NASA ADS]
- Benedettini, M., Viti, S., Codella, C., et al. 2007, MNRAS, 381, 1127 [CrossRef] [NASA ADS] (Paper I)
- Beltrán, M. T., Cesaroni, R., Neri, R., et al. 2004, ApJ, 601, L187 [CrossRef] [NASA ADS]
- Beltrán, M. T., Cesaroni, R., Neri, R., et al. 2005, A&A, 435, 901 [EDP Sciences] [CrossRef] [NASA ADS]
- Bisschop, S. E., Jörgensen, J. K., van Dishoeck, E. F., & de Wachter, E. B. M. 2007, A&A, 465, 913 [EDP Sciences] [CrossRef] [NASA ADS]
- Bisschop, S. E., Jörgensen, J. K., Bourke, T. L., Bottinelli, S., & van Dishoeck, E. F. 2008, A&A, 488, 959 [EDP Sciences] [CrossRef] [NASA ADS]
- Bottinelli, S., Ceccarelli, C., Lefloch, B., et al. 2004a, ApJ, 615, 354 [CrossRef] [NASA ADS]
- Bottinelli, S., Ceccarelli, C., Neri, R., et al. 2004b, ApJ, 617, L69 [CrossRef] [NASA ADS]
- Bottinelli, S., Ceccarelli, C., Williams, J. P., & lefloch, B. 2007, A&A, 463, 601 [EDP Sciences] [CrossRef] [NASA ADS]
- Cazaux, S., Tielens, A. G. G. M., Ceccarelli, C., et al. 2003, ApJ, 593, L51 [CrossRef] [NASA ADS]
- Churchwell, E., Walmsley, C. M., & Wood, D. O. S. 1992, A&A, 253, 541 [NASA ADS]
- Davis, C. J., & Eislöffel, J. 1995, A&A, 300, 851 [NASA ADS]
- Garay, G., Köhnenkamp, I., Bourke, T. L., Rodríguez, L. F., & Lehtinen, K. K. 1998, ApJ, 509, 768 [CrossRef] [NASA ADS]
- Garrod, R. T., & Herbst, E. 2006, A&A, 457, 927 [EDP Sciences] [CrossRef] [NASA ADS]
- Gueth, F., Guilloteau, S., & Bachiller, R. 1996, A&A, 307, 891 [NASA ADS]
- Gueth, F., Guilloteau, S., & Bachiller, R. 1998, A&A, 333, 287 [NASA ADS]
- Gusdorf, A., Cabrit, S., Flower, D. R., & Pineau Des Forêts, G. 2008a, A&A, 482, 809 [EDP Sciences] [CrossRef] [NASA ADS]
- Gusdorf, A., Pineau Des Forêts, G., Cabrit, S., & Flower, D. R. 2008b, A&A, 490, 695 [EDP Sciences] [CrossRef] [NASA ADS]
- Hirano, N., & Taniguchi, Y. 2001, ApJ, 550, L219 [CrossRef] [NASA ADS]
- Jörgensen, J. K., Bourke, T. L., & Myers, P. C. 2007, ApJ, 659, 479 [CrossRef] [NASA ADS]
- Loren, R. B., & Mundy, L. G. 1984, ApJ, 286, 232 [CrossRef] [NASA ADS]
- Müller, H. S. P., Schlóder, F., Stutzki, J., & Winnewisser, G. 2005, J. Mol. Struct., 742, 215 [CrossRef] [NASA ADS]
- Nisini, B., Codella, C., Giannini, T., et al. 2007, A&A, 462, 163 [EDP Sciences] [CrossRef] [NASA ADS]
- Olmi, L., Cesaroni, R., & Walmsley, C. M. 1993, A&A, 276, 489 [NASA ADS]
- Olmi, L., Cesaroni, R., & Walmsley, C. M. 1996a, A&A, 307, 599 [NASA ADS]
- Olmi, L., Cesaroni, R., Neri, R., & Walmsley, C. M. 1996b, A&A, 315, 565 [NASA ADS]
- Tafalla, M., & Bachiller, R. 1995, ApJ, 443, L37 [CrossRef] [NASA ADS]
- Ummoto, T., Mikami, H., Yamamoto, S., & Hirano, N. 1999, ApJ, 525, L105 [CrossRef] [NASA ADS]
- van Dishoeck, E. F., & Blake, G. A. 1998, ARA&A, 36, 317 [CrossRef] [NASA ADS]
- Viti, S., Collings, M. P., Dever, J. W., McCoustra, M. R. S., & Williams, D. A. 2004, MNRAS, 354, 1141 [CrossRef] [NASA ADS]
All Tables
Table 1: The L1157-B1 clumps: CH3CN rotational temperatures, and CH3CN and CH3OH column densities.
All Figures
![]() |
Figure 1:
Chemical differentiation in L1157-B1: the maps are centred at the following coordinates:
|
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Left panel: CH3CN J = 8K-7K spectra (in |
Open with DEXTER | |
In the text |
Copyright ESO 2009
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.