A&A 380, L1-L4 (2001)
DOI: 10.1051/0004-6361:20011503
N. Grosso1 - J. Alves2 - R. Neuhäuser1 - T. Montmerle3
1 - Max-Planck-Institut für extraterrestrische Physik,
PO Box 1312, 85741 Garching bei München, Germany
2 - European Southern Observatory
Karl-Schwarzschild-Str. 2,
85748 Garching bei München, Germany
3 - Service d'Astrophysique,
CEA Saclay,
91191 Gif-sur-Yvette, France
Received 31 August 2001 / Accepted 9 October 2001
Abstract
We report here the discovery of a 30
-chain of embedded Herbig-Haro (HH) objects in the
Ophiuchi dark cloud. These HH objects were first detected during a deep
-band observation (completeness magnitude for point source
19) made with NTT/SOFI. We confirm their nature with follow-up observations made with H2 v=1-0 S(1) narrow-band filter.
We argue that they belong to two different jets emanating from two Class I protostars: the main component of the recently resolved subarcsecond radio binary YLW15 (also called IRS43), and IRS54. We propose also to identify the [SII] knot HH224NW1 (Gómez et al. 1998) as emanating from a counterjet of YLW15. The alignment between these HH objects and the thermal jet candidate found in YLW15 by Girart et al. (2000) implies that this jet is not precessing at least on timescale
(2-4)
104yr.
Key words: open clusters and associations:
Ophiuchi dark cloud
- infrared: stars
- infrared: ISM
- stars: pre-main sequence
- Herbig-Haro objects
- ISM: jets and outflows
In the 1950's, objective prism surveys of dark clouds revealed optical small nebulae, associated with young stars, showing emission line spectra with very weak continua. Today, these Herbig-Haro (HH) objects are interpreted as shocks produced by the interaction of outflows from Young Stellar Objects (YSOs) and the interstellar medium (see review by Reipurth & Raga 1999).
The
Ophiuchi dark cloud is one of the nearest (
pc)
active site of low-mass star formation, displaying a rich embedded
cluster of
200 YSOs (see the updated census by the ISOCAM
survey; Bontemps et al. 2001). It is thus one of the most suitable locations
for the search of HH objects. Optical emission-line surveys
(H
,
[S II]) have led to the detection of only 10 bona
fide HH objects (Reipurth & Graham 1988; Wilking et al. 1997; Gómez et al. 1998;
Reipurth 1999). As outlined by Wilking et al.
and Gómez et al., most of these HH objects are located at the periphery
of the cloud, in the lowest extinction area: optical surveys do not
have access to embedded HH objects. Wilking et al. (1997) noted for all these
nebulae a relatively high [S II]/H
ratio, which is a tracer
of low excitation conditions, and thus indicates that H2 is not
destroyed after the bow shock; they suggested to use
the shocked-H2 transition at
m to probe deeply
embedded HH objects.
Up to now in the
Ophiuchi dark cloud, shocked-H2 imaging was only performed
to study molecular shocks in the CO outflow of the Class 0 protostar prototype VLA1623
(e.g., Davis & Eisloeffel 1995), which led to the detection of several H2-knots;
one of them was later detected in [S II] (HH313; Gómez et al. 1998).
We report here deep near-IR observations of the Ophiuchi dark cloud
unveiling bow-shape structures and knots, and complementary observation
with narrow-band filter confirming their nature as H2 shocks.
As these objects are not optically visible, we called them
embedded HH objects, and discuss the possible exciting sources.
During the period April 4-7 2001, deep
-band observations
of the
Ophiuchi dark cloud were made by N. G. with NTT/SOFI, as follow-up
of Chandra and XMM-Newton X-ray observations.
The complete report of this follow-up, consisting of 5 pointings,
will be published in an upcoming paper (Grosso et al., in preparation).
We will focus here only on the pointing related to the detection of
new embedded HH objects, and located SE of the millimetric
dense core Oph-B2 (Motte et al. 1998; see Fig. 1).
This area was already surveyed by optical emission-line surveys,
but no HH objects were found.
![]() |
Figure 1:
![]() |
Open with DEXTER |
We used the auto jitter imaging mode to obtain a total integration time of 40min,
with 6s detector integration time (DIT) and 1min-exposure frames.
The near-IR standard AS29-1 (Hunt et al. 1998) was observed several times during the night.
We took five 6s-darks at the end of the night.
The different stages of the usual IR data reduction, namely sky estimation and subtraction,
frame recentering and stacking, including also dark substraction and flat-field division
(we used the flat-field provided by the NTT team for our observational period),
were performed using the jitter routine of the ESO's eclipse package
(version 3.8-1).
Zero order astrometric correction was applied using the 2MASS field stars
(Cutri et al. 2000) as reference frame, reducing the offset residual to 0.7
.
The magnitude zero-point was obtained from our calibration star.
Figure 2a is an enlargement of the observation obtained with the
broad-band filter (
m,
m),
unveiling a complex object showing bow-shape structures and possible emission knots.
To detect weak structures in the J and H-band images,
we filtered it with the wavelet transform (MR/1 package; Starck et al. 1998)
to a 5
significance threshold.
In this manner, a few of knot candidates are also detected in the H-band,
and only one is detected in the J-band.
The shape of this complex object is reminiscent of the ones found in HH objects.
If this HH object interpretation is correct, this
-band detection
must be mainly due to the 2.121
m-line of the H2 shocked emission.
To confirm it, J. A. made a 20min-exposure observation (
s, 1min-exposure frames) with NTT/SOFI using the H2 v = 1-0 S(1) narrow-band filter (
m,
m) during the night of June 8 2001 (Fig. 2b).
The flat-field was computed from the median of the dark substracted frames.
The good seeing conditions (
)
of this shorter observation gives
much more details than the
-band image.
To differentiate between pure line emission knots and continuum emission from stars
(including scattered light) we follow Davis & Eisloeffel (1995) who noted that the narrow-band filter
reduces the intensity of the field stars by the ratio of the filter bandpasses (10),
whereas the intensity of the H2 knots is only reduced by
2
(indeed spectra of H2 molecular shocks have other emission lines in the
-band;
see Smith 1995). We scale down the
-band image by the ratio of the filter
bandpasses and the ratio of the seeing to have continuum features appearing with the same
brightness in Figs. 2a and b.
On the other hand, H2-shocks appear
5 times brighter
in Fig. 2b than in Fig. 2a and can thus be easily identified.
The H2-knots detected are labelled on the H2 contour map (Fig. 2c).
Table 1 gives the position of these knots
with photometry for a 1
-aperture from the 5
-filtered images.
Before applying wavelet filtering the H2 image was convolved with a Gaussian filter
to have the same seeing as in the
-band image.
To obtain the H2 magnitudes, we scaled up the H2 intensities
by the filter-band width ratio and applied the magnitude zero-point of
the
image.
The knots A display a bow-shape structure, characteristic of low-speed bow-shocks which produce such arcs and limb-brightened structures of conical appearance (Smith 1991). We also note in the H2 image some diffuse emission downstream of the bow-shock which may be related to the excitation of H2 in the pre-shock region of a J-shock. We thus propose the knots A as the leading part of a jet coming from the SW direction, and outlined by the upstream knots B. By contrast the knots C are displaced from the direction of this jet, and moreover the knots C1-C4 are clearly elongated, roughly along a NS direction. We propose to explain these features by shocks coming from the East and producing the arc C2-C1-C5, with the downstream knots C3 and C4. We will discuss in the following section the possible exciting sources of these two jets.
Knot |
![]() |
![]() |
J | H | ![]() |
H2 |
![]() |
names | 16![]() ![]() |
-24
![]() |
[mag] | [mag] | [mag] | [mag] |
![]() |
A1 | 43
![]() |
31
![]() ![]() |
- | 23.7 | 19.3 | 17.1 | 1.4 |
A2 | 43
![]() |
31
![]() ![]() |
- | - | 20.1 | 17.9 | 1.3 |
A3 | 43
![]() |
31
![]() ![]() |
- | - | 20.3 | 18.1 | 1.2 |
A4 | 43
![]() |
31
![]() ![]() |
- | - | 20.2 | 18.1 | 1.4 |
A5 | 43
![]() |
31
![]() ![]() |
- | - | 20.9 | 18.6 | 1.1 |
B1 | 42
![]() |
31
![]() ![]() |
- | - | 19.7 | 17.6 | 1.4 |
B2 | 43
![]() |
31
![]() ![]() |
- | - | 20.0 | 17.7 | 1.2 |
B3 | 42
![]() |
32
![]() ![]() |
- | - | 19.9 | 17.8 | 1.5 |
B4 | 42
![]() |
31
![]() ![]() |
- | - | 20.0 | 17.9 | 1.4 |
B5 | 42
![]() |
32
![]() ![]() |
- | - | 20.6 | 18.4 | 1.3 |
B6 | 42
![]() |
31
![]() ![]() |
- | - | 21.8 | 18.8 | 0.6 |
B7 | 43
![]() |
31
![]() ![]() |
- | - | 21.8 | 19.0 | 0.7 |
C1 | 43
![]() |
32
![]() ![]() |
24.5 | 21.2 | 18.5 | 16.6 | 1.7 |
C2 | 43
![]() |
32
![]() ![]() |
- | - | 19.6 | 17.5 | 1.4 |
C3 | 42
![]() |
32
![]() ![]() |
- | - | 19.6 | 17.8 | 2.0 |
C4 | 42
![]() |
32
![]() ![]() |
- | 22.2 | 19.7 | 17.8 | 1.8 |
C5 | 43
![]() |
32
![]() ![]() |
- | - | 20.0 | 17.9 | 1.4 |
![]() |
Figure 2:
Deep near-IR images of the embedded HH objects discovered with NTT/SOFI.
a) ![]() ![]() ![]() |
Open with DEXTER |
To obtain an estimate of the H2-shock velocity we compare the observed near-IR
photometry with the predictions for planar molecular shocks (Smith 1995).
We update first these results for the NTT set of filters.
We extract the filter transmission profiles from the plots provided on the SOFI web
page, and convolve them to the unfiltered molecular shock line fluxes (Michael D. Smith,
private communication). Zero points for J-H and
colors are derived
by requiring the colors of A0 spectrum template (Pickles 1998) to be zero.
We observe
)/
and
2.0 for the leading part
of the two jets, which must be compared to
)/
-3.3
(resp. 1.9-3.3) for J-shocks (resp. C-shocks) with velocity range
8-22kms-1 (resp. 20-45kms-1).
This implies J-shock (resp. C-shock) velocity
10kms-1
(resp. 20kms-1), well below the H2 dissociation velocity
in molecular cloud (
22,
47kms-1 for resp. J, C-shock;
see Smith 1995).
Hence, the H2-line emission is really mapping the bow-shock;
we can exclude the scenario where this emission would come only from the low-velocity wing
of a fast dissociative bow-shock located downstream, with strong iron lines in the J and
H bands tracing the apex.
These low-velocity embedded H2-shocks are reminiscent of the low excitation HH objects
already observed in the optical in this dark cloud.
From the intrinsic color of low-velocity H2-shocks,
,
and assuming the reddening law quoted by Cohen et al. (1981), we estimate the extinction
for knot C1:
[(J-H)-(J-H)0]
.
This large visual extinction explains the non-detection by previous optical surveys.
It will be possible to measure directly
by J, H-band spectroscopy using
[Fe II] lines at 1.25
m, 1.64
m, which arise from the same atomic upper level.
The optical survey of HH objects in Taurus has shown that their frequency decreases rapidly with the age of the YSO, from Class I sources to Class II sources (Gómez et al. 1997), i.e. from evolved protostars to classical T Tauri stars. This result is consistent with the decrease of the CO outflows observed from the Class 0 sources to the Class I sources (Bontemps et al. 1996). To find the exciting source candidates of the HH objects reported here, it is thus reasonable to look for YSOs with IR excesses. On the basis of the HH feature morphology, we proposed in the previous section shocks coming from the East (resp. SW) to explain the shape of knots C (resp. A-B). We checked by constructing a color-color diagram of the sources detected in our deep observation, that there is no new embebbed YSO with IR excess in these directions. Two known YSOs are Eastward (see Fig. 1): the Class II source GY350, and the Class I protostar IRS54. We propose to associate knots C with IRS54, the agreement with the knots C shape looking better (see Fig. 2c). Strong H2 v = 1-0 S(1) line emission was detected in the IRS54 spectrum (Greene & Lada 1996), it is unresolved in our H2 image, which shows only scattered light Eastward. To our knowledge this source has never been included in a CO outflow survey.
![]() |
Figure 3:
Optical finding chart (DSS2-red) of the HH objects in the
![]() ![]() |
Open with DEXTER |
Figure 2c shows that the Class I protostar YLW15,
also called IRS43, is on-axis of the knots A-B,
10
away (0.4pc for d = 140pc).
This evolved protostar is a radio source, and was recently announced to be a
subarcsecond radio binary (Girart et al. 2000).
The main radio component, YLW15-VLA1, is spatially extended with a position
angle of 25
(one sigma error), and Girart et al. proposed it
to be a thermal radio jet candidate.
The position angle (PA) of the knot chain, 22.3
,
is compatible with the PA
of this jet candidate within errors (Fig. 3).
As the probability for a coincidence by chance between these PAs
is very low (
/180
),
the identification of YLW15 as the exciting source of knots A-B is likely.
Figure 3 displays the position of the HH objects of this dark cloud.
The [S II] knot HH224NW1 (Gómez et al. 1998) was associated with
the complex HH object HH224. We note however that HH224NW1 is displaced
from the axis of HH224, and that
the PA of YLW15 is 23.4
,
thus also compatible with the orientation of
the thermal jet candidate. HH224NW1 may then be related to a counterjet of YLW15.
This association between YLW15 and these HH objects strengthens the proposal of
Girart et al. for a thermal jet in this object.
Moreover the alignment with the HH objects moving at
-20kms-1
implies that this jet is not precessing at least on timescale
(2-4)
104yr.
These observations confirm that HH objects excited by Class I protostars may be hidden by large extinctions in the optical, but can be easily unveiled and studied by deep near-IR observations. Embedded HH objects are probably very common, and it is essential now to obtain a reliable census, for instance to quantify the role of outflows in maintaining a high level of turbulence in molecular clouds and regulating star formation (Matzner & McKee 2000).
Acknowledgements
We would like to thank the referee P.T.P. Ho for his useful comments, M.D. Smith who provided us the shocked-H2 line fluxes partly published in Smith (1995), and the NTT-team for its efficient support during the observations. N.G. is supported by the European Union (HPMF-CT-1999-00228). R.N. acknowledges financial support from the BMBF through DLR grant 50 OR 0003.