Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://star.arm.ac.uk/~ambn/372.pdf Äàòà èçìåíåíèÿ: Thu Dec 13 15:36:00 2001 Äàòà èíäåêñèðîâàíèÿ: Mon Oct 1 20:26:39 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: m 106

Astronomy & Astrophysics manuscript no. (DOI: will be inserted by hand later)

December 3, 2001

Active star formation in the large Bok globule CB 34
T. Khanzadyan1 , M. D. Smith1 , R. Gredel2 , T. Stanke3 , and C. J. Davis
1 2 3 4

4

Armagh Observatory, Armagh BT61 9DG, Northern Ireland, UK e-mail: mds@star.arm.ac.uk Centro Astronomico Hispano-Aleman, 04080 Almeria, Spain, e-mail: gredel@caha.es MPI fur Radioastronomie, Auf dem Hugel 69, 53121, Bonn, Germany ¨ ¨ e-mail: tstanke@mpifr-bonn.mpg.de Joint Astronomy Centre 660 N. A'Ohoku Place, University Park, Hilo, Hawaii 96720 e-mail: c.davis@jach.hawaii.edu

Received 20 June 2001 / Accepted 29 Octob er 2001

Abstract. We present near-infrared and millimetre observations of the large Bok globule CB 34. Two long parallel trails of H2 knots are discovered on wide-field images in the 1­0 S(1) 2.12 µm emission line. These parsec scale H2 jets extend to the edge of the dark globule where they disapp ear without the trace of b ow shocks. This suggests that the outflows physically extend into a lower density ambient medium where their terminating b ows are b eyond present detection limits. The two outflows are extremely well collimated and parallel to within 3 . The outflow mechanical luminosity, derived from CO measurements, and the shocked luminosity, estimated from the H2 emission, are similar, consistent with jet-driven non-evolving outflow structure. The jets app ear to originate from the densest cores, as observed in H13 CO+ line emission. A central concentration of reddened stars and a lower density halo of less reddened stars within the globule are revealed by JHK photometry. Disordered motions are observed in the CO J = 2 - 1 line velocity channel maps and can b e driven by the p ower of the outflows emanating from dense cores. We sketch a picture for the star formation history of the globule in which two star phases have b een formed. A weak diffuse emission halo is detected in the near infrared with colours consistent with either scattered light or a ro-vibrational H2 cascade. We prop ose that the halo is produced by ongoing H2 formation. Cloud evolution and halo H2 formation timescales are then b oth a few â 105 yr. Thus, we may b e witnessing the formation of a molecular cloud out of diffuse atomic gas. This supp orts a scheme in which this Bok globule has formed indep endently rather than through dislocation from a nearby molecular cloud. Key words. infrared: ISM ­ stars: formation ­ ISM: jets and outflows ­ ISM: clouds

1. Intro duction
Bok globules are isolated clouds of dense molecular gas (Bok & Reilly 1947). They are recognised in the optical as sharply outlined dark patches against a background of stars. Bok & Reilly's (1947) statement that the globules were probably sites where stars would form has been confirmed (e.g., Yun & Clemens 1990). We have begun a pro ject to study Bok globules, to determine how they
Send offprint requests to : T. Khanzadyan, e-mail: tig@star.arm.ac.uk Based on observations collected at the German-Spanish Astronomical Centre, Calar Alto, op erated by the Max-PlanckInstitut fur Astronomie, Heidelb erg, jointly with the Spanish ¨ National Commission for Astronomy. Includes observations made at the 30 m IRAM telescop e (Pico Veleta, Spain).

form, evolve and disperse. Their significance lies not only in their illustration of a distinct mechanism leading to stars but also that a considerable proportion of gas in the Galaxy is pro cessed through them. A total of 1010 M is processed in the lifetime of the Galaxy, given a globule lifetime of 106 years (Launhardt & Henning 1997), a number of 105 (Clemens et al. 1991), along with a mass of 10 M . Hence, many field stars including our Sun, as well as the enigmatic isolated T Tauri stars (Yun et al. 1997; Launhardt & Henning 1997), could find their origin in Bok globules. Globules are not a homogeneous group of ob jects. They span a range in masses from 1­100 M and their angular sizes are typically a few arcminutes which corresponds to sizes of 0.1­1 parsec, according to their distance. Therefore, a single low-mass star could form or a small


2

T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34

group including high-mass stars could result. Moreover, we cannot expect to uncover a single evolutionary scheme for all Bok globules, if any scheme exists. Bok globules could (1) form spontaneously from small growing perturbations, (2) be triggered by a large pressure wave or (3) evolve after fragmenting from a larger molecular cloud. They are not perfectly isolated. Usually, some connection with a distant molecular cloud complex is found. We present here near-infrared and millimetre data for the Bok globule CB 34 (Clemens & Barvanis 1988), in order to determine its evolutionary stage. CB 34 belongs to a class of distant massive globules. We thus wish to examine if a few powerful outflows, corresponding to those expected from high-mass stars, are being generated. The Gemini OB1 cloud is 2 east of CB 34 at an estimated distance of 1.5 kpc from us. Supposing that the two formed at the same time, with the early CB 34 being stripped off from the main outskirts of the cloud and separating at 10 km s-1 , would then yield an age of at least 5 â 106 yr. This is uncomfortably old, being equal to the average age for globule dispersal derived by Launhardt & Henning (1997). CB 34 is certainly in an advanced stage. It has successfully formed dense molecular clumps and protostars (Launhardt & Henning 1997; Huard et al. 2000) as well as being asso ciated with young stars (Alves et al. 1997). The chemical age of the globule exceeds 105 yr (Scappini et al 1998) (although this result is sensitive to the assumed density of 104 cm-3 ). Collimated outflows also testify to active star formation. A bipolar CO outflow from co ol gas (Yun & Clemens 1994), sho cked optical knots of atomic emission lines from radiative sho cks (Herbig-Haro ob jects), as well as an infrared outflow from vibrationallyexcited molecular hydrogen (Yun & Clemens 1994) have been reported. Protostellar outflows can be several parsecs long. Therefore, to search for evidence for evolved outflows, we chose CB 34 for wide-field infrared observations in the 2.12 µm H2 line. We also present JH K photometry for CB 34 and compare it to previous results of Moreira & Yun (1995). We discuss the depth in the globule and the protostellar stage of individual stars. Finally, we have obtained velo citychannel maps in CO, H13 CO+ and SiO, in order to determine the lo cation of the densest cores and to search for further outflows.

2. Observations and data reduction 2.1. Near­infrared observations
The near infrared observations of CB 34 were carried out during the nights of December 8­9, 2000 using Omega Prime (Bizenberger et al. 1998) at the Calar Alto 3.5 m telescope. Omega Prime is equipped with a Rockwell 1024 â 1024 pixel HAWAII array detector. It provides a pixel scale of 0.4 arcsec/pixel scale and a total field of view of 6.8 â 6.8 arcmin2 on the sky.

Images were obtained in the J , H and Ks broadband filters, with filter passbands centered at = 1.275, 1.645, and 2.196 µm, respectively, as well as in the narrowband H2 1­0 S(1) filter centered at = 2.125 µm and a width of 0.0206 µm. The detector integration time (DIT) was 1.677 s, and the total sum per "frame" of 47 â 1.677 s for all four filters was achieved by using the new Full MPIA read mo de (Bizenberger et al. 1998). The total integration times were 1500 s in J and H , 2000 s in Ks and finally 3800 s in the H2 narrow band. Images were obtained by positioning the globule in the center and in each of the four quadrants of the detector. The method avoids the explicit acquisition of sky frames and maximises the time spent on the ob ject ­ see Cowie et al. (1990) and Stanke (2000) for a detailed description. Difference images obtained on a flat field screen with and without illumination by a halogen lamp were used to define the pixel-to-pixel response pattern of the detector. We identified hot and cold pixels using the flatfield frames taken each night. The bad pixels were flagged and excluded from processing during the following steps of data reduction. Individual frames in each filter were passed through two stages of the DIMSUM1 package: First Pass and Mask Pass. During the First Pass sky subtraction was performed automatically by median filtering and ob ject rejection using a 5 iterated rejection mechanism (Eisenhardt et al. 1996). This functions very well under the above observational strategy. Then, by applying a threshold algorithm to the ratio of the image and median filtered image, a cosmic ray mask was created and utilised. In the Mask Pass, the images were mosaiced and a wide field image obtained by choosing the reference star on the central frame and shifting and orientating the frames towards it. The thresholding for background detection and ob ject masking were performed for the mosaiced picture. Final mosaics of 8.8 â 8.8 arcmin2 are obtained after the 2 runs. The results for J , H , Ks and H2 are displayed in Figs. 14a-d, respectively. The astrometric calibration was performed by use of stars from the Hubble Space Telescope guide star catalog (GSC) and the SIMBAD2 data base which were identified on our mosaics. We compare our images with those from 2MASS3 . There is a go o d positional accuracy of the order of 1 , which is similar to the positional errors in the SIMBAD and GSC data bases. In the final stage of data reduction, the H2 image was continuum subtracted. We to ok H2 (sub) = (H2 â C ) ­ Ks . Where the constant C was determined from the width of the Ks and H2 filters and by eye examination of
DIMSUM: Deep Infrared Mosaicing Software, Um ... IRAF external package develop ed by Eisenhardt et al. (1996). 2 SIMBAD is online database, op erated at CDS, Strasb ourg, France. 3 2MASS two micron all sky survey pro ject is collab oration b etween The University of Massachusetts and the Infrared Processing and Analysis Center (JPL/Caltech).
1


T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34

3

02

CB34FU

N1 N2

SMM3

N3 G
HH291Y HH290 & HH290N1

01

T Q5
Q2 Q1 SMM4 Q HH291X

N5

N4
Q3 Q4

CB034YC1 SMM2 CB34 VLA1 SMM1

N7
Declination (ddd:mm:ss) 21:00
SMM5

N6
HH290N2

H1
HH291

H2 H3 N8
59

CB034YC1-I

H4
HH290IRS HH290S

58

H5 H6
15 10 05 Right ascension (hh:mm:ss.s) 5:47:00 46:55

Fig. 1. The central part of CB 34 in H2 . Total integration et al. (2000) and the "â" are HH ob jects and H2 knots "diamond" (Yun et al. 1996). The YSOs NIR counterparts (1994). Our discovered ob jects are indicated as chains of found ­ Q5.

of this frame is 3800 s. Crosses "+" are SMM ob jects from Huard discovered by Moreira & Yun (1995). CB34 VLA1 is marked by a CB034YC1 and CB034YC1-I ("triangles") are from Yun & Clemens H and N letters as well as G and T. One more Q knot has b een

the subtracted image (in order to subtract the ma jority of field stars). From the subtracted image, H2 excitation knots were revealed. Finally, an H2 excitation mask was made (defined as 1 where H2 excitation exist and 0 where not) which was multipled by the actual H2 image in order to overlay it on the unsubtracted H2 image. This distinguishes the H2 excitation ob jects in Fig. 1.

IRAM 30m telescope on Pico Veleta, Spain, on April 28th, 29th and 30th of 2001. The observations were done using 4 receivers simultaneously, with the backends configured to observe the following lines: the 12 CO(2­1) line (at 230.538 GHz) with the A230 receiver and the autocorrelator as backend with a frequency resolution of 320 kHz (velocity resolution of 0.42 km s-1 ); the 13 CO(2­1) line (at 220.399 GHz) with the B230 receiver and the autocorrelator (same resolution as for the 12 CO(2­1) line); the 12 CO(1­0) line (at 115.271 GHz) with the B100 receiver and the auto correlator (frequency resolution 320 kHz corresponding to

2.2. H13 CO+ and

12

CO (2-1) observations

We obtained millimetre line maps in several molecular tracers of the central regions of CB 34 using the


4

T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34 Table 1. New ob ject p ositions and fluxes.
RA (2000.0) 05h 05h 05h 05h 05h 05h 05h 05h 05h 05h 05h 05h 05h 05h 47m 47m 47m 47m 47m 47m 47m 47m 47m 47m 47m 47m 47m 47m 00s . 03s . 04s . 04s . 10s . 10s . 00s . 01s . 02s . 05s . 05s . 05s . 06s . 09s . 9 5 5 4 1 7 5 4 4 0 2 6 3 1 Decl. (2000.0) 21 20 20 20 20 20 21 21 21 21 21 21 21 20 00 59 59 59 57 57 01 01 01 00 00 00 00 59 08 41 26 23 46 25 35 21 12 35 30 16 06 17 H2 v = 1 - 0 S(1) (â 10-18 ) (Wm-2 arcsec-2 ) 2. 3. 1. 1. 6. 2. 1. 1. 1. 2. 2. 3. 0. 3. 45 29 85 56 93 86 92 21 11 84 08 06 76 58

0.83 km s-1 ); the SiO(2­1 v = 0) line (at 86.847 GHz) with the A100 receiver and the 1 MHz backend (velo city resolution 3.5 km s-1 ); and finally the H13 CO+ (1­0) line (at 86.754 GHz) also with the A100 receiver and the 100 kHz backend (velo city resolution 0.35 km s-1 ). The final maps consist of a number of smaller raster maps. Of those, the first map centred roughly on the north-eastern core containing SMM4 was taken on a 4 grid covering 40 â 40 (i.e., fully sampled even for the high-frequency CO(2­1) lines; the beamsize at 230 GHz is about 11 ), the remaining maps were taken on 8 or 10 grids, i.e., fully sampled for the low-frequency data (the beamsize is about 27 at 86 GHz and 22 at 115 GHz), but undersampling the high-frequency data. The data were flux calibrated using the standard hot/cold load calibration pro cedure at the 30 m. However, it should be noted that part of the data were taken in cloudy or foggy weather, thus the calibration (particularly for the high frequency CO(2­1) data) might be severely inaccurate for some of the data. We will therefore only use general morphological information rather than interpret the brightness of the features found in the data. We present here only velo city channel maps of the 12 CO(2­1) line and the H13 CO+ line. No emission was detected in the SiO line (a sho ck tracer, e.g., Schilke et al.1997); apparently the sho cks in CB 34 are not strong enough and to o far away to pro duce detectable SiO emission.

Ob ject H1 H2 H3 H4 H5 H6 N1 N2 N3 N4 N5 N6 N7 N8 G T Q5

05h 47m 05s 1 21 00 44 . . 05h 47m 05s 7 21 00 49 . 05h 47m 07s 9 21 00 43

3.78 0.85 1.69

3. Results 3.1. Newly discovered objects
The central part of CB 34 in the continuum subtracted (1, 0) S(1) line of molecular hydrogen is shown in Fig. 1. The lowest H2 contours are two standard deviations above the noise in the sky level. The marked ob jects are taken from Huard et al. (2000; crosses "+"), Moreira & Yun (1995; "â"s), Yun et al. (1996; "diamonds") and Yun & Clemens (1994; "triangles") which are overlaid using the SIMBAD online data base. Newly found H2 ob jects are marked with the letters H1­H6, N1­N6, G, T and Q5. The ob jects Q1­Q4 are H 2 emission ob jects discovered by Moreira & Yun (1995). Our knot Q5 appears to be asso ciated with Q1­Q4. An estimate of the surface brightness averaged over a 2 aperture of the newly discovered emission ob jects is given in Col. 4 of Table 1. Columns 2 and 3 of Table 1 contain the right ascension and declination of the ob jects.

3.2. Infrared photometry
The sky conditions during the observations were less than optimal and non-photometric in particular, which do es not allow us to derive independent absolute photometry. Here, we estimate fluxes of the emission ob jects in CB 34 from a

relative photometry obtained from a comparison with the same knots which were found by Alves & Yun (1995). For the JH Ks photometry, we use the IRAF DAOPHOT routine (Stetson 1987) to identify stars which are simultaneously present in the J , H , and Ks frames and to obtain instrumental stellar magnitudes. Then, instrumental magnitudes were compared with the same objects' absolute magnitudes obtained from the work of Alves & Yun (1995) in order to detect magnitude transformation constants. Then, using the same constants for all our detected stars, the instrumental magnitudes were transformed to absolute magnitudes. This transformation metho d yields errors of no more than 15% in comparison with the results from Alves & Yun (1995), although in some cases the values are identical (sf. ob jects A, C and D). Figure 2 identifies the 14 most reddened ob jects found in this way, labelled with the same letters A­S as used by Alves & Yun 1995. These ob jects are represented by filled circles in the (J - H ) and (H - K ) colour-colour diagrams of Fig. 3. The solid line in the figure represents the lo cation of the un-reddened main­sequence stars (Ko ornneef 1983). The two parallel dashed lines are the reddening vectors which define the reddening band for normal stellar photospheres. Ob jects with colours that fall outside and to the right of this band are sources with intrinsic infrared excess emission (Lada & Adams 1992). Figure 4 shows the spatial distribution of all sources detected simultaneously in the J , H and Ks images. Those with (H - K ) > 1.0 are represented by filled circles in Figs. 2, 4 and 5. Crosses mark ob jects with (H - K ) 1.0. From Figs. 4 and 5 one can derive a core radius of the Bok globule of approximately 50 ­60 which contains all those ob jects with (H - K ) > 1.6.


T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34
02:00

5

30

01:00

M K

30

Q E F C B D N

A

21:00:00 Declination

P
30

G

S

59:00

30

58:00

57:30

10

05

5:47:00

46:55

Right ascension

Fig. 2. Isophotes of the central region of CB 34 in Ks covering an area of ab out 2.5 â 2.5 arcmin2 . East is left, north is up. Letters A, B, C, D, E, F, G, K, L, M, N, P, Q and S are the 14 reddest ob jects. Notice Ob ject S which has large [H­ K] and is located at the edge of the cloud, suggesting that it is a background star.

A number of ob jects with large (H - K ) are lo cated at the edge of the cloud (cf. Fig. 3 in comparison with Figs. 2, 4 and 5). Those ob jects (sf. S , N and K ) are probably background stars. Because of the large distance of 1500 pc to the globule (Carpenter et al. 1995), 1­2 of the red stars lo cated within the globule may be red foreground ob jects and not embedded. Ob jects C and G are lo cated in the region of the [(J - H )vs.(H - K )]-plane which has been explained by disk emission (Strom et al. 1989) and are class II ob jects (Adams et al. 1987). Ob jects A, B, D, E, and L have large (H - K ) and are lo cated in the region of the diagram characterized by large extinction corresponding mostly to class I sources (Lada & Adams 1992). The near­infrared colours of the objects A, B, C, D, E, G and L give us confidence that they are young stellar ob jects (YSOs), surrounded by varying amounts of circumstellar material, and that they formed recently as a part of ongoing star formation in the CB 34 Bok globule. This is in go o d agreement with results from Alves & Yun (1995).

Fig. 3. Near­infrared colour­colour diagram for all sources detected in the J , H and Ks bands simultaneously. The solid line is the location of the un-reddened main­sequence stars (Koornneef 1983). Dashed lines are two parallels to the reddening vector which define the reddening band for normal stellar photospheres. Ob jects with (H - K ) > 1.0 are marked with filled circles (see Figs. 2, 4 and 5). Crosses are ob jects with (H - K ) 1.0. Letters represent the same ob jects as in Fig. 2. The size of the error bars is shown in the lower right­hand corner.

4. Discussion 4.1. The outflows
We identify three powerful near-infrared H2 outflows which we will call "H2 jets", as shown in Fig. 6. These are designated the H, N and Q jets. The Q H2 jet has been previously reported in the NIR by Moreira & Yun (1995).

Note the proximity of optical HH ob jects to the H jet. Note also that there are several knots of H2 emission not directly identified with these jets (e.g. T and G Fig. 15) which may belong to further jets. Here, however, we remain conservative, requiring the presence of at least four aligned knots before we declare an H2 jet. The 2.12 µm emission line is from the 1­0 S(1) transition of molecular hydrogen. This arises from vibrationally excited H2 , probably with a temperature in the range 1000­3000 K, excited in the remote lo cations by sho ck waves. No lo cal source of UV radiation to suggest fluorescent excitation is present near the ma jority of H2 knots. The shocks are being driven, however, from distant protostars through jets or accumulated gas in bullets. The H and N jets have required wide-field cameras to be revealed. They alter our view of the evolution of the Bok globule. These two jets are considerably longer than the Q jet. We assume here a distance to CB 34 of 1500 pc (Carpenter et al. 1995) to yield pro jected sizes of 1.59 pc (219 ) for H1­H6, 1.39 pc (191 ) N1­N8 and 1.17 pc (140 ) Q5-Q4-HH290N1. We find no evidence for further extensions to these jets out to a scale of 5 .


6

T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34

Fig. 4. Plot of the [H - K ] colour as a function of pro jected angular distance to an estimated centre p oint ob ject­E, which was chosen by visual examination of the cloud (s.f. Figs. 1 and 2. The symb ols are the same as describ ed in Fig. 3.

Fig. 5. Spatial distribution of all the sources seen in the 3 bands simultaneously. The offset centre is the ob ject E.

Supposing a total expansion speed of 100 km s-1 yields dynamical ages of just 1.2­1.4â104 yr. This would indicate that the driving protostars are in the late Class 0 regime. This age is misleading, however, since the outer parts of the outflows decelerate in the Class 1 stage. Class 1 near infrared outflows are not significantly longer than Class 0 outflows (Stanke 2000), consistent with mo dels which account for outflow deceleration (Smith 2000). One of the most striking features is the high collimation of the H2 jets. Angles subtended by individual components, as well as the overall collimation, of the H and N jets are less than three degrees. We define a collimation angle which is independent of the lo cation of the driving source. This is: collimation angle = arctan(maximum transverse extent/total linear extent). This yields collimation angles of 2.7 (H-flow) and 2.3 (Nflow). Recent three dimensional simulations of hydro dynamic molecular jets show very high collimation, with jet material being fo cussed while ambient gas is swept aside (Rosen & Smith 2001). The energy of the CO J = 1 - 0 outflow, scaled to the distance of 1500 pc, is 1.8â1045 erg (Yun & Clemens 1994) and the CO momentum is 26.3 M km s-1 . Yun & Clemens (1994) derived a mechanical luminosity Lmech = 0.37 L and a dynamical age of just 3.8â104 years for the CO. This, however, assumed the existence of a high velocity wind rather than the measured CO speed of 5km s-1 . The latter speed, as usually assumed to calculate CO dy-

namical ages, is confirmed here in our CO J = 2 - 1 velo city channel maps (see below). The dynamical age is then 7.6â105 years. The mechanical luminosity is reduced to Lmech = 0.019 L and the momentum flow rate becomes FCO = 3.5 â 10-5 M km s-1 yr-1 . This is a mo derately high momentum flux, and would typically be driven from a Class 1 protostar of bolometric luminosity 10 L or a Class 0 protostar of luminosity 1 L (e.g. Smith 2000). The submillimetre cores could certainly contain either object, since their bolometric luminosities are 40 L (SMM1) and 20 L (SMM4). The total luminosities in the H2 1­0 S(1) line are 2.2 â 10-4 L (H flow) and 1.9 â 10-4 L (N flow), on summing up the individual knot luminosities. Applying factors based on general models, these values convert to total sho cked H2 radiation of 2-5 â 10-3 L per outflow, for typical C-type or J-type sho cks (Smith 1995). We then estimate that the total sho ck power from all molecular and neutral atomic species (CO, H2 O, O[I], etc.) is 2 - 5 â 10-2 L per outflow (Smith 1991). We assume here that the outflows are predominantly located on the edges of the globule where there will be no significant infrared extinction. Thus, the mechanical and sho cked luminosities are comparable given an outflow age of 7.6â105 yr. Such a relationship is expected when an outflow is jet-driven for long extended perio ds. At this age, however, we could expect that the outflow has had time to amass a reservoir of CO gas whereas only feeble sho cked H2 gas would be detected as the jets' power diminishes. Hence, we suggest that the present sho cked H2 is being generated from a younger outflow while previous


T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34
02

7

N1
SMM3

N2 N3

01

Q5
Q1 SMM4
Declination (ddd:mm:ss)

T
Q2

G
Q4

CB034YC1 SMM2 SMM1

N5 N4 N6 N7
HH290IRS

H1

21:00

SMM5

HH291

H2 H4

H3 N8
59

HH290S

58

H5 H6
15 57 10 05 Right ascension (hh:mm:ss.s) 5:47:00 46:55

Fig. 6. Schematic map of the H, N and Q near-infrared H2 outflows in CB 34. The squares are H2 knots discovered here. All other symb ols are the same as in Fig. 1.

bipolar outflows dominate the CO emission. Such a conclusion would apply to any region where there are multiple outflows. We remark that the age estimated from the CO outflow would be sufficient for jets to penetrate out to many tens of parsecs, given a propagation speed of 40 km s-1 . This suggests that molecular flows reach a limiting size of 1­5 parsecs during the Class 0/1 phase, and then shrink in apparent size as the largest scale structure becomes disconnected from the source and rapidly fades. Then, during the Class 2 phase the outflows are further reduced to of order 0.1 pc, as measured (Hogerheijde et al. 1998). Evidence is presented below that the observed infrared outflow is indeed of a similar size to the globule. The new CO data displayed in Fig. 7 provide information on the co ol molecular gas which has been disturbed or swept up by the protostellar outflows. Each map covers an interval of 1 km s-1 around the inner 2 . The overall structure of the high velo city lobes is very similar to the maps presented by Yun & Clemens (1994) in the 12 CO J = 1 - 0 line. However, the south-western lobe of high velo city gas (which is visible both in redshifted and blueshifted emission at relatively low velo cities) displays a clear bipolar morphology at higher velo cities, with a redshifted lobe extending to the south-east and a blueshifted lobe to the north-west of the position of SMM1. This bipolar flow seems to be fairly well collimated. Regarding the high-

velocity gas in the north-eastern part of the map, mainly red-shifted emission is detected with this emission showing considerable substructure. Blue-shifted emission is not or only marginally detected; better data might be useful to further constrain the presence and distribution of highvelo city gas in the north-eastern part of the cloud.

4.2. Individual H2 structures
As opposed to aligned knots, a true jet-type structure N6 is detected within the SMM4 core (Fig. 15). The compact knot located at the northern edge of N6 suggests that the source of the N jet is located within this core. In fact, the Q and N jets appear to cross at the position of SMM4. We note that other knots are extended in the direction of the N6 jet, including H5 (Fig. 6). The extension of H5 towards H6 is evidence that these outer knots, aligned with the inner knots H1­H4, are indeed part of the H outflow rather than a chance conjunction. Knot H6 (Fig. 16a) is at the end of the observed H jet but does not possess a bow or arc structure. It would appear to be a structure within the supposed jet or along a cavity wall, with the jet actually extending beyond H6, out of the globule, although no further infrared structure south of the globule edge is found (see Figs. 16a, 14d and 1). Similarly, the knots extending to the north of the N jet possess no bow


8

T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34

Fig. 7. 12 CO(2­1) velocity channel maps of CB 34. The data have b een resampled from the original velocity resolution of 0.42 km s-1 to 1.0 km s-1 . The central velocity of each channel is indicated. The pixel size is 4 , the b eam size at 230 GHz is ab out 11 (note however that only a 40 â 40 area centered at the small cross has spatially fully sampled data; in the remaining area, the data are undersampled).

or cometary appearance, again suggesting that this jet has also exited from the globule. These structures suggest an explanation for the strong asymmetry in the jet extensions from the probable sources. If the cloud were spheroidal with the N outflow near the rear surface, then N8 would demarcate the location where the jet exits from the rear surface of the globule, while N1 (Fig. 16b) would travel further through the globule to reach and exit on the near surface. On the other hand, if the SMM1 core is nearer to the front of the globule, then the receding southern H jet (as suggested by the CO lobes) could exit from near the middle section of the globule (see Fig. 1), pro ducing a longer jet as it cuts its way through more globule material. The forwardmoving northern part of the jet would then exit in a short length from the front of the cloud. Hence, the asymme-

tries can be pro duced by the source lo cations within the globule. Against this picture is the lo cation of the optical knot HH 290S which is located in the receding N-flow (see Fig. 6), thus deep within the globule. The optical extinction is 5­15 mag. for the embedded stars within the inner arcminute. HH 290S is, however, 85 south of the central star E (Fig. 2, and lies in a region associated with less reddened stars (Fig. 2).

4.3. The protostars and cores
Matching protostars to outflows and jets is not unique without spectroscopic information or proper motions. Furthermore, at the distance of CB34, each submillimetre source could consist of a group of protostars, as explicitly demonstrated by Huard et al. (2000).


T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34

9

a distance of 100 from the globule centre, (ob ject E). Slightly reddened stars (1.0 > H - K > 0.7) extend out to 200 . Furthermore, Fig. 4 demonstrates that the number of embedded stars (H - K > 0.7) per unit radius is roughly constant, with 18, 19, 17 and 12 ob jects lying within the first four 50 radius intervals. Hence the observed number per unit surface area is inversely proportional to the radial distance. This indicates that the reddened stars are internal to the globule, not background. The 200 radius halo of stars can be generated from a 40 radius core (i.e. the region containing SMM1 and SMM4) in 4 â 106 yr, if stars formed at a constant rate with an average radial expansion speed of 0.3 km s-1 . This speed would appear plausible, given the relative motions observed in the core in the H13 CO+ line (Fig. 8). Note that there are fewer highly reddened stars to the south of ob ject E, suggesting that the star formation has shifted to the north by a distance of 20 (Fig. 5). Clearly, stars have been forming for a long time within CB 34 and the centre of activity may have shifted northwards from HH 290IRS to the vicinity of Ob ject E, while more dense gas still remains further north. The channel maps in H13 CO+ , a tracer of high density material in dense cores, demonstrate some evidence for large­scale rotation of 1.5 km s-1 over 90 (Fig. 8). With this interpretation, a mass of 90/(cos2 ) M within a radius of 45 is derived, where is the inclination of the rotation axis to the line of sight. This is consistent with the mass derived from CS measurements of 170 M (Launhardt & Henning 1998). SMM4 is redshifted by 0.5 km s-1 relative to SMM1. SMM2 is also identified in this molecule. We also find complex internal velocity structure on the scale of 10 , in the 12 CO J = 1 - 0 line (Fig. 7). This structure suggests that the inner 2 of the globule is being stirred up by many outflows. The numerous H2 knots spread about the inner globule may thus be part of several cavity walls. The turbulent motions of v 1.0 km s-1 combined with a total globule mass of Mc 170 M and radius of Rc = 0.75 pc (100 ), imply an energy decay rate of order Mc (v 3 /(2Rc ) = 0.04 L (Mac Low 1999). Therefore, quantitatively, the mechanical power of the outflows is sufficient to maintain the turbulence. The decay timescale from the largest scales is Rc /v = 7.5 â 105 yr. A further feature to note is the gap in CO emission at all velo cities directly between SMM1 and SMM4, with stronger blueshifted emission to the south. This suggests that star formation, which first o ccured near HH290IRS, was later centered in this gap, dispersing the material and compressing the surrounding cores. This has now triggered a third phase of star formation in the adjacent cores. Note also the appearance of weak blue-shifted emission at all velocities in the range -2-9 km s-1 from SMM2. This is the third most massive core in the globule, and may thus be just beginning to form a low-mass star as star formation gradually moves onwards in the north-west direction.

Fig. 8. H13 CO+ velocity channel maps of CB 34. The pixel size is 11 , the b eam size at 86 GHz is ab out 27 . H13 CO+ is a tracer of high density gas. Two ma jor condensations are seen, corresp onding roughly to SMM1 to the south-west, centered at v 0.4 km s-1 , and SMM4 to the north-east, centered at v 0.9 km s-1 .

Protostars in the SMM1 core may be responsible for both the H and Q jets (Fig. 6). This core, of size 10 , contains the IRAS source IRAS 0540+2059 but was previously not asso ciated with near-infrared emission. It is flattened at 450 µm with a minor axis of 70 (Huard et al. 2000), almost orthogonal to the Q-outflow of 68 . In contrast, the large-scale globule, on scales of 1­2 arcmin, is flattened almost in the perpendicular direction to the core. This is evident in the HCN molecule (148 , Afonso et al. 1998) and the far-infrared continuum (149 , Huard et al. 2000). In comparison, the orientations of the H jet and N jet are both 150. Hence, the ma jor outflows are aligned with either their core or the globule minor axis. It is also interesting to note that the parallel outflows are also parallel to the galactic plane (152 ). In contrast, the CO J = 2 - 1 outflow from SMM1 has a rather badly defined orientation of 115­130, which supports our assertion above that the CO lobes are produced by several outflows with a wide spread in ages while the H2 emission originates only from the latest outflows.

4.4. The globule
The extent of the globule is found from the photometric results (Fig. 4). Analysis of our wide field reveals that the the reddened stars (i.e. H - K > 1.0) lie predominantly within


10

T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34

Table 2. The diffuse halo brightness in mag/arcsec2 , found on insp ection of large-scale regions not confused by p oint sources in individual J , H and Ks images. Radial distances are from ob ject E. Name Extended halo Main halo Inner halo Radius( ) 65 50 30 J H K
s

21.0 19.8 18.7 20.2 19.0 18.2 19.5 18.2 17.0

CB 34 is not fully isolated. Its neighbour, CB 33, is at a pro jected distance of just 5 pc and connects on to a wide ridge of material as shown in the IRAS HIRES4 image requested from IPAC5 (Fig. 9). Filamentary structures are also present between the two globules, with clumps within 2 pc. For material which now forms CB 34 to have become disconnected from the ridge about 106 yr ago would require a separation speed of 2 km s-1 , about equal to the gravitational escape speed from an average molecular cloud. This material would still, however, have to condense into the Bok globule.

4.5. A diffuse near infrared halo
We detect a weak halo of diffuse emission in the J , H and Ks filters centered on the submillimetre and far infrared cores. Fig. 10 displays the diffuse halo in the H band. Figure 11 presents the radially-averaged halo brightness for all three bands. The numerous lo cations where stars have been subtracted hamper the extraction of radial profiles. Furthermore, the halo consists of a largescale diffuse component plus numerous lo calised regions of scattered light from young stars within the globule, as discussed below. We estimate a main halo of radius 50 (0.37 pc at 1500 pc). An inner halo which covers just the submillimetre cores SMM1 and SMM4 is also defined in Table 2. An extended halo, down to the Ks band noise level, is of radius of 65 (0.48 pc). The extended halo has an intensity of 18.7 mag/arcsec2 in the Ks band (Table 2). This converts to 4 â 10-4 erg s-1 cm-2 sr-1 . To determine the physical pro cess pro ducing the halo, we generated maps which lo cate the pixels which possess H - K in specific ranges. Quite remarkably, for the range of values (H - K ) = 0.5 --0.9 mag, the halo closely silhouettes the dense far-infrared cores mapped by Huard et al. (2000), shown here in Figs. 12 and 13.
4 HIRES uses the Maximum Correlation Method (MCM, Aumann et al. 1990) to produce images with b etter than the nominal resolution of the Infrared Astronomical Satellite (IRAS) data. 5 IPAC is Infrared Processing and Analysis Center established as the NASA archival center for the IRAS data and op erated by the JPL and California Institute of Technology (Caltech) http://www.ipac.caltech.edu

The halo is not a result of incomplete sky subtraction. Firstly, the halo region was well within all the frames taken. Secondly, the halo extends far beyond the core where a number of small halos around bright stars could combine to produce weak general diffuse emission. Thirdly, the J - H and H - K colours of the diffuse halo are distinct from the embedded stars. Magnitudes for the diffuse halo, estimated directly from contour maps in the three bands, are given in Table 2. These values are not consistent with thermal H2 emission from sho cks, which predicts (H - K ) = 1.5 - -3.4 mag (Smith 1995), far to o red. Scattering from an extended distribution of reddened stars within the globule would appear relevant to specific locations but the observed diffuse halo is generally weaker in the H band, the halo possessing H - K = 0.8 - 1.2 and J - H = 1.2 - 1.3. These values are, therefore, also not consistent with unreddened or reddened main-sequence stellar photospheres. We also expect that CB 34 is shielded from the UV photons of the ambient interstellar radiation field. Scattering of radiation from the deeply embedded infrared-excess stars could produce the halo. These stars are redder than the halo (J - H exceeding 1.7), which would be consistent with the stronger scattering of shorter wavelength radiation although the extra path length of this radiation would tend to have the opposite effect, reddening the emission. We find 23 embedded protostars with estimated luminosities of 10 L within a core radius of 65 . In comparison, the total bolometric luminosity of the cold dust in the submillimetre cores is 240 L , consistent with the stellar content. Therefore, given about two magnitudes of K -band extinction, it is clear that the infrared radiation from the protostars is largely absorbed and re-radiated by the cold dust with an average intensity of 2.5 â 10-3 erg s-1 cm-2 sr-1 . The fraction of scattered near-infrared radiation which escapes the cloud depends on several factors, including the clumpiness and albedo. With two magnitudes of extinction and an effective albedo of 0.25, the halo intensity level would be 1 â 10-4 erg s-1 cm-2 sr-1 . This is an estimate for the total scattered light, and is a factor of 4 lower than the observed K -band halo. This tentatively suggests, given the uncertainties involved, that scattered protostellar light is a minor contribution to the diffuse halo. The H - K and J - H colours correspond to that of reddened H2 fluorescence or any other mechanism which pro duces a cascade down the rotational and vibrational energy ladders of H2 . The inner halo displays more reddening, thus accounting for the anticorrelation of low H - K with the dust maps. We have calculated the fluorescent colours from the line fluxes measured for NGC 2023 (Burton et al. 1998) and predicted from detailed calculations (Black & van Disho eck 1987). These are: H - K = 0.5, J - H = 0.6 and H - K = 0.7, J - H = 0.5, respectively. These colours are indeed consistent with stronger H -band emission than from stellar photospheres, as found for the diffuse halo. Fluoresence implies a halo reddening corresponding to (J - H ) 0.6 - 0.8. This is consistent


T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34

11

30

Declination

CB34
21:00

CB33

20:30

5:50

49

48

47

46

45

44

Right ascension
Fig. 9. IRAS HIRES 100 µm map of the CB 34 region after 20 iterations. Image size is 2 â 2 . The displayed contour levels are 0.1, 1, 2, 4 and 8 Jy. Courtesy of the IPAC CALTECH http://www.ipac.caltech.edu.

with the reddening of the stars pro jected onto the diffuse halo of (J - H ) 0.5­0.9, despite the general difference in absolute J - H values.

4.6. Origin of the diffuse halo
We now discuss the various physical pro cesses which could pro duce diffuse near-infrared H2 emission: H2 -formation on grains; impact excitation of H2 by secondary electrons produced by the absorption of X-rays from embedded young ob jects; impact excitation of electronic states of H2 and direct impact excitation of vibrational levels in the H2 electronic ground state by secondary electrons pro duced by cosmic-ray ionisations; fluorescent excitation of H2 by internal UV photons emitted from embedded young stars. Interstellar molecular hydrogen is formed on the surface of grains (Hollenbach & Salpeter 1971; Jura 1975). Duley & Williams (1993) calculated the total H2 flux for canonical cloud parameters. For their assumptions, including steady state cloud parameters, formation rate of 10-17 /ns-1 and five per cent of the reformations leading to 1­0 S(1) photons, their calculation yields a 1­0 S(1)

intensity of just 10-8 erg cm-3 s-1 , for a cloud of mean depth 0.6 pc. This is far below detectable limits. Duley & Williams noted, however, that mixing between gas phases would raise the atomic density from the steady state value of n = 1 cm-3 to 100 cm-3 , increasing the 1­0 S(1) intensity to 10-6 erg cm-3 s-1 . This is still well below present measurement capabilities. Le Bourlot et al. (1995) calculated the resulting near-infrared emission spectra of H2 leaving the grains. The authors presented H and K band spectra assuming a) equipartition between translational, ro-vibrational, and vibrational degrees of freedom, b) H2 formation in the highest vibrational level, and c) H2 formation in v = 6 but at a low rotational excitation temperature of 65 K. The latter is suggested from catalysis experiments (Duley & Williams 1986). They found total K -band fluxes below 10-7 erg cm-3 s-1 , which confirm the above conclusion. Molecule formation, on the other hand, would provide a full explanation for the measured diffuse halo. In a formation scenario, as opposed to reformation, the entire cloud halo forms dynamically on the same time scale


12

T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34

Declination (2000)

20:59:30

21:00:00

30

01:00

30

08

06

04
Right Ascension (2000)

02

5:47:00

46:58

Fig. 10. Near infrared H-band image of CB 34 with the Calar Alto 3.5 m and the Omega Prime camera. The three lowest contour levels, corresp onding to the values in Table 2, are 19.8, 19.0 and 18.2 mag/arcsec2 .

as molecule formation. Hence the atomic density is high because equilibrium has not been reached. This scenario will be explored separately (Smith et al. 2001). Here, we note that with an atomic density of 104 cm-3 in the halo, the predicted 1­0 S(1) intensity is 10-4 erg cm-3 s-1 . The integrated K s intensity is then predicted to be 10-3 erg cm-3 s-1 , close to that observed. Smith et al. (2001) demonstrate more generally that molecule and cloud formation can both o ccur on time scales of 106 yr, much shorter than expected in many previous mo dels, yet consistent with Bok globule statistics. This is due to the revision in our ideas of the interstellar medium away from the concept of pressure-matched phases to that of supersonic turbulence.

In a gas of low fractional ionisation, the emission from v = 0, 1, 2 is dominated by direct impact excitation from non-thermal secondary electrons produced by cosmic ray absorptions in dark clouds (Gredel & Dalgarno 1995). Cosmic ray ionisations are always present, and provide in fact the main ionisation source which drives the chemistry in the dense molecular parts. The resulting H2 emission spectrum, no matter what is assumed for the H2 formation pro cess, is extremely red, with most of the emission o ccurring in the K -band (see e.g. Figs. 3a-d of Le Bourlot et al. 1995). Colours are dominated by values of (H - K ) 1 mag, which is in conflict with the observations presented above.


T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34

13

Fig. 11. The brightness in the globule (p er square arcsecond) averaged over each 2.5 for all three NIR bands after extraction of the p oint sources. The radial distance is measured from source E (s.f. Figs. 4 and 5). The background values are reached at a radius of 65 for J and Ks and 80 for the H -band.

20:59:30

21:00:00

Declination 30

01:00

30

08

06

02 04 Right ascension

5:47:00

46:58

Fig. 13. The 450 µm contour map of Huard et al. (2000) overlaid on the greyscale image of the locations associated with the H - K colour in the range 0.5­0.9.

01:00

30

08

06

02 04 Right ascension

5:47:00

46:58

Fig. 12. The 850 µm contour map, kindly provided by Huard et al. (2000), overlaid on the greyscale surface brightness map of the locations associated only with the H - K colour in the range 0.5­0.9. This is the range of values corresp onding to cascade H2 mechanisms with little or no reddening. Note that there is a quite detailed, although not p erfect, anti-correlation.

Apart from penetrating cosmic rays, absorptions of Xrays in dense molecular material resulting from young objects embedded in CB 34 may pro duce non-thermal, secondary electrons as well. The response of H2 to X-rays was calculated by Gredel & Dalgarno (1995) and Tin´ et al. e (1997). Again, the preferential excitation of the v = 1, 2 levels by electron impact results in high intensities of the transitions out of these levels. Consequently, the expected near-infrared emission from X-ray absorptions is very red,

with (H - K ) 1 mag. This is also in conflict with the observed colour of the halo in CB 34. The near-infrared emission spectrum expected from embedded sources of UV photons is significantly different from that expected from X-rays. Opposite to the direct excitation of ro-vibrational levels in all electronic states of H2 , including its ground state, UV photons populate electronic states allowed by the selection rules of electric dipole radiation only. The fluorescence process mainly in the Lyman and Werner bands populates the ro-vibrational levels of the ground state which is followed by cascading. (Black & van Disho eck 1987). The expected color of the pure fluorescence pro cess is (H - K ) = 0.5-0.7 mag. UV fluorescent excitation, however, requires the presence of massive stars. Even NGC 2023, the brightest source of fluorescent H2 , powered by a B star generates a 1­0 S(1) intensity of 7 â 10-5 erg s-1 cm-2 sr-1 . Here, however, we have no such massive star. We conclude that fast molecule formation within an interstellar medium which is constantly evolving provides the most plausible explanation for the diffuse halo. This scenario will be explored elsewhere (Smith et al. 2001). Tests for the model could be based on the low H2 /CO abundance, the effects of reformation heating and the H2 cascade spectrum. Near-infrared polarimetry would test the extent of the halo as a reflection nebula with internal sources of illumination.

20:59:30

21:00:00

Declination 30

4.7. Origin of CB 34
How did the globule form? If the diffuse halo is identified with hydrogen molecules, they form on the same time


14

T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34

04

CB34|J

04

CB34|H

03

03

02

02

01 Declination Declination

01

21:00

21:00

59

59

58

58

20

15
56

10

05

5:47:00

55

50

20

15
56

10

05

5:47:00

55

50

Right ascension

Right ascension

(a)

(b)

04

CB34|Ks

04

CB34|H2

03

03

02

02

01 Declination Declination

01

21:00

21:00

59

59

58

58

20

15
56

10

05

5:47:00

55

50

20

15
56

10

05

5:47:00

55

50

Right ascension

Right ascension

(c)

(d)

Fig. 14. Mosaics in J(a), H(b), Ks (c) and H2 (d): X-axes are Right Ascension and Y-axes are Declination (J2000). Images have b een obtained after reducing and mosaicing data by DIMSUM package under IRAF. Sizes of images are 8.78 â 8.75 arcmin2 .

scale of 106 yr as the globule age. This favours a mo del in which the cloud has formed by co oling and collapse of an atomic cloud. This cloud would have originally been, perhaps, a few parsecs in size and 10­100 cm-3 in density. Cloud formation is, however, a continuous pro cess and the interpretation of the halo in terms of molecule formation suggests that the cloud is still growing in mass.

The original warm atomic cloud may have been associated with the end of a ridge of gas seen in the IRAS 100 µm map. An alternative is that the cloud was torn off from a distant molecular cloud at the edge of the Gemini OB1 association. This would, however, imply an age of at least 107 yr which do es not agree with other estimates. We thus conclude that CB 34 has formed, along with


T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34
Right ascension

15

T Q5
SMM3 Q1 Q2 Q

G
CB034YC1

N4 N5
Q3 Q4 HH290 HH290N1 HH291Y
21:00:00

SMM2
30

SMM4 HH291X

CB34VLA1

N7
Declination

N6

H1
SMM1

SMM5

HH290N2

HH291 CB034YC1-I

HH290IRS

H2
30

H3 N8 H4

HH290S
09 08 07 06

05

04

03

02

01

5:47:00

Fig. 15. The central part of CB 34 as describ ed in Fig. 1.

Right ascension 50

Right ascension

N1

35

H5

45

30

40
Declination Declination

N2

25

35

20

20:57:30

N3

15

H6
11.6 11.2 10.8 10.4 5:47:10.0 09.6

25

21:01:10

02.4

5:47:02.0

01.6

01.2

00.8

00.4

09.2

Fig. 16. Close-up views of the H2 ob jects in the continuum-subtracted 1­0 S(1) filter. a) Details of H5 and H6 knots; b) details of N1, N2 and N3 knots.

CB 33, near the end of a dust ridge. Triggering by compression may have instigated the collapse and shortened the molecule formation time scale. A triggering agent is, however, not evident and slow spontaneous formation, followed by the present rapid evolution, cannot be excluded.

How has the globule evolved? One possibility is that a roundish globule has undergone an inside-out collapse. This formed a central flattened core about 106 yr ago. Collapse proceeded fast 40 south of this core, between SMM1 and SMM4, and these protostars evacuated a cen-


16

T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34

tral cavity in the flattened core. Turbulent motions, driven by the outflows, of 1­2 km s-1 led to members of the young group expanding to their present lo cations, up to 1 pc from the original core. The protostellar winds have driven into the annulus, forming a core SMM1 on the front side and a core SMM4 towards the rear. This has led to predominantly red-shifted and blue-shifted molecular H13 CO+ lines, respectively, with accelerated material further away from the initial core. This may also explain the main CO bipolar outflow, which appears to be aligned along the cores rather than along the H2 flows (Fig. 7), with only a weaker blue lobe extending along the N outflow. The second phase of star formation begun about 5â105 yr ago, within the two cores which have (in this scenario) a common rotation axis, which is ultimately responsible for generating parallel jets.

A diffuse halo is detected in J , H and K -band images, at a level consistent with that asso ciated with the rapid formation of molecular hydrogen at a density of 104 cm-3 in a time of order 106 yr. This suggests that the globule has formed directly out of an atomic cloud asso ciated with a dust ridge, from which it continues to accrete mass. The alternative remains that the diffuse halo is purely scattered light from infrared-excess protostars. There is, however, a clear observed deficit in the number of protostars as well as the dust bolometric luminosity, which makes this interpretation problematic. Two phases of star formation are suggested. The first phase generated the wide distribution of young stars. This activity triggered further core formation in the north of the globule, leading to the present generation of dense cores, outflows and protostars.
Acknow ledgements. We would like to thank Tracy Huard for supplying the far infrared images. We are extremely grateful to Martin Murphy for his supp ort in establishing the data reduction software.

5. Summary
We have obtained and analysed wide-field near-infrared images of the Bok Globule CB 34 in the JH K continua and the H2 2.12 µm, which reveal some of the properties of the cloud mass, distribution of young stars and their jets and sho cked outflows. We have also obtained velocity channel maps in molecular lines, CO maps revealing the co ol outflows and turbulent motions, and H13 CO+ maps displaying the dense cores. Two new long extended infrared jets were discovered. These jets are detected over 1.2 pc and 1.6 pc to the outskirts of the globule. Since no terminating bow sho ck is found and the tracing H2 knots are confined to the dark cloud, it is proposed that the jets extend much further into a lower density halo. The two jets are extremely well collimated and are parallel to within 3 . We estimate that the total sho cked luminosity is comparable to the mechanical luminosity, derived from CO measurement, consistent with jet-driven non-evolving outflow structure. With an estimated dynamical age of 7.6â105 years, the mechanical luminosity and sho cked luminosities are both Lmech = 0.02 L . These outflows may be driven by low-mass very young stars. We note, however, that a more viable explanation is that the CO outflows are dominated by the accumulated outflowing mass from older outflows whereas the H2 is from sho cks generated by the latest generation of outflows. JH K photometry yields the presence of a central concentration of reddened stars, and a surrounding lower density region of less reddened stars within the globule. The spatial spread suggests that the stars have been generated for 106 yr within the cores. The presence of a pre-mainsequence star, CB34FU (see Fig. 1), with an age of 106 yr supports this as a minimum age (Alves et al. 1997). The two outflows appear to be driven from the two densest submillimetre cores which we detect here in the H13 CO+ molecule, a tracer of dense gas. The axis separating these cores is roughly transverse to the outflow direction.

References
Adams, F. C., Lada, C. J., & Shu, F. H. 1987, ApJ, 312, 788 Afonso, J. M., Yun, J. L., & Clemens, D. P. 1998, AJ, 115, 1111 Alves, J. F., & Yun, J. L. 1995, ApJ, 438, L107 Alves, J. F., Hartmann, L., Brieno, C., & Lada, C. J. 1997, ~ AJ, 113, 1395 Alves, J. F., Lada, C. J., & Lada, E. 2001, in preparation Aumann, H. H., Fowler, J. W., & Melnyk, M. 1990, AJ, 99, 1674 Bizenb erger, P., McCaughrean, M. J., Birk, C., Thompson, D., & Storz, C. 1998, Proc. SPIE, 3354, 825 Black, J. H., & van Dishoeck, E. F. 1987, ApJ, 322, 412 Bok, B. J., & Reilly, E. F. 1947, ApJ, 105, 255 Burton, M. G., Howe, J. E., Geballe, T. R., & Brand, P. W. J. L. 1998, PASA, 15, 194 Carp enter, J. M., Snell, R. L., & Schloerb, F. P. 1995, ApJ, 445, 246 Clemens, D. P., & Barvanis, R. E. 1988, ApJS, 68, 257 Clemens, D. P., Yun, J. L., & Heyer, M. H. 1991, ApJS, 75, 877 Cowie, L. L., Gardner, J. P., Lilly, S. J., & McLean, I. S. 1990, ApJ, 360, L1 Duley, W. W., & Williams, D. A. 1986, MNRAS, 223, 177 Duley, W. W., & Williams, D. A. 1993, MNRAS, 260, 37 Eisenhardt, P., Dickinson, M., Stanford, S. A., Ward, J., & Valdes, F. 1996, May, available from ftp>iraf.noao.edu/iraf/contrib/dimsumV2/ Gredel, R., & Dalgarno, A. 1995, ApJ, 446, 852 Hogerheijde, M. R., van Dishoeck, E. F., Blake, G. A., & van Lngevelde, H. J. 1998, ApJ, 502, 315 Hollenbach, D., & Salp eter, E. E. 1971, ApJ, 163, 155 Huard, T. L., Weintraub, D. A., & Sendell, G. 2000, A&A, 362, 635 Jura, M. 1975, ApJ, 197, 575


T. Khanzadyan et al.: Active star formation in the large Bok globule CB 34 Koornneef, J. 1983, A&A, 128, 84 Lada, C. J., & Adams, F. C. 1992, ApJ, 393, 278 Launhardt, R., Evans I I, N. J., & Wang, Y. 1998, ApJS, 119, 59 Le Bourlot, J., Pineau des For^ s, G., Roueff, E., Dalgarno, A., et & Gredel, R. 1995, ApJ, 449, 178 Launhardt, R., & Henning, Th. 1997, A&A, 326, 329 Mac Low, M.-M. 1999, ApJ, 524, 169 Moreira, M. C., & Yun, J. L. 1995, ApJ, 454, 850 Rosen, A., & Smith, M. D. 2001, in preparation Scappini, F., Cecchi-Pestellini, C., Olberg, M., Casolari, A., & Fanti, C. 1998, ApJ, 504, 866 Schilke, P., Walmsley, C. M., des For^ G. P., & Flower, D. R. ets, 1997, A&A, 321, 293 Smith, M. D. 1991, MNRAS, 253, 175 Smith, M. D. 1995, A&A, 296, 789

17

Smith, M. D. 1996, in Disks and Outflows around Young Stars, ed. S. Beckwith, et al., Springer Heidelberg, CD-ROM, 809815 Smith, M. D. 2000, IAJ, 27, 25 Stanke, T. 2000, Ph.D. Thesis, Univ. Potsdam Stetson, P. B. 1987, PASP, 99, 191 Strom, K. M., Strom, S. E., Edwards, S., Cabrit, S., & Skrutskie, M. F. 1989, AJ, 97, 1451 Tin´ S., Lepp, S., Gredel, R., & Dalgarno, A. 1997, ApJ, 481, e, 282 Yun, J. L., & Clemens, D. P. 1994, ApJ, 365, L73 Yun, J. L., & Clemens, D. P. 1994, ApJS, 92, 145 Yun, J. L., Moreira, M. C., Torrelles, J. M., Afonso, J. M., & Santos, N. C. 1996, AJ, 111, 841 Yun, J. L., Moreira, M. C., Alves, J. F., & Storm, J. 1997, A&A, 320, 167