A&A 376, L5-L8 (2001)
DOI: 10.1051/0004-6361:20011008
J. V. Keane1 - A. M. S. Boonman2 - A. G. G. M. Tielens1,3 - E. F. van Dishoeck2
1 - Kapteyn Institute, PO Box 800,
9700 AV Groningen, The Netherlands
2 -
Leiden Observatory, PO Box 9513, 2300 RA Leiden,
The Netherlands
3 -
SRON, PO Box 800, 9700 AV Groningen, The Netherlands
Received 18 June 2001 / Accepted 11 July 2001
Abstract
We present the first detection of the
ro-vibrational
band of gas-phase
in absorption in the mid-infrared spectral region around
7.3
of a sample of deeply embedded massive protostars. Comparison with
model spectra shows that the derived excitation temperatures correlate with
previous C2H2 and HCN studies, indicating that the same warm gas
component is probed. The
column densities are similar along all lines
of sight suggesting that the
formation has saturated, but not destroyed,
and the absolute abundances of
are high (
10-7). Both the
high temperature and the high abundance of the detected
are not
easily explained by standard hot core chemistry models. Likewise,
indicators of shock induced chemistry are lacking.
Key words: star-formation: gas-phase molecules - ISM: abundances - ISM: molecules
Observations with the Infrared Space Observatory (ISO) have dramatically increased our knowledge of the active chemistry occurring within star-forming regions. Extensive studies have revealed a vast richness of solid-state molecules embedded in icy grain mantles (Ehrenfreund & Schutte 2000) which highlight the crucial role of grain surface chemistry in molecule formation. Paralleling this, the infrared and submillimeter (Boonman et al. 2000; Lahuis & van Dishoeck 2000; van der Tak et al. 2000a) observations of gas-phase molecules directly probe the chemical and physical conditions of the star-formation process. By combining the solid-state and gas-phase observations, a detailed picture of the evolving chemistry emerges which can serve as a stringent test of proposed chemical models of star-forming regions.
Sulphur-bearing species are particularly interesting to study as they
were originally proposed as tracers of shocks since the
increased availability of OH radicals will lead to enhanced abundances of
specific molecular species (Hartquist et al. 1980).
The chemistry of sulphur-bearing molecules in warm gas is essentially governed
by neutral-neutral reactions involving
formed on grain surfaces and
subsequently evaporated into the gas-phase (cf. Charnley 1997).
The destruction of
frees atomic sulphur, which can then readily
react with OH and O2 to produce SO.
is easily formed through the
conversion of SO by OH. Above
200-300 K, the OH radicals are
driven into H2O and the formation of
is halted. Except for
the detection of the
band of gas-phase
in emission towards
Orion (van Dishoeck et al. 1998), only purely rotational lines of
in the submillimeter have been detected toward massive star-forming
regions. Abundances of
10-9 (Schreyer et al. 1997) are
typically derived which are much lower than model predictions
(Charnley 1997).
abundances up to
are only found in
Orion-KL in the so-called plateau gas associated with the low-velocity
outflow (e.g. Blake et al. 1987; Sutton et al. 1995). This gas is
known to contain high abundances of OH (Melnick et al. 1987), and
hence, the formation of
is intimately connected with the
availability of reactive OH.
Here we present the first detection of infrared gas-phase
in absorption
in the
ro-vibrational band towards a sample of embedded massive
protostars.
![]() |
Figure 1: ISO-SWS AOT 6 spectra towards six massive protostars. The thin solid lines indicate the locally defined 4th order polynomials adopted as the continua. Some of the sources were offset for clarity by a constant factor indicated in the brackets. |
Open with DEXTER |
High resolution AOT6 (
)
grating mode observations of the massive protostars
presented in this article were made with the Short Wavelength
Spectrometer (SWS) on-board ISO.
The
ro-vibrational mode of
lies within Band 2C which suffers from
instrumental fringing of varying severity. In order to extract unimpeded
absorption profiles, the fringes were corrected for by dividing the
observed fluxes by cosine functions fitted to the data in wavenumber space
(Lahuis & van Dishoeck 2000). The data were flat-fielded to the
average level and then rebinned to the wavelength grid with a constant binsize
of 0.003
which corresponds to
The fully reduced 7-8
spectra are shown in Fig. 1, where the noise
level is approximately 4-5% for 3
significance.
The 7-8
spectra (Fig. 1), towards all lines of sight,
display a richness of broad and narrow absorption features
attributable to solid-state and gas-phase molecular species.
The region is dominated by the red wing of the 6.85
feature and the
blue wing of the 10
silicate feature (Keane et al. 2001). The feature
near 7.6
is well fitted by gas-phase and/or solid CH4 (Boogert
et al. 1997; Dartois et al. 1998). The spectra in Fig. 1 show
evidence for weak features between 7.2
and 7.4
m. The spectrum
of MonR2:IRS3 is particularly revealing in that it shows
a narrow absorption feature at 7.342
flanked by broader red- and
blue-shifted bands. This structure is reminiscent of the P, Q, and R branch
structure of gaseous molecules. The other sources show similar structure
albeit less pronounced due to the presence of gas-phase H2O absorption
lines (Boonman et al. 2000). A weak broad feature has been seen toward
W33A centered at 7.25
and has been attributed to solid HCOOH
(Keane et al. in prep.). However, this feature is easily distinguished from
the spectral structure observed here as it is shifted to the blue and
cannot explain the observed Q- and P-branch structure.
We attribute the spectral structure between 7.2
and 7.4
to
gas-phase
.
![]() |
Figure 2:
Synthetic gas-phase
![]() ![]() |
Open with DEXTER |
![]() |
Figure 3: The continuum divided spectra upon which the best fitting models (grey) are superimposed. Also shown are the H2O model spectra (offset) used for the modeling except in the case of NGC7538:IRS1. The position of the gas-phase CH4 band is indicated for the sources where it was included in the modeling along with the SO2 Q, R, and P branches (thick solid lines). |
Open with DEXTER |
The modeling of the spectra has been performed using synthetic spectra from
Helmich (1996) combined with the molecular line data from the HITRAN 2000
database (http://www.hitran.com). Following the same analysis as in Lahuis
& van Dishoeck (2000) and Boonman et al. (2000) a homogeneous source has
been assumed with a single temperature
and column density N.
Since the SO2 models are not sensitive to the linewidth, a Doppler b
parameter of 3 kms-1 is adopted here, corresponding to the mean value
of the submillimeter SO2 lines. Figure 2 illustrates the
expected spectral structure of the
ro-vibrational band of gas-phase
for different column densities and excitation temperatures. A global
comparison with the observations shows that the observed features imply
typically a column of a few times 1016 cm-2 of warm (
200 K)
.
Using these models, we have made detailed fits to the observed absorption
features. Half of the sources show the presence of strong gas-phase H2O
absorption in the
ro-vibrational band extending well into the
7.2-7.5
region. Therefore the H2O model fits of Boonman et al.
(2000) have been included in the modeling of gas-phase SO2. In the sources
GL2591, NGC7538:IRS1, and GL2136 gas-phase CH4 is also present and
this has been included in the models, although it affects the SO2 band only
moderately (
2%). The best fitting models have been determined using
the reduced
-method and are shown in Fig. 2. The
corresponding excitation temperatures and column densities are listed in
Table 1.
Source |
![]() |
N(1016 cm-2)
![]() |
![]() ![]() ![]() ![]() |
![]() ![]() ![]() ![]() ![]() |
|
MonR2:IRS3 | 225
![]() |
4 ![]() |
W3:IRS5 | 450
![]() |
5 ![]() |
GL2136 | 350
![]() |
6 ![]() |
GL2591 | 750
![]() |
6 ![]() |
GL4176 | 350
![]() ![]() |
4 ![]() |
NGC7538:IRS1 | 700
![]() |
4 ![]() |
a Boonman et al. (2000) unless otherwise noted; b Mitchell et al. (1990) unless otherwise noted, using 13CO and assuming | ||
12CO/13CO = 60; c Lahuis & van Dishoeck (2000); d Giannakopoulou et al. (1997). |
Molecular abundances can serve as a direct means of probing the chemical
history of star-formation. The derived
excitation temperatures
range from 200-700 K and are in good agreement with those inferred
for HCN and
(Table 1), which are good tracers of warm gas
(Lahuis & van Dishoeck 2000). The
column densities show little
variation from source to source with typical abundances of
4-
relative to the total H2. The infrared
abundances are roughly consistent with the
abundances of
10-7 observed in the submillimeter towards the Plateau, the
Compact Ridge, and the Hot Core in Orion (Sutton et al. 1995). The
relative constancy of the Orion
abundances is striking given the
physical differences that exist between the afore mentioned regions in Orion.
On the other hand, the derived
abundances in hot cores are much
higher than
abundances in dark clouds (Irvine et al. 1983).
More recently, Hatchell et al. (1998) have found gas-phase
abundances of
to
in hot core regions,
which are a factor of
10 less than what is derived here. Some of this
difference may well reflect the beam dilution suffered by the submillimeter
observations. Thus, given the Orion template, there are two possible
origins for the high abundance of gaseous
:-hot core chemistry or
shock induced chemistry.
In hot core chemistry, the
originates from oxidation of sulphur bearing
species by OH (Charnley et al. 1997). This limits the
to gas with
temperatures in the range
100 to
200 K. The
observed
temperature is well above this, though it is possible that
this high temperature reflects radiative pumping by the
dust. Moreover, HCN and C2H2 appear in gas with similar
temperatures to that of the
.
This is difficult to reconcile since
the formation of
becomes very inefficient for temperature above
230 K, whereas the route to HCN is greatly enhanced (Charnley
1997; Boonman et al. 2001; Rodgers & Charnley 2001). Thus,
cannot be abundant in gas which has become enriched in HCN through the
removal of OH. The
and HCN abundances must therefore peak at
different radii from the protostar in order for standard hot core
models to be compatible. Alternatively, the
may be formed on
grains surfaces and then be released into the gas by evaporation.
Grain surface chemistry models predict an abundance of
relative to H2O on the ice (Tielens &
Hagen 1982; Tielens private communication), which is a factor
10
less than what is derived here. Another aspect is that the
column
density does not vary between the sources, whereas H2O shows large
variation. The constancy of the
column density can be explained
by the fact that the
is only abundant in a narrow zone between
90 K (the ice evaporation temperature) and
230-300 K
(the OH
H2O transition), whose mass does not vary much
in spite of the different total masses of the sources (Doty et al., in
prep.).
The presence of SO and
has often been quoted as evidence for shocks
(Hartquist et al. 1980). The degree to which the
abundance
is enhanced depends on whether most of the sulphur is initially in atomic
form (
10-7, Pineau des Forêts et al. 1993) or locked up in
stable molecules (
10-8, Leen & Graff 1988). A good aspect
of the shock induced chemistry hypothesis is that the
and HCN may be
colocated in gas which contains freshly (i.e.,
3.
yr) released grain mantle
molecules. However, the lack of increased
line widths at
submillimeter wavelengths would seem to indicate the decay of
shock activity within the region sampled by the submillimeter
observations. In addition, the presence of fragile molecules
sensitive to destruction by shocks (e.g. H2CO; van der Tak et al. 2000b) makes shock-induced chemistry less likely as the source
of
for these sources.
In general, the gaseous
/
ratio serves
as a sensitive chemical clock for the star formation process and
searches for these molecules at high spectral resolution are needed to
help resolve the issue of the origin of the gas-phase
.
Acknowledgements
The authors are grateful to F. Helmich for setting up thesynthetic spectra and to D. Kester for insightful discussions on fringe removal from ISO data. This works was supported by the Netherlands Organization for Scientific Research (NWO) through grant 614-041-003.