X. Chen1 - Y. Yao2 - J. Yang1, 2 - Z. Jiang1, 2 - M. Ishii3
1 - Purple Mountain Observatory, Chinese Academy of Sciences,
2 Beijing Westroad, Nanjing 210008, PR China
2 - National Astronomical Observatory, Chinese Academy of
Sciences, Beijing 100012, PR China
3 -
National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka 181, Japan
Received 10 May 2004 / Accepted 30 July 2004
Abstract
We present JHK' imaging polarimetry for the massive
star-forming region S87, and K-band spectrometry from 2.0 to
2.35 m for several infrared sources in the region. The
polarimetric patterns of S87E reveal a deeply embedded source
(DES-S87) associated with reflection nebulosities. Exciting OH
maser and 1.3 mm dust continuum emission, the DES-S87 is suggested
to be the youngest star in S87 and the driving source of bipolar
CO outflow in the region. A blue cavity is discovered in the
[J-K'] color image, to the southwest of DES. Radial
dependence of polarization degrees in the cavity could be fitted
to the slab model of outflows at an inclination of
60
out of the plane of the sky (Minchin et al. 1991b).
The nebular complex of S87E can be divided into three area. The
North nebula, exhibiting a bipolar structure, may represent a
bipolar outflow in the northwest-southeast direction. The
southwest nebula traces the blueshifted CO outflow and corresponds
to the cavity structure discovered in the [J-K'] color
image. The nebular structure of S87E suggests a quadrupolar
outflow driven by DES-S87. A well-defined centrosymmetric pattern
is found in the nebula of S87W with high degrees of polarization,
which indicates single scattering off small dust grains in the
reflection nebula. The polarimetric vectors around NIRS A show a
polarization disk feature oriented in the southeast-northwest
direction. Strong Br
emission and H2 emissions are
found in the K-band spectra, indicating the presences of strong
stellar wind, envelope/disk, shocked gas, and high-density PDRs in
the circumstellar environment of S87E and S87W. According to the
K-band spectrum, the mass of DES-S87E is estimated to be
20
.
Key words: ISM: individual: objects S87 - ISM: lines and bands - polarization - stars: formation - instrumentation: polarimeters
S87 is well-known as a site of massive star formation at a
distance of 2.1 kpc (Clemens 1985). In the optical, S87 is a faint
and diffuse nebulosity, extending over 10 arcmin. IRAS observations
suggest it to be a strong far-infrared source at all
observed wavelengths with the luminosity of
(Barsony 1989). Near-infrared imaging presents two
infrared reflection nebulae in the east and west of the region,
denoted as S87E and S87W, respectively (Chen et al. 2003,
hereafter Paper I). A young cluster was discovered within the
nebula S87E (Paper I), where active star formation is indicated by
an IRAS source, a compact HII region (Sharpless 1959), highly
variable H2O masers (Blair et al. 1980; Henkel
et al. 1986), and a massive bipolar outflow (Barsony
1989). On the other hand, the nebula S87W is relatively simple and
evolved, in which two bright sources are obvious (Paper I).
VLA maps of the S87 HII region reveal a radio continuum source
consisting of a compact core about 1'' (0.6 pc for a distance of
2.1 kpc) in size and a fan-shaped tail approximately 7''
6'' extending to the southeast (Bally & Predmore 1983;
Barsony 1989). Both Barsony and Bally & Predmore concluded that a
single B0 star could be responsible for the ionization of the HII
region. CO observations show that the HII region coincides
spatially with the peak CO emission from a 3.6 pc diameter
molecular cloud (Bally 1981). A molecular cloud core, which is
2'' in diameter and
in mass (Bally &
Lada 1983), is found in the center of the cloud. Traced by high
spatial resolution CS (2-1) and 13CO (1-0) lines, a
northeast-southwest oriented molecular outflow has been discovered
by Barsony (1989). The millimeter-wavelength data indicate that
the bipolar outflow has a wide opening angle (>60
)
and a substantial inclination angle (
50
)
(Barsony 1989).
Infrared imaging polarimetry serves as a powerful technique to study the morphology of star formation regions (e.g., Tamura et al. 1991). In infrared cluster regions, the polarization pattern can provide unique information on the structure of reflection nebulae and their relationship to the illuminating sources. Specifically, polarimetric images can be used to identify deeply embedded sources (DESs), even though they cannot be detected directly in the near-infrared (e.g., Weintraub & Kastner 1993; Yao et al. 2000), and to disentangle the geometry and structure of outflow cavities (Yamashita et al. 1987; Minchin et al. 1991a).
In this paper, we present for the first time the ' imaging
polarimetry for S87 and K-band spectroscopy from 2.0
to 2.35
m for the infrared sources in it. Our polarimetry
observations reveal a deeply embedded source associated with
infrared nebulae, and allow us to investigate the properties of
the illuminating source, nebular structure, and outflow activity
in S87. The moderate-resolution spectrometry offers physical
information about the young stellar objects (YSOs) in S87. We
describe in Sect. 2 the observations and present the results in
Sect. 3. The observational results are discussed in Sect. 4
and summarized in Sect. 5.
Observations were carried out on 1997 October 21 with the 1.88 m
telescope at Okayama Astrophysical Observatory, using the infrared
camera OASIS (Okayama Astrophysical System for Infrared Imaging
and Spectroscopy; Okumura et al. 2000). OASIS, equipped with a
NICMOS3 HgCdTe array, provides a field of view of
arcmin2 with a plate scale of 0
97 pixel-1. A
polarimeter, which consists of a rotating achromatic half-wave
plate and a fixed cold polarizer, was attached to the OASIS. Such
a configuration provides
4' field of view. The
instrumental polarization was measured to be less than 1%.
S87 was observed in the J, H, and K' (2.16 m) bands
in step-and-integrate mode. Five dithered sets of images were
obtained through each filter. Each set consists of four images
taken consecutively at four positions of the wave plate stepped by
22
5. At each position of the wave plate, integration
time of 10, 20, and 30 s were taken for the J, H, and K'bands, respectively. The images were reduced using the standard
IRAF routines with the same procedures as described in Yao et al.
(1997). The Stokes parameter images were established from the
reduced images at the four position angles. The zero position
angle of the polarization images was calibrated through
observations toward the polarized source AFGL 2591
(
J = 171
,
H = 167
,
K = 171
;
Minchin et al. 1991b). The
FWHM of the seeing disk was measured to be
1
8.
The 5
limiting magnitudes at J, H, and K' bands
are 18.0, 16.5, and 15.2 mag pixel-1, respectively.
The spectra of S87E and S87W were obtained with OASIS on 1999
November 18 and 2000 November 24, respectively. The spectral
resolution of OASIS was set to be /
500 with a 300 mm-1 grating and a 2
4 wide slit.
The wavelength coverage was from 2.0 to 2.35
m. The slit was
about 230'' long and was aligned roughly along the east-west
direction. In the observations, the slit was centered on the
K-band peak of S87E and S87W. The objects were observed 10 times
with the telescope dithered along the slit to reduce the effects
of bad pixels and cosmic rays. Each exposure time was 3 min,
and the telescope was guided during the exposure by monitoring
nearby optically visible stars with the slit-viewer of the camera.
Standard stars were also observed for atmospheric correction and
flux calibration. The seeing condition during the observations was
1
8.
Each spectrum frame was dark-subtracted and flat-fielded. The flat
field was constructed by two sets of dome flat frames taken with
an illuminating lamp on and off. The spectra were extracted using
the IRAF APALL task. The extracted spectra were then
level-adjusted to each other and median combined to produce a
final spectrum. The spectrum was further divided by A0 V standard
star (HR 9019) reduced in the same manner to remove atmospheric
absorption features; Br
absorption in the A0 V standard
had been removed with the SPLOT task after being divided by the
late-type standard (HR 9079), which shows little Br
absorption (Wallace & Hinkle 1997). The spectrum was finally
multiplied by a Planck function of an 9790 K blackbody,
representative of the A0 V star. The accuracy of flux level was
estimated to be
20% from the change of signals obtained at
different slit positions.
![]() |
Figure 1:
K' image
of S87 in an area of
![]() |
Open with DEXTER |
![]() |
Figure 2:
a) K'-band brightness contour map of S87E. The
offsets are in arcsecs from the S87E-peak
[
![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 3: a) J-band polarization map of S87E with the same axis offsets as Fig. 2a. The length of the vectors represents the percentage polarization, and the orientation of the vectors represents the position angle of polarization with respect to the equatorial coordinate system. The 100% polarization vector is shown at bottom right of the frame. The superposed contours are shown in linear scale to mark the relative position of vectors. b) The same as a, but at H. c) The same as a, but at K'. The positions of the DES are marked by plus signs in three frames. The reflection knots are denoted as A-C in Fig. 3b. |
Open with DEXTER |
Table 1: Software aperture polarimetry of the reflection knots in S87E and S87W.
Figure 1 shows the K' image of S87, in which two bright
nebulae presented in Paper I are labelled. S87E is an M-shaped
nebula associated with a young cluster (
yr; Paper I). The symmetry axis, passing through the
nebula, is in the direction from northeast to southwest. Located
1.5 arcmin to the west is the S87W nebula with a fairly
round shape. K' brightness contour maps of the two nebulae are
presented in Fig. 2. The offsets are in arcsec from the
astrometry origin (
,
)
and
(
,
), respectively, for which we
assume the K peak of S87E and S87W.
The polarization vector maps for S87E at J, H, and K' are
presented in Figs. 3a-3c, respectively. It should be noted that
the signal-to-noise (S/N) at K' is relatively low compared to
that at J and H. A centrosymmetric pattern is exhibited in the
S87E nebula at all three wavelengths, especially the in H band.
Referring to the procedure described by Weintraub & Kastner
(1993), we located a polarization centroid at
,
(B1950), with positional
uncertainties of
1'' in the north-south direction and
1
5 in the east-west direction. In the discussion
below, this polarization centroid is demonstrated to be a deeply
embedded source in S87E (Sect. 4.1.1). To the southeast and
northwest of the polarization centroid, a pair of polarized
nebular knots are found, which are labelled as "A'' and "B'' in
Fig. 3b. Simultaneously, a moderate polarized elongation
extending from the polarization centroid to the southwest is
labelled as "C'' in Fig. 3b. Table 1 gives the degree of
polarization and position angle (PA) of polarization for the
centroid and nebular knots "A-C'' in the J, H and
K' bands. Directly to the west of the polarization centroid, the
polarization vectors at K' deviate from a pure centrosymmetric
pattern, and polarization degrees are particularly low, generally
less than 10% (Fig. 3c). This disordered region includes the Kpeak of S87E and some red sources (see Fig. 2a), which may
interfere with the arrangement of K' polarimetric vectors in
this region.
Figure 4 presents the [J-K'] color image of the S87E nebula. In the image, the nebula shows different colors in the eastern and western area. The reddest color is seen in the northern part of the S87E nebula, while the bluest color is seen in the southwest part of S87E with an elongation to the southwest, corresponding to the polarized elongation found on polarization maps. This elongation also coincides with the blueshifted lobe of CO outflow. If the bluer colors imply less extinction, this elongation could suggest a cavity structure.
Figures 5a-5c display polarization vector maps for S87W in J, H and K', respectively. Offsets are in arcsec from NIRS A, which is saturated in the polarimetric images at all three wavelengths and is identifiable as a circular area of zero polarization. The polarization vectors in J and H bands roughly show a centrosymmetric pattern around NIRS A. In the K' band, for the low S/N, the alignment of polarization vector is a little irregular. Two highly polarized knots are found to the north and south of NIRS A, which are labelled as "D'' and "E'' in Fig. 5a. Both the polarization degree and position angle in J, H, K' for knots "D'' and "E'' are also listed in Table 1.
![]() |
Figure 4: J-K' color image of S87E. The J-K' color is linearly scaled from 0 mag (white) to 5 mag (black). The superposed contour map is same as the K'-contour in Fig. 2a. Deeply embedded source is labelled "DES''. The CO molecular outflow is marked by a solid arrow (blueshifted) and dashed arrow (redshited). A dashed circle delineates the HII region. |
Open with DEXTER |
Table 2: Equivalent widths and fluxes of the observed objects.
The spectra of S87E and S87W are shown in Figs. 6a and 6b,
respectively. In our observations, the K-band peaks of S87E and
S87W were centered on the slits. We also searched for emission
from nebula or nearby stars along the long slit of S87E. When such
emission was found, the spectrum was extracted in the same manner
as for the peak of S87E. The positions of these extracted sources
are given in Table 2 and marked by "NS 1-5'' in Fig. 2. Some
K-band spectral features of young stars can be found in the
spectra, such as Br
(2.166
m), H2 (2.122, 2.223,
2.248
m) and Na I (2.206
m). We regard the features as
"detected'' if they are more than 3
above the local
continuum. Equivalent widths (EW) were measured over a 0.01
m
wide interval centered at the detected lines (1
of the
equivalent widths is typically
0.5 Å; see Ishii et al.
2001). The equivalent widths and fluxes measured for the
Br
,
H2 v=1-0 S(1), and H2 v=1-0 S(0) lines are
also listed in Table 2.
It can be found in Fig. 6 that the continuum of the "S87E-peak''
is the reddest and rises towards longer wavelengths, and those of
other objects are relatively flat or blue. Both H2 1-0 S(1)
and 1-0 S(0) emission lines are detected in all spectra. In the
S87E nebula (Fig. 6a), the spectrum of the "S87E-peak'' is the
only example that does not show Br
emission and instead
shows relatively strong H2 emission. About 5
1 east of
"S87E-peak'', strong Br
emission (
Å),
H2 emissions (
Å for 1-0 S(1),
for 1-0 S(0), and
for 2-1 S(1)), and
CIV emission (
Å) are found in the spectrum
of "NS 1''. A spectrum of a star located
12
8 east was
also taken, which shows Br
(
Å) and
H2 emission (
Å for 1-0 S(1),
for 1-0 S(0), and
for 2-1 S(1)). Going
farther east, relatively weak Br
and H2 emissions are
found (NS 3). In the S87W nebula (Fig. 6b), strong Br
emission (
Å), H2 emissions (
Å for 1-0 S(1)), and Na I emission (
Å) are detected in the spectrum of the "S87W-peak''.
Through the orientation of polarimetry vectors, one polarization centroid is inferred in the polarimetric maps of S87E. However, no point sources in the near-infrared images were found to coincide with this polarization centroid. Analogous to many observations (e.g., Weintraub & Kastner 1993, 1996; Yao et al. 2000), the centroid implies a deeply embedded source in the nebula of S87E (hereafter DES-S87). A low degree of polarization is observed toward this source (see Table 1). Small but rather centrosymmetric polarimetric patterns are found to the northwest, southeast and southwest of DES-S87 (Fig. 3).
Several observations, referred to earlier, strongly support the
existence of DES-S87. The first evidence is the OH maser spot
detected in S87E (
,
,
B1950), which coincides
well with the position of DES-S87 (Braz & Epchtein 1983; Braz &
Sivagnanam 1987; Wouterloot et al. 1993). Since there is no other
infrared source closer to the maser spot, we suggest that the OH
maser is powered by the DES-S87. The second evidence is the 1.3 mm
dust continuum emission at the position of DES-S87 (Chini et al.
1986), which generally comes from a dust envelope around the YSOs.
Thirdly, the far-infrared source IRAS 19442+2428
(19
,
24
28'00'', B1950)
detected at this position provided another piece of evidence.
Additionally, the spectral feature of 3.1
m H2O ice
absorption, which is ubiquitously found in spectra of YSOs deeply
embedded in molecular clouds, was detected in S87E in the 3
m
spectra survey of YSOs by Ishii et al. (1998). Although the
aperture size of their observations was too large (5 arcsec or 9
arcsec) to decide an accurate observed position, the 3.1
m
H2O ice feature indicates the presence of deeply embedded
sources in the S87E nebulae, which is consistent with our
polarimetry result.
![]() |
Figure 5: a) J-band polarization map of S87W with the same axis offsets as Fig. 2b. The length of the vectors represents the percentage polarization, and the orientation of the vectors represents the position angle of polarization with respect to the equatorial coordinate system; The 100% polarization vector is shown at top right of the frame. The superposed contours are shown in linear scale to mark the relative position of vectors. b) The same as a, but at H. c) The same as a, but at K'. The reflection knots are denoted as "D'' and "E'' in Fig. 5a. |
Open with DEXTER |
To investigate the wavelength dependence of polarization in the
reflection nebulosities around DES-S87, we calculated the ratios
of H to J and K' to H percentage polarization
(PH/PJ and PK'/PH, respectively)
along radial cuts starting from DES-S87 as described by Minchin
et al. (1991a). The cuts begin at a radial distance 3 arcsec from
DES-S87 to avoid contamination. The ratios for each pixel along
the cuts were averaged. The results are given in Table 3. The
average values of both PH/PJ and
PK'/PH along cut 1 (
and
respectively), cut 2 (
and
respectively), and cut 3 (
and 0.72
0.16 respectively) are clearly less than unity. Thus there is a
general trend for the observed polarization around DES-S87 to
decrease with wavelength, which indicates single Rayleigh
scattering in the reflection nebulosities. At the same time, this
decreasing percent polarization as a function of wavelength
provides further evidence that DES-S87 is a thermal source. If the
DES-S87 was a purely scattering center, then the degree of
polarization should increase or remain constant with increasing
near-infrared wavelengths regardless of grain size distribution
(Pendleton et al. 1990), which is contrary to what we observe.
Table 3: The wavelength dependence of scattered radiation from DES.
![]() |
Figure 6:
K-band moderate
resolution (![]() ![]() ![]() ![]() |
Open with DEXTER |
Yao et al. (2000) have proposed that DESs could be a manifestation
in the near-infrared of the earliest phase of massive star
formation and possess characteristics similar to hot molecular
cores. The lines of evidence above indicate that DES-S87 is an
extremely young massive star, probably in the earliest stage of
protostellar evolution of massive stars. In light of these facts,
we suggest that it is the youngest source in S87. The technique of
using a polarization centroid to locate an embedded protostar has
successfully been employed by Weintraub & Kastner (1993, 1996)
and Yao et al. (2000). However, the DESs unveiled by near-infrared
polarimetric imaging need further confirmation with longer
wavelength or high-resolution (sub)millimeter-wave observations.
The DESs revealed by Weintraub et al. in near-infrared
polarimetric studies of molecular outflow regions, such as L1287
(IRAS 00338+6312; Weintraub & Kastner 1993), LkH 234
(PS1; Weintraub et al. 1994), and AFGL 437 (WK 34;
Weintraub & Kastner 1996), all have been confirmed by the
follow-up L'-band (3.8
m) imaging (Weintraub et al. 1996).
There are also increasing reports of confirmations based on
(sub)millimeter-wave observations, such as IRAS 20126+4104 (Hofner
et al. 1999) and S233B (Yao et al. 2000; Jiang et al. 2001;
Beuther et al. 2002). We believe that high resolution observations
at (sub)millimeter-wave will be helpful to further understand the
properties of the DES-S87.
The bipolar outflow centered on S87E cluster is a massive outflow,
whose blueshifted and redshifted CO flows run roughly in the
southwest and northeast directions, respectively. A close-up view
of the outflow by Barsony (1989) suggested a massive
(20
)
pre-main-sequence driving source, which is still
embedded in its parent molecular cloud. Nevertheless, the
identification of a driving source for the outflow is not
straightforward in the near-infrared images, since some bright
infrared sources are located in the center of cluster and the
resolution of CO observations are not high enough to point out the
driving source.
Weintraub & Kastner (1993, 1996) have argued that the DESs are
most likely the drivers of molecular outflows, which could be
consistent with the more recent concept that the youngest
protostars tend to drive the most energetic and most highly
collimated outflows (e.g., Richer et al. 2000). As discussed
above, DES-S87 is proposed as the youngest source in S87 powering
OH maser and 1.3 mm dust continuum emission. At the same time, the
strong Br
and H2 emissions found in the K-band
spectrum of DES-S87 implies it powers strong wind or outflow. The
mass of DES-S87 is estimated to be
20
based the
K-band spectrum (see Sect. 4.3), which is in close agreement with
the mass suggested by Barsony (1989). According to its young and
massive nature, we suggest that DES-S87 is the driving source of
bipolar CO outflow in S87.
The [J-K'] color image in Fig. 4 exhibits an elongated
cavity structure to the southwest of DES-S87. This cavity is
spatially coincident with the blueshifted lobe of the CO outflow,
which could indicate that the cavity is excavated by the blue lobe
of outflow approaching the observer and sweeping away the material
in the molecular cloud. Along the axis of the cavity to the
southwest, the polarization degrees increase with the distance
from DES-S87 (Fig. 3). Yamashita et al. (1987) and Minchin et al.
(1991b) had calculated the radial dependence of polarization for
parabolic outflows with varying inclination of the outflow axis.
Comparing our data to their models, the radial dependence of
polarization degrees in the cavity could be fitted to the slab
model of outflows at an inclination of 60
out of
the plane of the sky (Minchin et al. 1991b, Fig. 10). It is
approximately consistent with the inclination angle proposed by
Barsony (
50
;
1989).
The nebulae of S87E show complex morphology in the polarization
images. This morphology could be roughly divided into three areas
by polarimetry vectors: the first is the centrosymmetric pattern
predominating over most of the northern part of the nebulae
(hereafter referred to as the northern nebula); the second is the
bright nebula in the southwest (hereafter the southwest nebula);
the third is the disordered polarimetry area in the southeast that
coincides with the HII region in S87 (hereafter the ionized
nebula). This feature is consistent with the [J-K'] color
of S87E (Fig. 4) that the reddest color is seen around the
northern nebula, the bluest color is seen in the southwest nebula,
and an intermediate color is seen in the ionized nebula. The
polarimetric vectors in the ionized nebula is arranged
irregularly, especially at K' (Fig. 3c). In our near-infrared
images, this region is relatively weak (Figs. 1 and 2a). The
degree of polarization in the region is in the range of
.
The polarization vectors in this region have large
errors because the unpolarized intensity of the nebula is low. We
suggest that the most polarimetric vectors in this region are
spurious given the weakness of nebula.
The overall morphology of the northern nebula shows a bipolar
structure. The whole nebula is centered on DES-S87 and extends
roughly in the southeast-northwest direction. This structure is
confirmed by the [J-K'] image of S87E, in which the
northern nebula also presents a bipolar structure with DES-S87 on
the apex. Two polarized nebular knots, "A'' and "B''
(Figs. 3 and 4), have nearly the same projected distance from DES-S87
and display good alignment with DES-S87. As previously noted, strong
Br emission and H2 emissions are detected at the
position of DES-S87, which suggests strong stellar wind and
outflow activity. Furthermore, H2 emission may imply the
presence of a circumstellar envelope around DES-S87 (see below).
All the factors fit the general picture of star formation
involving a central source, a flattened disk/envelope, and a
bipolar outflow. The two polarized nebular knots ("A'' and "B'')
thus most likely indicate a pair of reflection lobes of bipolar
outflow from DES-S87. The large degrees of polarization (
)
in the lobes, as well as the similar infrared colors seen to
the east and west of DES-S87, further indicate that the axis of
the bipolar lobes is not far from the plane of the sky.
The southwest nebula traces the cavity evacuated by the blue lobe of bipolar CO outflow. Suggested by the bipolar structure of the CO outflow, it is possible that the southwest nebula is also one reflection lobe of the bipolar nebula. However, the other lobe has not been detected in our images. A scenario for the production of this feature is that the other lobe traces the redshifted lobe of CO outflow that is moving away from us and obscured by dense cloud material of the northern nebula. Thus, there are maybe two bipolar nebulae detected around DES-S87: one lies nearly in the plane of the sky with a pair of lobes projecting in the northwest-southeast direction; the other extends in the northeast-southwest direction with a pair of lobes tracing the redshifted and blueshifted CO outflow, respectively. This nebular structure suggests that DES-S87 may drive two distinct bipolar outflows nearly perpendicular to each other. The configuration of 13CO outflows in the S87 could be support for the quadrupolar outflow hypothesis. In Fig. 12 of Barsony (1989), red and blue 13CO wings expand to the northeast and southwest, respectively; furthermore, both red and blue line centers extend in the southeast-northwest direction. This configuration is in close agreement with the quadrupolar structure suggested by our polarimetric images.
Such quadrupolar (multipolar) outflow structures have been
observed in a few cases in millimeter-wavelength observations
(e.g., IRAS 20050+2720; Bachiller et al. 1995) and
near-infrared H2 observations (e.g., Hodapp & Ladd 1995;
Ladd & Hodapp 1997). Yao et al. (1998) reported a quadrupolar
outflow in S140 based on polarimetry imaging observations. It can
be found that these quadrupolar (multipolar) outflow structures
are closely associated with extremely young and deeply embedded
sources. However, the physics behind these multipolar outflows
remains unknown. Two plausible explanations have been proposed: a
single time- and angle-variable driving source has generated the
multipolar outflows, or multiple unresolved driving sources are
independently responsible for each bipolar outflow (e.g., Anglada
et al. 1996; Ladd & Hoddap 1997). In S87, there
is currently no observational evidence to exclude either
possibility. In preliminary discussions, we tend to adopt the
latter explanation for S87, because: (1) there are some red
sources in the center of bipolar CO outflow, thus we cannot
exclude the possibility that other sources, not DES-S87, drive the
bipolar CO outflow; (2) within the positional errors of the
polarimetric centroid (1
5), there are possibly multiple
embedded sources which are unresolved by our polarimetric vectors.
As suggested previously, further high resolution observations at
longer wavelengths will be helpful to clarify the properties of
S87.
The well-defined centrosymmetric pattern with high degrees of
polarization (
)
found in the S87W nebula indicates
single Rayleigh scattering off small dust grain in the nebula
(Fischer et al. 1996; Lucas & Roche 1998). The polarization
vectors in a
10'' region around the illuminating source
NIRS A form an elliptical pattern (Figs. 5a, 5b). Small degrees
of polarization are seen along the southeast-northwest direction.
Such a polarization pattern, commonly referred to as a
"polarization disk'', has been widely observed in reflection
nebulae associated with young stellar objects, such as AFGL 2591
(Tamura et al. 1991; Minchin et al. 1991b), HL Tau (Weintraub
et al. 1995) and Mon R2 (Yao et al. 1997). According to the dust
scattering models of Bastien & Ménard (1988), Whitney
& Hartmann (1993) and Fischer et al. (1994, 1996), a
polarization disk can be expected for scattered light emerging
from an optically thick disk around the illuminating source. The
elongated direction of the polarization disk indicated the
orientation of the presumed disk plane. The fluorescent 2.206
m Na I emission line found in the spectrum of NIRS A
provides further evidence for the presence of a circumstellar
disk, where the Na I region can be shielded from the direct
stellar radiation since sodium has a low first ionization
potential (5.1 eV) (McGregor et al. 1988).
The morphology of the S87W nebula is roughly round, compared to the complex morphology of S87E. It is generally recognized that nebulae evolve as the YSOs in them do, while the relationship between the evolution of nebulae and YSOs still remains unclear. Jiang et al. (2001) classified the structure of the nebulae into three groups: (1) compact, referring to those nebulae with small and roughly round geometry; (2) bipolar, referring to those showing bipolar geometry and centrosymmetric polarization patterns; (3) complex, referring to those of irregular geometry but being disentangled by polarization patterns, and proposed that while YSOs evolved, the associated nebulae evolved morphologically from complex to bipolar and then to compact. According to their definition, the structures of nebulae S87E and S87W are "complex'' and "compact'', respectively. Since the YSOs in S87W are relatively more evolved than those in S87E (Paper I), our polarimetric results confirm the conjecture of Jiang et al. (2001). Furthermore, the polarimetric imaging of the S87E nebula is also consistent with the statistical result found in Table 2 of Jiang et al. (2001) that in the cluster regions, when the YSOs are DESs, the associated nebulae all display complex structure.
The K-band spectra of S87E and S87W are presented in Figs. 6a
and 6b, respectively. The Br
and H2 emission
lines found in our spectra are frequently seen in the
near-infrared spectra of YSOs. The 2.166
m Br
line is
the only hydrogen line expected in this wavelength region. This
line emission is generally thought to occur in stellar winds close
to the stars (Ishii et al. 2001 and references therein). In the
K-band spectroscopy study of YSOs, Br
emission is
correlated with the source luminosities rather than the spectral
energy distributions (SEDs) (Persson et al. 1984; Greene & Lada
1996; Ishii et al. 2001). It indicates that the emission is
closely associated with the mass of YSOs. On the contrary, H2emission clearly depends on the SEDs from class I to class II
(Greene & Lada 1996; Ishii et al. 2001). It is rarely detected in
class III SEDs, but mostly detected in class I SEDs. The fact that
H2 emission is found for those with class I SEDs indicates
that the emission is closely associated with the circumstellar
envelopes.
It has been known that collisional excitation in shocked gas is
the dominant H2 emission mechanism in protostellar flows.
However, it should be kept in mind that UV fluorescence is also an
important mechanism of H2 excitation. In general, the line
ratio of H2 2-1 S(1)/1-0 S(1) lines is calculated to infer
the emission mechanism such that 0.1 are observed in
shocked regions and
0.5 in the photodissociation regions
(PDRs) (Shull & Bechwith 1982; Martini et al. 1999). The small
ratios of 2-1 S(1)/1-0 S(1) found in our spectra could support
excitation by shocks. However, there is another possibility, that
the H2 emission occurs in dense PDRs with densities of
105 cm-3 where the 2-1 S(1)/1-0 S(1) ratios could be
similar to those in shocked regions (Draine & Bertoldi 1996;
Luman et al. 1998). To differentiate between
excitation in shocks and in dense PDRs, the lines ratios of 2-1
S(1)/1-0 S(1) and 1-0 S(0)/1-0 S(1) are shown in Fig. 7 compared
with several theoretical models. It is noted that although the
wavelengths of these transitions are similar, the high extinction
toward S87 regions (
;
Paper I) can significantly
alter the observed line ratios from their intrinsic values. For
modification, each line ratio would be multiplied by a factor of
0.8, causing the positions in Fig. 7 to move toward the
origin. In Fig. 7, our spectral data are away from the region of
"low density PDR'' (
cm-3,
,
T0 = 500 K) and mainly distributed in the zone
between the "shock'' (
T0 = 2000 K) and "high density PDR''
(
cm-3,
,
T0 =
1000 K). Considering the small line ratio of 2-1 S(1)/1-0 S(1) in our
spectra, we suggest that the H2 emissions in S87 are
generated by the cooperative action of shocks and dense PDR. It
possibly implies that H2 emission in S87 originates in the
inner high-density regions of the shocked stellar wind. In short,
the Br
and H2 emission lines in our K-band spectra
show properties of strong stellar wind, envelope/disk, shocked
gas, and dense PDRs in the circumstellar environment of S87E and
S87W.
![]() |
Figure 7:
2-1 S(1)/1-0 S(1) vs. 1-0 S(0)/1-0 S(1) diagram for the
H2 emissions. The data (filled circles) are compared to the
theoretical predictions (open circle) of Draine &
Bertoldi (1996) for pure fluorescence in a low density PDR
(
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
No absorption feature is detected in our spectra. The absence of near-infrared absorption features is likely consistent with continuum veiling of photospheric features by the large near-infrared excess (Paper I). This excess is probably caused by warm circumstellar material and may indicate active star-forming activity. Higher resolution with high signal-to-noise spectroscopy at J or H will be helpful to further understand the properties of these objects.
It can be found in Fig. 2a that the extracted position of
"NS 1'' coincides well with the position of DES-S87. Because our
spectra were extracted from 4'' aperture along 2
4
wide slit, we consider that the spectrum of "NS 1'' covers most
emission from DES-S87. Thus, the spectrum of "NS 1'' in fact
represents the K-band spectrum of DES-S87, in which strong
Br
emission (
Å), H2 emission
(
Å for 1-0 S(1),
for 1-0
S(0), and
for 2-1 S(1)), and 2.078
m CIV
emission (
Å) are found (Fig. 6a
and Table 2). It is possibly the first K-band spectroscopy
observation of a deeply embedded source.
The 2.078 m CIV line found in the spectrum is usually
apparent in mid-O and late-O type stars (Hanson et al.
1996). Both CIV and Br
in emission without NIII and HeI
lines may imply kO7-kO8 spectral type of YSOs (cf. Table 6 of
Hanson et al. 1996). (The "k'' in this notation is an indication
that these are K-band spectral types, not optical MK spectral
types.). Spectral classes of kO7-kO8 correspond closely to MK
spectral classes of O8-O9 (Hanson et al. 1996). This range of
spectral types (O8-O9) corresponds to effective temperatures of
K and a luminosity of
assuming a ZAMS star (cf. Table 1 of Panagia
1973). On the other hand, it has been shown that there is a
correlation between the source luminosity and the luminosity of
Br
line (Persson et al. 1984; Greene & Lada 1996; Ishii
et al. 2001). The Br
flux of DES-S87 corresponds to a
luminosity of
104
(cf., Fig. 2 in Ishii et al.
2001). If an extinction of
is taken into
account (Paper I), the luminosity would be corrected to
.
Both luminosities derived from the
CIV line and Br
line are a little larger than the
bolometric luminosity of
determined from energy distribution from 0.001 to 100 mm of S87
(Barsony 1989). Barsony (1989) suggested the B0 ZAMS luminosity be
an underestimate of the true luminosity of an embedded source. In
this paper we suggest O8 ZAMS luminosity to be an upper limit of
the true luminosity of the DES-S87. Therefore, DES-S87 could be
reconciled as a deeply embedded pre-main-sequence source with a
mass of 20
.
This mass is equal to the mass suggested by
Barsony (1989) based on radio observations.
Near-infrared polarimetric imaging has been carried out for the
massive star-forming region S87, and moderate-resolution
spectrometry from 2.0 to 2.35 m has been performed for the
infrared sources in the region. The main results can be summarized
as follows:
Acknowledgements
The authors wish to thank the OAO staff for their assistance in the observations. X. P. Chen thanks Dr. Y. F. Chen for his kind help in the manuscript revision. This work is supported by NSFC grants 10133020 and 10273022, and Ministry of Science and Technology G19990754.