Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.mrao.cam.ac.uk/~eic22/Papers/curtisricher2010b.pdf Äàòà èçìåíåíèÿ: Tue Jul 27 17:42:25 2010 Äàòà èíäåêñèðîâàíèÿ: Tue Oct 2 07:45:35 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: massive stars

Mon. Not. R. Astron. Soc. 000, 1­17 (2010)

Printed 27 July 2010

A (MN L TEX style file v2.2)

A submillimetre survey of the kinematics of the Perseus molecular cloud ­ III. Clump kinematics
Emily I. Curtis1,2 and John S. Richer1,2 1
2

Astrophysics Group, Cavendish Laboratory, J. J. Thomson Avenue, Cambridge, CB3 0HE Kavli Institute for Cosmology, c/o Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 0HA

Accepted 2010 July 26; Received 2010 July 16; in original form 2010 May 20

ABSTRACT

We explore the kinematic properties of dense continuum clumps in the Perseus molecular cloud, derived from our wide-field C18 O J = 3 2 data across four regions ­ NGC 1333, IC348/HH211, L1448 and L1455. Two distinct populations are examined, identified using the automated algorithms C L FI N D (85 clumps) and G AU S S C L U M P S (122 clumps) on existing SCUBA 850 µm data. These kinematic signatures are compared to the clumps' dust continuum properties. We calculate each clump's non-thermal linewidth and virial mass from the associated C18 O J = 3 2 spectrum. The clumps have supersonic linewidths, NT /cs = 1.76 ± 0.09 (C L FI N D population) and 1.71 ± 0.05 (with G AU S S C L U M P S). The linewidth distributions suggest the C18 O line probes a lower-density `envelope' rather than a dense inner core. Similar linewidth distributions for protostellar and starless clumps implies protostars do not have a significant impact on their immediate environment. The proximity to an active young stellar cluster seems to affect the linewidths: those in NGC 1333 are greater than elsewhere. In IC348 the proximity to the old IR cluster has little influence, with the linewidths being the smallest of all. The virial analysis suggests that the clumps are bound and close to equipartition, with virial masses similar to the masses derived from the continuum emission. In particular, the starless clumps occupy the same parameter space as the protostars, suggesting they are true stellar precursors and will go on to form stars. We also search for ordered C18 O velocity gradients across the face of each core. Approximately one third have significant detections, which we mainly interpret in terms of rotation. However, we note a correlation between the directions of the identified gradients and outflows across the protostars, indicating we may not have a purely rotational signature. The fitted gradients are in the range G = 1 to 16 km s-1 pc-1 , larger than found in previous work, probably as a result of the higher resolution of our data and/or outflow contamination. These gradients, if interpreted solely in terms of rotation, suggest that the rotation is not dynamically significant: the ratios of clump rotational to gravitational energy are rot 0.02. Furthermore, derived specific angular momenta are smaller than observed in previous studies, centred around j 10-3 km s-1 pc, which indicates we have identified lower levels of rotation, or that the C18 O J = 3 2 line probes conditions significantly denser and/or colder than n 105 cm-3 and T 10 K. Key words: submillimetre ­ stars: formation ­ stars: evolution ­ ISM: kinematics and dynamics ­ ISM: individual: Perseus.

1 INTRODUCTION Stars form inside dense cores deep within molecular clouds. The similarity of the core mass function (CMF) to the initial mass function of stars (IMF, see Motte, Andre & Neri 1998; Alves, Lombardi & Lada 2007) has sometimes been used (e.g. Enoch et al. 2008) as evidence against theories of core formation where no correspon-

dence is anticipated, e.g. in the competitive accretion picture (Bonnell et al. 2001; Bate & Bonnell 2005). However, the mapping of a given CMF on to the resultant IMF is complicated by many factors (see e.g. Hatchell & Fuller 2008; Curtis & Richer 2010), such as the cores' varying multiplicity (Goodwin et al. 2008), star-forming efficiencies (Swift & Williams 2008) and/or mass-dependent lifetimes (Clark, Klessen & Bonnell 2007). In fact, as shown by Swift & Williams (2008) and Goodwin & Kouwenhoven (2009), diverse evolutionary schemes can map the observed CMFs on to the IMF,

E-mail:e.curtis@mrao.cam.ac.uk E-mail: jsr10@cam.ac.uk c 2010 RAS


2

E. I. Curtis and J. S. Richer
1.1 Nomenclature and clump populations

implying the current mass function data are not sufficient to discriminate between different theoretical models by themselves. The kinematics of star-forming cores, probed with spectralline observations of dense molecular tracers, are arguably the best discriminator between different models of core formation (e.g. ´ Andre et al. 2007; Kirk, Johnstone & Tafalla 2007). For example, cores created by shocks in large-scale flows exhibit large velocity gradients and are located at local maxima in the line-of-sight ´ velocity dispersion (e.g. Ballesteros-Paredes, Klessen & VazquezSemadeni 2003; Klessen et al. 2005). Alternatively, the magnetically controlled scenario has more quiescent velocity fields (e.g. Nakamura & Li 2005), well-matched to observations of isolated starless cores with subsonic levels of turbulence (e.g. Myers 1983; Caselli et al. 2002). Such quiescent cores are in opposition to purely hydrodynamic models of gravoturbulent fragmentation (see Mac Low & Klessen 2004) which produce a majority of cores with supersonic velocity dispersions. Offner, Klein & McKee (2008a) and Krumholz, McKee & Klein (2005) maintain that the physical mechanism of star formation depends on the much debated presence (or absence) of turbulent feedback. In one view, molecular clouds are short-lived (on timescales of order one dynamical time), non-equilibrium structures (e.g. Elmegreen et al. 2000; Hartmann 2001; Dib et al. 2007) characterized by transient turbulence, which dissipates quickly. In this scenario, stars could form by the collapse of discrete cores. Conversely, molecular clouds might form slowly and be quasiequilibrium objects (see e.g. Shu, Adams & Lizano 1987; McKee 1999; Krumholz & Tan 2007; Nakamura & Li 2007). This would require turbulence to be constantly injected into the clouds, either externally (from e.g. supernova blast waves or HI I regions) or internally (from e.g. protostellar outflows) and cores could form through competitive accretion (Bonnell et al. 2001). Offner et al. (2008b) demonstrate that these two disparate views of star formation give rise to distinguishable kinematics in dense stellar precursors and thus the study of core kinematics potentially offers a way to probe the turbulent state of the surrounding molecular cloud. This paper is concerned with the kinematical properties of dense clumps in the Perseus molecular cloud (hereafter simply Perseus) and compares them to theoretical models. We examine the linewidths, virial masses and ordered velocity gradients of two populations of clumps (Curtis & Richer 2010) identified across SCUBA dust continuum maps (Hatchell et al. 2005). These kinematic signatures are derived for each clump from our wide-field survey (Curtis, Richer & Buckle 2010a, Paper I) of four clusters of star-forming cores towards NGC 1333, IC348/HH211 (simply IC348 hereafter), L1448 and L1455 in the J = 3 2 rotational lines of CO and its common isotopologues 13 CO and C18 O. In the preceding paper of this series (Curtis et al. 2010b, Paper II), we explored the evolution of molecular outflows across these regions. We organize this paper as follows: the remainder of the introduction explains the naming conventions we have adopted to describe star-forming clumps/cores. §2 briefly describes the continuum data and associated clump catalogues alongside our spectralline datasets, from which we derive the clump kinematics. Various clump kinematic signatures are investigated in §3, including the clumps' linewidths (§3.1) and virial masses (§3.2), before we search for ordered velocity gradients across the face of the clumps in §3.3, which we interpret with a solid-body-rotation model. Finally, we summarize this work in §4.

We adhere to the clump hierarchy explained in Curtis & Richer (2010), which follows Williams, Blitz & McKee (2000). Briefly, molecular clouds contain clumps which in turn harbour cores, the direct precursors of individual or multiple stars. Clumps without embedded objects are starless, unless they are gravitationally bound when we refer to them as prestellar.

2 2.1

OBSERVATIONAL DATA Dust continuum

The SCUBA 850 µm data analysed have a beam size of 14 arcsec (0.017 pc at 250 pc, our adopted distance to Perseus, see Paper I) and are sampled on a 3 arcsec grid with a typical rms noise level, rms = 35 mJy beam-1 (Hatchell et al. 2005).

2.2

Clump catalogues

We investigate the kinematics of two populations of dust continuum clumps, which we identified in the SCUBA 850 µm data in a previous paper (Curtis & Richer 2010), using the two most popular automated algorithms: C L FI N D (Williams, de Geus & Blitz 1994) ¨ and G AU S S C L U M P S (Stutzki & Gusten 1990). Our motivation in that study was to determine which clump properties can be robustly measured and which depend on the extraction technique. We continue that approach in this work, analysing the kinematic signatures for each population in turn. In Curtis & Richer (2010), we located 85 and 122 clumps using C L FI N D and G AU S S C L U M P S respectively with peak flux densities 4rms = 140 mJy beam-1 , across the same areas where we also have C18 O J = 3 2 data (Paper I). The clumps were identified as starless, Class 0 or Class I protostars using the classifications of Hatchell et al. (2007, hereafter HFR07), which are based on source SEDs incorporating Spitzer data. We associate a clump with a HFR07 source, if the separation between the clump and HFR07 source peak positions is less than the clump's diameter (see Curtis & Richer 2010). Given the similarities in the extraction method and resulting population, most of the C L FI N D clumps are directly on top of the HFR07 sources and are unambiguously identified. However, three G AU S S C L U M P S sources are large enough to encapsulate two HFR07 sources and their results are correspondingly included in the population statistics twice.

2.3

Spectral-line

The C18 O J = 3 2 observations used in this analysis were undertaken with HARP (the Heterodyne Array Receiver Project; Buckle et al. 2009) on the James Clerk Maxwell Telescope and have been described in detail previously (Paper I). Striping artefacts resulting from systematic differences in calibration between HARP channels were removed using a `flatfield' procedure (fully described in Paper I). The final datasets are sampled on a 3 arcsec grid, distributed using a Gaussian gridding kernel with a 9 arcsec full-width half maximum (FWHM), yielding an equivalent beam size of 17.7 arcsec. The median rms spectral noise values, averaged across each region's map (measured in 0.15 km s-1 channels), are 0.18, 0.15, 0.15 and 0.14 K for NGC 1333, IC348, L1448 and L1455 respectively.
c 2010 RAS, MNRAS 000, 1­17


The kinematics of SCUBA clumps in Perseus
3 THE KINEMATICS OF THE SCUBA CORES In molecular clouds, the C18 O J = 3 2 line is expected to be optically thin and excited in dense regions (ncrit 104 cm-3 ). It should therefore trace material in the vicinity of or within the dense agglomerations surveyed with SCUBA. At the peak of every clump in our two catalogues (see §2.2), we extracted the C18 O J = 3 2 spectrum and fitted a Gaussian profile to it using S P L AT1 if the peak of the C18 O line was 3sp , where sp is the rms spectral noise. We did not detect C18 O at the 3 level towards 4 (5 per cent) C L FI N D and 6 (5 per cent) G AU S S C L U M P S positions. Of the detections, 13 (15 per cent) C L FI N D and 18 (15 per cent) G AU S S C L U M P S spectra have double-peaked line profiles. For comparison, Kirk et al. (2007, hereafter KJT07) made pointed N2 H+ J = 1 0 and C18 O J = 2 1 observations on the Institut de Radio Astronomie ´ Millimetric (IRAM) 30-m telescope (with beam sizes of 25 and 11 arcsec respectively) towards 157 candidate cores in Perseus ­ 44 pointings in our fields ­ and found 66 (42 per cent) C18 O spectra required two component Gaussian fits with three in one case. Generally, J = 2 1 data trace lower density (ncrit 103 cm-3 ) and/or colder material than the 3 2 at lower optical depths by a factor of 2­3 in the optically thin limit. Therefore the J = 2 1 line is unlikely to saturate if the J = 3 2 data are optically thin. This suggests there are separate emitting sources along the lineof-sight with different temperatures and/or densities which require multiple-component fits for one transition and not the other. An alternative is simply that at our velocity resolution (0.15 km s-1 compared to KJT07's 0.05 km s-1 ) many of KJT07's separate components are blended into one, which might partly explain our larger linewidths (see §3.1). 3.1 Linewidths The linewidths in molecular clouds are substantially larger than the sound speed of their constituent gas (e.g. Zuckerman & Evans 1974; Solomon et al. 1987). Such supersonic linewidths are thought to arise from chaotic and/or turbulent motions of unknown fundamental origin (see Elmegreen & Scalo 2004 and references therein). The velocity dispersion in a region scales with its size (e.g. Larson 1981) akin to the Kolmogorov law for subsonic incompressible turbulence, although turbulence in the interstellar medium is actually highly compressible and often supersonic causing shocks. However, dense star-forming cores harboured in cloud interiors possess near-thermal linewidths almost devoid of the turbulent movements that cause the broadening in their outer envelopes (e.g. Benson & Myers 1989; Barranco & Goodman 1998; Goodman et al. 1998; Pineda et al. 2010). The C18 O J = 3 2 transition is expected to probe a region intermediate in scale between the dense core and its environment, although it is not the tracer of choice for young prestellar cores. In the cold ( 10 K), high-density ( 105 cm-3 ) centres of such objects, the conditions may result in depletion, where many molecules (including CO and its isotopologues) freeze-out on to dust grains. Such zones are better probed by molecules such as N2 H+ and H2 D+ , but in any case the CO J = 3 2 transitions should probe higher densities than lower J ones. We take a core's total one-dimensional velocity dispersion, C18 O to be the quadrature sum of its thermal and non-thermal dispersions, T and NT respectively:
Part of the Starlink software collection, see http://starlink.jach. hawaii.edu. c 2010 RAS, MNRAS 000, 1­17
1

3

Table 1. Summary of the turbulent fractions, fturb , derived from Gaussian fits to the C18 O J = 3 2 line at the peak of every SCUBA clump identified with C L FI N D (Curtis & Richer 2010). The results are broken into either different regions in Perseus or clump classifications (see §2.2). The quoted errors ( ) are errors on the mean ( / N ) not sample deviations. Population All NGC 1333 IC348 L1448 L1455 Starless Protostars Class 0 Class I


Number 81 37 23 14 7 17 33 19 14

f

turb



f

1.76 2.05 1.36 1.8 1.5 1.8 1.89 2.04 1.7

0.09 0.12 0.09 0.2 0.4 0.2 0.12 0.16 0.2

The Class 0 and I populations combined.

2 2 2 C18 O = T + NT .

(1)

To find NT (plotted in Fig. 1), we calculate T :
2 T =

kT , mC18 O mp

(2)

where mp is the mass of a proton and mC18 O the relative molecular mass of C18 O (mC18 O = 30). The clump temperature, T , is taken from the kinetic temperatures of ammonia gas (Rosolowsky et al. 2008), where available, or 10 K and 15 K for identified starless and protostellar cores respectively where not. For those clumps without a HFR07 association we assume a temperature of 12 K. We calculate the turbulent fraction, plotted in Figs. 2 to 5, fturb = NT /cs , where cs = kT /mmp is the thermal sound speed in the bulk ¯ of the gas of mean molecular mass m = 2.33mp (assuming 1 ¯ He for every 5 H2 ) and the adiabatic index ( 7/5 for a diatomic molecule). For T = 10 K, cs = 0.22 km s-1 and at 15 K cs = 0.27 km s-1 . Summaries of fturb for various populations are listed in Tabs. 1 and 2. A majority of the linewidths are supersonic with means fturb = 1.76 ± 0.092 and 1.71 ± 0.05 for the C L FI N D and G AU S S C L U M P S populations respectively. The level of non-thermal broadening should take into account motions from both: (a) the source itself i.e. infall, rotation and outflows and (b) the source's natal environment. Young stars can inject significant energy into their surroundings, increasing the non-thermal motions within nearby star-forming cores. For example Caselli & Myers (1995) found an inverse relationship between a core's NH3 linewidth and the distance to the nearest young stellar cluster in Orion B. There are three such young clusters in Perseus (see Hatchell et al. 2005): IC348, NGC 1333 and the Per OB2 association. Therefore, of our fields, NGC 1333 and IC348 are very near to young clusters, while L1448 and L1455 are further away (L1455 being the furthest). The variations in fturb region-to-region (see Fig. 2), which we might attribute to each region's distance from a young stellar cluster, are slightly different for the two populations. For the C L FI N D clumps (see Tab. 1), NGC 1333 and L1448 have statistically similar fractions: 2.05 ± 0.12 and 1.8 ± 0.2 respectively. We ignore the

2

Errors quoted on averages throughout this paper are errors on the mean ( / N ) and not simple standard deviations ( ).


4

E. I. Curtis and J. S. Richer
1600
24

1400 1200

06

Declination (J2000)

21

03

1000
18

800

+32 : 00

600 +31 : 57
15

400
+31 : 12

45 : 00

40

20

44 : 00
-1

3 : 43 : 40

200
24 12 29 : 00 48 36 3 : 28 : 24

Non-thermal linewidth (km s 0.1 0.2 0.3 0.4 0.5 0.6

) 0.8 0.9

0

0.7

Right ascension (J2000)

48 18 46 16 44 14 12 +30 : 10 28 : 00 44 28 27 : 12 3 : 26 : 56 42

+30 : 40 26 : 00 48 36 24 12 3 : 25 : 00

Figure 1. Spatial distribution of C18 O non-thermal linewidth across NGC 1333 (upper left), IC348 (upper right), L1448 (lower right) and L1455 (lower left). The grey-scale is SCUBA 850 µm flux density in mJy beam-1 (Hatchell et al. 2005) with contours at 100, 200, 400, 800, 1600 and 3200 mJy beam-1 . Circles are positioned at the peak of every SCUBA clump in our catalogue identified using C L FI N D (Curtis & Richer 2010) where there is a 3 C18 O detection. The circles' size and colour reflect the linewidth as shown in the key.

Table 2. Summary of the turbulent fractions, f clump population as in Tab. 1. Population All NGC 1333 IC348 L1448 L1455 Starless Protostars Class 0 Class I Number 119 64 25 23 7 24 40 27 13 f
turb

turb

, for the G AU S S C L U M P S



f

1.71 1.91 1.44 1.52 1.5 1.61 1.80 1.73 2.0

0.05 0.07 0.07 0.13 0.3 0.11 0.10 0.11 0.2

2) has clumps in NGC 1333 with much larger turbulent fractions than the similar fturb means in IC348 and L1448: 1.91 ± 0.07 compared to 1.44 ± 0.07 and 1.52 ± 0.13. Therefore, for both populations, the IC348 linewidths are smaller than in NGC 1333, even though both regions are close to stellar clusters. This is possibly (as noted in Paper I) because the IC348 IR cluster is old and no longer actively forming stars, so it affects its environs much less than the NGC 1333 cluster. There is a wide range of clump fturb in L1448, which could be a result of a significant contribution to the C18 O linewidth from molecular outflows. Many of the protostellar clumps in L1448 drive particularly strong flows (Paper II). On dividing the clumps into protostellar and starless subsets (see Fig. 3), there appears to little difference between the distributions of turbulent fraction for both source types and clump populations. At the peak positions of the C L FI N D objects, we find fturb = 1.8 ± 0.2 for starless clumps compared to 1.89 ± 0.12 for protostars. Similarly, fturb = 1.61 ± 0.11 and 1.80 ± 0.10 for starless and protostellar G AU S S C L U M P S sources respectively. Further separation of the protostars into Class 0 and I clumps does not yield
c 2010 RAS, MNRAS 000, 1­17

L1455 populations in these comparisons as they are so few in number. IC348, on the other hand, has much smaller linewidths overall: fturb = 1.36 ± 0.09. The G AU S S C L U M P S population (see Tab.


The kinematics of SCUBA clumps in Perseus
0.5 0.4 Fraction of cores 0.3 0.2 0.1 0 0 0.5 1 1.5 0.5 0.4 Fraction of cores 0.3 0.2 0.1 0 0 0.5 1 1.5
NT NT

5

(a) clfind

All NGC1333 IC348

(b) clfind

All L1448 L1455

2 /c
s

2.5

3

3.5

0.5

1

1.5
NT

2 /c
s

2.5

3

3.5

(c) gaussclumps

All NGC1333 IC348

(d) gaussclumps

All L1448 L1455

2 /c
s

2.5

3

3.5

0.5

1

1.5
NT

2 /c
s

2.5

3

3.5

Figure 2. Distribution of clump turbulent fractions, fturb = NT /cs , measured from the HARP C18 O J = 3 2 data (see text), broken down by region for the C L FI N D (upper panels) and G AU S S C L U M P S populations (lower panels).

0.4 (a) clfind Fraction of cores 0.3

Protostars Starless

(b) clfind

Class 0 Class I Starless

0.2

0.1

0 0 0.5 1 1.5 0.4 (c) gaussclumps Fraction of cores 0.3
NT

2 /c
s

2.5

3

3.5

0.5

1

1.5
NT

2 /c
s

2.5

3

3.5

Protostars Starless

(d) gaussclumps

Class 0 Class I Starless

0.2

0.1

tions of KJT07 than their N2 H+ results. The high levels of nonthermal motions, shown by some of our sources, are simply not present in N2 H+ data. A striking contrast between the C18 O distributions is that many of our cores have high fturb ( 3.5), which are not seen for either protostars or starless cores by KJT07. In addition, many KJT07 starless cores have sub-thermal linewidths ( fturb < 1). These differences may be partly caused by the spectral resolution of our observations (v = 0.15 km s-1 compared to KJT07's 0.05 km s-1 ), rather than intrinsically different results. Our slightly poorer resolution means we cannot probe as narrow linewidths and we possibly blend many of the separate components, noted in KJT07's C18 O spectra, into one. Given the similarities it seems likely our C18 O J = 3 2 data probe the same regions as the J = 2 1 observations of KJT07, namely the outer parts of star-forming cores, referred to by KJT07 as their envelopes. This is probably because of C18 O freeze-out in the dense, cold interiors of star-forming cores. Such regions should see an enhancement in the abundance of say N2 H+ relative to C18 O. KJT07 investigated this possibility for their similar population of candidate cores in Perseus by measuring the variation in the integrated C18 O-to-N2 H+ ratio with peak SCUBA flux density (an approximate proxy for central H2 volume density). They found high C18 O-to-N2 H+ ratios (i.e. little C18 O depletion) occurred mainly in starless cores, which have the smallest SCUBA fluxes (i.e. lowest densities), whereas highflux cores (mainly protostars) have low ratios suggesting high levels of C18 O depletion. Given that our starless and protostellar C18 O distributions are quite similar, we may also conclude (like KJT07) that protostars do not affect their environment significantly. The starless distribution may differ from the true prestellar one as starless clumps with the largest linewidths are likely to be unbound, and thus the true prestellar turbulent fractions will be smaller. Finally, we also note little correspondence with the simulations of Klessen et al. (2005), whose results for a large-scale driven (LSD) gravoturbulent model are also plotted in Figs. 4 and 5. KJT07 claim to have a good fit to the Klessen et al. data with their C18 O starless population, although the model was designed to match N2 H+ observations. 3.2 Virial theorem

0 0 0.5 1 1.5
NT

2 /c
s

2.5

3

3.5

0.5

1

1.5
NT

2 /c
s

2.5

3

3.5

Figure 3. Distribution of clump turbulent fractions, fturb = NT /cs , measured from the HARP C18 O J = 3 2 data (see text), broken down by the different source types for the C L FI N D (upper panels) and G AU S S C L U M P S populations (lower panels).

statistical differences between their average fturb (see Tabs. 1 and 2). Unfortunately, such simple averages do not provide the whole picture we can see in the distributions (in Fig. 3). For example, if the small but significant population of high fturb starless cores were eliminated it would seem that protostars have marginally higher fturb on average. However, the large uncertainties on the average fturb do emphasize the need for larger samples of objects. For the C L FI N D population (but not the G AU S S C L U M P S sources) most of the high fturb protostars are Class 0, and all of the low fturb Class I, which is perhaps what we might expect if outflow power decreases from the Class 0 to I stage (Bontemps et al. 1996; Paper II). Overall, there are few significant differences between the starless and protostellar distributions. We can draw interesting comparisons with the C18 O (J = 2 1) and N2 H+ fturb distributions presented by KJT07 (see Figs. 4 and 5). Our distributions are closer to the C18 O distribuc 2010 RAS, MNRAS 000, 1­17

If the C18 O data trace material that will go on to form the final star and not just the ambient gas in which a YSO is embedded, the core linewidth can inform us about its stability, i.e. whether the cores are gravitationally bound. We can then distinguish bound, starless i.e. prestellar cores from unbound ones that will dissipate without forming a star. The virial theorem is often used to examine this stability and we will apply it to our two populations of clumps identified in SCUBA 850 µm data. The clump virial mass, Mvir , is used as an estimate of the internal energy of the clump whilst its dust mass at 850 µm, M850 , is as an estimate of the potential energy. If we assume that the linewidths only reflect gravity and we have spherical clumps the virial theorem is Mvir = a 2 R , (3) G where is the average three-dimensional velocity dispersion, R the clump radius and a a constant that depends on the form of the density profile (a derivation can be found in e.g. Binney & Tremaine 2008; Rohlfs & Wilson 2004). If we assume a power law density profile, (r) r-n , then (e.g. MacLaren, Richardson & Wolfendale 1988):


6
0.4

E. I. Curtis and J. S. Richer
(a) clfind

Fraction of cores

0.3

This work KJT07 C18 O KJT07 N2 H+

(b) gaussclumps

This work KJT07 C18 O KJT07 N2 H+

2 =

3 8 ln 2kT v2 18 O + C 8 ln 2 mp

1 1 - m mC18 ¯

,
O

(5)

0.2

0.1

0 0 0.5 1 1.5 0.4 (c) clfind Fraction of cores 0.3
NT

where mp is the proton mass, m = 2.33 to factor in the abundance of ¯ He relative to H2 and mC18 O = 30. This yields for a shallow power law clump, n = 1.5, consistent with the shapes found by Enoch et al. (2008)
2 /c
s

2.5

3

3.5

0.5

1

1.5
NT

2 /c
s

2.5

3

3.5

Mvir
This work Klessen et al.

=

3.4

This work Klessen et al.

(d) gaussclumps

a 4/3 â

vC18 O 1 km s-1 M.

2

+ 0.22

T 12 K (6)

0.2

R 0.02 pc

0.1

0 0 0.5 1 1.5
NT

2 /c
s

2.5

3

3.5

0.5

1

1.5
NT

2 /c
s

2.5

3

3.5

Figure 4. Distribution of fturb , for the starless clump populations identified with C L FI N D (left panels) or G AU S S C L U M P S (right panels), measured from the HARP C18 O J = 3 2 data. Overlaid for comparison are the distributions found by KJT07 for 157 candidate cores using C18 O J = 2 1 (upper panels, pink) and N2 H+ J = 1 0 (upper panels, blue). A prediction for cores from the gravoturbulent simulation of Klessen et al. (2005) is also shown for turbulence driven on large scales (lower panels, brown). For that simulation the turbulent fraction extends beyond 3.5 to 4.3 (not plotted), though there are small numbers of clumps in those bins ( 2 per cent).

0.4 (a) clfind Fraction of cores 0.3

This work KJT07 C18 O KJT07 N2 H+

(b) gaussclumps

This work KJT07 C18 O KJT07 N2 H+

0.2

0.1

0 0 0.5 1 1.5 0.4 (c) clfind Fraction of cores 0.3
NT

2 /c
s

2.5

3

3.5

0.5

1

1.5
NT

2 /c
s

2.5

3

3.5

This work Klessen et al.

(d) gaussclumps

This work Klessen et al.

0.2

0.1

0 0 0.5 1 1.5
NT

2 /c
s

2.5

3

3.5

0.5

1

1.5
NT

2 /c
s

2.5

3

3.5

Figure 5. Distribution of fturb , for the protostellar clump populations identified with C L FI N D (left panels) or G AU S S C L U M P S (right panels), measured from the HARP C18 O J = 3 2 data. Overlaid for comparison are the distributions found by KJT07 for 157 candidate cores using C18 O J = 2 1 (upper panels, pink) and N2 H+ J = 1 0 (upper panels, blue). A prediction for cores from the gravoturbulent simulation of Klessen et al. (2005) is also shown for turbulence driven on large scales (lower panels, brown). For that simulation the turbulent fraction extends beyond 3.5 to 4.3 (not plotted), though there are small numbers of clumps in those bins ( 2 per cent).

a=

5 - 2n . 3-n

(4)

Provided the C18 O line of FWHM, vC18 O , traces the bulk of the gas with a Gaussian velocity distribution and the same non-thermal linewidth as H2 , we can estimate

We expect objects in equipartition to have Mvir M850 whilst those that are self-gravitating should have Mvir 2M850 . We take data on the clumps' dust properties from Curtis & Richer (2010). Each clump's 850 µm mass (M850 ) assumes the dust is optically thin, has an opacity, = 0.012 cm2 g-1 and is at a single temperature, TD . Again these temperatures are taken as the NH3 kinetic temperatures (Rosolowsky et al. 2008), where available, or 10 K (for starless cores) and 15 K (for protostars) where not. The core radius, Rdec , is the geometric mean of the two core semi-major and -minor axis `sizes' each deconvolved with the beam size. These `sizes' are the standard deviation of the pixel coordinates about the core centroid, weighted by the pixel values. There are considerable uncertainties in any dust and virial mass estimates. The errors are very difficult to quantify for individual clumps without detailed modelling. Thus, we follow complementary studies (e.g. Buckle et al. 2010; Enoch et al. 2008) and do not attempt to account for the uncertainties in our analysis, except for a discussion of their magnitude, which follows. The dust masses depend on the assumed distance to Perseus and the dust properties (temperature and opacity). We try to minimize the effects of dust temperature by using the NH3 kinetic temperature as an estimate of TD , which should be an accurate measurement at high volume densities ( 104 cm-3 , e.g. Galli, Walmsley & Goncalves 2002), ¸ where the gas and dust are thermally coupled. Nevertheless our dust temperature estimates still do not account for variations in the dust temperature across a clump. A range of distances have been used in the literature for Perseus (220 to 350 pc see e.g. Paper I) and indeed it may not be a contiguous cloud at a single distance. These result in an uncertainty of a factor of 5 in the dust mass estimates. The virial masses depend on the assumed clump profile, with steeper profiles producing smaller masses, however it only varies by a factor of 1.7 between a constant density and 1/r2 profile. This combined with the uncertainties in the distance and linewidth produce again around a factor of 5 in uncertainty. As we previously noted (§3.1), it is likely that the C18 O J = 3 2 line only traces the envelope of a clump, which could lead to an over-estimate of the non-thermal H2 linewidth (compared to estimates from say NH3 or N2 H+ ) and correspondingly the virial mass. Even if our estimates of the dust and virial masses were infallible, assuming that a clump with Mvir M850 is unbound may be misleading; external pressure and/or magnetic fields may contain such a clump. For instance KJT07 calculate that cores with external pressures consistent with their previous Bonnor-Ebert sphere modelling (Kirk, Johnstone & Di Francesco 2006), should be considered in equipartition not merely self-gravitating if Mvir 2M850 . The C18 O virial masses, which we plot in Fig. 6 against the 850 µm mass, lie in the range 0.8 to 22.9 M for the C L FI N D sources and 0.4 to 17.8 M for the G AU S S C L U M P S sources, reflecting their smaller radii. Most of the clumps lie scattered near the `equipartition' line, Mvir = M850 . The results are similar for
c 2010 RAS, MNRAS 000, 1­17


The kinematics of SCUBA clumps in Perseus
1.5

7

1 /M )

(a) clfind NGC 1333 IC348 L1448 L1455

(b) gaussclumps NGC 1333 IC348 L1448 L1455

Table 3. Summary of virial mass statistics for the identified in Curtis & Richer (2010), which have All the quoted errors ( ) are errors on the mean ( deviations. Population No 50 24 16 5 5 17 33 19 14 Mvir (M ) 7.8 8.9 4.0 12.6 7.9 4.6 9.6 10.1 9 M (M ) 0.9 1.2 0.8 3.9 1.6 0.8 1.2 1.4 2

SCUBA clumps a HFR07 classification. / N ) not simple sample
C L FI N D

vir

log(M

0.5

M850 /Mvir 1.32 0.96 1.38 2.7 1.4 1.4 1.3 1.3 1.2

ratio 0.13 0.15 0.22 0.5 0.6 0.2 0.2 0.2 0.3

0

All
-0.5

(c) clfind Starless Protostars

(d) gaussclumps Starless Protostars

1 /M )

NGC 1333 IC348 L1448 L1455 Starless Protostars Class 0 Class I

vir

log(M

0.5

0

-0.5

1 /M )

(e) clfind Starless Class 0 Class I

(f ) gaussclumps Starless Class 0 Class I

Table 4. Summary of virial mass statistics for the G AU S S C L U M P S SCUBA clumps identified in Curtis & Richer (2010), which have a HFR07 classifi cation. All the quoted errors ( ) are errors on the mean ( / N ) not simple sample deviations. Population No 64 30 17 12 5 24 40 27 13 Mvir (M ) 7.1 8.9 5.0 6.5 4.3 6.0 7.8 8.4 7 M (M ) 0.7 1.1 0.8 1.5 1.1 0.6 1.0 1.1 2 M850 /Mvir 0.9 0.68 0.75 2.2 0.8 0.87 0.91 1.06 0.67 ratio 0.1 0.08 0.14 0.4 0.3 0.16 0.13 0.18 0.14

vir

log(M

0.5

0

All
-0.5 -1 - 0. 5 0 0.5 log(M
850

1 /M )

1.5

2

-0.5

0

0.5 log(M
850

1 /M )

1.5

2

Figure 6. Virial mass versus clump dust mass for the C L FI N D (left panels, (a), (c) and (e)) and G AU S S C L U M P S (right panels, (b), (d) and (f)) populations, differentiated by region (top panels, (a) and (b)) or the classifications of HFR07 (the rest). Virial masses are estimated from the C18 O J = 3 2 linewidth at the catalogue positions of Curtis & Richer (2010). The masses are calculated using the temperatures from Rosolowsky et al. (2008) where possible or 10 and 15 K for starless and protostellar cores respectively where not. Lines denote where we expect cores to be in equipartition, i.e. Mvir = M850 (dashed) or the limit where they are self-gravitating, Mvir = 2M850 (dot-dashed), although the reader should note the considerable uncertainties in both masses. Only clumps which have a HFR07 source association are plotted.

NGC 1333 IC348 L1448 L1455 Starless Protostars Class 0 Class I

the different algorithms, regions and evolutionary types (see Tables 3 and 4 for summaries of the population statistics). Least-squares straight-line fitting results in poorly-constrained power law exponents of 0.8­0.9 for almost all the correlations. Given the uncertainties in both the virial and dust mass estimates, it is not possible to draw unequivocal conclusions about the stability of individiual clumps. However, we will compare distinct populations. The different regions show similar behaviour with a lot of scatter. The L1448 population is outlying for both algorithms with a high M850 /Mvir ratio. L1448 is interesting as it contains large SCUBA flux densities and massive outflows but has relatively weaker C18 O emission. The NGC 1333 population has a smaller M850 /Mvir ratio than the others, which we might interpret as being more unbound and therefore transient, with a larger number of cores above the Mvir = M850 threshold. This could be a result of its more perturbed environment. Protostellar clumps, by definition should be gravitationally bound. If starless clumps occupy a similar Mvir - M850 sample space, then this suggests that they too are a gravitationally bound
c 2010 RAS, MNRAS 000, 1­17

population. In Fig. 6, the starless and protostellar clumps occupy the same regions of the plot and their mass ratios are very similar in Tabs. 3 and 4. This implies that the starless population is gravitationally bound as well and truly prestellar in nature (as found by Enoch et al. 2008). The clumps also seem spread around the `equipartition' line Mvir = M850 . Approximate equipartition was also found by Caselli et al. (2002) in a sample of 60 starless cores using N2 H+ linewidths. Other studies e.g. the C18 O cores looked at by Tachihara et al. (2002) found the virial mass considerably larger than the core mass estimate. The latter result is successfully explained by the gravoturbulent models of Klessen et al. (2005), who find a large majority of starless cores with Mvir M for both large and small driving-scale turbulence models, with their protostellar sources non-coincident in parameter space at M > Mvir . Our results would not favour such a model for Perseus, even in its clustered environment. 3.3 Localized velocity gradients

Systematic variations in the C18 O line centre velocity are apparent across the face of many of the identified SCUBA clumps, e.g. the velocity gradually increases along a particular direction. There are a number of plausible explanations for such gradients: (i) rotation, (ii) outflows, (iii) motions between smaller unresolved constituent


8

E. I. Curtis and J. S. Richer
possibly freeze-out at cores. 3.3.1 5000 AU in starless and young protostellar

clumps or (iv) if gravoturbulent models are true, cores form at the stagnation points in convergent flows (e.g Padoan et al. 2001) and velocity gradients may arise from colliding gas streams. It is difficult to distinguish between each scenario so we follow the majority of studies and focus on analysing the rotational properties of the clumps from measured velocity gradients. A strong motivation to understand the details of rotation in star-forming cores is provided by the discs out of which most, if not all, stars are born (e.g Shu et al. 1987). Planets are thought to originate inside such protoplanetary discs and the details of their formation crucially depend on various disc parameters, such as surface density, controlled by the detailed evolution of the angular momentum of the parent star-forming core (e.g. Lissauer 1993; Ruden 1999). Additionally, there is the classical `angular momentum problem' of star formation (e.g. Spitzer 1978): the angular momentum of prestellar cores is orders of magnitude larger than that which can be contained within a single star, even though cores are observed to be rotating much less than originally predicted (Goodman et al. 1993; Caselli et al. 2002)3 . A plausible solution is provided by magnetic braking at the early low-density phases, where the field lines are strongly coupled to the gas and transfer angular momentum from the contracting core to the surrounding medium (e.g. Mouschovias 1987). Recent results suggest that gravitational interactions may dominate over magnetic braking. MHD models of self-gravitating, decaying (Gammie et al. 2003) and driven (Li et al. 2004) turbulence have angular momenta consistent with each other and observations but Jappsen & Klessen (2004, hereafter JK04) and Tilley & Pudritz (2004) find similar results in purely hydrodynamic frameworks. Even if magnetic fields do dominate, they cannot indefinitely strip away angular momentum or support a core against collapse because ambipolar diffusion will eventually lead to dynamic collapse. Finally, for cores to fragment into multiple systems some angular momentum must be present and perhaps cores with the largest angular momenta will go on to form binaries (e.g. Larson 2003; Goodwin et al. 2007). The pre-eminent observational study in this area is the search for solid-body-rotation in 40 NH3 cores by Goodman et al. (1993, hereafter GBF93). In all their objects the rotational energy is at most a few per cent of the gravitational and cannot provide support. An evolutionary sequence can be built using their results with others (e.g. Caselli et al. 2002), demonstrating the specific angular momentum ( j = J /M ) decreases with decreasing scale (JK04). Rotational signatures become more complicated at higher resolution and once collapse has started in protostars. Belloche et al. (2002) studied a young Class 0 protostar, IRAM 04191+1522, finding two distinct regimes of collapse: the inner (r 2000 - 4000 AU), rapidly collapsing and rotating whilst the outer (4000 r 11000 AU) has only moderate infall and rotation. The fall in rotational velocity beyond the 4000 AU boundary and the flat inner profile suggests that the inner region's angular momentum is conserved whilst it is dissipated in the outer perhaps by magnetic braking. At the distance to Perseus the r 4000 AU boundary is bigger than the JCMT beam (32 arcsec diameter compared to 15 arcsec) so we may be able to probe the inner regions, although C18 O will

Velocity gradient fitting

A clump undergoing solid-body rotation will display a linear velocity gradient across its face perpendicular to the rotation axis. we follow GBF93 and fit a linear gradient, vLSR , across each clump using: vLSR = v0 + a + b (7)

where v0 is the systemic clump velocity, and the angular ascension and declination offsets from the clump centre and a and b the projections of the gradient per radian on to the and axes. The gradient has a magnitude, G = |vLSR | = (a2 + b2 )1/2 /D, where D is the distance to the object, at an angle (east of north), G = arctan(a/b). We explore the C18 O J = 3 2 velocity field of our two populations of SCUBA 850 µm clumps. At each SCUBA map pixel the corresponding spectrum in the HARP data has been fitted with a Gaussian profile to extract its centre velocity and linewidth. Fits were only performed where the spectral peak was greater than three times the rms noise, estimated on a line-free portion of the spectrum for every spatial position. L1455 was omitted, as its weak C18 O emission had too few fits. We then used the Levenberg-Marquardt algorithm implemented in the scientific python module, SciPy4 , to perform non-linear least-squares fitting of Equation 7 to every clump about its SCUBA peak. GBF93 performed the same analysis on sets of randomly generated maps with no gradient and report that approximately 10 per cent of those random maps were found to have significant gradients if a "3 " criterion is used i.e. G 3G . Using the same significance criterion we expect similar levels of reliability. 3.3.2 A signature of rotation?

Furthermore, Dib et al. (2010) point out that measurements of core angular momenta from global velocity gradient fitting tend to overestimate the intrinsic (three-dimensional) angular momenta by a factor of 810, as complicated fluctuations in the three-dimensional velocity field are smoothed out.

3

We noted there are four plausible causes of velocity gradients. In this section, we look at some of the suggestions other than rotation in more detail. First, some of our clumps may be composed of smaller constituents. The gradient then might measure the velocity dispersion of the multiple cores (GBF93). If the SCUBA clumps are gravitationally bound, as suggested in §3.2, such dispersion in a multi-core clump may signal its rotation about a common centre. HFR07 estimate within their similar core catalogue, 10 per cent of clumps will break into multiple sources, by comparing the most luminous sources to higher resolution observations. Second, the exact region that is assigned to each clump can change the fitted gradient's magnitude, direction and level of significance. One of GBF93's original sample, TMC-1C, a starless core in Taurus, has been mapped over a larger area at higher resolution by Schnee et al. (2007). They found a more complicated pattern of local gradients no longer consistent with solid-body or differential rotation, which had previously been noted by Caselli et al. (2002). Finally, if the C18 O line centre is affected by outflows, then this could be misinterpreted as a sign of rotation. GBF93 also fitted velocity gradients to C18 O J = 1 0 data in a sub-sample of their cores ­ both with and without embedded objects ­ and found they were in the same direction as the gradients from NH3 . However, they did not formally examine any outflows. Towards starless
4

www.scipy.org. c 2010 RAS, MNRAS 000, 1­17


The kinematics of SCUBA clumps in Perseus
Table 5. Approximate outflow orientation compared to the velocity gradient direction for C L FI N D clumps with clear bipolar outflows. Classifications are from HFR07. Sub-region NGC 1333 IC348 IC348 IC348 L1448 L1448 L1448
C L FI N D

9

Table 6. Number of clumps with significant velocity gradients by region, source type and algorithm. Percentages give the detection rate for the particular type of source in each region or overall. Source Type NGC 1333
C L FI N D

Class 0 0 0 0 0 0 0

ID 6 1 2 4 1 2 4

G (deg E of N) 126 -51 50 4 139 78 154

Outflow orientation (deg E of N) -113 -64 68 0 124 165 168

Region IC348 Population 2(17 %) 3(100 %) 1(50 %) 1(17 %) 7(30 %)

L1448 0(0 3(75 0(0 1(9 4(25 %) %) %) %) %) %) %) %) %) %)

Total 6(32 13(72 5(45 3(11 27(35 10(48 7(35 3(35 15(24 35(31 %) %) %) %) %) %) %) %) %) %)

Starless Class 0 Class I No HFR07 ID All Starless Class 0 Class I No HFR07 ID All

4(67 7(64 4(44 1(8 16(41 4(50 6(50 3(33 11(31 24(37

%) %) %) %) %) %) %) %) %) %)

G AU S S C L U M P S

cores, outflows should not be a problem. To explore the effect on our gradients, we compare the direction of the velocity gradient to the orientation of the outflow from every protostellar source with a significant gradient in the C L FI N D catalogue. We estimated the outflow position angle from our CO J = 3 2 datacubes (a full outflow analysis was presented in Paper II). For some of the protostellar clumps, the outflows are complicated and we have ignored such cases, only comparing clumps where a characteristic bipolar structure is apparent (see Table 5). In a number of these cases the velocity gradient is in the same direction as the outflow, implying we are tracing the outflow rather than solid-body rotation with the gradient. If we were viewing solid-body rotation of the core or a disc around it, we would expect the gradient and outflow to lie perpendicular to one another. The data themselves present a more complicated picture. In Fig. 7, we overlay the outflows on top of the C18 O centre velocities for the sources in Table 5. Often the higher C18 O line centres (i.e. red-shifted gas) correspond closely to the extent of the redshifted CO outflow lobe and similarly for the blue-shifted gas. For instance, the complex outflow structure in L1448-1 is closely followed by the line centre data. Furthermore, the blue-shifted outflow lobe in NGC 1333-6 is coincident with blue-shifted C18 O line centres. The outflow itself extends away from the clump to the west, resulting in orientations that seem very disparate in Table 5. In two cores there seems to be little outflow-gradient correspondence: IC348-2 and L1448-4. Perhaps it is more convincing to argue that none of the cores exhibit gradients perpendicular to their outflows. Therefore, there appears to be a link between the outflow orientation and the velocity gradient, although how much this influences the gradient found in every protostar is difficult to quantify. Thus, we must view the interpretation of a linear velocity gradient as unequivocal evidence of solid-body rotation with some scepticism. However, in the absence of a strong outflow i.e. for starless clumps and assuming we can consider each clump as a single object, rotation would seem to be a likely cause of a coherent velocity gradient and the only one we consider in the following analysis.

Population 6(55 %) 0(0 1(33 %) 0(0 0(0 %) 0(0 3(38 %) 1(6 10(42 %) 1(4

3.3.3 Fitted gradients Table 6 lists the rates of detecting significant velocity gradients towards various core populations, separated by their different clumpfinding algorithms, regions or types of source. They are similar for both algorithms with approximately a third of clumps displaying significant velocity gradients. The detection rate for Class 0 sources in the C L FI N D catalogue is nearly twice as large as for the other types and Class 0s with G AU S S C L U M P S. This is probably because Class 0 sources have typically stronger outflows than Class Is
c 2010 RAS, MNRAS 000, 1­17

(e.g. Bontemps et al. 1996; Paper II). It is clear in Fig. 7 that many of the outflows extend beyond their associated clump. However, the clumps found by G AU S S C L U M P S are smaller (Curtis & Richer 2010) with a definite elliptical shape, which does not necessarily follow the orientation of the outflow. The lower detection rate for Class 0 sources with G AU S S C L U M P S is perhaps another example of how their velocity gradients are dominated by outflows. If the smaller G AU S S C L U M P S cores trace their shape less effectively, we would expect a lower detection rate. The fitted gradients are in the range 1.5 to 13.0 km s-1 pc-1 and 2.5 to 16.0 km s-1 pc-1 for C L FI N D and G AU S S C L U M P S respectively. We provide a summary of all the different rotational properties for different categories of clumps in Tables 7 and 8. Fig. 8 is a histogram of their magnitude by region. The distributions are qualitatively similar for the different regions with the usual uncertainty due to the small sample sizes; Kolmogorov-Smirnov (K-S) tests do not conclusively reject or confirm the hypothesis that each sample is drawn from the same population. On average across all the regions, the gradients are G = (5.7 ± 0.5) km s-1 pc-1 for C L FI N D and G = (6.9 ± 0.6) km s-1 pc-1 for G AU S S C L U M P S . The fitted G AU S S C L U M P S gradients are larger than those for C L FI N D but may be misleading as there are only 16 clumps with detections in both. We plot the gradient found with the C L FI N D designation versus that with G AU S S C L U M P S for these common sources in Fig. 9. There is fairly good agreement on the gradient with the protostellar sources deviating the most from the line of equal gradients. This might again be because the C18 O line centres are affected by outflowing gas which is not well traced by the G AU S S C L U M P S outlines. There are too few sources of each classification to draw firm conclusions about any trends with age. It does seem in the distributions of Fig. 10 that the starless cores have G 10 km s-1 pc-1 and larger gradients are protostellar. The protostellar average is larger than the starless one for the G AU S S C L U M P S population: G = (7.0 ± 1.0) compared to G = (5.0 ± 0.6) km s-1 pc-1 . Although they are similar for C L FI N D G = (5.8 ± 0.7) compared to G = (5.5 ± 1.0) km s-1 pc-1 . This might point to an outflow contribution to the gradient. The gradients we find in Perseus are much greater than those originally found by GBF93 or Caselli et al. (2002) in their N2 H+ J = 1 0 survey of 60 low mass cores: 0.3 to 3.9 and 0.5 to 6.0 km s-1 pc-1 respectively. Outflows cannot explain all of this


10

E. I. Curtis and J. S. Richer
8.9 (a) NGC 1333 Clump 6
20

30

(b) IC348 Clump 1

9.1 9.0 8.9 8.8 8.7 8.6 8.5 8.4 8.3

8.8

19 : 00

8.7
01 : 00

40

Declination (J2000)

8.6 Declination (J2000)

8.5

30

20

8.4

+32 : 00 : 00
18 : 00

8.3

8.2
+31 : 17 : 40

+31 : 59 : 30

18

16

Right ascension (J2000)

14

12

3 : 29 : 10

8.1

44 : 00

56

Right ascension (J2000)
(d) IC348 Clump 4

52

3 : 43 : 48

8.2

45

9.00 (c) IC348 Clump 2 8.94
30

9.4

9.3

30

8.88 8.82

9.2

Declination (J2000)

04 : 00

9.1

15

Declination (J2000)

8.76 8.70

9.0

03 : 00

30

8.64 8.58

8.9

8.8
03 : 00

45

8.52 8.46

8.7

+32 : 02 : 30 02 44 : 00

8.6
+32 : 02 : 30

Right ascension (J2000)
4.9

58

56

54

3 : 43 : 52

8.40

54

52

Right ascension (J2000)

50

48

3 : 43 : 46

8.5

30

(e) L1448 Clump 1

(g) L1448 Clump 4

4.48

5.40
4.8
47 : 00

46 : 00

30

4.40

5.25
4.7
30
30

4.32

44 : 00

5.10
4.24 Declination (J2000)

4.6 Declination (J2000)

Declination (J2000)

46 : 00

4.95
30

45 : 00

4.5

4.16

30

4.80

4.4
43 : 00

30

4.08

4.65

4.00

45 : 00

4.3
44 : 00

4.50
+30 : 44 : 30

3.92

4.2

+30 : 42 : 30

(f) L1448 Clump 2
40 36 32 3 : 25 : 28

4.35
+30 : 43 : 30

3.84
31 29 28 27 25 3 : 25 : 24

4.1

42

40

38

Right ascension (J2000)

Right ascension (J2000)

36

34

32

3 : 25 : 30

Right ascension (J2000)

Figure 7. Outflow orientation compared to the velocity gradient direction for clumps in Table 5. The colour-scale is the C18 O J = 3 2 line centre in km s-1 for map pixels allocated to each clump. Overlaid are contours of CO J = 3 integrated antenna temperature TA dv for the blue- and red-shifted line wings, integrated from: -5 to 3 (blue) and 12 to 18 km s-1 (red) for NGC 1333, -5 to 6 (blue) and 11 to 15 km s-1 (red) for IC348 and -5 to 0 (blue) and 8 to 15 km s-1 (red) for L1448. Contours are at various heights in K km s-1 (n is an integer from zero to nmax ): (a) (2+4n), nmax = 4 (blue) and (4+4n), nmax = 6 (red). (b) (6+2n), nmax = 6 (blue) and (6+2n), nmax = 3 (red). (c) (8+1.5n), nmax = 6 (blue) and (4+1.5n), nmax = 8 (red). (d) (5+n), nmax = 4 (blue) and (2+n), nmax = 2 (red). (e) (2+n), nmax = 8 (blue) and (2+2n), nmax = 7 (red). (f) (2+4n), nmax = 5 (both blue and red). (g) (2+4n), nmax = 5 (blue and red).

Table 7. Average rotational properties derived from the velocity-gradient fitting for the C L FI N D clump population. All the errors listed ( ) are errors on the mean and not sample deviations. Population All NGC 1333 IC348 L1448 Starless Protostars Class 0 Class I Number 27 16 7 4 6 19 14 5 G G (km s-1 pc-1 ) 5.7 5.7 5.0 7 5.5 5.8 6.1 4.9 0.5 0.7 0.6 2 1.0 0.7 0.8 1.5
rot

0.003 0.0010 0.006 0.016 0.007 0.003 0.004 0.001

j (km s 0.0019 0.0017 0.0017 0.0032 0.0021 0.0019 0.0021 0.0011


-1

j

pc)

0.014 0.0040 0.015 0.032 0.013 0.008 0.010 0.003

0.0002 0.0002 0.0003 0.0009 0.0005 0.0003 0.0003 0.0001

c 2010 RAS, MNRAS 000, 1­17


The kinematics of SCUBA clumps in Perseus

11

Table 8. Average rotational properties derived from the velocity-gradient fitting for the G AU S S C L U M P S clump population. All the errors listed ( ) are errors on the mean and not sample deviations. Population All NGC 1333 IC348 Starless Protostars Class 0 Class I Number 35 24 10 10 10 7 3 G G (km s-1 pc-1 ) 6.9 7.4 5.6 5.0 7.0 6.4 8 0.6 0.7 0.8 0.6 1.0 1.0 3
rot

0.004 0.006 0.001 0.001 0.004 0.005 0.010

j (km s 0.00120 0.0012 0.0011 0.0010 0.0012 0.0012 0.0012

-1

j pc) 0.00014 0.0002 0.0002 0.0002 0.0002 0.0002 0.0003

0.013 0.015 0.007 0.006 0.009 0.008 0.012

0.5 0.4 Fraction 0.3 0.2 0.1 0 0

0.7
(a) clfind

All NGC1333 IC348

(b)

All L1448
Fraction

0.6 0.5 0.4 0.3 0.2 0.1

(a) clfind

Protostars Starless

(b)

Class 0 Class I Starless

2

4

6

8
-1

10 pc
-1

12

0

2

4

6

8
-1

10 pc
-1

12

0 0 2 4 6 8
-1

10 pc
-1

12

0

2

4

6

8
-1

10 pc
-1

12

Gradient/km s
0.5 0.4 Fraction 0.3 0.2 0.1 0 0

Gradient/km s

Gradient/km s
0.8 (c) gaussclumps 0.6 Fraction 0.4 0.2

Gradient/km s
(d)

(c) gaussclumps

All NGC1333 IC348

Protostars Starless

Class 0 Class I Starless

4

8 Gradient/km s
-1

12 pc
-1

16

0 0 4 8 Gradient/km s
-1

12 pc
-1

16

4

8 Gradient/km s
-1

12 pc
-1

16

Figure 8. Distribution of significant (G /G 3) velocity gradients for the C L FI N D (top) and G AU S S C L U M P S (bottom) clump populations. The L1448 distribution is not plotted for G AU S S C L U M P S as it has only one significant gradient.

Figure 10. Distribution of significant (G /G 3) velocity gradients for the classifications of HFR07 for the C L FI N D (top) and G AU S S C L U M P S (bottom) clump populations.

10

(a)

(b)

gaussclumps gradient/km s

-1

pc

8

NGC1333 IC348 L1448

Starless Class 0 Class I

6

4

2

0 0 5 clfind gradient/km s
-1

10 pc
-1

0

5 clfind gradient/km s
-1

10 pc
-1

Figure 9. Fitted velocity gradient for G AU S S C L U M P S versus C L FI N D designations for the 16 clumps which have significant detections with both algorithms, by region (left) or by the classifications of HFR07 (right). The line marks where both gradients are equal.

difference because our starless clumps have significantly larger gradients as well. The discrepancy is reminiscent of that found by GBF93 towards B361, which was also investigated by Arquilla & Goldsmith (1985) with 13 CO J = 1 0 data. GBF93 estimated a gradient five times smaller than the earlier study, 0.7 versus 3.3 km s-1 pc-1 , which had led Arquilla & Goldsmith (1985) to conclude that rotation was dynamically significant in dark clouds.
c 2010 RAS, MNRAS 000, 1­17

GBF93 explain the difference by emphasizing that tracers of lower critical densities can produce complicated velocity patterns dominated by outflows or clump-to-clump motions that mimic solid body rotation. However, Olmi, Testi & Sargent (2005) fitted gradients to a number of cores in Perseus, finding similar magnitudes for the two common cores also in GBF93 and generally in the same range 0.14 to 2.32 km s-1 pc-1 . Their study used C18 O J = 1 0, CS J = 2 1 and N2 H+ J = 1 0 data with their C18 O fitting not displaying markedly different gradients. A critical factor will undoubtedly be resolution. The angular resolution of previous surveys is worse than the 17.7 arcsec of our maps: 88, 54 and 46 arcsec for GBF93, Caselli et al. (2002) and Olmi et al. (2005) respectively. Hence even for their closest sources in Taurus (at 140 pc), their best linear resolution (0.037 pc) is nearly twice as large as ours (0.021 pc). A larger beam will tend to smooth out differences between regions and reduce the overall gradient magnitude. This effect may be large enough to explain the differences for cores perhaps three times as distant as Perseus5 , in a beam four times as large as ours.

-1

5

The most distant sources examined by Caselli et al. (2002) are L1031B at 900 pc and L1389 at 600 pc.


12
15 (a)

E. I. Curtis and J. S. Richer
NGC1333 IC348 L1448 Starless Class 0 Class I
26 24
(b)

-1

pc

10

Gradient/km s

-1

47 27

22

5

39 12 8 5 21 35 6 25 34 17 18 9

31:20

33 58

64

4 41 29 17 18 16

0

(c)

(d)

NGC1333 IC348 L1448
-1

Starless Class 0 Class I

Declination (J2000)

clfind

-1

pc

10

18

Gradient/km s

16

63 16

8

2 9 26 3 50 26 44

10 7 20
11

5

14

13
gaussclumps 1 1.5 Axis Ratio 2 2.5 1.5 Axis Ratio 2 2.5

1

Figure 11. Velocity gradient versus axis ratio for C L FI N D (top panels) and G AU S S C L U M P S sources (bottom panels) broken down by region (left panels) or by the classifications of HFR07 (right panels).

10

12

0

15

30
08

20

10

3:29:00

50

40

28:30

Right ascension (J2000)

0.35 0. 3 0.25 Fraction 0. 2 0.15 0. 1 0.05 0 0 10 20 30 40 50 60 70 80 90 10 20 30 40 50 (a)

Starlink GAIA::Skycat SURF:remdbm

06

3:28:56.673

04

All NGC1333 IC348

(b)

Emble

Protostars Starless
Declination (J2000) 7 6 02 3 22 18 21 10 58 10 23 5 17 14 1 7 2 4 4

Aug 20,

60

70

80

90

Angle between ma jor axis and gradient/deg

Angle between ma jor axis and gradient/deg

3.3.4 Gradient orientation Centrifugal stresses due to energetic rotation will flatten cores along the rotation axis. This tends to produced oblate cores (GBF93), currently the favoured shape (e.g. Jones, Basu & Dubinski 2001), which could also be caused by strong magnetic fields. If rotation dominates, we might expect the the degree of core elongation to increase with increasing velocity gradient and a core's major axis to lie parallel to the gradient i.e. perpendicular to the rotation axis. Fig. 11 shows no relationship between the magnitude of the velocity gradient and the elongation of a clump, quantified through the axis ratio (taken from Curtis & Richer 2010). It comes as no surprise then that the angle between the clump major axis and the velocity gradient in Fig. 12 (only computed for the G AU S S C L U M P S sources), seems to be distributed at random as well. Therefore, clump rotation is unlikely to be energetically significant. Is there a correlation between the orientation of the rotation axes of neighbouring clumps or an anti-correlation to keep the net angular momentum small? In Fig. 13, we overlay the velocity gradients on dust emission maps of the various regions. Reassuringly

46

Starlink GAIA::Skycat SURF:remdbm Emble

47

48

Figure 12. Distribution of the angle between the velocity gradient and the major axis of the G AU S S C L U M P S clumps by region (left) or by the classifications of HFR07 (right). If rotational motions are elongating the clumps we would expect the gradient and axis to be in the same direction making the angle zero.

31:56 45:00

32:00

50

40

30

20 10 3:44:00 Right ascension (J2000)

50

40

43:30

7 1 4

ic348_850.sdf 3:44:13.995 32:01:45.80 J2000

18
Aug 20, 2009 at 11:02:18

Declination (J2000)

44

45

2 43 30:40 3:26:00 41 42

50

40

30 Right ascension (J2000)

20

10

25:00

Starlink GAIA::Skycat SURF:remdbm Emble

Figure 13. Spatial distribution of the velocity gradient detections in NGC 1333 (top), IC348/HH211 (middle) and L1448 (bottom). Arrows are drawn from the peak position of the associated clump for clumps where there are significant detections in either the C L FI N D or G AU S S C L U M P S catalogues. The direction and length of the arrow represents the direction and magnitude of the gradient and the detections are labelled with the reference from our SCUBA catalogues (Curtis & Richer 2010). The colour-scale is SCUBA 850 µm emission from 0 to 1600 mJy beam-1 . c 2010 RAS, MNRAS 000, 1­17

l1448_850.sdf 3:25:30.406 30:43:54.02 J2000 Aug 20, 2009 at 11:03:30


The kinematics of SCUBA clumps in Perseus
the directions and magnitudes of the gradients measured using both algorithms' allocations match very well. In small patches the arrows of neighbouring clumps line up but on large scales there seems to be little correlation. This is similar to the gravoturbulent models of JK04, where shortly after collapse starts, the angular momenta in small regions are well aligned over a moderate correlation length with cores further away spinning in random directions. Their explanation is that neighbouring cores are formed from the same reservoir of material and therefore can be expected to have similar angular momenta. During subsequent accretion and evolution, the associated correlation length decreases and neighbours lose their alignment, as embedded cores are ejected and turbulence disrupts the current material or brings in new gas.
0.8

13

0.6 Fraction

clfind Goodman et al. Caselli et al.

gaussclumps Goodman et al. Caselli et al.

0.4

0.2

0 0 0.05 0. 1
rot

0.15

0

0.05

0.1
rot

0.15

0.2

Figure 14. Distribution of the ratio of rotational to gravitational energies, rot . For comparison the distributions found by GBF93 (green) and (Caselli et al. 2002, orange) are overlaid.

0.4 0.3

3.3.5 Dynamic support? To quantify the level of support that the velocity gradients might provide against gravitational collapse we calculate rot , the ratio of rotational to gravitational energy. GBF93 define this ratio for a uniform density sphere: rot = (1/2)I 2 R3 G 2 = 2 /5R 3GM 3GM sin2 i (8)

Fraction

All clfind NGC1333 IC348 L1448

All gaussclumps NGC1333 IC348

0.2 0.1 0 -4 -3.5 -3 log(j /km s
-1

- 2. 5 pc)

-2

-3.5

-3 log(j /km s
-1

-2.5 pc)

-2

Figure 15. Distribution of j by region and algorithm: C L FI N D on the left and G AU S S C L U M P S on the right.

where the moment of inertia is I = pM R2 with p = 2/5 for a uniform density sphere and the angular momentum is = G / sin i with i the angle of inclination to the line of sight. Assuming values representative of these clumps and sin i = 1, this becomes: rot = 6.2 â 10
-4

3.3.6

Angular momenta

The core rotation speed can be quantified via the specific angular momentum, j = J /M , i.e. the angular momentum per unit mass: 2 j = R2 5 (10)

Rdec 0.02 pc

3

M M

-1

G 1 km s-1 pc-1

2

. (9)

Across all the clumps the average values of rot are: rot = (0.014 ± 0.003) and (0.013 ± 0.004) for C L FI N D and G AU S S C L U M P S respectively. At most the rotational energy is just six per cent for the C L FI N D sample and 14 per cent for G AU S S C L U M P S of the gravitational energy in the clumps. Fig. 14 shows the distribution of rot for the two clump populations. There is a strong preference for low ratios, with the fraction of clumps rapidly falling off as rot increases. Additionally, both populations closely match each other and the ratios found by Caselli et al. (2002), while the GBF93 distribution is wider, spreading to higher values. K-S tests could not confirm the samples are drawn from the same population as each other or either of the examples from the literature to any degree of significance. Even though we found larger velocity gradients than the previous studies, rot is not any higher. This is probably because our clumps are significantly smaller: on average the clump radii are 0.03 and 0.02 pc for C L FI N D and G AU S S C L U M P S respectively while the Caselli et al. (2002) cores have an average of 0.06 pc. As rot depends on R3 this will reduce our values by 9 relative to theirs for an equally massive source, compensating for our larger derived gradients. Thus, rotation is not dynamically significant in star-forming cores and should not support the core against collapse. It would be interesting to examine how the level of support varies region-to-region and as a clump evolves. However, the numbers here are too small to draw any firm conclusions, with one or two objects heavily distorting overall averages. This is further compounded by the outflow ambiguities of the protostellar sources. Nevertheless, there are no vast differences in the populations by region or age with almost all the clumps having rot 0.02.
c 2010 RAS, MNRAS 000, 1­17

where = G / sin i is again the angular velocity of the clump, R its radius and the factor 2/5 is for a constant density sphere. This yields for representative values: j = 1.6 â 10
-4

G 1 km s-1 pc-1

Rdec 0.02 pc

2

km s

-1

pc

(11)

The distribution of j has been important, with modellers using the one derived by GBF93 either as an input to their models or a target (e.g. Burkert & Bodenheimer 2000; JK04). Furthermore, the distribution of angular momentum across a core and its evolution may prove whether magnetic fields are necessary to strip off angular momentum as collapse ensues. In Fig. 15, we plot the distribution of j by region. There is little difference between the different algorithms or clumps from one region to the next. For C L FI N D the overall average is j = (1.9 ± 0.2) â 10-3 km s-1 pc whilst j = (1.7 ± 0.2) â 10-3 , (1.7 ± 0.3) â 10-3 and (3.2 ± 0.9) â 10-3 in NGC 1333, IC348 and L1448 respectively. The G AU S S C L U M P S data are skewed to lower j due to their smaller clump radii but there are still few differences across the regions: j = (1.20 ± 0.14) â 10-3 , (1.2 ± 0.2) â 10-3 , (1.1 ± 0.2) â 10-3 km s-1 pc overall and in NGC 1333 and IC348 respectively. JK04 distinguish between starless and protostellar sources in their gravoturbulent simulations which focus particularly on rotational properties. They find wide continuous distributions for j, stretching over two orders of magnitude, as we do. Their prestellar cores have an order of magnitude larger j than the protostars. Fig. 16 compares our distributions of j (overall and for the starless clumps) to previously observed distributions and JK04. To match a


14
0.5

E. I. Curtis and J. S. Richer
(a) (b)

0.4

Caselli et al. JK04, scaling 1

clfind Starless

Fraction

0.3

0.2

0.1

0 (c) (d)

GBF93+ JK04, scaling 2
0.4

gaussclumps Starless

Fraction

0.3

0.2

0.1

0 (e)

0.4

Pirogov et al. JK04, scaling 3

-4

-3

-2 log(j /km s
-1

-1 pc)

0

1

Fraction

0.3

0.2

0.1

0 -4 -3 -2 log(j /km s
-1

-1 pc)

0

1

Figure 16. Distribution of specific angular momenta, j for starless cores. (a) N2 H+ J = 1 0 cores of Caselli et al. (2002, grey) compared to the JK04 model (brown) scaled by n = 105 cm-3 and T = 10 K. (b) All our C L FI N D detections (red) and just those in starless clumps (black). (c) NH3 cores of GBF93, Barranco & Goodman (1998) and Jijina et al. (1999) (green) compared to the JK04 model (brown) scaled by n = 104 cm-3 and T = 10 K. (d) all our G AU S S C L U M P S detections (blue) and just those in starless clumps (black). (e) N2 H+ J = 1 0 high-mass cores of Pirogov et al. (2003) compared to the JK04 model (brown) scaled by n = 104 c -3 and T = 50 K.

particular set of observations, JK04 scale their code output according to the mean density, n, and temperature, T , with j depending on these as j T / n. With appropriate scalings their data match three different observational regimes: low-mass cores mapped in the high-density tracer N2 H+ (Caselli et al. 2002), low-mass cores in the slightly lower density tracer NH3 (GBF93 with Barranco & Goodman 1998 and Jijina, Myers & Adams 1999) and high-mass cores in the high-density tracer N2 H+ (Pirogov et al. 2003). Our observed distribution has no overlap with the high-mass cores of Pirogov et al. (2003) and little with GBF93, having much lower j. For the starless clumps only, we find j = (2.1 ± 0.5) â 10-3 and (1.0 ± 0.2) â 10-3 km s-1 pc for C L FI N D and G AU S S C L U M P S respectively. The most similar distribution is that of Caselli et al. (2002), which is over twice as wide with a much larger average, j = 2.2 â 10-2 km s-1 pc, but does overlap with ours at low j. For our average to match that of Caselli et al. (2002), presumably our clumps would have to be considerably denser and/or colder than the values assumed by JK04 of n = 105 cm-3 and T = 10 K. In §3 we noted that the critical density is expected to be ncrit = 104 cm-3 .

However, in radiative transfer models of a free-falling protostellar envelope, we found peak H2 number densities of greater than a few times 105 cm-3 were required to reproduce the C18 O J = 3 2 line strengths seen in our data (Curtis 2009). Alternatively or additionally, we are seeing considerably smaller angular momenta then those measured by previous authors. The narrow width of the distribution compared to previous work is likely to reflect the uniform environment of our cores in a single cloud rather then spread over many different ones as in Caselli et al. (2002) and GBF93. On the other hand, the smaller angular momenta observed in our starless clumps may not indicate lower degrees of rotation in the clump core but may reflect a more slowly rotating outer envelope, making comparisons with tracers such as NH3 inappropriate. Redman et al. (2004) examined the asymmetric HCO+ line profiles across L1689B, a prestellar core in Ophiuchus, which they modelled as a rapidly-rotating inner core contained within a static envelope. Again, the key diagnostics are the densities and temperatures traced by the C18 O J = 3 2 line, which can only be properly answered with radiative transfer modelling. We would expect to probe high densities, ncrit 104 cm-3 , but earlier we saw the similar non-thermal linewidths for the J = 3 2 line compared to the J = 2 1, suggesting we might not be detecting very different material. A further consideration in the low temperature environments of starless cores is CO depletion. Estimating the degree of depletion in these cores is difficult without similar observations in "late-depleting" molecules, such as N2 H+ . However, models ^ (e.g. Walmsley, Flower & Pineau des Forets 2004) suggest that CO freezes-out in the central 7000 AU of a prestellar core (this is really a temperature and density dependence). Furthermore, we noted in §3.1 that there is some existing evidence for considerable C18 O depletion in candidate star-forming cores in Perseus (KJT07) Thus, our C18 O signal in evolved starless clumps and young protostars at least should have little contribution from the dense inner core and this reasoning supports the idea we are probing a more slowly rotating envelope. In JK04, protostars have a distribution of j that narrows with time, staying around the same average. Their distribution is also similar to that observed for binaries (see Fig. 17). Our distributions for all the sources are as narrow as those for the protostars alone in JK04's models. For C L FI N D there is a trend of decreasing j with protostellar age: j = (2.1 ± 0.3) â 10-3 and (1.1 ± 0.1) â 10-3 km s-1 pc for Class 0 and I clumps, which is not reproduced with G AU S S C L U M P S: j = (1.2 ± 0.2) â 10-3 (Class 0) and (1.2 ± 0.3) â 10-3 km s-1 pc (Class I). Of course in these sources there is potential confusion of the rotational signatures with outflows which is perhaps why the expected decrease in j from starless to protostellar clumps is not observed. 3.3.7 Scalings with clump size

GBF93 found that many of the derived rotation parameters scale with radius, including G ( R-0.4 ) and j ( R1.6 ) but not rot . We plot the same parameters as a function of clump size in Fig. 18. Most of the relations found by GBF93 hold for these clumps, although there is considerable scatter. GBF93 show that the relation for the gradient is implied from a linewidth-size relation (e.g. Larson 1981), v = R0.6 , for a core in virial equilibrium. we find rot is independent of R as well, implying M R3 for solid body rotation ( = constant) or for differential rotation ( R-1 ) that M R. The SCUBA clumps follow M R2 (Curtis & Richer 2010), so we can infer a more complicated form of rotation than pure solid body or differential to ensure rot is independent of size.
c 2010 RAS, MNRAS 000, 1­17


The kinematics of SCUBA clumps in Perseus
0. 8 (a) (b) (c) (d)

15

JK04 15 %
0. 6

JK04 30 %

JK04 45 %

JK04 60 %

Fraction

0. 4

0. 2

0 (e) (f ) (g) (h)

Duquennoy & Mayor
0. 6

Simon

All clfind Protostars Class 0 Class I

All gaussclumps Protostars Class 0 Class I

Fraction

0. 4

0. 2

0 -4 -3 log(j /km s
-1

-2 pc)

-1

-4

-3 log(j /km s
-1

-2 pc)

-1

-4

-3 log(j /km s
-1

-2 pc)

-1

-4

-3 log(j /km s
-1

-2 pc)

-1

Figure 17. Distribution of the specific angular momenta, j, for protostars. (a) to (d) show JK04's model M6k2 with different star formation efficiencies: (a) 15, (b) 30, (c) 45 and (d) 60 per cent. (e) Binaries among nearby G stars of Duquennoy & Mayor (1991). (f) Binaries in Taurus (Simon 1992). (g) and (h) show the distributions from paper for the C L FI N D and G AU S S C L U M P S designations respectively, decomposed according to the designations of HFR07.

20
-1

(a) clfind

(b)

gaussclumps Gradient/km s-1 pc 10 6 4

Gradient/km s-1 pc

10 5

2 1

Starless Class 0 Class I
(c)

Starless Class 0 Class I
(d)

3 2 0.1

10 6
rot

0.01
rot

4 3 2 0.001

Starless Class 0 Class I
(e)

Starless Class 0 Class I
(f )

subcritical environment, where j scales with R. In our data a horizontal line is almost as convincing as the slope of GBF93, indeed the spread of clump sizes is not really wide enough to constrain the trend unambiguously. Our results are therefore consistent with this magnetically-controlled scenario: there is a lot of spread about j 10-3 km s-1 pc but no j much below this value. How the rotational velocity changes across a clump is less clear and is a key discriminator between different theories. Such measurements would require high-resolution, interferometric observations of many cores in different stages of collapse.

-1



10

-4

0. 1

0.01

0.003
-1

pc

0.01 0.001 0.001

4

SUMMARY

j /km s

j /km s

10

-4

Starless Class 0 Class I
0.01 0.02 R
dec

Starless Class 0 Class I
0.03 /pc 0.04 0.06 0.02 R
dec

3 â 10 10 0.05 0.06
-4

-4

0.03 /pc

0.04

Figure 18. Rotational parameters as functions of clump size for C L FI N D (left) and G AU S S C L U M P S (right). G versus clump radius; rot versus Rdec (middle panels); j versus Rdec (bottom panels). Lines mark the GBF93 relations: G R-0.4 and j R1.6 .

We have analysed the kinematic signatures of two populations of dust continuum clumps identified in SCUBA 850 µm data using either C L FI N D or G AU S S C L U M P S (Curtis & Richer 2010), across four regions (NGC 1333, IC348/HH211, L1448 and L1455) of active star formation in the Perseus molecular cloud. The kinematics were derived from our large-scale (600 arcmin2 ) C18 O J = 3 2 survey (Paper I). First, the non-thermal contribution to the C18 O linewidth was calculated and compared to the results of KJT07 and the gravoturbulent models of Klessen et al. (2005): (i) Most clumps have supersonic non-thermal linewidths, NT /cs = 1.76 ± 0.09 (C L FI N D population) or 1.71 ± 0.05 (G AU S S C L U M P S), with overall distributions similar to the C18 O J = 2 1 data of KJT07 towards dense core candidates in Perseus. This implies we are also probing the `envelopes' of the star-forming clumps rather than deep within their interiors. (ii) There is little difference in the level of the non-thermal linewidth contribution between protostars and starless clumps, implying that protostars do no affect their environment significantly.

The scaling of j with size is of particular interest. Ohashi et al. (1997) looked at rotation in protostars in Taurus using C18 O J = 1 0 interferometric spectra. They found at distances 3500 AU, that j was at a constant value of 10-3 km s-1 pc. These findings are consistent with magnetically-controlled collapse elucidated by Belloche et al. (2002), where in a dense inner core, the angular momentum is locked at a minimum value. This region could be interpreted as a magnetically supercritical core, decoupling from its
c 2010 RAS, MNRAS 000, 1­17

-1

pc


16

E. I. Curtis and J. S. Richer
found in our clumps, suggestively hovers around the background value, 10-3 km s-1 pc-1 , found by Ohashi et al. (1997) at radii 0.03 pc in Taurus's protostellar cores. To see how exactly j varies with radius we need to probe across the clumps themselves, using higher-density tracers or higher-resolution observations.

Clumps in NGC 1333 and IC348 are in the vicinity of young stellar clusters, explaining the wider linewidths in NGC 1333 ( NT /cs = 1.91 ± 0.07) than L1448 (1.52 ± 0.13) for the G AU S S C L U M P S objects (although their linewidths are similar for C L FI N D population: 2.05 ± 0.12 and 1.8 ± 0.2 respectively). IC348 however, has the narrowest lines of all (1.36 ± 0.09 with C L FI N D and with 1.44 ± 0.07 G AU S S C L U M P S ), possibly related to its statistically younger population (Hatchell et al. 2005), containing more starless clumps. Second, the C18 O linewidths allowed an estimate of each clump's virial mass: (iii) Most clumps in Perseus appear bound and close to equipartition. The ratio of the SCUBA to virial mass for clumps based on flux allocation using C L FI N D is M850 /Mvir = 1.32 ± 0.13, compared to 0.9 ± 0.1 for the G AU S S C L U M P S decomposition. (iv) Starless clumps occupy a similar section of M850 -Mvir parameter space to the protostars suggesting that they are gravitationally bound and therefore truly prestellar. This has implications for our previous work investigating the starless CMF in Perseus (Curtis & Richer 2010). These results are in contrast to the models of Klessen et al. (2005) where a majority of starless core have Mvir M and occupy a different parameter space to the protostars. Finally, we fitted a linear velocity gradient to the C18 O line centres, across the face of each SCUBA clump, again identified either with C L FI N D or G AU S S C L U M P S. Significant detections were found in approximately a third of the clumps, which we interpreted in terms of solid-body rotation. (v) A correlation between CO J = 3 2 outflows and the C18 O line centres is observed, implying any gradients in outflow sources, i.e. protostars, may reflect outflows rather than rotation. (vi) The fitted gradients, in the range 1 to 16 km s-1 pc-1 , are larger than those previously observed. This probably results from the higher angular resolution of our observations and/or outflow contamination. (vii) There is no correlation between the gradient direction and the clump orientation nor between the gradient and clump axis ratio, implying that the rotation is not energetically significant. Furthermore, the ratio of the clump rotational to gravitational energy is typically 0.02, demonstrating that rotation is not dynamically important as well, with a distribution very similar to that of the N2 H+ cores of Caselli et al. (2002). (viii) The distribution of specific angular momentum is narrower and centred around lower values, j 10-3 km s-1 pc, than previous studies. Interpreting the results of JK04 this suggests a denser and/or colder environment for the Perseus clumps than seen in the Caselli et al. (2002) sample (taken to have n 105 cm-3 and T 10 K). This would seem unlikely given the large linewidths we observed towards the clumps, although radiative transfer modelling does suggest we are probing moderately high densities ­ peak H2 number densities of a few â105 cm-3 for a free-falling protostellar envelope (Curtis 2009). Thus, we are probably seeing lower levels of rotation in our clumps. (ix) There are no strong trends in the rotational parameters with radius. The somewhat inconclusive findings of this paper motivate further work on the star-forming cores in Perseus. First, the construction of accurate radiative transfer models to determine precisely the conditions traced by the CO J = 3 2 transitions. Second, higher-resolution or higher-tracer-density observations. j

5

ACKNOWLEDGMENTS

EIC thanks the Science and Technology Facilities Council (STFC) for studentship support while carrying out this work. We thank Jane Buckle and Gary Fuller for reading carefully an early version of this paper. We are grateful to the referee, for comments and suggestions which improved the clarity of this paper. The JCMT is operated by The Joint Astronomy Centre (JAC) on behalf of the STFC of the United Kingdom, the Netherlands Organisation for Scientific Research and the National Research Council (NRC) of Canada. We have also made extensive use of the SIMBAD data base, operated at CDS, Strasbourg, France. We acknowledge the data analysis facilities provided by the Starlink Project which is maintained by JAC with support from STFC. This research used the facilities of the Canadian Astronomy Data Centre operated by the NRC with the support of the Canadian Space Agency.

REFERENCES Alves J., Lombardi M., Lada C. J., 2007, A&A, 462, L17 ´ Andre P., Belloche A., Motte F., Peretto N., 2007, A&A, 472, 519 Arquilla R., Goldsmith P. F., 1985, ApJ, 297, 436 ´ Ballesteros-Paredes J., Klessen R. S., Vazquez-Semadeni E., 2003, ApJ, 592, 188 Barranco J. A., Goodman A. A., 1998, ApJ, 504, 207 Bate M. R., Bonnell I. A., 2005, MNRAS, 356, 1201 ´ Belloche A., Andre P., Despois D., Blinder S., 2002, A&A, 393, 927 Benson P. J., Myers P. C., 1989, ApJS, 71, 89 Bertoldi F., McKee C. F., 1992, ApJ, 395, 140 Binney J., Tremaine S., 2008, Galactic Dynamics. Princeton Univ. Press, Princeton, NJ Bonnell I. A., Bate M. R., Clarke C. J., Pringle J. E., 2001, MNRAS, 323, 785 Bontemps S., Andre P., Terebey S., Cabrit S., 1996, A&A, 311, 858 Buckle J. V., et al., 2009, MNRAS, 399, 1026 Buckle J. V., et al., 2010, MNRAS, 401, 204 Burkert A., Bodenheimer P., 2000, ApJ, 543, 822 Caselli P., Myers P. C., 1995, ApJ, 446, 665 Caselli P., Benson P. J., Myers P. C., Tafalla M., 2002, ApJ, 572, 238 Clark P. C., Klessen R. S., Bonnell I. A., 2007, MNRAS, 379, 57 Curtis E. I., 2009, PhD thesis, Univ. of Cambridge Curtis E. I., Richer J. S., 2010, MNRAS, 402, 603 Curtis E. I., Richer J. S., Buckle J. V., 2010a, MNRAS, 401, 455 (Paper I) Curtis E. I., Richer J. S., Swift J. J., Williams J. P., 2010b, MNRAS, in press (Paper II) ´ Dib S., Kim J., Vazquez-Semadeni E., Burkert A., Shadmehri M., 2007, ApJ, 661, 262 Dib S., Hennebelle P., Pineda J. E., Csengeri T., Bontemps S., Audit E., Goodman A. A., 2010, ApJ, in press Duquennoy A., Mayor M., 1991, A&A, 248, 485
c 2010 RAS, MNRAS 000, 1­17


The kinematics of SCUBA clumps in Perseus
Elmegreen B. G., Efremov Y., Pudritz R. E., Zinnecker H., 2000, in Mannings V., Boss A. P., Russell S. S., eds, Protostars and Planets IV. Univ. of Arizona Press, Tucson, p. 179 Elmegreen B. G., Scalo J., 2004, ARA&A, 42, 211 Enoch M. L., et al., 2006, ApJ, 638, 293 Enoch M. L., Evans II N. J., Sargent A. I., Glenn J., Rosolowsky E., Myers P., 2008, ApJ, 684, 1240 Galli D., Walmsley M., Goncalves J., 2002, A&A, 394, 275 ¸ Gammie C. F., Lin Y.-T., Stone J. M., Ostriker E. C., 2003, ApJ, 592, 203 Goodman A. A., Benson P. J., Fuller G. A., Myers P. C., 1993, ApJ, 406, 528 (GBF93) Goodman A. A., Barranco J. A., Wilner D. J., Heyer M. H., 1998, ApJ, 504, 223 Goodwin S. P., Kouwenhoven M. B. N., 2009, MNRAS, 397, L36 Goodwin S. P., Kroupa P., Goodman A., Burkert A., 2007, in Reipurth B., Jewitt D., Keil K., eds, Protostars and Planets V. Univ. of Arizona Press, Tucson, p. 133 Goodwin S. P., Nutter D., Kroupa P., Ward-Thompson D., Whitworth A. P., 2008, A&A, 477, 823 Hartmann L., 2001, AJ, 121, 1030 Hatchell J., Fuller G. A., 2008, A&A, 482, 855 Hatchell J., Richer J. S., Fuller G. A., Qualtrough C. J., Ladd E. F., Chandler C. J., 2005, A&A, 440, 151 Hatchell J., Fuller G. A., Richer J. S., Harries T. J., Ladd E. F., 2007, A&A, 468, 1009 (HFR07) Jappsen A.-K., Klessen R. S., 2004, A&A, 423, 1 (JK04) Jijina J., Myers P. C., Adams F. C., 1999, ApJS, 125, 161 Jones C. E., Basu S., Dubinski J., 2001, ApJ, 551, 387 Kirk H., Johnstone D., Di Francesco J., 2006, ApJ, 646, 1009 Kirk H., Johnstone D., Tafalla M., 2007, ApJ, 668, 1042 (KJT07) ´ Klessen R. S., Ballesteros-Paredes J., Vazquez-Semadeni E., ´ Duran-Rojas C., 2005, ApJ, 620, 786 Krumholz M. R., Tan J. C., 2007, ApJ, 654, 304 Krumholz M. R., McKee C. F., Klein R. I., 2005, Nat, 438, 332 Larson R. B., 1981, MNRAS, 194, 809 Larson R. B., 2003, Rep. of Progress in Phys., 66, 1651 Li P. S., Norman M. L., Mac Low M.-M., Heitsch F., 2004, ApJ, 605, 800 Lissauer J. J., 1993, ARA&A, 31, 129 Mac Low M.-M., Klessen R. S., 2004, Rev. of Modern Phys., 76, 125 MacLaren I., Richardson K. M., Wolfendale A. W., 1988, ApJ, 333, 821 McKee C. F., 1999, in Lada C. J., Kylafis N. D., eds, The Origin of Stars and Planetary Systems. Kluwer, Dordrecht, p. 29 Motte F., Andre P., Neri R., 1998, A&A, 336, 150 Mouschovias T. C., 1987, in Morfill G. E., Scholer M., eds, Physical Processes in Interstellar Clouds. Reidel, Dordrecht, p. 453 Myers P. C., 1983, ApJ, 270, 105 Nakamura F., Li Z.-Y., 2005, ApJ, 631, 411 Nakamura F., Li Z.-Y., 2007, ApJ, 662, 395 Offner S. S. R., Klein R. I., McKee C. F., 2008a, ApJ, 686, 1174 Offner S. S. R., Krumholz M. R., Klein R. I., McKee C. F., 2008b, AJ, 136, 404 Ohashi N., Hayashi M., Ho P. T. P., Momose M., Tamura M., Hirano N., Sargent A. I., 1997, ApJ, 488, 317 Olmi L., Testi L., Sargent A. I., 2005, A&A, 431, 253 ° Padoan P., Juvela M., Goodman A. A., Nordlund A., 2001, ApJ, 553, 227 Pineda J. E., Goodman A. A., Arce H. G., Caselli P., Foster J. B., Myers P. C., Rosolowsky E. W., 2010, ApJ, 712, L116
c 2010 RAS, MNRAS 000, 1­17

17

Pirogov L., Zinchenko I., Caselli P., Johansson L. E. B., Myers P. C., 2003, A&A, 405, 639 Redman M. P., Keto E., Rawlings J. M. C., Williams D. A., 2004, MNRAS, 352, 1365 Rohlfs K., Wilson T. L., 2004, Tools of Radio Astronomy. Springer, Berlin Rosolowsky E. W., Pineda J. E., Foster J. B., Borkin M. A., Kauffmann J., Caselli P., Myers P. C., Goodman A. A., 2008, ApJS, 175, 509 Ruden S. P., 1999, in Lada C. J., Kylafis N. D., eds, The Origin of Stars and Planetary Systems. Kluwer, Dordrecht, p. 643 Schnee S., Caselli P., Goodman A., Arce H. G., BallesterosParedes J., Kuchibhotla K., 2007, ApJ, 671, 1839 Shu F. H., Adams F. C., Lizano S., 1987, ARA&A, 25, 23 Simon M., 1992, in McAlister H. A., Hartkopf W. I., eds, IAU Colloq. 135, Complementary Approaches to Double and Multiple Star Research. Astron. Soc. Pac., San Francisco, p. 41 Solomon P. M., Rivolo A. R., Barrett J., Yahil A., 1987, ApJ, 319, 730 Spitzer L., 1978, Physical processes in the interstellar medium. Wiley-Interscience, New York, NY ¨ Stutzki J., Gusten R., 1990, ApJ, 356, 513 Swift J. J., Williams J. P., 2008, ApJ, 679, 552 Tachihara K., Onishi T., Mizuno A., Fukui Y., 2002, A&A, 385, 909 Tilley D. A., Pudritz R. E., 2004, MNRAS, 353, 769 ^ Walmsley C. M., Flower D. R., Pineau des Forets G., 2004, A&A, 418, 1035 Williams J. P., de Geus E. J., Blitz L., 1994, ApJ, 428, 693 Williams J. P., Blitz L., McKee C. F., 2000, in Mannings V., Boss A. P., Russell S. S., eds, Protostars and Planets IV. Univ. of Arizona Press, Tucson, p. 97 Zuckerman B., Evans II N. J., 1974, ApJ, 192, L149

APPENDIX A: DERIVED CORE DYNAMICAL AND ROTATIONAL PROPERTIES We list linewidths, virial masses and rotational parameters alongside other properties for every clump examined in Tables A1, A2, A3 and A4. In Figs. A1 and A2, we plot the C18 O J = 3 2 line centre velocity and best-fitting velocity gradient for every SCUBA clump with a significant detection.


18

E. I. Curtis and J. S. Richer

Table A1. Dynamical properties of SCUBA clumps found with C L FI N D. The full version of this table is available as Supporting Information to the online version of this article. Sub-region
a

Clump ID 1 2 3 4 5

b

Hatchell ID 42 43 44 41 45

c

TNH3 (K) 13.5 16.3 16.5 15.0 16.4

d

vc (C18 O) (km s-1 ) 7.36 8.24 7.62 7.48 7.79

e

vC18 (km s

f O -1

)

NT g (km s-1 ) 0.81 0.64 0.50 0.77 0.59

Mvir h (M ) 11.3 11.9 12.3 13.7 22.9

NGC 1333 NGC 1333 NGC 1333 NGC 1333 NGC 1333 ...
a b

1.91 1.53 1.19 1.81 1.40

Name of sub-region map. Clump number from Curtis & Richer (2010). c H07 identifier. d Rosolowsky et al. (2008) NH kinetic temperature where they exist or 10.0, 15.0 and 12.0 K for starless clumps, protostars and clumps with no H07 3 identification where not respectively. e, f Measured C18 O J = 3 2 line centre velocity, v and FWHM, v c C18 O from a Gaussian fit to the line profile at the clump peak. g Non-thermal contribution to the linewidth estimated using the temperatures of column (d). h Virial masses calculated using the temperatures of column (d).

Table A2. Dynamical properties of SCUBA clumps found with G AU S S C L U M P S. The full version of this table is available as Supporting Information to the online version of this article. Sub-region
a

Clump ID 1 2 3 4 5

b

Hatchell ID 41,42 43 ­ 45 46

c

TNH3 (K) 15.0 16.3 16.5 16.4 14.3

d

vc (C18 O) (km s-1 ) 6.77 8.18 7.62 7.82 8.40

e

vC18 (km s

f O -1

)

NT g (km s-1 ) 0.44 0.64 0.50 0.57 0.54

Mvir h (M ) 11.9 5.8 8.7 4.0

NGC 1333 NGC 1333 NGC 1333 NGC 1333 NGC 1333 ...
a b

1.04 1.51 1.19 1.35 1.27

Name of sub-region map. Clump number from Curtis & Richer (2010). c H07 identifier. d Rosolowsky et al. (2008) NH kinetic temperature where they exist or 10.0, 15.0 and 12.0 K for starless clumps, protostars and clumps with no H07 3 identification where not respectively. e, f Measured C18 O J = 3 2 line centre velocity, v and FWHM, v c C18 O from a Gaussian fit to the line profile at the clump peak. g Non-thermal contribution to the linewidth estimated using the temperatures of column (d). h Virial masses calculated using the temperatures of column (d).
30

17:00

30

15:00

30

Declination 31:14:00 30

Declination

Declination

31:16:00

30

15:00

06

04

02 3:29:00 Right ascension

58

28:56

13:00

18:30

19:00

30

30

31:20:00

30

21:00

3:29:00

58

56

54 52 Right ascension

28:50

48

46

06

04

02

3:29:00

58

56

54

28:52

Right ascension

NGC 1333 2

NGC 1333 3

NGC 1333 5

Figure A1. SCUBA 850 µm clumps identified with C L FI N D containing significant velocity gradients, G /G 3 in NGC 1333, IC348 and L1448. The colourscale denotes the C18 O J = 3 2 line centre velocity, linearly scaled from the clump minimum (blue/black) to maximum (red/white). Contours of 850 µ m flux density are overlaid (black) at 100, 200, 400, 800, 1600 and 3200 mJy beam-1 . A white/black circle marks the position of the clump peak emission with a white line positioned to show the direction of the fitted velocity gradient. The full version of this figure is available as Supporting Information to the online version of this article.

c 2010 RAS, MNRAS 000, 1­17


The kinematics of SCUBA clumps in Perseus

19

Table A3. Rotational properties of SCUBA clumps located with C L FI N D. These properties have been found by non-linear least squares fitting of a linear velocity gradient to the first moment of the C18 O J = 3 2 data. Only significant detections are shown with G /G 3. The full version of this table is available as Supporting Information to the online version of this article. Sub-region NGC1333 NGC1333 NGC1333 NGC1333 NGC1333 ...
a b c e f g h i a

Clump ID 2 3 5 6 8

b

Hatchell class I 0 I 0 I

c

Number d of points 1117 2106 2229 978 330

G ± G e (km s-1 pc-1 ) 3.16 ± 0.76 4.85 ± 0.35 1.54 ± 0.41 4.66 ± 0.84 10.63 ± 2.67

G f (deg E of N) 155.44 26.88 21.10 125.79 -58.59

G /G 4 14 4 6 4

g



rot

h

ji (km s-1 pc) 1.05E-03 3.51E-03 1.63E-03 1.44E-03 9.83E-04

6.98E-04 5.39E-03 8.91E-04 3.35E-03 4.67E-03

Name of sub-region map. Clump identification from Curtis & Richer (2010). H07 clump class, S=starless, 0=Class 0 protostar and I=Class I protostar. d Number of points fitted in the map. Fitted velocity gradient and its error. Angle of the fitted velocity gradient, east of north. Level of significance of the fit, G /G . Ratio of the kinetic to gravitational energy. Specific angular momentum.

Table A4. Rotational properties of SCUBA clumps located with G AU S S C L U M P S. These properties have been found by non-linear least squares fitting of a linear velocity gradient to the first moment of the C18 O J = 3 2 data. Only significant detections are shown with G /G 3. The full version of this table is available as Supporting Information to the online version of this article. Sub-region NGC1333 NGC1333 NGC1333 NGC1333 NGC1333 ...
a b c e f g h i a

Clump ID 1 4 7 8 9

b

Hatchell class 0 I I 0 0

c

Number d of points 800 922 1222 812 1374

G ± G e (km s-1 pc-1 ) 5.33 4.93 6.73 4.25 5.37 ± ± ± ± ± 1 1 0 1 0 .33 .16 .84 .11 .65

G f (deg E of N) 135.67 -79.61 -77.95 26.44 -38.35

G /G 4 4 8 4 8

g



rot

h

ji (km s-1 pc) 7.86E-04 7.41E-04 1.56E-03 5.54E-04 1.63E-03

6.90E-04 1.22E-03 4.27E-03 1.07E-03 3.70E-03

Name of sub-region map. Clump identification from Curtis & Richer (2010). H07 clump class, S=starless, 0=Class 0 protostar and I=Class I protostar. d Number of points fitted in the map. Fitted velocity gradient and its error. Angle of the fitted velocity gradient, east of north. Level of significance of the fit, G /G . Ratio of the kinetic to gravitational energy. Specific angular momentum.

c 2010 RAS, MNRAS 000, 1­17


20

E. I. Curtis and J. S. Richer
10 14:00

21:00

10

40

30

Declination

40

50

Declination

Declination

20

31:13:00

10

20

30

12:30

19:50

31:20:00

40

15

14

13

12 11 Right ascension

3:29:10

09

08

06

05

04

03

02

01

3:29:00

59

58

28:57

15:00

50

10

30

31:16:00

50

30

11

3:29:10

09

08

Right ascension

07 06 Right ascension

05

04

03

02

NGC 1333 1

NGC 1333 4

NGC 1333 7

Figure A2. SCUBA 850 µm clumps identified with G AU S S C L U M P S containing significant velocity gradients, G /G 3 in NGC1333, IC348 and L1448. The colour-scale denotes the C18 O J = 3 2 line centre velocity, linearly scaled from the clump minimum (blue/black) to maximum (red/white). Contours of 850 µ m flux density are overlaid (black) at 100, 200, 400, 800, 1600 and 3200 mJy beam-1 . A white/black circle marks the position of the clump peak emission with a white line positioned to show the direction of the fitted velocity gradient. The full version of this figure is available as Supporting Information to the online version of this article.

c 2010 RAS, MNRAS 000, 1­17