Issue |
A&A
Volume 535, November 2011
|
|
---|---|---|
Article Number | A77 | |
Number of page(s) | 10 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/201117132 | |
Published online | 11 November 2011 |
The formation of active protoclusters in the Aquila rift: a millimeter continuum view⋆
1
ESO, Karl Schwarzschild Strasse 2, 85748 Garching bei München, Germany
e-mail: amaury@eso.org
2
Laboratoire AIM, CEA/DSM-CNRS-Université Paris Diderot, IRFU/SAp, CEA Saclay, 91191 Gif-sur-Yvette, France
3
CNRS/INSU, Laboratoire d’Astrophysique de Bordeaux, UMR5804, BP 89, 33271 Floirac Cedex, France
Received: 25 April 2011
Accepted: 30 July 2011
While most stars are believed to form in stellar clusters, the formation and early evolution of young stellar clusters is still largely unknown. Improving our knowledge of the earliest phases of clustered star formation is crucial for understanding the origin of the stellar initial mass function and the efficiency of the star formation process, which both play a key role in the evolution of galaxies. Here, we present an analysis of the Aquila rift complex which addresses the questions of the star formation rate (SFR), star formation efficiency (SFE) and typical lifetime of the Class 0 protostellar phase in two nearby cluster-forming clumps: the Serpens South and W40 protoclusters. We carried out a 1.2 mm dust continuum mapping of the Aquila rift complex with the MAMBO bolometer array on the IRAM 30 m telescope. Using a multi-scale source extraction method, we perform a systematic source extraction in our millimeter continuum map. Based on complementary data from the Herschel Gould Belt survey and Spitzer maps, we characterize the spectral energy distributions (SEDs) of the 77 mm continuum sources detected with MAMBO and estimate their evolutionary stages. Taking advantage of the comprehensive dataset available for the Serpens South region, spanning wavelengths from 2 μm to 1.2 mm, we estimate the numbers of young stellar objects (YSOs) at different evolutionary stages and find a ratio of Class 0 to Class I protostars N(0)/N(I) = 0.19−0.27. This low ratio supports a scenario of relatively fast accretion at the beginning of the protostellar phase, and leads to a Class 0 lifetime of ~4−9 × 104 yr. We also show that both the Serpens South and W40 protoclusters are characterized by large fractions of protostars and high SFRs ~ 20−50M⊙ Myr-1 pc-2, in agreement with the idea that these two nearby clumps are active sites of clustered star formation currently undergoing bursts of star formation, and have the potential ability to form bound star clusters. While the formation of these two protoclusters is likely to have been initiated in a very different manner, the resulting protostellar populations are observed to be very similar. This suggests that after the onset of gravitational collapse, the detailed manner in which the collapse has been initiated does not affect much the ability of a clump to form stars.
Key words: stars: formation / circumstellar matter / ISM: clouds / ISM: structure / dust, extinction / stars: protostars
© ESO, 2011
1. Introduction
The formation and early evolution of protoclusters While most of the young pre-main sequence (PMS) stars in molecular clouds are found in groups and clusters rather than being uniformly distributed throughout the clouds (Testi et al. 1999; Carpenter 2000; Lada & Lada 2003, see also Bressert et al. 2010 for a discussion of clustered star formation), the formation and early evolution of stellar clusters are still poorly known. Improving our knowledge of the earliest phases of clustered star formation holds the key to understanding the distribution of stellar masses at birth since the similarity of the stellar initial mass function (IMF) of young stellar clusters with the field stars IMF may imply it is linked with the cluster formation process itself (Palla & Stahler 1999; Hillenbrand & Carpenter 2000; Meyer et al. 2000).
Protoclusters are parsec-scale gas-dominated cluster-forming clumps, which contain a large amount of dense gas material (several Jeans masses), and are actively forming protostars, as testified by the significant population of prestellar cores and embedded Class 0 objects detected toward them at (sub-)millimeter wavelengths (see Motte et al. 1998; Hurt & Barsony 1996; Myers 2001; Testi & Sargent 1998). While protoclusters are the direct progenitors of infrared embedded clusters, only a small fraction of them eventually result in bound star clusters, as it has been suggested (Lada & Lada 2003) that most embedded clusters are not able to remain bound when they rapidly lose their gas mass and evolve into revealed clusters. Their stellar component therefore disperse to become field stars. This high infant mortality has been proposed to be due to insufficient star formation rates (SFR) to fight against fast gas loss (Goodwin 2009; Baumgardt & Kroupa 2007), but may also rely on the dynamical state of the young stars and left-over material during the later evolution of the embedded cluster (Smith et al. 2011).
With ages of only a few ~105 yr in some cases (Peretto et al. 2007), protoclusters are therefore the youngest precursors of star clusters. Their youth makes them ideal targets to investigate the formation and early evolution mechanisms of clustered star formation, since their properties are likely to still reflect the physics of their formation. Deep (sub-)millimeter continuum observations of cluster-forming clumps can set strong constraints on the early properties and conditions for clustered star formation (cf. Motte et al. 1998; Testi & Sargent 1998), the efficiency of the conversion of gas into stars from well-defined mass reservoirs (cf. Enoch et al. 2007; Belloche et al. 2011), and the typical lifetimes of the youngest embedded objects, i.e., Class 0 protostars (André et al. 2000; Visser et al. 2002; Enoch et al. 2009). Moreover, only (sub-)millimeter line observations can probe the dynamics of cluster-forming clumps (Peretto et al. 2006; André et al. 2007; Kirk et al. 2007), as well as the role of feedback from embedded young stellar objects (YSOs) on the early evolution of young protoclusters (Knee & Sandell 2000; Ginsburg et al. 2009; Maury et al. 2009; Nakamura et al. 2011a).
Several observational studies targeting both low-mass and intermediate-mass protoclusters have shown that embedded clusters form stars more efficiently than global molecular clouds. Indeed, typical star formation efficiencies (SFEs) of 10–40% (Wilking & Lada 1983; Olmi & Testi 2002; Lada & Lada 2003) have been estimated for embedded clusters, while molecular clouds have SFEs of only a few percent (Myers et al. 1986; Evans et al. 2009; Lada et al. 2010). The SFE and SFR are known to vary from region to region, but few studies have focused on comparing these properties for protoclusters formed from different initial conditions. Moreover, past estimations of SFEs and SFRs have generally relied on only partial population counts and limited source classification, as it is very difficult to detect and characterize complete protostellar populations with limited wavelength coverage. Indeed, in order to probe complete YSO samples from Class 0 protostars to Class II objects and properly determine their evolutionary stages, deep imaging observations from ~2 μm to millimeter wavelengths, including the submillimeter domain, are required.
The Herschel Space Observatory now makes it possible to bridge the submillimeter gap by providing deep, wide-field dust continuum imaging at 6 wavelengths from 70 μm to 500 μm. In particular, the Herschel Gould Belt survey (André et al. 2010) is imaging the bulk of the nearest (d ≤ 0.5 kpc) molecular cloud complexes with the SPIRE and PACS cameras and will eventually provide a complete census of Class 0 protostars and prestellar cores in these clouds.
Here, we present a millimeter continuum study of the Aquila rift complex which addresses the question of the lifetime of the Class 0 phase and that of the SFR in two nearby cluster-forming clumps: the Serpens South and W40 protoclusters.
The Aquila rift complex The Aquila rift molecular cloud complex corresponds to a large extinction feature (see Bontemps et al. 2010), located above the Galactic plane at l = 28°, only ~3° south of the well-known Serpens Main star-forming region. The Aquila rift extinction-defined area (see Bontemps et al. 2010) mainly harbours two known sites of star formation: Serpens South is the western young embedded cluster revealed by infrared Spitzer observations (Gutermuth et al. 2008), and W40 is the eastern cluster associated to the eponymous HII region, also known as Sharpless 2-64 (Smith et al. 1985; Vallee 1987).
A recent Herschel mapping of ~a 3.3° × 3.3° field was carried out toward this complex as part of the Gould Belt survey (André et al. 2010). The corresponding far-infrared maps have unprecedented spatial resolution and sensitivity in 5 bands between 70 and 500 μm. The first results of these Herschel observations of the Aquila complex are presented in Bontemps et al. (2010), and Könyves et al. (2010). Briefly, the Herschel maps revealed a large number of new protostars and prestellar cores mainly spread over three currently active star-forming sites: the filamentary molecular cloud harboring the Serpens South infrared embedded cluster in its center, the eastern HII region W40/Sh2-64, and the southernmost MWC297 / Sh2-62 region.
Distance of the star forming clouds in the Aquila rift The distances of the Serpens South and W40 protoclusters, both located inside the extinction feature of the Aquila rift, are still uncertain.
Whether or not the Aquila rift complex and the Serpens Main cloud are part of the same complex and therefore at the same distance is still a matter of debate (Gutermuth et al. 2008; Dzib et al. 2010; Loinard et al. 2011), however from the extinction maps obtained by Bontemps et al. (2010), their respective extinction features are seen as clearly distinguishable separate regions. Indeed, assuming that Serpens Main and Aquila rift are different entities also has the advantage of explaining the discrepancy found between their extinction-based distance estimates (Straižys et al. 1996, 2003; Bontemps et al. 2010; Knude 2011): ~ 200–280 pc for the Aquila rift complex (due to the known extinction wall at a distance of 225 ± 55pc), and ~415 pc for Serpens Main (from the recent VLBI-based distance by Dzib et al. 2010), suggesting that Serpens Main is located behind the extinction wall associated to Serpens South.
Previous studies targeting the Serpens South protocluster (Gutermuth et al. 2008; Bontemps et al. 2010; Könyves et al. 2010) used the extinction-related distance of ~260 pc for this star forming region, Serpens South being the highest extinction region found inside the Aquila rift complex. Also based on extinction measurements (Smith et al. 1985) and radio spectral line data (Vallee 1987; Vallee et al. 1992), the nearby HII region W40 has usually been considered to be at a distance ranging from 300 to 700 pc. A recent Chandra study by Kuhn et al. (2010) emphasizes, on the basis of X-ray luminosity function arguments, a most-likely distance of ~600 pc for the W40 cluster, although not excluding the possibility of a smaller (~300 pc) distance.
We argue that both Serpens South and W40 are likely to be at the same distance because (i) they belong to the same continuous extinction feature (see Bontemps et al. 2010; Schneider et al. 2011) associated to the Aquila rift; and (ii) they are populating a continuum of LSR velocities between 4 and 10 km s-1.
For consistency with the Herschel Gould Belt survey (André et al. 2010) results for that region (Bontemps et al. 2010; Könyves et al. 2010) and easier comparison with previous studies, we adopt a distance of 260 pc for the Aquila rift complex, including both Serpens South and W40, throughout this paper.
2. MAMBO 1.2 mm dust continuum observations and data reduction
Fig. 1 1.2 mm dust continuum MAMBO mosaic of the Aquila rift complex, smoothed to an effective angular resolution of 13. Due to the high dynamic range of fluxes in the map, and to enhance visibility of the structure of emission in the whole area mapped, the color scale chosen here has a maximum of 400 mJy/beam. The highest flux density of order ~1400 mJy/beam is found in the center of the Serpens South protocluster (see the blow-up in Fig. 2). In this smoothed map, the mean rms noise level is ~13 mJy/13-beam. |
We used the MAx-Planck Millimeter BOlometer array (MAMBO) on the IRAM 30-m telescope near Pico Veleta, (Spain) during March 2007 to carry out a millimeter continuum mapping of the Aquila rift complex. Ten adjacent sub-fields were selected thanks to the near-IR extinction map derived from 2MASS data (Bontemps et al. 2010; Schneider et al. 2011), in order to image the dust continuum associated with the main extinction features of the Aquila rift complex, at 1.2 mm (~240 GHz). These sub-fields were mapped using the on-the-fly technique in the double-beam switching mode, i.e. chopping was performed by switching the wobbling secondary with a throw ranging from 62 to 160. The scanning speed was 8/s. The atmospheric optical depth at zenith was estimated between 0.09 and 0.41 at λ = 1.2 mm. Pointing and focus checks were made every 3 h, resulting in a pointing accuracy better than 3, and the calibration obtained by mapping Mars (primary calibrator) results in a flux accuracy better than ~15%. The angular resolution with IRAM 30-m antenna is ~11 (HPBW) at λ = 1.2 mm, and the average rms achieved is ~15 mJy/11-beam. We reduced and co-added all the MAMBO sub-fields using standard methods in MOPSIC, part of the GILDAS1 software, and produced a global mosaic covering a total area of ~60 × 50 (see Fig. 1).
3. Dust continuum mapping: results and analysis
3.1. Main features of the millimeter continuum maps
The large-scale 1.2 mm mosaic obtained with the MAMBO camera is shown in Fig. 1. Four main parsec-scale regions of interest can be identified in this mosaic, which are labelled 1, 2, 3, and 4. Indeed, these regions show both extended emission and compact structures in dust continuum, and are therefore likely tracing active star-forming areas.
Region labelled 1 is the previously known Serpens South protocluster, showing up as an elongated ~1 pc filament in the 1.2 mm dust continuum map (see blow-up in Fig. 2), coincident with the Spitzer near-infrared absorption feature discussed by Gutermuth et al. (2008). Very condensed in its center, the filament sees its density decreasing while going along the two arms that have a cross-shape in more diffuse emission. Compact sources are found up to 15 away from the center of the filament. This high-density region (column densities up to 3 × 1023 cm-2) at the convergence zone of the filament shows several dust continuum peaks at the MAMBO resolution: it is the center of a protocluster, with several compact millimeter sources tracing very young protostars at the beginning of their accretion phase (separated by less than 5000 AU), as well as numerous more evolved infrared YSOs.
Region labelled 2 is shown with much greater details in Fig. 3. It mainly features a ring-like structure following the outer border of the W40 HII region (Smith et al. 1985; Vallee 1987; Kuhn et al. 2010), and spatially coincident with the CO emission mapped by Zhu et al. (2006). In addition to this western dense shell, the map shows elongated structures (also containing compact sources) in the eastern part of the photo-dissociated region, while the central part of the infrared cluster containing most of the pre-main sequence objects (Kuhn et al. 2010) does not show any evidence for cold dust emission at 1.2 mm.
Finally, regions labelled 3 and 4 are more diffuse, less rich star-formation sites, for which no previous literature was found.
The ~11 angular resolution of MAMBO at 1.2 mm allowed us to resolve individual condensations within the dust cores at the assumed distance of Aquila rift, and we therefore used our millimeter data to complement the Herschel dataset and characterize the protostellar population in the mapped area.
Fig. 2 Blow-up of the Region 1 from the MAMBO 1.2 mm dust continuum map shown in Fig. 1, around the Serpens South filamentary cloud. The mean rms noise level is ~12 mJy/13-beam. Black contours show the levels of 5σ, 8σ, 15σ, 30σ, 70σ, 90σ, and 120σ (1.44 Jy/beam). |
Fig. 3 Blow-up of the Region 2 from the 1.2 mm continuum map shown in Fig. 1, around the HII region W40. The mean rms noise level is ~14 mJy/13-beam. Black contours show the levels of 5σ, 8σ, 15σ, 20σ, 25σ, and 35σ. One of the candidate exciting stars, Sh64-OS2a, is detected in our MAMBO mapping and indicated here. |
3.2. Source extraction and properties
In order to improve the detection of compact sources in our field, we smoothed the original 1.2 mm continuum map to a 13-beam, and performed a systematic source detection on the MAMBO map using the new source extraction method getsources2 (Men’shchikov et al. 2010, and in prep.). This analysis led to the extraction of a total of 77 dust continuum sources in the whole region shown in Fig. 1, among which 25 are located in the Serpens South filament, while 36 lie in the W40 region. The resulting MAMBO compact sources catalog is reported in Table 1. Starting from this 1.2 mm dust continuum set of sources, parallel measurements were performed at other wavelengths with getsources in order to estimate the counterparts of the MAMBO sources in both the Herschel (from 70 μm to 500 μm), the Spitzer (8 μm and 24 μm) and the 2MASS K-band maps, when data was available. All these sources were also examined individually by eye, ensuring association to the right counterparts, which is not always straightforward while working with different sizes of beams in crowded regions, such as Serpens South or W40.
Characteristics of the MAMBO millimeter sources in the Aquila rift region.
continued.
We used these flux measurements to establish spectral energy distributions (SEDs) for all the 77 MAMBO sources, and derived other source properties (dust temperature, core/envelope mass, bolometric temperature and luminosity) when a good fit was found to reproduce the SED obtained between 2 μm and 1.2 mm with grey-body functions. The condensation masses (core masses for starless sources and envelope masses – Menv – for protostellar objects) and dust temperatures Tdust were estimated from the fitting of these SEDs by grey-body functions (Könyves et al. 2010). We adopted the recommended values (Ossenkopf & Henning 1994; Henning et al. 1995) for the dust opacity per unit mass column density (of dust and gas): κ1.2mm = 0.01cm2g-1 for protostellar envelopes, and κ1.2mm = 0.005cm2g-1 for starless clumps. The condensations have masses ranging from ~0.03M⊙ to ~6M⊙, with an average factor ~2 absolute uncertainty (see Table 1). Note also that for the most evolved sources (mainly Class I protostars), these masses might be contaminated or even dominated by disk material, while they represent the core masses for the starless sources.
The MAMBO sources showing no counterpart shortward of 160 μm with Herschel-PACS, but detected in at least two Herschel-SPIRE bands (250, 350, or 500 μm) were classified as starless sources. The MAMBO sources with compact counterparts at Herschel wavelengths <160μm were classified as candidate YSOs. The bolometric luminosities Lbol, submillimeter luminosities , and bolometric temperatures Tbol of the candidate YSOs were estimated by integrating their SEDs, and are reported in Table 1. Note that the values of and Lbol derived here are expected to be more accurate (estimated uncertainty ≤ 40%) than previous estimates from the literature, thanks to the Herschel data providing a good coverage of the far-infrared/submillimeter domain.
3.3. Evolutionary classification of protostellar sources
Fig. 4 Envelope mass versus bolometric luminosity diagram comparing the locations of the protostars detected with MAMBO in the Aquila rift region (open triangles for the Serpens South sub-field, open squares for the other sub-fields) with the positions of known low-mass Class I (filled circles) and Class 0 objects (open circles) in Ophiuchus, Perseus, and Orion (André et al. 2000; Motte & André 2001). Numbers reported inside triangles are the reference numbers of the sources reported in the last column of Table 1. Model protostellar evolutionary tracks, computed for various final stellar masses (indicated above the Lbol axis) assuming the accretion/luminosity history described in the text, are superimposed (4 blue tracks representing protostellar evolution for objects with final masses 0.1, 0.3, 1 and 3M⊙). Evolution proceeds from the upper left to the lower right as indicated by two arrows on each track, plotted when 50% and 90% of the final stellar mass has been accreted, respectively. The straight lines show two Menv–Lbol relations marking the conceptual border zone between the Class 0 (Menv > M ⋆ /ϵ) and the Class I (Menv < M ⋆ /ϵ) stage, where ϵ is the local star formation efficiency taken here to be 50%. The dashed line is such that Menv ∝ Lbol (cf. André et al. 1993; André & Montmerle 1994), while the dotted relation follows as suggested by the accretion scenario adopted in the tracks (see also Bontemps et al. 1996). The typical uncertainty in the plotted values of Lbol is ~30–40%, while the uncertainty in the Menv values is a factor of ~2 (see Sect. 3.2). The black arrow in the upper right indicates how the positions of the W40 protostars would move in the diagram if a distance of 600 pc was adopted for W40, instead of 260 pc. |
All of the candidate embedded protostellar objects detected with MAMBO were then placed in the envelope mass versus bolometric luminosity (Menv–Lbol) diagram shown in Fig. 4, in order to distinguish between envelope-dominated (Menv ≫ M ⋆ ) Class 0 protostars and more evolved (Class I) protostars. In this diagram, the bolometric luminosity of a protostellar object can be used as an indirect tracer of the stellar mass M ⋆ it has already accreted (André et al. 1993; André & Montmerle 1994), based on comparison with simple accretion model tracks (cf. Bontemps et al. 1996; Saraceno et al. 1996; André et al. 2000).
The specific evolutionary tracks shown in Fig. 4 are similar to the tracks used by André et al. (2008) and were computed assuming each protostar forms from a prestellar core of finite initial mass Mcore ≡ Menv(0) and has Lbol = GM ⋆ (t)Ṁacc(t)/R ⋆ (t) + L ⋆ (t), where R ⋆ (M ⋆ ) is the protostellar radius (mass) and L ⋆ the PMS birthline luminosity (Stahler 1988). The mass accretion rate and the envelope mass were assumed to be related by Ṁacc(t) = ϵMenv(t)/τ, where ϵ = 50% is the typical SFE for individual cores (cf. Matzner & McKee 2000) and τ = 105 yr is the characteristic timescale of protostellar evolution, leading to Ṁacc(t) and Menv(t) functions declining exponentially with time as suggested by the observed decline of outflow power from Class 0 to Class I (see Bontemps et al. 1996). Note that we adopt ϵ = 50% here as in André et al. (2008) instead of ϵ = 100% as in André et al. (2000), since recent comparisons of the prestellar core mass function with the stellar IMF suggest ϵ50% (e.g., Alves et al. 2007; Nutter & Ward-Thompson 2007; Könyves et al. 2010; André et al. 2010).
A Class 0 status was attributed to sources with Menv > M ⋆ /ϵ according to the diagram of Fig. 4, i.e., protostars which have accreted less than half of their final stellar mass. Conceptually, protostars at the border line between the Class 0 and the Class I stage are objects which have accreted exactly half of their final stellar mass, and which thus have Menv = M ⋆ /ϵ according to the tracks. In practice, sources falling into the “Border Zone” region located between the dashed and dotted lines in Fig. 4 have been classified as border-line Class 0/I objects in Table 1. Sources lying below this border zone in the diagram were classified as Class I or Class II, depending on their infrared spectral index αIR ≡ dlog λFλ/dlog λ between 2.2 μm and 24 μm (as estimated by us from independent getsources measurements at 2.2 μm, 8 μm, 24 μm when available – see Sect. 4).
The use of this evolutionary diagram (see Fig. 4) leads to the identification of 9 Class 0 (and 3 border-line Class 0/I) protostars in the Serpens South filament. Likewise, 8 Class 0 (and 4 border-line Class 0/I) protostars are found in the W40 region. The other sub-fields labelled 3 and 4 gather 8 Class 0 (and 2 border-line Class 0/I) protostars.
Two other criteria, []3 (e.g. André et al. 2000) and [Tbol < 70 K] (e.g. Myers et al. 1998) respectively, have been commonly used in the literature to select candidate Class 0 protostars among embedded YSOs. In order to assess the importance of uncertainties in our classification, we have applied these two alternative methods of classifying embedded YSOs to the Aquila rift sample and compared the results to those obtained above on the basis of the Menv – Lbol diagram.
We find (see Table 1) that using the criterion (which traces a ratio of envelope to stellar mass larger than unity – see André et al. 1993) to select Class 0 protostars leads to similar results to those obtained from the Menv – Lbol diagram. Only in a few cases does the use of this criterion lead to the selection of border-line Class 0/I sources (according to Menv – Lbol diagram) as candidate Class 0 protostars. For example, using the criterion leads to the selection of 12 Class 0 protostars in the Serpens South filament (compared to 9 Class 0 and 3 border-line sources when using the Menv – Lbol diagram), and 11 Class 0 protostars in the W40 region (compared to 8 Class 0 and 4 border-line sources when using the Menv – Lbol diagram). Based on the present study, we suggest that setting the border line at may better match the results obtained with the more conceptual classification based on the tracks in the Menv – Lbol diagram that we use here.
We find that using the Tbol < 70 K criterion for distinguishing Class 0 protostars from Class I objects (Chen et al. 1995; Myers et al. 1998) tends to systematically overestimate the number of candidate Class 0 objects, compared to the numbers obtained from both the Menv – Lbol diagram and the criterion. This is especially true in the “warm” W40 region where the estimates of Tbol are more uncertain due to confusion at mid-infrared wavelengths from the photo-dissociation region. For example, using the Tbol < 70 K criterion leads to the identification of 12 Class 0 protostars in the Serpens South filament, and 19 Class 0 protostars in the W40 protocluster (twice the number of Class 0 protostars found in this region if using the two other criteria presented above).
We therefore conclude, in agreement with Young & Evans (2005), that the observational criterion to distinguish Class 0 protostars making use of the ratio seems to be more reliable than the Tbol criterion, also because it is less sensitive to environment and geometry effects. Comparing the various methods of classifying protostars, which is roughly equivalent to taking into account the border-line Class 0/I protostars from the Menv – Lbol diagram, we find that the typical uncertainty in the number of Class 0 protostars is ~50% in W40, and ~30% in Serpens South.
4. Discussion and conclusions
4.1. Class 0 phase lifetime and star formation rates
Here we propose a statistical analysis of the properties of the protostellar populations in the two main regions showing active star formation from our MAMBO map: the Serpens South filamentary cloud and the W40 HII region.
4.1.1. The Serpens South protocluster
First, we take advantage of having a complete population of Class 0 protostars in a spatially well-defined star-forming cloud to derive an estimate of the lifetime of the Class 0 phase.
The Serpens South protocluster was observed with IRAC and MIPS onboard the Spitzer telescope, and the analysis of IRAC data by Gutermuth et al. (2008) revealed 37 Class II YSOs and 54 Class I protostars (including flat spectrum sources) in a region of 14 × 10 around the filament. Based on the publicly available MIPS-24 μm, IRAC-8 μm, and 2MASS/K-band data, we independently performed an extraction of MIPS-24 μm point-like sources with getsources in the 0.8 × 1.2 pc area around the Serpens South filament shown in Fig. 2, which revealed a total of 112 objects with an infrared spectral index αIR ≡ dlog λFλ/dlog λ > 0.3 between 2 μm and 24 μm, that we considered as candidate Class I protostars.
Fig. 5 Same as Fig. 4, but showing only the sources from the Serpens South cluster forming clump. Triangles show the locations of the MAMBO protostars (green triangles are Class 0 and red triangles are Class I protostars, while black triangles represent the border-line objects). Red circles show the locations of the Class I protostars detected by Spitzer and without MAMBO detection. The Class I sources without any MAMBO counterpart are indicated as upper-limit values for Menv, around the detection limit of 2.5 × 10-2M⊙. The green double arrow along the Menv axis shows the range of masses for the prestellar cores detected by Herschel, based on a preliminary analysis of the Gould Belt survey data, which are 85% complete down to 0.3M⊙ in this region (Könyves et al. 2010). The respective locations of the prestellar cores, Class 0 and Class I protostars in the Serpens South protocluster are roughly consistent with what is expected from a single population of objects evolving along the model evolutionary tracks shown in blue. |
As discussed in Sect. 3.2 above, our MAMBO mapping finds 9 candidate Class 0 protostars in the same area (not counting border-line Class 0/Class I sources). Our census of Class 0 protostars in Serpens South is sensitive to protostars with final masses down to ~0.1M⊙ (see evolutionary tracks in Fig. 5). According to the grid of protostellar models by Robitaille et al. (2007), Class I protostars with final masses are expected to have 24μm flux densities F24μm > 35 mJy at d = 260 pc (assuming a random distribution of inclination angles, i.e. computed using a mean value of ~57°). We will therefore consider here only the 48 ± 2 candidate Class I sources which have F24μm > 30−40 mJy from our extracted sample. Note that at such a 24μm flux density level, our Class I sample should not be contaminated by extragalactic sources (cf. Harvey et al. 2007). Among these 48 objects, 16 coincide with 21 IRAC Class I sources reported by Gutermuth et al. (2008), while 21 are new detections made possible thanks to the use of the 24 μm data. The Spitzer, Herschel, and MAMBO maps are very likely wide enough to probe complete samples of YSOs associated with the Serpens South filamentary clump, since there are no YSOs, neither high column-density gas, immediately outside this 0.8 × 1.2 pc area (shown in Fig. 2) in available larger Herschel, Spitzer, and MAMBO maps of the region.
Assuming that the SFR in Serpens-South has been constant over the duration of the protostellar phase, we can estimate the lifetime of the Class 0 phase by comparing the relative numbers of Class 0 and Class I sources in this region. Combining the MAMBO and Spitzer samples discussed here, the Serpens South cluster-forming filament harbors a total of 57 protostellar objects with expected final stellar masses >0.1 M⊙, and more precisely 84% of Class I and 16% of Class 0 objects. Adopting a lifetime of 0.2–0.5 Myr for the Class I phase (Greene et al. 1994; Evans et al. 2009), the ratio of Class 0 to Class I protostars we find, N(0)/N(I) = 0.19, suggests a typical lifetime of ~4−9 × 104 yr for the Class 0 phase. This result is intermediate between the early estimate of 1−3 × 104 yr by André et al. (2000) and the more recent estimate of 0.10–0.16 Myr by Evans et al. (2009).
Our analysis makes use of an unprecedented census of YSOs in Serpens South, using wavelengths from 2 μm to 1.2 mm and including the Herschel submillimeter domain. It should therefore provide a more complete population of embedded YSOs in Serpens-South and a more reliable assessment of their evolutionary states (via better determination of their envelope masses and bolometric luminosities) than earlier studies.
However, the lifetime we derive for the Class 0 phase is still subject to large uncertainties, partly due to remaining uncertainties on the lifetime of the Class I phase. Estimates of the latter are themselves highly dependent on assumptions concerning the lifetime of the Class II (or T Tauri) phase (see Evans et al. 2009).
The method chosen to distinguish Class 0 protostars from more evolved YSOs also impacts the determination of this lifetime, because different methods lead to slightly different source counts (see Sect. 3.3). However, the two favored methods to classify YSOs were compared in Sect. 3.3, and the differences they show (typical ~30% uncertainty on the number of Class 0 protostars detected, i.e. N(0)/N(I) ~ 0.19−0.27) are expected to result in a typical uncertainty of ~0.02 Myr in the Class 0 lifetime we find here, which is a factor of 2 lower than the main source of error due to the uncertainty in the Class I lifetime (leading to a ~0.05 Myr uncertainty on the Class 0 lifetime as shown above). Moreover, we do not claim our result to be universal, as our study is focused on only one star-forming region and does not address the question of environmental effects.
Even taking these caveats into account, we stress that our results tend to support a scenario of relatively fast accretion at the beginning of the protostellar phase in protoclusters (cf. Henriksen et al. 1997), since they suggest that a typical YSO accretes the first half of its final stellar mass in <0.1 Myr (Class 0 phase) and the second half in ~0.4 Myr (Class I phase – cf. Evans et al. 2009).
Finally, we remark, in Fig. 5, an over-abundance of low-luminosity Class I protostars if compared to what is expected from the population of Class 0 protostars and associated model evolutionary tracks in the Menv–Lbol diagram for Serpens South, which is likely to be an illustration of the well-known “luminosity problem” first noted by Kenyon et al. (1990) and more recently stressed by Evans et al. (2009). However, we point out that the locations of prestellar cores, Class 0 and Class I protostars in the Menv–Lbol diagram for Serpens South shown in Fig. 5, are roughly consistent with what expected for a single population of objects evolving along the model evolutionary tracks. This illustrates the validity of our method for determining the evolutionary status of protostars.
Second, using the column-density map of the Aquila rift derived from Herschel data (Könyves et al. 2010), we can estimate a mass for the Serpens South filamentary cluster-forming clump shown in Fig. 2, and derive a value of the SFE for this region. We find a total mass of ~610M⊙ in the 15′ × 25′ area containing the protostellar population discussed here, with a mean column density of NH2 ~ 3.1 × 1022cm-2. Considering that this region has formed 9 Class 0, 48 Class I, and 37 Class II objects, i.e. a total of 94 YSOs of mean mass 0.5 M⊙, we find: SFE = Mstar/(Mstar + Mgas) = 47/(47 + 610) ~7%. Assuming a standard duration of 2 Myr for going from the Class 0 to the Class II phase (Evans et al. 2009) and a projected area of ~1 pc2 for the region harboring YSOs, we derive the following SFR per unit surface for the Serpens South cluster-forming clump: SFR ~ 23M⊙ Myr-1 pc-2.
This value is significantly higher than the typical SFR values reported by Evans et al. (2009) for entire molecular clouds (SFR ~ 1.6M⊙ Myr-1 pc-2), which is consistent with the fact that we are focusing here on a compact cluster-forming clump with higher column densities (up to NH2 ~ 1023cm-2) than the clouds targeted by Evans et al. (2009), also showing evidence of a recent burst of star formation. A recent burst of star formation is indicated by the unusually large number of Class 0 and Class I protostars found in the Serpens South star-forming clump compared to the relatively small number of Class II YSOs, which was already pointed out by Gutermuth et al. (2008) and that we confirm here. Note that if Serpens-South has indeed undergone a recent burst of star formation, then the observed ratio of Class 0 to Class I protostars (N(0)/N(I)) may overestimate the actual ratio of Class 0 to Class I lifetimes.
Finally, the SFE value derived here for Serpens South is of the same order as that reported by Evans et al. (2009) for molecular clouds (i.e. SFE ~ 3−6%). The fact that Serpens South shows both a high SFR and a relatively low SFE may be interpreted as an additional sign of its very early evolutionary stage (Lada & Lada 2003).
4.1.2. The W40 HII region
The present MAMBO 1.2 mm continuum study allows us to identify 8 candidate Class 0 protostars in the W40 cluster region, where copious photo-dissociation processes (HII region probably powered by the massive young star Sh64-OS2a, see Fig. 3) make it difficult to obtain a census of YSOs at shorter wavelengths. Indeed, the Lyman-continuum luminosity towards the HII region of W40 (~1.5 × 1048 photons s-1) inferred from the radio continuum emission by Altenhoff et al. (1970) could be produced by a single O9 V star (Smith et al. 1985). In a recent study, Kuhn et al. (2010) found 93 YSOs with NIR K-band excess belonging to the W40 region, based on their Chandra observations toward the young cluster. However, their census did not distinguish between Class I, Class II, and later evolutionary stages. It is therefore difficult to derive accurate lifetimes for the Class 0 and/or Class I phases in the W40 region.
We estimate a gas mass of 310 M⊙ from the Herschel column density map (Könyves et al. 2010) of the region shown in Fig. 3, similar to what was found (320M⊙) from the 13CO(2–1) map by Zhu et al. (2006). The mean column density over this region is NH2 ~ 1.8 × 1022cm-2. Considering that this cluster presently contains a total of 20 protostellar sources detected with MAMBO and 93 Class II sources showing NIR K-band excess and detected with Chandra (Kuhn et al. 2010), we estimate that the current SFE of the protocluster is SFE = 56/(310 + 56) ~ 15% assuming a mean stellar mass of 0.5 M⊙4. We stress here that the mass derived from the Herschel submillimeter continuum maps may be significantly lower than the original clump mass reservoir that formed the PMS objects presently observed, as a fairly large fraction of the gas may have already been expelled from the region. This is suggested for example by the lack of any submillimeter/millimeter continuum emission toward the center of the cluster where most of the PMS stars are detected.
Assuming a distance of 260 pc for W40, the projected area of the W40 region harboring YSOs (area shown in Fig. 3) is 0.9 pc2. Assuming a period for star formation of 2 Myr (up to the Class III population), we can derive a lower-limit value for the SFR of 34M⊙ Myr-1 pc-2, which is similar to what is found for the Serpens South protocluster (see Sect. 4.1.1) or more massive star forming regions (see for example Nguyen Luong et al. 2011). However, Kuhn et al. (2010) argue that the high disk fraction they find thanks to K-band excess measurements would support an age for the W40 cluster of only ~1 Myr, which would lead to a upper-limit SFR of ~57M⊙ Myr-1 pc-2. For comparison, the SFR in the Serpens Main cluster, i.e. 56 M⊙ Myr-1 pc-2 (Harvey et al. 2007), is of the same order of magnitude as the values we find here.
We note that the classification we adopted here to determine the evolutionary status of the MAMBO sources, based on the Menv – Lbol diagram shown in Fig. 4, is independent of the distance used for W40 as the envelope mass and bolometric luminosity have the same dependency with distance. At first sight, we may expect the SFR and SFE values estimated for the W40 complex to scale as 1/d2, since both the gas mass and the cloud area scale as d2, while the number of YSOs detected stays the same. However, assuming a standard IMF, the total stellar mass is roughly proportional to the total bolometric luminosity of the YSO population, which also scales as d2. Therefore, we argue that both the SFR and the SFE we estimate for W40 do not depend much on the adopted distance.
4.2. Comparing evolutionary models and IMF predictions in our Serpens South sample
Considering that we detected 57 protostellar objects with final masses M ⋆ > 0.1M⊙ in the Serpens South cluster-forming clump, and following the IMF (Kroupa 2001; Chabrier 2003), one would expect the most massive object from our sample to have (taking into account the ~ 30% uncertainty on the number of protostellar objects in Serpens South). Similarly, applying the IMF to our sample of YSOs predicts a median final mass , for our sample.
The masses of most of the Class I sources are not constrained by our MAMBO observations (43 out of 48 are Spitzer sources not detected with MAMBO). Moreover, variable accretion rates during the embedded phases (only scenario likely to solve the luminosity problem, see Sect. 3.3, seen in our data by an over-abundance of low-luminosity Class I protostars in Fig. 5), do not allow us to use luminosities to derive accurate masses for Class I protostars not being detected at sub-millimeter or millimeter wavelengths. Therefore, we use the properties from our sample of Class 0 protostars to compute the median stellar mass and maximum stellar mass, expected from the model evolutionary tracks of Fig. 5. We then compare these values with values expected from the standard IMF. We stress however that the low number of Class 0 protostars we use here is not allowing to draw robust conclusions.
We find a median envelope mass . If we follow the evolutionary model tracks shown in Fig. 5, we predict a final median mass of . Moreover, the most massive Class 0 object from our sample (SerpS-MM18) has , and evolutionary models therefore suggest a maximum final stellar mass .
Thus, the use of the model evolutionary tracks drawn in Fig. 5 predicts both a median and maximum stellar mass in agreement with a random sampling from the standard IMF (Kroupa 2001; Chabrier 2003).
4.3. Two modes of efficient clustered star formation?
The high values of SFRs we derive for both Serpens South and W40 are consistent with the fact that we consider dense and compact regions actively forming large numbers of protostars from large mass reservoirs, i.e. protoclusters (cf. Hurt & Barsony 1996; Myers 2001; André et al. 2000; Lada & Lada 2003). These two adjacent star-forming regions are likely to share the same distance and have similar projected areas with ongoing star formation as well as similar gas reservoirs (gas masses of 610 M⊙ and 310 M⊙, respectively). However, the formation of each of these two young protostellar populations is likely to have been initiated in a different way:
-
(1)
Large-scale collapse of a gravitationally unsta-ble filament for the Serpens South cluster-formingclump, similar to what was observed in the intermediate-mass protocluster NGC 2264-C(Peretto et al. 2006).Indeed, it has been proposed based on the firstresults from Herschel (Andréet al. 2010) that a dense fil-ament with a highly supercritical mass per unit length(cf. Inutsuka & Miyama 1997)such as Serpens South may become globally gravitationallyunstable and undergo overall collapse, eventually leadingto the formation of a rich protocluster. Our MAMBO 1.2 mm map confirms that the high-column-densityfilamentary clump surrounding the central mid-infrared embed-ded cluster detected by Spitzer (Gutermuthet al. 2008) in Serpens-South is actively forming young Class 0protostars. The high fraction of Class 0 andClass I protostars compared to the total number(30–40) of Class II YSOs, found in the centralcore of the protocluster but also in its filamentary dusty arms,is suggestive of a recent burst of star formation in this cluster-forming clump. This view is also supported by the recent resultsof Nakamura et al. (2011b), whodetected numerous protostellar outflows toward this clump.
-
(2)
Pressure-triggered star cluster formation (see Elmegreen 1998; Deharveng et al. 2005) around the HII region W40. The winds and energetic radiation from the massive young stars of W40 have created an expanding HII region and compressed the gas clump at the origin of the cluster, forming a ring-like structure of dense gas now harboring a second generation of young protostars. Indeed, all of the youngest (Class 0) YSOs are found in the high column-density regions surrounding the infrared cluster containing most of the more evolved Class II and Class III objects, while not many of these more evolved objects are found in the outer dense regions harboring the Class 0 protostars. If, as suggested by the high disk fraction reported by Kuhn et al. (2010), the young infrared cluster is indeed ~1 Myr old, then W40 is currently undergoing a small burst of star formation, and may end up forming a rich bound cluster.
We have shown that these two protoclusters, although resulting from different collapse initial conditions, are presently harboring similar numbers of Class 0 protostars. This suggests that after the onset of gravitational collapse, the detailed manner in which the collapse has been initiated does not affect much the ability of the clump to form stars. Moreover, both W40 and Serpens South are characterized by large SFRs, an order of magnitude higher than the typical SFRs observed for embedded infrared clusters (Lada & Lada 2003). This suggests that (sub)millimeter protoclusters are star forming clumps caught in the very act of efficiently forming stars, while embedded infrared clusters may already be in the decay phase following the initial burst of star formation. If they last long enough, such high SFRs are likely to produce a large number of stars before gas dispersal processes (outflows, radiation, etc.) from the growing cluster membership terminate the process (see the conclusions of Nakamura et al. 2011b regarding Serpens South), i.e. reach SFEs > 30% which would allow them to remain bound (Lada & Lada 2003; Baumgardt & Kroupa 2007; Goodwin 2009). W40 and Serpens South are therefore different but good candidates to form rich bound clusters.
Grenoble Image and Line Data Analysis System, software provided and actively developed by IRAM (http://www.iram.fr/IRAMFR/GILDAS).
This new algorithm analyzes fine spatial decompositions of original images across a wide range of scales and across all wavelengths, cleans those decomposed images of noise and background fluctuations separately at each wavelength, and constructs a set of wavelength-independent detection images that preserve information in both spatial and wavelength dimensions as well as possible. Sources are detected in the set of decomposed combined detection images by following the evolution of their segmentation masks across all spatial scales. Measurements of the sources’ properties are performed in the original images at each wavelength after the background has been subtracted by interpolation under the sources “footprints” and overlapping sources have been deblended, see Men’shchikov et al. (in prep.).
The original criterion introduced by André et al. (1993) was , but only YSOs with were considered bona-fide Class 0 objects by André et al. (2000), while sources with between ~ 0.5% and ~ 1% were classified as border-line objects (see Table 1 of André et al. 2000).
Note that if we take into account the additional Class III objects identified in X-rays by Kuhn et al. (2010) (~109 Chandra YSOs without K-band excess), then the SFE goes up to ~26% and the upper-limit of the SFR to ~120 M⊙ Myr-1 pc-2 (assuming a lower-limit age of ~1 Myr, see Kuhn et al. (2010)).
Acknowledgments
We are grateful to the specialist astronomy group dedicated to star formation within the SPIRE Consortium (“SAG 3”) for letting us use a preliminary version of their Herschel photometric results in advance of publication. A first-generation paper including the full Herschel sources catalog in the Aquila rift complex will be published shortly (Könyves et al., in prep.). We thank R. Gutermuth for providing us with its list of Spitzer IRAC Class I sources in the Serpens South field prior to publication. We thank M. Kuhn, E. Feigelson and K. Getman for useful comments on the Chandra PMS population in W40.
References
- Altenhoff, W. J., Downes, D., Goad, L., Maxwell, A., & Rinehart, R. 1970, A&AS, 1, 319 [NASA ADS] [Google Scholar]
- Alves, J., Lombardi, M., & Lada, C. J. 2007, A&A, 462, L17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- André, P., & Montmerle, T. 1994, ApJ, 420, 837 [NASA ADS] [CrossRef] [Google Scholar]
- André, P., Ward-Thompson, D., & Barsony, M. 1993, ApJ, 406, 122 [NASA ADS] [CrossRef] [Google Scholar]
- André, P., Ward-Thompson, D., & Barsony, M. 2000, Protostars and Planets IV, 59 [Google Scholar]
- André, P., Belloche, A., Motte, F., & Peretto, N. 2007, A&A, 472, 519 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- André, P., Minier, V., Gallais, P., et al. 2008, A&A, 490, L27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- André, P., Men’shchikov, A.,Bontemps, S., et al. 2010, A&A, 518, L102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Baumgardt, H., & Kroupa, P. 2007, MNRAS, 380, 1589 [NASA ADS] [CrossRef] [Google Scholar]
- Belloche, A., Schuller, F., Parise, B., et al. 2011, A&A, 527, A145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bontemps, S., André, P., Terebey, S., & Cabrit, S. 1996, A&A, 311, 858 [Google Scholar]
- Bontemps, S., André, P., Könyves, V., et al. 2010, A&A, 518, L85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bressert, E., Bastian, N., Gutermuth, R., et al. 2010, MNRAS, 409, L54 [NASA ADS] [Google Scholar]
- Carpenter, J. M. 2000, AJ, 120, 3139 [NASA ADS] [CrossRef] [Google Scholar]
- Chabrier, G. 2003, PASP, 115, 763 [NASA ADS] [CrossRef] [Google Scholar]
- Chen, H., Myers, P., Ladd, E., & Wood, D. 1995, ApJ, 445, 377 [NASA ADS] [CrossRef] [Google Scholar]
- Deharveng, L., Zavagno, A., & Caplan, J. 2005, A&A, 433, 565 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dzib, S., Loinard, L., Mioduszewski, A. J., et al. 2010, ApJ, 718, 610 [NASA ADS] [CrossRef] [Google Scholar]
- Elmegreen, B. G. 1998, in Origins, ed. C. E. Woodward, J. M. Shull, & H. A. Thronson Jr., ASP Conf. Ser., 148, 150 [Google Scholar]
- Enoch, M. L., Glenn, J., Evans, II, N. J., et al. 2007, ApJ, 666, 982 [NASA ADS] [CrossRef] [Google Scholar]
- Enoch, M. L., Evans, N. J., Sargent, A. I., & Glenn, J. 2009, ApJ, 692, 973 [NASA ADS] [CrossRef] [Google Scholar]
- Evans, N. J., Dunham, M. M., Jørgensen, J. K., et al. 2009, ApJS, 181, 321 [NASA ADS] [CrossRef] [Google Scholar]
- Ginsburg, A. G., Bally, J., Yan, C., & Williams, J. P. 2009, ApJ, 707, 310 [NASA ADS] [CrossRef] [Google Scholar]
- Goodwin, S. P. 2009, Ap&SS, 324, 259 [NASA ADS] [CrossRef] [Google Scholar]
- Greene, T., Wilking, B., André, P., Young, E., & Lada, C. 1994, ApJ, 434, 614 [NASA ADS] [CrossRef] [Google Scholar]
- Gutermuth, R. A., Bourke, T. L., Allen, L. E., et al. 2008, ApJ, 673, L151 [NASA ADS] [CrossRef] [Google Scholar]
- Harvey, P., Merín, B., Huard, T. L., et al. 2007, ApJ, 663, 1149 [NASA ADS] [CrossRef] [Google Scholar]
- Henning, T., Michel, B., & Stognienko, R. 1995, Planet. Space Sci., 43, 1333 [NASA ADS] [CrossRef] [Google Scholar]
- Henriksen, R., Andre, P., & Bontemps, S. 1997, A&A, 323, 549 [NASA ADS] [Google Scholar]
- Hillenbrand, L. A., & Carpenter, J. M. 2000, ApJ, 540, 236 [NASA ADS] [CrossRef] [Google Scholar]
- Hurt, R., & Barsony, M. 1996, ApJ, 460, L45 [NASA ADS] [CrossRef] [Google Scholar]
- Inutsuka, S., & Miyama, S. M. 1997, ApJ, 480, 681 [NASA ADS] [CrossRef] [Google Scholar]
- Kenyon, S. J., Hartmann, L. W., Strom, K. M., & Strom, S. E. 1990, AJ, 99, 869 [NASA ADS] [CrossRef] [Google Scholar]
- Kirk, H., Johnstone, D., & Tafalla, M. 2007, ApJ, 668, 1042 [NASA ADS] [CrossRef] [Google Scholar]
- Knee, L. B. G., & Sandell, G. 2000, A&A, 361, 671 [NASA ADS] [Google Scholar]
- Knude, J. 2011 [arXiv:1103.0455] [Google Scholar]
- Könyves, V., André, P., Men’shchikov, A., et al. 2010, A&A, 518, L106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kroupa, P. 2001, MNRAS, 322, 231 [NASA ADS] [CrossRef] [Google Scholar]
- Kuhn, M. A., Getman, K. V., Feigelson, E. D., et al. 2010, ApJ, 725, 2485 [NASA ADS] [CrossRef] [Google Scholar]
- Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57 [NASA ADS] [CrossRef] [Google Scholar]
- Lada, C. J., Lombardi, M., & Alves, J. F. 2010, ApJ, 724, 687 [NASA ADS] [CrossRef] [Google Scholar]
- Loinard, L., Mioduszewski, A. J., Torres, R. M., et al. 2011, XIII Latin American Regional IAU Meeting, ed. W. J. Henney & S. Torres-Peimbert, 40, 205 [Google Scholar]
- Matzner, C., & McKee, C. 2000, ApJ, 545, 364 [NASA ADS] [CrossRef] [Google Scholar]
- Maury, A., André, P., & Li, Z.-Y. 2009, A&A, 499, 175 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Men’shchikov, A., André, P., Didelon, P., et al. 2010, A&A, 518, L103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Meyer, M. R., Adams, F. C., Hillenbrand, L. A., Carpenter, J. M., & Larson, R. B. 2000, Protostars and Planets IV, 121 [Google Scholar]
- Motte, F., & André, P. 2001, A&A, 365, 440 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Motte, F., André, P., & Neri, R. 1998, A&A, 336, 150 [Google Scholar]
- Myers, P. C. 2001, in From Darkness to Light: Origin and Evolution of Young Stellar Clusters, ed. T. Montmerle, & P. André, ASP Conf. Ser., 243, 131 [Google Scholar]
- Myers, P. C., Dame, T. M., Thaddeus, P., et al. 1986, ApJ, 301, 398 [NASA ADS] [CrossRef] [Google Scholar]
- Myers, P., Adams, F., Chen, H., & Schaff, E. 1998, ApJ, 492, 703 [NASA ADS] [CrossRef] [Google Scholar]
- Nakamura, F., Kamada, Y., Kamazaki, T., et al. 2011a, ApJ, 726, 46 [NASA ADS] [CrossRef] [Google Scholar]
- Nakamura, F., Sugitani, K., Shimajiri, Y., et al. 2011b, ApJ, 737, 56 [NASA ADS] [CrossRef] [Google Scholar]
- Nguyen Luong, Q., Motte, F., Hennemann, M., et al. 2011, A&A, 535, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nutter, D., & Ward-Thompson, D. 2007, MNRAS, 374, 1413 [NASA ADS] [CrossRef] [Google Scholar]
- Olmi, L., & Testi, L. 2002, A&A, 392, 1053 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943 [NASA ADS] [Google Scholar]
- Palla, F., & Stahler, S. W. 1999, ApJ, 525, 772 [NASA ADS] [CrossRef] [Google Scholar]
- Peretto, N., André, P., & Belloche, A. 2006, A&A, 445, 979 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Peretto, N., Hennebelle, P., & André, P. 2007, A&A, 464, 983 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Robitaille, T. P., Whitney, B. A., Indebetouw, R., & Wood, K. 2007, ApJS, 169, 328 [NASA ADS] [CrossRef] [Google Scholar]
- Saraceno, P., Andre, P., Ceccarelli, C., Griffin, M., & Molinari, S. 1996, A&A, 309, 827 [NASA ADS] [Google Scholar]
- Schneider, N., Bontemps, S., Simon, R., et al. 2011, A&A, 529, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Smith, J., Bentley, A., Castelaz, M., et al. 1985, ApJ, 291, 571 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, R., Fellhauer, M., Goodwin, S., & Assmann, P. 2011, MNRAS, 414, 3036 [NASA ADS] [CrossRef] [Google Scholar]
- Stahler, S. 1988, ApJ, 332, 804 [NASA ADS] [CrossRef] [Google Scholar]
- Straižys, V., Černis, K., & Bartašiūtė, S. 1996, Baltic Astron., 5, 125 [NASA ADS] [Google Scholar]
- Straižys, V., Černis, K., & Bartašiūtė, S. 2003, A&A, 405, 585 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Testi, L., & Sargent, A. 1998, ApJ, 508, L91 [NASA ADS] [CrossRef] [Google Scholar]
- Testi, L., Palla, F., & Natta, A. 1999, A&A, 342, 515 [NASA ADS] [Google Scholar]
- Vallee, J. P. 1987, A&A, 178, 237 [NASA ADS] [Google Scholar]
- Vallee, J. P., Guilloteau, S., & MacLeod, J. M. 1992, A&A, 266, 520 [NASA ADS] [Google Scholar]
- Visser, A., Richer, J., & Chandler, C. 2002, AJ, 124, 2756 [NASA ADS] [CrossRef] [Google Scholar]
- Wilking, B., & Lada, C. 1983, ApJ, 274, 698 [NASA ADS] [CrossRef] [Google Scholar]
- Young, C. H., & Evans, II, N. J. 2005, ApJ, 627, 293 [NASA ADS] [CrossRef] [Google Scholar]
- Zhu, L., Wu, Y., & Wei, Y. 2006, Chinese J. Astron. Astrophys., 6, 61 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
Fig. 1 1.2 mm dust continuum MAMBO mosaic of the Aquila rift complex, smoothed to an effective angular resolution of 13. Due to the high dynamic range of fluxes in the map, and to enhance visibility of the structure of emission in the whole area mapped, the color scale chosen here has a maximum of 400 mJy/beam. The highest flux density of order ~1400 mJy/beam is found in the center of the Serpens South protocluster (see the blow-up in Fig. 2). In this smoothed map, the mean rms noise level is ~13 mJy/13-beam. |
|
In the text |
Fig. 2 Blow-up of the Region 1 from the MAMBO 1.2 mm dust continuum map shown in Fig. 1, around the Serpens South filamentary cloud. The mean rms noise level is ~12 mJy/13-beam. Black contours show the levels of 5σ, 8σ, 15σ, 30σ, 70σ, 90σ, and 120σ (1.44 Jy/beam). |
|
In the text |
Fig. 3 Blow-up of the Region 2 from the 1.2 mm continuum map shown in Fig. 1, around the HII region W40. The mean rms noise level is ~14 mJy/13-beam. Black contours show the levels of 5σ, 8σ, 15σ, 20σ, 25σ, and 35σ. One of the candidate exciting stars, Sh64-OS2a, is detected in our MAMBO mapping and indicated here. |
|
In the text |
Fig. 4 Envelope mass versus bolometric luminosity diagram comparing the locations of the protostars detected with MAMBO in the Aquila rift region (open triangles for the Serpens South sub-field, open squares for the other sub-fields) with the positions of known low-mass Class I (filled circles) and Class 0 objects (open circles) in Ophiuchus, Perseus, and Orion (André et al. 2000; Motte & André 2001). Numbers reported inside triangles are the reference numbers of the sources reported in the last column of Table 1. Model protostellar evolutionary tracks, computed for various final stellar masses (indicated above the Lbol axis) assuming the accretion/luminosity history described in the text, are superimposed (4 blue tracks representing protostellar evolution for objects with final masses 0.1, 0.3, 1 and 3M⊙). Evolution proceeds from the upper left to the lower right as indicated by two arrows on each track, plotted when 50% and 90% of the final stellar mass has been accreted, respectively. The straight lines show two Menv–Lbol relations marking the conceptual border zone between the Class 0 (Menv > M ⋆ /ϵ) and the Class I (Menv < M ⋆ /ϵ) stage, where ϵ is the local star formation efficiency taken here to be 50%. The dashed line is such that Menv ∝ Lbol (cf. André et al. 1993; André & Montmerle 1994), while the dotted relation follows as suggested by the accretion scenario adopted in the tracks (see also Bontemps et al. 1996). The typical uncertainty in the plotted values of Lbol is ~30–40%, while the uncertainty in the Menv values is a factor of ~2 (see Sect. 3.2). The black arrow in the upper right indicates how the positions of the W40 protostars would move in the diagram if a distance of 600 pc was adopted for W40, instead of 260 pc. |
|
In the text |
Fig. 5 Same as Fig. 4, but showing only the sources from the Serpens South cluster forming clump. Triangles show the locations of the MAMBO protostars (green triangles are Class 0 and red triangles are Class I protostars, while black triangles represent the border-line objects). Red circles show the locations of the Class I protostars detected by Spitzer and without MAMBO detection. The Class I sources without any MAMBO counterpart are indicated as upper-limit values for Menv, around the detection limit of 2.5 × 10-2M⊙. The green double arrow along the Menv axis shows the range of masses for the prestellar cores detected by Herschel, based on a preliminary analysis of the Gould Belt survey data, which are 85% complete down to 0.3M⊙ in this region (Könyves et al. 2010). The respective locations of the prestellar cores, Class 0 and Class I protostars in the Serpens South protocluster are roughly consistent with what is expected from a single population of objects evolving along the model evolutionary tracks shown in blue. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.