Water Masers in the Circumstellar Environments of Young Stellar ...

2 downloads 0 Views 372KB Size Report
Aug 4, 1994 - continuum observations are used to locate precisely the young stellar object(s) responsible for the far- infrared emission and to investigate their ...
THE ASTRONOMICAL JOURNAL, 115 : 1599È1609, 1998 April ( 1998. The American Astronomical Society. All rights reserved. Printed in U.S.A.

WATER MASERS IN THE CIRCUMSTELLAR ENVIRONMENTS OF YOUNG STELLAR OBJECTS LEBE E S. GRISSOM MEEHAN AND BRUCE A. WILKING Department of Physics and Astronomy, University of Missouri at St. Louis, 8001 Natural Bridge Road, St. Louis, MO 63121 ; lebeem=newton.umsl.edu, brucew=newton.umsl.edu

MARK J. CLAUSSEN National Radio Astronomy Observatory, Array Operations Center, P.O. Box O, Socorro, NM 87801 ; mclausse=nrao.edu

LEE G. MUNDY Department of Astronomy, University of Maryland, College Park, MD 20742 ; lgm=astro.umd.edu

AND ALWYN WOOTTEN National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903 ; awootten=nrao.edu Received 1997 August 25 ; revised 1997 December 10

ABSTRACT We present a high-resolution radio and millimeter-wavelength study of Ðve low-mass young stellar objects with known water maser emission : RNO 15 FIR, Orion AÈW, L1157, B361, and L1251A. These objects are cold IRAS sources with far-infrared luminosities ranging from less than 6 to 40 L . Radio continuum observations are used to locate precisely the young stellar object(s) responsible for_ the farinfrared emission and to investigate their relationship to the water masers and tracers of their stellar winds. Compact radio continuum emission was detected within the IRAS error ellipse for Orion AÈW, L1157, and L1251A ; the spectral indices of their radio emission are consistent with thermal ionized winds. High-resolution VLA H O observations located pairs of masers associated with these radio 2 of 50, 235, and 150 AU, respectively, clearly placing these masers in sources within projected distances the circumstellar environments of the young stellar objects. In Orion AÈW, the strongest maser feature was used to self-calibrate the line and continuum data, resulting in the detection of a j \ 1.3 cm continuum source o†set 50 ^ 17 AU from the strongest maser. In part because of the small separations of the masers and continuum sources, none of the masers could be identiÐed with gravitationally unbound motions expected for a stellar wind origin. Key words : circumstellar matter È ISM : jets and outÑows È masers È radio continuum È stars : preÈmain-sequence 1.

INTRODUCTION

dissociative shocks (J-type shocks) ; nondissociative shock models are more efficient in exciting maser emission and better at explaining higher frequency water maser emission (Kaufman & Neufeld 1996). Alternatively, models of infalling gas onto protostars have shown that water maser emission in low- to intermediate-luminosity YSOs may be excited in accretion (Ceccarelli, Hollenbach, & Tielens 1996). Ultimately, high-resolution studies that reveal the positions and space velocities of masers relative to the YSOs can distinguish between outÑow or infall models for maser emission. For example, VLA observations of the water masers associated with the cold YSO L1448C show that the masers are well aligned with the highly collimated molecular outÑow and that their velocities are too high to arise in a Keplerian disk (Chernin 1995). OutÑow excitation is also suspected for water masers found toward two YSOs in the Ophiuchus cloud : YLW 16A and the southern component of the protoÈbinary system IRAS 16293[2422 (Terebey et al. 1992 ; Wootten 1989). While evidence for maser excitation via a disk/wind interaction or infall has not been found in low- to intermediate-luminosity YSOs, they have been suggested for higher luminosity YSOs. For example, SiO masers in Orion IRc2 have been modeled as arising in a rotating and expanding disk (Plambeck, Wright, & Carlstrom 1990). A disk origin for water masers has also been proposed for the luminous YSO IRAS 0038]6312 (Fiebig et al. 1996).

Water maser emission is a common feature of deeply embedded young stellar objects (YSOs). These YSOs have Class 0 or Class I spectral energy distributions that rise steeply at far-infrared wavelengths with little or no emission at near-infrared wavelengths, indicative of extensive circumstellar dust (Lada 1991 ; Andre 1995). A recent survey of a sample of Class 0/Class I sources identiÐed by the Infrared Astronomical Satellite (IRAS) has revealed that H O maser 2 24 of activity is common in YSOs with L \ 120 L , with _ the studyÏs original 42 objects displaying highly variable H O maser emission (Wilking et al. 1994 ; Claussen et al. 2 1996). The low to intermediate luminosity of these YSOs implies a low mass for the central object (\3 M ). For low-mass YSOs, unlike massive stars, maser _ activity appears to be concentrated within several hundred AU of the YSO (Wootten 1989 ; Terebey, Vogel, & Myers 1992 ; Chernin 1995). The presence of water maser emission from the 6 ] 5 16 23 transition at 22 GHz is well correlated with the presence of high-velocity molecular gas and appears to be another indicator of mass loss in YSOs (Felli, Palagi, & Tofani 1992 ; Wilking et al. 1994). This suggests that water masers arise in dense shock fronts created by the interaction of a strong stellar wind with circumstellar or ambient cloud material (Elitzur, Hollenbach, & McKee 1989). This maser emission may be produced by a combination of low-velocity, nondissociative shocks (C-type shocks, v \ 40 km s~1) or fast, 1599

1600

MEEHAN ET AL. TABLE 1 PROPERTIES OF IRAS SOURCES

IRAS NAMEa

OTHER NAME

R.A. (s)

03245]3002 . . . . . . 05302[0537 . . . . . . 20386]6751 . . . . . . 21106]4712 . . . . . . 22343]7501 . . . . . .

RNO 15 FIR Orion AÈW L1157 B361 L1251A

34.9 14.5 39.6 40.9 22.0

S (Jy) l

DECL. (arcsec)

12 km

25 km

60 km

100 km

DISTANCE (pc)

36 52 33 01 32

\0.25 4.3 \0.25 \0.05 5.0

3.4 19 0.29 0.95 26

47 55 11 17 66

92 82 43 38 79

300 450 450 350 200

L b (LFIR) _ \12 40 \8.5 \5.9 9.3

v (kmLSR s~1) 5.5 7.0 9.0 2.7 [4.0

a B1950.0 positions. b IRAS in-band luminosity from 7 to 135 km, computed using the technique described by Emerson 1988.

In this study, we present high-resolution observations for Ðve Class 0/Class I YSOs taken from the Claussen et al. (1996) study. These observations are part of a larger VLA study to investigate water maser emission from low- to intermediate-luminosity YSOs (Wootten, Wilking, & Claussen 1998). The positions and properties of the Ðve sources are presented in Table 1. All sources are known to display H O maser emission, and all but one is associated 2 with a molecular outÑow. Prior to our study, the locations of the YSOs inside the IRAS error ellipse were largely unknown. Our main motivation is to locate precisely the YSO responsible for the far-infrared emission by detecting the radio continuum emission and associated H O maser activity. The radio continuum observations detect2 free-free emission due to either the stellar wind or an accretion shock, thus providing the location of the YSO. In the past, single-dish studies had insufficient resolution to locate the radio sources or masers relative to the IRAS position, nearinfrared sources, or other signs of mass loss. Our ultimate goal, after locating the YSO and its surrounding features, is to investigate the origin of the water maser emission. 2.

OBSERVATIONS

High-resolution observations in the millimeter and centimeter continuum, the J \ 1 ] 0 transition of CO, and the 6 ] 5 transition of H O were made for Ðve YSOs with 16 23 O maser emission. 2 The frequencies of the continknown H 2 uum observations were 4.86, 8.46, 22.255, and 107.1 GHz, hereafter referred to as j \ 6 cm, 3.6 cm, 1.3 cm, and 2.8 mm. Centimeter continuum observations of L1251A, L1157, Orion AÈW, RNO 15 FIR, and B361, and H O 2 6 ] 5 observations of L1251A, L1157, and Orion AÈW, 16 23 were obtained using the NRAO Very Large Array (VLA).1 The millimeter continuum observations for L1157 and L1251A and observations of the CO J \ 1 ] 0 transition for L1251A were made with the BIMA array.2 The observations are summarized in Table 2. The details of the observing mode are given in columns (3)È(9), including the frequency (col. [3]), the instrument and conÐguration (col. [4]), the total bandwidth (col. [5]), the velocity resolution for spectral line mode (col. [6]), the Ðeld of view (col. [7]), and the synthesized beam size and position angle east of north (cols. [8] and [9]). The amplitude and phase calibrators and their respective Ñux densities are given in columns (10) and (11). The rms noises for continuum and ÈÈÈÈÈÈÈÈÈÈÈÈÈÈÈ 1 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. 2 The BIMA array is operated by the Berkeley-Illinois-Maryland Association under funding from the National Science Foundation.

spectral line observations are found in Tables 3 and 4, respectively. 2.1. Data Reduction Calibrations of VLA line and continuum data were performed using the CLCAL task in the AIPS3 package. Standard amplitude calibrators were used to bootstrap the Ñuxes for phase calibrators for all VLA observations (with the exception of those taken in 1994 August). An error in the observing conÐguration for the 1994 August program led to a spectrum of L1251A with only the channels redshifted from the ambient cloud velocity. Also, in lieu of an amplitude calibrator for L1251A, a 1994 May observation of the radio source 0333]321 was used with its bootstrapped Ñux of 1.25 Jy. For all continuum data, phase errors of more than 20¡ and amplitude errors of more than 20% were Ñagged and excluded from the calibration. The strong (210 Jy beam~1) water maser in Orion AÈW was used to self-calibrate both the spectral line and j \ 1.3 cm continuum data for this source. Using the CALIB task in AIPS, several iterations were performed, with the last solution solving for both amplitude and phase. As a result of applying the new solutions to the simultaneous observations of the j \ 1.3 cm continuum and water maser lines, the rms noise improved in the spectral channel map from 0.96 to 0.12 mJy beam~1, and in the continuum channel from 1.1 to 0.52 mJy beam~1. All data sets from the BIMA array were calibrated using the MIRIAD4 software. Analysis of the calibrated data was performed using AIPS software. The Ñuxes of the phase calibrators were bootstrapped from the primary calibrator W3(OH). 2.2. Positional Uncertainties Absolute positional uncertainties were estimated for each set of observations in the following way : For the VLA data, we used the larger of ^0A. 1 and ^0.1 times the synthesized beam size projected onto the right ascension and declination axes. The error of ^0A. 1 is set by uncertainties in the reference frame of the phase calibrators. An error of ^0.1 times the synthesized beam is valid if phase errors greater than 36¡ are excluded, as is the case for our continuum data. For the BIMA data, absolute positions are estimated to be ^0.2 times the synthesized beam. Absolute positional errors for source positions in right ascension and declination are presented in columns (9) and (11) of Tables 3 and 4. ÈÈÈÈÈÈÈÈÈÈÈÈÈÈÈ 3 Astronomical Image Processing Software. 4 Multichannel Image Reconstruction, Image Analysis, and Display software, developed by the Berkeley-Illinois-Maryland Association.

22.235 22.235 22.235 22.235 22.235 22.215 22.255 22.255 107.1 107.1 8.46 8.46 8.46 8.46 8.46 4.86 4.86 4.86 4.86

Orion AÈW L1251A L1157 L1157 L1251A L1251A L1157 Orion AÈW RNO 15 FIR B361 L1251A L1157 RNO 15 FIR B361

115.271

Frequency (GHz) (3)

Orion AÈW L1251A L1157 L1157 L1251A

L1251A

Source (2)

(B) (B) (B) (A) (B)

VLA (B) VLA (B) VLA (B) BIMA (C) BIMA (C) VLA (C) VLA (C) VLA (C) VLA (C) VLA (C) VLA (C) VLA (C) VLA (C) VLA (C)

VLA VLA VLA VLA VLA

BIMA (C)

Instrument (4)

50 50 50 800 800 100 100 100 100 100 100 100 100 100

3.125 3.125 3.125 3.125 1.54

24

Bandwidtha (MHz) (5)

0.66 0.66 0.66 0.33 0.33

0.51

Velocity Resolutionb (km s~1 channel~1) (6)

2 2 2 1.9 1.9 5.4 5.4 5.4 5.4 5.4 9 9 9 9

2 2 2 2 2

1.7

Field of View (arcmin) (7) 79 2 70 64 90 1 2 70 64 88 58 [75 [37 1 [4 79 [83 [57 0 [87

0.43 ] 0.30 0.40 ] 0.21 0.37 ] 0.21 0.20 ] 0.11 0.33 ] 0.23 0.43 ] 0.30 0.41 ] 0.22 0.37 ] 0.21 6.0 ] 4.1 6.0 ] 4.5 4.4 ] 2.6 3.5 ] 2.6 3.7 ] 2.6 2.8 ] 2.7 4.7 ] 2.7 8.0 ] 4.5 6.6 ] 4.6 4.9 ] 4.9 6.5 ] 4.8

P.A. (deg) (9)

5.5 ] 3.9

HPBW (arcsec2) (8)

3C 48 (1.19) 3C 286 (2.52) 3C 286 (2.52) W3(OH) (3.85) W3(OH) (3.85) 3C 48 (3.23) 3C 48 (3.23) 3C 48 (3.23) 3C 48 (3.23) 3C 48 (3.23) 3C 48 (5.49) 3C 48 (5.49) 3C 48 (5.49) 3C 48 (5.49)

3C 48 (1.19) 3C 286 (2.52) 3C 286 (2.52) 3C 48 (1.17) 0333]321 (1.25)

W3(OH) (3.85)

Flux Referencec (10)

0539[057 2229]695 2021]614 1928]738 0059]581 2229]695 2021]614 0605[085 0333]321 2200]420 2229]695 2021]614 0333]321 2200]420

0539[057 2229]695 2021]614 1928]738 0016]731

(1.20) (0.38) (2.03) (2.0) (3.6) (0.37) (3.0) (3.2) (1.5) (5.8) (0.41) (2.9) (1.7) (6.6)

(1.20) (0.38) (2.03) (4.4) (0.95)

0059]581 (3.2)

Phase Referencec (11)

a Bandwidths for continuum observations are the total for two intermediate frequencies. b The frequency resolution listed is the theoretical resolution. Because of Hanning smoothing, the VLA data taken on 1994 August 4 and 1997 April 6 and 7 have an e†ective resolution of 2.0 times the values listed. The rest of the VLA spectral line data are not Hanning smoothed, leaving an e†ective frequency resolution of 1.2 times the values listed. c Values in parentheses are Ñux densities (Jy). Phase reference sources are listed in B1950.0.

CO J \ 1 ] 0 : 1996 May 24 . . . . . . H O 6 ]5 : 2199716Apr 723. . . . . . . . 1997 Apr 6 . . . . . . . . 1997 Apr 6 . . . . . . . . 1995 Aug 21 . . . . . . 1994 Aug 4 . . . . . . . . Continuum : 1997 Apr 7 . . . . . . . . 1997 Apr 6 . . . . . . . . 1997 Apr 6 . . . . . . . . 1996 Jun 12 . . . . . . . 1996 May 25 . . . . . . 1996 Feb 26 . . . . . . 1996 Feb 26 . . . . . . 1996 Feb 26 . . . . . . 1996 Feb 26 . . . . . . 1996 Feb 26 . . . . . . 1996 Feb 26 . . . . . . 1996 Feb 26 . . . . . . 1996 Feb 26 . . . . . . 1996 Feb 26 . . . . . .

Date (1)

TABLE 2 SUMMARY OF OBSERVATIONS

1602

C1

E

D

C

C

B

C

B

A

B

A

A

A

COMPONENT (2)

107.1 107.1

8.46 4.86 8.46 4.86 8.46 4.86 8.46 4.86 8.46 4.86 8.46 4.86 8.46 4.86 8.46 4.86

8.46 4.86 22.215 8.46 22.255 8.46 4.86 8.46 4.86 22.255 8.46 4.86 8.46 4.86

36.3d \11.7

0.52 1.81 4.46 10.17 2.78 7.79 0.14 0.24 0.78 2.68 \0.08 0.21 0.19 0.21 1.39 0.80

\0.12 \0.24 2.8 0.92 \1.2 0.18 0.15 \0.06 \0.09 \1.2 0.31 0.28 0.19 \0.11

PEAKa (mJy beam~1) (4) a (6) Error (7)

SPECTRAL INDEX

[0.99

0.19

0.33

1.15

0.31

0.42

0.19

0.03 0.35 0.05

[1.86 [0.98 [2.24

1.02

[0.18

0.03

0.45

0.02

[1.49

\ [1.74

0.16

[2.26

3.9 3.9

...

BIMA Sources near IRAS Objects

0.04 0.08 0.04 0.08 0.04 0.08 0.02 0.03 0.02 0.03 0.03 0.04 0.03 0.04 0.03 0.04

VLA Sources outside IRAS Error Ellipse

0.04 0.08 0.5 0.02 0.4 0.02 0.03 0.02 0.03 0.4 0.03 0.04 0.03 0.04

VLA Sources near IRAS Objects

RMS (mJy beam~1) (5)

20 38 39.34

22 34 18.78 (1)

22 34 23.6 22 34 26.88 (8)

20 38 47.664 (6)

20 38 14.11 (4)

03 24 38.744 (1)

03 24 37.424 (1)

03 24 26.052 (7)

22 34 22.54 (8)

22 34 21.10 (5)

20 38 39.23 (3)

05 30 14.415 (2)

R.A.b (8)

R.A. (B1950.0)

0.2

0.1

0.2 0.1

0.05

0.05

0.02

0.02

0.02

0.1

0.1

0.05

0.007

Errorc (9)

67 51 34.6

75 04 10.13 (3)

75 01 45.2 75 01 34.1 (2)

67 49 05.34 (4)

67 50 27.8 (2)

30 01 11.02 (2)

30 01 23.71 (1)

30 04 32.2 (1)

75 01 37.6 (2)

75 01 34.1 (1)

67 51 35.8 (2)

[05 37 50.82 (4)

Decl.b (10)

0.8

0.3

0.5 0.3

0.3

0.3

0.3

0.3

0.3

0.3

0.3

0.3

0.1

Errorc (11)

DECL. (B1950.0)

NOTE.ÈUnits of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. a Upper limits are 3 p values. b The astrometric resolution of the last digit follows in parentheses. c Absolute error calculated from the prescription in ° 2.2. Errors in right ascension are in seconds, errors in declination in arcseconds. d Integrated emission at a level of 64 ^ 18 mJy is observed from an area of 7A. 3 ] 5A. 9, P.A. \ 105¡. REFERENCES.È(1) See also Tofani et al. 1995 ; (2) RD-D, Rosvick & Davidge 1995 ; (3) RD-A, Rosvick & Davidge 1995 ; (4) see also Schwartz et al. 1985 ; (5) see also Gueth et al. 1997.

L1157 . . . . . . . . . . . . . . L1251A . . . . . . . . . . . .

L1251A . . . . . . . . . . . .

L1157 . . . . . . . . . . . . . .

RNO 15 FIR . . . . . .

L1251A . . . . . . . . . . . .

B361 . . . . . . . . . . . . . . .

L1157 . . . . . . . . . . . . . .

Orion AÈW . . . . . . .

RNO 15 FIR . . . . . .

SOURCE (1)

FREQUENCY (GHz) (3)

VLA AND BIMA CONTINUUM RESULTS

TABLE 3

5

4

4

3

2

1

REFERENCE (12)

1603

Apr 7 Apr 7 Apr 6 Aug 21 Aug 4 Aug 4 Apr 6

1996 May 24 1996 May 24 1996 May 24

1997 1997 1997 1995 1994 1994 1997

DATE (2)

3.96 2.64 0.66 1.98 0.99 2.31 ...

... ... ...

3.70 9.64 5.0 8.67 [3.34 5.55 ...

[21.02 [20.58 [20.01

*va (km s~1) (4)

6.1 ^ 0.6 4.0 ^ 0.7 2.8 ^ 0.8

210.2 ^ 0.3 15.44 ^ 0.05 0.027 ^ 0.006 19.7 ^ 0.7 0.22 ^ 0.04 13.5 ^ 0.4 ...

PEAKb (Jy beam~1) (5)

INTEGRATEDb (Jy) (7)

0.53 0.53 0.53

BIMA CO Sources 21.2 ^ 3.9 10.3 ^ 3.3 4.9 ^ 2.6

VLA H O Sources 2 0.007 212.4 ^ 0.4 0.007 15.63 ^ 0.08 0.003 0.036 ^ 0.015 0.06 26.5 ^ 1.6 0.02 0.24 ^ 0.07 0.02 14.9 ^ 0.7 0.009 ...

RMSc (Jy beam~1) (6)

30 30 38 38 34 34

14.41206 (1) 14.41600 (2) 39.30 39.29 21.0853 (25) 21.0901 (5) ...

R.A.d (8)

R.A. (B1950.0)

22 34 21.71 (7) 22 34 22.02 (10) 22 34 21.91 (14)

05 05 20 20 22 22

NOTE.ÈUnits of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. a Full velocity width at the detection limit. b Flux and error are from AIPS task JMFIT. The error listed is that given in JMFIT based on the rms noise. c The rms noise listed is the noise found from channels free of any strong source. d Astrometric resolution of the last digit(s) follows in parentheses. e Absolute error calculated from the prescription in ° 2.2. Errors in right ascension are in seconds, errors in declination in arcseconds.

L1251A . . . . . . . . . .

L1251A . . . . . . . . . .

L1157 . . . . . . . . . . . .

Orion AÈW . . . . . .

SOURCE (1)

v (km lsrs~1) (3)

TABLE 4 VLA AND BIMA SPECTRAL LINE RESULTS

0.3 0.3 0.3

0.007 0.007 0.02 0.02 0.03 0.03 ...

Errore (9)

37 37 51 51 01 01

50.7213 (1) 50.8061 (3) 35.6 35.4 34.382 (14) 34.360 (3) ...

75 01 34.5 (2) 75 01 34.8 (3) 75 01 34.9 (4)

[05 [05 67 67 75 75

Decl.d (10)

DECL. (B1950.0)

0.8 0.8 0.8

0.1 0.1 0.1 0.1 0.1 0.1 ...

Errore (11)

1604

MEEHAN ET AL.

In cases where multiple sources were detected in a single observation, the relative positional error can be estimated from the astrometric resolution of each source. The astrometric resolution of a source is given by h /2(S/N), where h b b is the synthesized beam diameter and S/N is the signal-tonoise ratio. The relative positional error between two sources is determined largely by the astrometric resolution of the source with the lowest S/N. When relevant, the astrometric resolution is given in parentheses for the last digits in the right ascension seconds and the declination arcseconds in columns (8) and (10) of Tables 3 and 4. 3.

RESULTS AND ANALYSIS

A summary of the IRAS positions and properties for our sample is given in Table 1. Continuum detections near the IRAS position are listed in Table 3, as well as detections outside of the IRAS error ellipse but still within the j \ 3.6 cm primary beam. All continuum sources were unresolved, and hence only their peak Ñux densities are given in Table 3. For sources detected in two centimeter continuum bands, the spectral index a is calculated, assuming the emission varies as la. Spectral line detections are listed in Table 4. Each of our sources, except B361, has been imaged at j \ 2.1 km by Hodapp (1994). The following sections discuss the sources individually. 3.1. IRAS 03245]3002 (RNO 15 FIR or L 1455 IRS 1) RNO 15 FIR lies in the L1455 dark cloud. A j \ 2 km source has been found displaying positional and spectral properties that suggest it corresponds to the far-infrared source (Ladd, Lada, & Myers 1993 ; Hodapp 1994 ; Persi, Palagi, & Fall 1994 ; Tapia et al. 1997). Using these combined data, Tapia et al. derived a spectral energy distribution typical of an extreme Class I object with a dust temperature of T \ 42 K. Strong H O maser activity (as seen in single-dishD observations), [S II]2 emission knots, and a collimated molecular outÑow and jet trace the energetic wind from RNO 15 FIR (Bally et al. 1997 ; Davis et al. 1997a, 1997b ; Tapia et al. 1997). Despite its association with signs of mass loss, we detected no centimeter emission at the IRAS position. The 3 p upper limits for the emission are 0.12 and 0.24 mJy beam~1 at j \ 3.6 and 6 cm, respectively. Our j \ 6 cm upper limit is consistent with the upper limit of 0.5 mJy beam~1 established by Schwartz, Frerking, & Smith (1985) at the same wavelength. Three nonthermal sources were detected in the Ðeld, which are probably extragalactic in origin (Table 3). The j \ 3.6 and 6 cm detections of our sources B and C conÐrm the double nonthermal source seen by Schwartz et al. (1985). 3.2. IRAS 05302[0537 (Orion AÈW ) Orion AÈW is a Class I YSO associated with a molecular outÑow (Fukui 1989). Fairly weak H O maser activity, 2 been detected by which appears to be highly variable, has several groups. Felli et al. (1992), Tofani et al. (1995), and Claussen et al. (1996) have reported detections over several years of observations, with velocities of 0, 8, and 14 km s~1. Tofani et al. (1995) also have reported the detection of an extended source (0A. 53 ] 0A. 43) at the IRAS position at j \ 3.6 cm using the VLA in the A conÐguration. As seen in Table 3, we have detected an unresolved source at the IRAS position in the j \ 1.3 and 3.6 cm bands. Our j \ 3.6 cm Ñux (0.92 mJy beam~1) is very com-

Vol. 115

parable to the integrated Ñux measured by Tofani et al. The spectral index is 1.15 ^ 0.19, consistent with thermal emission. Unusually strong maser emission was observed in 1997 April with a peak Ñux of 210 Jy beam~1 from a blueshifted component at 3.7 km s~1. A second, redshifted component was observed with a peak Ñux of 15.4 Jy beam~1 at 9.6 km s~1. Within our absolute positional uncertainties (^0A. 1), the masers and continuum source are coincident and at the same position reported by Tofani et al. However, the relative positions of the j \ 1.3 cm continuum and water maser sources are much better determined and, as shown in Figure 1, have signiÐcant o†sets. In particular, the strongest maser feature is o†set from the j \ 1.3 cm continuum source at the 3 p level (0A. 04 ^ 0A. 03 west and 0A. 10 ^ 0A. 04 north) and from the weaker maser component at the greater than 100 p level (0A. 0588 ^ 0A. 0002 west and 0A. 0848 ^ 0A. 0003 north). At the distance of Orion AÈW, the strong maser is o†set by about 50 AU from the continuum source and weaker maser. Figure 1 shows the locations of the masers and the j \ 1.3 cm continuum source relative to the IRAS position. 3.3. IRAS 20386]6751 (L 1157) This deeply embedded YSO in the L1157 cloud shows several signs of mass loss consistent with being a very young, possibly Class 0, object. It has only recently displayed H O maser activity (since 1996), as reported by 2 al. (1996). It lies at the center of a highly colliClaussen et mated CO outÑow with associated clumpy SiO and molecular hydrogen emission (Zhang et al. 1995 ; Davis & EisloŽ†el 1995). The outÑow lies nearly in the plane of the sky (i D 80¡ ; Gueth, Guilloteau, & Bachiller 1996). The

FIG. 1.ÈPlot of the Orion AÈW region. The (0, 0) position corresponds to the location of the source IRAS 05302[0537 given in Table 1. Only a portion of the IRAS error ellipse, which is 26A ] 5A along P.A. \ 87¡, is shown. In addition to the IRAS source, we show the relative locations of our j \ 1.3 cm source [““cm(A)ÏÏ] and the two H O masers. The error bar representing the astrometric resolution of the 2centimeter source (0A. 04 ] 0A. 03) is shown ; comparable error bars for the masers are not shown, since they are smaller than the symbols marking their locations. The absolute positional uncertainty of the centimeter and maser source positions is ^0A. 1.

No. 4, 1998

WATER MASERS IN YSOs

blueshifted CO outÑow shows evidence for strong interaction between the stellar wind and ambient gas. CH OH 3 emission is also present in the bipolar outÑow, displaying proÐles similar to other lines tracing the shocked molecular material (Bachiller et al. 1995). Gueth et al. (1996) have reported high-resolution 12CO observations of two limbbrightened cavities, which they propose may be created by large bow shocks due to a precessing, highly collimated jet. Further evidence of stellar winds comes from observed ammonia emission, which may also trace bow shocks (Tafalla & Bachiller 1995). As seen in Table 3, we have detected a point source at j \ 3.6 and 6 cm within the IRAS error ellipse (Fig. 2). The spectral index of this radio source, 0.3 ^ 0.4, is consistent with thermal emission. Two nonthermal sources, most likely extragalactic in origin, were also observed in the Ðeld and are listed in Table 3. A resolved continuum source was detected at j \ 2.8 mm with the BIMA array (Table 3). The source is within the IRAS error ellipse and o†set at the 1.5 p level, east and south of the j \ 3.6 cm source. However, we do not regard this o†set as statistically signiÐcant. The peak j \ 2.8 mm Ñux of 36 mJy beam~1 is in excellent agreement with the j \ 2.7 mm measurement by Gueth et al. (1997) using the Institut de Radio Astronomie Millimetrique (IRAM) interferometer. The emission appears extended with an integrated Ñux of 64 mJy and a deconvolved source size of 4A. 8 ] 3A. 5 along a position angle of 146¡. The elongation of the source is in the direction of the outÑow axis (P.A. D 155¡), as noted by Gueth et al. In Figure 2, the positions of the centimeter, millimeter, and maser detections are plotted, as well as the IRAS position and error ellipse. The H (B1) position corresponds to 2

FIG. 2.ÈPlot of the L1157 region. The (0, 0) position corresponds to the location of the source IRAS 20386]6751 given in Table 1. Only a portion of the IRAS error ellipse, which is 10A ] 5A along P.A. \ 12¡, is shown. In addition to the IRAS source and error ellipse, we show the relative locations of our VLA j \ 3.6 cm source [““cm(A)ÏÏ], BIMA j \ 2.8 continuum source (““3 mm ContÏÏ), and two VLA H O masers. Also displayed are an H knot found by Davis & EisloŽ†el (1995)2 [““H (B1)ÏÏ] and the outÑow axis 2 ; Davis & EisloŽ†el 1995). as2deÐned by the H knots A1 and C1 (dashed line 2

1605

the molecular hydrogen knot B1 from Davis & EisloŽ†el (1995). The dashed line shows the position angle of the outÑow as deÐned by H knots A1 and C1 (Davis & EisloŽf2 fel 1995). The H O masers and the centimeter source appear 2 coincident to within their 2 p positional uncertainties (\235 AU). The spectral energy distribution of IRAS 20386]6751 (Fig. 3) shows that the millimeter continuum emission arises from cold dust with a temperature of D45 K. This spectral energy distribution is characteristic of extreme Class I or Class 0 sources, with no detected emission at j \ 25 km. The bolometric luminosity of the source is estimated to be 10 L . We can estimate the mass of gas and dust from the _ j \ 2.8 mm emission by assuming a constant temperature of 45 K and a dust mass opacity of i \ 0.1(0.25/j )b. mm Using the integrated Ñux, the estimates for the total mass vary from 0.5 to 1.6 M as b varies from 1.0 to 1.5. The strength of the compact_ j \ 2.8 mm emission compared with the bolometric luminosity solidiÐes the classiÐcation of L1157 as a Class 0 source. This relative strength implies a high ratio of circumstellar to stellar mass, which is a signpost of objects in the accretion phase of evolution. Following Andre (1995), we can estimate the ratio Speak (mJy)(d/160)2/L for L1157 and several prototypi2.7Class mm 0 sources. For bolL1157, this quantity is 25 mJy L ~1 cal _ compared with 17 mJy L ~1 for the combined IRAS _ 16293[2422 protobinary (Mundy et al. 1992), 24 mJy L ~1 _ for B335 (Chandler & Sargent 1993), and 100 mJy L ~1 for _ VLA 1623 (Andre, Ward-Thompson, & Barsony 1993). 3.4. IRAS 21106]4712 (B361) B361 contains a low-luminosity YSO that has displayed a single, weak H O maser (Claussen et al. 1996). This source 2 included in this paper that has no known is the only region molecular outÑow. The lack of a molecular outÑow is based

FIG. 3.ÈSpectral energy distribution for L1157. Flux densities (and their 1 p uncertainties) are plotted from the IRAS Point Source Catalog (12, 25, 60, and 100 km), IRAM j \ 1.3 mm observations by Andre (1997), our BIMA j \ 2.8 mm observations, and our radio continuum observations. The best Ðt to the data from 25 km to 2.8 mm by a singletemperature blackbody of 45 K is shown, yielding an L D 10 L . The bol _ distribution is characteristic of a Class 0 source.

1606

MEEHAN ET AL.

on the absence of broad wings in the CO J \ 1 ] 0 proÐles observed toward the IRAS position (Wouterloot & Brand 1989). However, given the high correlation between masers and molecular outÑows, it is expected that high-sensitivity CO mapping of the source will reveal an outÑow. Our observations in the radio continuum yielded upper limits for both the j \ 3.6 and 6 cm bands, which are presented in Table 3. 3.5. IRAS 22343]7501 (L 1251A) This IRAS source in the L1251A cloud is associated with a cluster of near-infrared sources spread over a 10A ] 10A area (Rosvick & Davidge 1995). Millimeter continuum emission at j \ 1.3, 1.1, and 0.85 mm has been detected from the cluster using the James Clerk Maxwell Telescope (Rosvick & Davidge 1995). Signs of stellar outÑows in this area include an optical jet (Balazs et al. 1992) and a largescale bipolar molecular CO outÑow (Sato & Fukui 1989). The source has also displayed a strong, single-velocity H O maser at 4 km s~1, as seen by Claussen et al. (1996), and2 a weak H O maser at D3 km s~1, as seen by Toth & Walms2 and Claussen et al. (1996). In this paper, we will ley (1994) refer to the Ðve near-infrared sources labeled AÈE by Rosvick & Davidge as RD-A to RD-E. Based upon its red infrared color from J (j \ 1.25 km) to K@ (j \ 2.11 km), they associate the brightest K@ source, RD-A, with the IRAS source. As shown in Figure 4, only two of the near-infrared sources were detected in our radio continuum observations. Although Rosvick & Davidge do not provide accurate absolute positions, the relative positions of their sources RD-A and RD-D are identical to our VLA sources B and A,

FIG. 4.ÈContour plot of the L1251A region taken from our VLA j \ 3.6 cm observation, with the lowest contour at a level of 2 times the rms noise (0.052 mJy beam~1) and increasing in increments of 0.052 mJy beam~1. The centroids of the j \ 3.6 cm emission peaks correspond to our point sources VLA A and B. We have used these positions to register the near-infrared map of Rosvick & Davidge (1995) with our VLA map, assuming their sources RD-D and RD-A correspond to our VLA sources A and B. The plus signs represent the locations for Ðve near-infrared sources observed by Rosvick & Davidge. VLA A (RD-D) is coincident with our VLA H O maser position. 2

Vol. 115

respectively (Table 3). Both VLA sources lie within the IRAS error ellipse (Fig. 5) and have spectral indices consistent with thermal emission. However, our observations suggest that RD-A may not be the YSO primarily responsible for the far-infrared emission. Not only is the radio continuum emission from RD-D stronger, but the water maser emission is nearly coincident with RD-D (Table 4). Given these properties and its large value of J[K, RD-D may be the major contributor to the farinfrared Ñux. Ultimately, the relative contribution of these YSOs to the IRAS Ñux densities can be determined from high-resolution mid-infrared or submillimeter imaging. Unfortunately, our BIMA observations failed to detect either source ; the 3 p upper limit to the j \ 2.8 mm Ñux densities for VLA sources A and B was 12 mJy beam~1. In addition to the VLA sources corresponding to RD-A and RD-D, we detected two radio sources, D and E, outside of the IRAS error ellipse that are possible thermal sources associated with YSOs (Table 3). While they correspond to no known near-infrared, IRAS, or Ha emission-line sources (Kun & Prusti 1993), we cannot rule out that they are associated with low-luminosity YSOs in the L1251 cloud. VLA radio source C is a nonthermal source that lies just outside of the IRAS error ellipse and is most likely a background source (Table 3). Two water masers were detected in our A-conÐguration observation and associated with VLA source A. The masers are coincident with each other to within twice the astrometric resolution of the weaker source, which is 0A. 03 or 6 AU. They are coincident with VLA source A to within their mutual 2 p uncertainties (\150 AU). High-velocity CO emission was detected with the BIMA array about midway between RD-A and RD-D. The emission is blueshifted by 16È17 km s~1 from the ambient cloud velocity (Table 4). Given the association of VLA source A/RD-D with H O 2

FIG. 5.ÈPlot of the L1251A region. The (0, 0) position corresponds to the location of the source IRAS 22343]7501 given in Table 1. In addition to the IRAS source and error ellipse, we show the relative locations of our VLA j \ 3.6 cm sources AÈD [““cm(A)ÏÏÈ““cm(D)ÏÏ], VLA H O maser, and BIMA CO sources. The sources cm(B) and cm(A) correspond2 to RD-A and RD-D, respectively.

No. 4, 1998

WATER MASERS IN YSOs

1607

masers, one is tempted to assign the high-velocity gas to it as well. The position of the CO emission is shown in Figure 5 relative to the centimeter and H O detections, as well as 2 the IRAS position and error ellipse. 4.

NATURE OF THE RADIO CONTINUUM EMISSION

Radio continuum emission is common toward low- to intermediate-luminosity YSOs and characterized by a power-law spectral distribution, S P la. Several mechal nisms have been advanced to explain this emission. The Ðrst is thermal emission from an ionized, spherically symmetric stellar wind ; such emission has a spectral index of a \ 0.6 (Panagia & Felli 1975 ; Panagia 1991). Consideration of nonspherical outÑows, adiabatic cooling, velocity gradients, and recombination can a†ect stellar wind emission, and models of such e†ects predict spectral indices ranging from a \ 0.25 to a \ 1.1 (Reynolds 1986). Stellar wind models have been criticized because they predict too much radio Ñux from sources that are extended on scales of a few arcseconds (Curiel et al. 1989). A second mechanism has been proposed in which a largely neutral stellar wind collisionally ionizes dense circumstellar material (Torrelles et al. 1985). In this case, the extended radio emission arises from an optically thin, postshock H II region with an expected spectral index of a \ [0.1. This latter mechanism is consistent with the observed correlation between the radio continuum luminosity and the CO outÑow momentum Ñux (Rodri guez et al. 1989). Shock ionization models still have problems, however. The spectral index for the radio emission from many low-luminosity YSOs is much larger than [0.1, suggesting the emission is partially optically thick. In this case, shock ionization models only work if the emission regions are extremely compact, with radii on the order of 10 AU (Villuendas et al. 1996). A third mechanism, relevant to YSOs in the earliest phase of evolution, involves ionization from extreme-ultraviolet radiation and collisions produced in an accretion shock in the inner regions of a circumstellar disk (Neufeld & Hollenbach 1996). Such emission would be compact and partially optically thick with a spectral index near 2. To attain the required degree of ionization, stellar masses must be at least 3.5È4 M , which implies that the YSOs in our sample would be_underluminous relative to their eventual mainsequence luminosities. The spectral indices for the sources associated with IRAS emissionÈOrion AÈW (A), L1157 (A), and L1251A (A and B)Èvary between a \ 0.19 and a \ 1.15, consistent with predictions of ionized stellar wind models. Yet our radio continuum data are not able to rule out other mechanisms, because the weakness of the emission does not constrain precisely the spectral index. In addition, the resolution of the VLA C conÐguration ([2A. 6) is not sufficient to measure extended emission components with angular scales of an arcsecond that are best explained by shock ionization by a stellar wind. Shock ionization may be necessary to explain the extended nature of the j \ 3.6 cm emission from Orion AÈW (Tofani et al. 1995) and extended radio continuum emission reported for other maser sources, such as IRAS 16293[2422 (Wootten 1989). A plot of the radio continuum luminosity, 4nd2S , 6 cm versus far-infrared luminosity for YSOs with water maser emission is shown in Figure 6. In ionized-wind models, the radio continuum luminosity is proportional to the wind mass-loss rate (Panagia 1991). In most cases, the IRAS

FIG. 6.ÈLinear-linear plot of the radio luminosity in units of kpc2 mJy vs. the IRAS in-band luminosity (7È135 km) in L for 15 YSOs with water maser emission. Data for sources not included in_Table 1 were taken from data compiled in Table 3 of Claussen et al. (1996). The radio luminosity at j \ 6 cm is used except for Haro 4-255 FIR and L483, where the j \ 3.6 cm Ñux density was used. For Orion AÈW and IRAS 4B, the j \ 6 cm Ñux density was extrapolated from shorter radio wavelength data. The solid line shows an unweighted linear least-squares Ðt to the data. The 1 p error bars shown assume a 10% error in the distance to the source and add this quadratically (with appropriate proportionality constants) to the uncertainties in the radio Ñux (y error) and IRAS Ñux densities (x error). A similar plot on a log-log scale, minus data for Orion AÈW, L1157, and L1251A, has been presented in Fig. 23 of Claussen et al. (1996).

in-band luminosity is used to estimate the far-infrared luminosity, which varies between about 1 and 120 L . As noted _ between by Claussen et al. (1996), there is a linear relation these two properties. With the addition of the new detections in this study, the correlation is even tighter. An unweighted linear least-squares regression to the 15 known maser sources with radio emission yields the following relation and statistical error : 4nd2(kpc)S (mJy) \ L (L )/(32 ^ 3) 6 cm FIR _ with a correlation coefficient of r \ 0.89. Two maser sources do not Ðt this relation : T Tauri, which is a complicated system with two or three components (Ray et al. 1997), and GSS 30 IRS 1, which, despite its bolometric luminosity of 25 L , has no detectable j \ 6 cm emission (Leous et al. 1991). _ 13, in NGC 1333, was also excluded from the Ðt since SSV recent observations have revealed a second VLA source that could be responsible for the Herbig-Haro and H O 2 emission (Rodri guez, Anglada, & Curiel 1997). A similar relationship between radio luminosity and bolometric luminosity has been reported for a similarly selected sample of YSOs compiled by Rodri guez et al. (1989) and Cabrit & Bertout (1992), and for a higher luminosity sample of YSOs in the L1641 cloud (Morgan, Snell, & Strom 1990). The correlation between the radio continuum and farinfrared luminosity indicates that the radio luminosity is approximately a constant fraction of the bolometric luminosity for YSOs in this phase of evolution. This relationship suggests that a common mechanism generates the radio

1608

MEEHAN ET AL.

continuum Ñux in all of these objects. As noted by Cabrit & Bertout, a wind ionized by Balmer continuum photons can account for the level of the observed radio Ñux as a function of luminosity. However, this model requires that a Class 0/ Class I source have the same fraction of Balmer continuum photons as a zero-age main-sequence star of the same luminosity and a mechanism to populate the n \ 2 level of hydrogen. 5.

NATURE OF THE H2O MASER EMISSION

Our VLA observations of Orion AÈW, L1157, and L1251A (Figs. 1, 2, 5) have revealed three compact pairs of masers. The Orion AÈW masers are clearly resolved with a projected separation of 46 AU, while the L1157 and L1251A masers are unresolved at the 2 p level (upper limits of 95 and 6 AU, respectively). All of the water masers are very near compact radio continuum sources : within 50 AU for Orion AÈW, within 235 AU for L1157, and within 150 AU for L1251A. It appears that the masers lie in the circumstellar environments of their respective YSOs. In accordance with larger scale surveys (Wilking et al. 1994 ; Claussen et al. 1996), the maser features are all low in radial velocity relative to the ambient cloud (\10 km s~1). Given the spatial resolution of these data and the positional uncertainties, we cannot determine whether the masers are aligned with larger scale outÑows. However, if we assume the compact radio continuum source marks the position of the YSO, then we can investigate whether the radial velocities of the maser features indicate that they are gravitationally unbound. Unbound motions would be expected if the dense masing gas was carried by a stellar wind. The escape velocity from a solar-mass object is given by v \ 2.0(1A/h)0.5(450 pc/d)0.5 km s~1 , esc where h is the angular separation from the YSO in arcseconds and d is the distance to the source in parsecs (Table 1). The projected angular separation observed between the maser and continuum source provides a lower limit to the true value of h. If the radial velocity of the maser relative to the ambient cloud velocity exceeds v , then the maser is esc velocities less than gravitationally unbound. Relative radial v are inconclusive mainly because the relative radial esc velocity is a lower limit to the true relative space velocity. This could be a large e†ect, since many of the low- to intermediate-luminosity maser sources (such as L1157) have outÑows that lie nearly in the plane of the sky (Wilking et al. 1994). No gravitationally unbound masers can be identiÐed from our data. Only the strongest Orion AÈW maser has a resolved separation from the continuum source (0A. 108 ^ 0A. 038). However, its relative radial velocity of 3.3 km s~1 is less than the escape velocity of 6 km s~1 at its projected distance from the continuum source. The highest velocity maser features in L1157 and L1251A could be gravitationally unbound if h is greater than 0A. 25 and 0A. 1, respectively. Unfortunately, such separations are less than or equal to the positional uncertainties for these continuum sources. Future observations that simultaneously observe the maser and continuum source (such as for Orion AÈW in

Vol. 115

this study) are critical for measuring accurately the small projected separations between masers and continuum sources. Also of great importance to investigations of the wind origin of the maser emission are multiepoch studies that allow for estimates of the masersÏ tangential velocities.

6.

CONCLUSIONS

We have conducted a continuum (j \ 1.3, cm 3.6 cm, 6 cm, and 2.8 mm) and spectral line (H O and CO) study of 2 Ðve low-mass YSOs identiÐed by the IRAS survey and with previously known H O maser activity. The major results 2 are as follows : 1. For three of the Ðve YSOs, we are able to constrain the location of the protostar responsible for the known farinfrared emission by detecting thermal radio emission inside the IRAS error ellipse. Two thermal sources were found inside the L1251A IRAS ellipse, with one source coincident with H O masers. 2 thermal radio emission, millimeter con2. In addition to tinuum emission from cold dust was detected toward IRAS 20386]6751 (L1157). The high ratio of millimeter continuum luminosity to bolometric luminosity is characteristic of a Class 0 source. 3. Five VLA radio continuum sources (AÈE) were observed toward IRAS 22343]7501, two of which correspond to near-infrared sources observed by Rosvick & Davidge (1995). Using this correspondence, we were able to register the absolute positions of the cluster members in the near-infrared image by Rosvick & Davidge (1995) and identify RD-D as the dominant luminosity source. One or, perhaps, two sources found outside of the IRAS error ellipse have thermal radio emission and are probably associated with less luminous YSOs in the cloud. 4. Radio emission associated with the IRAS sources has spectral indices between 0.19 and 1.15, consistent with a thermal ionized wind. The radio luminosity of the IRAS sources is nearly a constant fraction of the far-infrared luminosity, following a trend observed for more luminous YSOs. This relationship suggests a common radio emission mechanism for YSOs in this evolutionary phase. 5. Water masers mapped toward the IRAS positions in Orion AÈW, L1157, and L1251A originate within the circumstellar environments of these YSOs, with projected separations within 50, 235, and 150 AU of their associated compact radio source. In part because of the proximity of the masers and YSOs, no gravitationally unbound maser motions could be established. We gratefully acknowledge Philippe Andre, David Devine, and Joanne Rosvick for communicating data in advance of publication. We also thank Pat Murphy for answering all of our AIPS questions and an anonymous referee who gave us many helpful suggestions that improved the paper. L. S. G. M. and B. A. W. acknowledge support from a University of MissouriÈSt. Louis Research Award and a Univerity of Missouri Research Board Award.

No. 4, 1998

WATER MASERS IN YSOs

1609

REFERENCES Andre, P. 1995, Ap&SS, 224, 29 Morgan, J. A., Snell, R. L., & Strom, K. M. 1990, ApJ, 362, 274 Andre, P., Ward-Thompson, D., & Barsony, M. 1993, ApJ, 406, 122 Mundy, L. G., Wilking, B. A., Blake, G. A., & Sargent, A. I. 1992, ApJ, 385, Bachiller, R., Liechti, S., Walmsley, C. M., & Colomer, F. 1995, A&A, 295, 306 L51 Neufeld, D. A., & Hollenbach, D. J. 1996, ApJ, 471, L45 Balazs, L. G., EisloŽ†el, J., Holl, A., Keleman, J., & Kun, M. 1992, A&A, Panagia, N. 1991, in The Physics of Star Formation and Early Stellar 255, 281 Evolution, ed. C. J. Lada & N. D. KylaÐs (NATO ASI Ser. C, 342) Bally, J., Devine, D., Alten, V., & Sutherland, R. S. 1997, ApJ, 478, 603 (Dordrecht : Kluwer), 565 Cabrit, S., & Bertout, C. 1992, A&A, 261, 274 Panagia, N., & Felli, M. 1975, A&A, 39, 1 Ceccarelli, C., Hollenbach, D. J., & Tielens, A. G. G. M. 1996, ApJ, 471, 400 Persi, P., Palagi. F., & Felli, M. 1994, A&A, 291, 577 Chandler, C. J., & Sargent, A. I. 1993, ApJ, 414, L29 Peter, W., & Eichler, D. 1996, ApJ, 466, 840 Chernin, L. 1995, ApJ, 440, L97 Plambeck, R. L., Wright, M. C. H., & Carlstrom, J. E. 1990, ApJ, 348, L65 Claussen, M. J., Wilking, B. A., Benson, P. J., Wootten, A., Myers, P. C., & Ray, T. P., Muxlow, T. W. B., Axon, D. J., Brown, A., Corcoran, D., Dyson, Terebey, S. 1996, ApJS, 106, 111 J., & Mundt, R. 1997, Nature, 385, 415 Curiel, S., Rodri guez, L. F., Bohigas, J., Roth, M., & Canto, J. 1989, AstroReynolds, S. P. 1986, ApJ, 304, 713 phys. Lett. Commun., 27, 299 Rodri guez, L. F., Anglada, G., & Curiel, S. 1997, ApJ, 480, L125 Davis, C. J., & EisloŽ†el, J. 1995, A&A, 300, 851 Rodri guez, L. F., Myers, P. C., Cruz-Gonzalez, I., & Terebey, S. 1989, ApJ, Davis, C. J., EisloŽ†el, J., Ray, T. P., & Jenness, T. 1997a, A&A, 324, 1013 347, 461 Davis, C. J., Ray, T. P., EisloŽ†el, J., & Corcoran, D. 1997b, A&A, 324, 263 Rosvick, J. M., & Davidge, T. J. 1995, PASP, 107, 49 Elitzur, M., Hollenbach, D. J., & McKee, C. F. 1989, ApJ, 346, 983 Sato, F., & Fukui, Y. 1989, ApJ, 343, 773 Emerson, J. P. 1988, in Formation and Evolution of Low Mass Stars, ed. Schwartz, P. R., Frerking, M. A., & Smith, H. A. 1985, ApJ, 295, 89 A. K. Dupree & M. T. V. T. Lago (NATO ASI Ser. C, 241) (Dordrecht : Tafalla, M., & Bachiller, R. 1995, ApJ, 443, L37 Kluwer), 193 Tapia, M., Persi, P., Bohigas, J., & Ferrari-Toniolo, M. 1997, AJ, 113, 1769 Felli, M., Palagi, F., & Tofani, G. 1992, A&A, 255, 293 Terebey, S., Vogel, S. N., & Myers, P. C. 1992, ApJ, 390, 181 Fiebig, D., Duschl, W. J., Menten, K. M., & Tscharnuter, W. M. 1996, Tofani, G., Felli, M., Taylor, G. B., & Hunter, T. R. 1995, A&AS, 112, 299 A&A, 310, 199 Torrelles, J. M., Ho, P. T. P., Rodri guez, L. F., & Canto, J. 1985, ApJ, 288, Fukui, Y. 1989, in ESO Workshop on Low Mass Star Formation and 595 PreÈMain Sequence Objects, ed. B. Reipurth (Garching : ESO), 95 Toth, L. V., & Walmsley, C. M. 1994, Inf. Bull. Variable Stars, No. 4107 Gueth, F., Guilloteau, S., & Bachiller, R. 1996, A&A, 307, 891 Villuendas, E., Anglada, G., Estalella, R., Rodri guez, L. F., Curiel, S., & Gueth, F., Guilloteau, S., Dutrey, A., & Bachiller, R. 1997, A&A, 323, 943 Torrelles, J. M. 1996, in ASP Conf. Ser. 93, Radio Emission from the Hodapp, K.-W. 1994, ApJS, 94, 615 Stars and the Sun, ed. A. R. Taylor & J. M. Paredes (San Francisco : Kaufman, M. J., & Neufeld, D. A. 1996, ApJ, 456, 250 ASP), 56 Kun, M., & Prusti, T. 1993, A&A, 272, 235 Wilking, B. A., Claussen, M. J., Benson, P. J., Myers, P. C., Terebey, S., & Lada, C. J. 1991, in The Physics of Star Formation and Early Stellar Wootten, H. A. 1994, ApJ, 431, L119 Evolution, ed. C. J. Lada & N. D. KylaÐs (NATO ASI Ser. C, 342) Wootten, A. 1989, ApJ, 337, 858 (Dordrecht : Kluwer), 329 Wootten, A., Wilking, B., & Claussen, M. 1998, in preparation Ladd, E. F., Lada, E. A., & Myers, P. C. 1993, ApJ, 410, 168 Wouterloot, J. G. A., & Brand, J. 1989, A&AS, 80, 149 Leous, J. A., Feigelson, E. D., Andre, P., & Montmerle, T. 1991, ApJ, 379, Zhang, Q., Ho, P. T. P., Wright, M. C. H., & Wilner, D. J. 1995, ApJ, 451, 683 L71