arxiv: v1 [astro-ph.ga] 18 Mar 2010

Transcription

1 Astronomy & Astrophysics manuscript no. CV AA V c ESO August, The fate of the pre-main sequence-rich clusters Collinder 97 and vdb 9: dissolution? C. Bonatto and E. Bica Universidade Federal do Rio Grande do Sul, Departamento de Astronomia CP 55, RS, Porto Alegre 95-97, Brazil charles@if.ufrgs.br, bica@if.ufrgs.br arxiv:.v [astro-ph.ga] Mar Received ; accepted ABSTRACT Context. Since most of the star clusters do not survive the embedded phase, their early dissolution appears to be a major source of field stars. However, catching a cluster in the act of dissolving is somewhat elusive. Aims. We investigate the nature and possible evolution of the young Galactic star clusters Collinder 97 (Cr 97) and vdb 9. Methods. Photometric and structural properties are derived with near-infrared photometry and field-star decontamination. Kinematical properties are inferred from proper motions of the main sequence (MS) and pre-ms (PMS) member stars. Results. The colour-magnitude diagrams (CMDs) are basically characterised by a poorly-populated MS and a dominant fraction (> 75%) of PMS stars, and the combined MS and PMS CMD morphology in both clusters consistently constrains the age to within 5 ± 4Myr, with a Myr spread in the star formation process. The MS+PMS stellar masses are + 59 M (Cr97) and M (vdb9). Cr97 and vdb9 appear to be abnormally large, when compared to clusters within the same age range. They have irregular stellar radial density distributions (RDPs) with a marked excess in the innermost region, a feature that, at less than Myr, is more likely related to the star formation and/or molecular cloud fragmentation than to age-dependent dynamical effects. The velocity dispersion of both clusters, derived from proper motions, is in the range 5 kms. Conclusions. Both clusters appear to be in a super-virial state, with velocity dispersions higher than those expected of nearly-virialised clusters of similar mass and size. A possible interpretation is that Cr 97 and vdb 9 deviate critically from dynamical equilibrium, and may dissolve into the field. We also conclude that early cluster dissolution leaves detectable imprints on RDPs of clusters as massive as several M. Cr97 and vdb9 may be the link between embedded clusters and young stellar associations. Key words. (Galaxy:) open clusters and associations: general; (Galaxy:) open clusters and associations: individual: Collinder97 and vdb9. Introduction It is now well-established that a significant fraction of the embedded star clusters dissolve into the field on a timescale of a few 7 Myrs. Basically, dissolution occurs mainly because the gravitational potential can be rapidly reduced by internal processes, such as the impulsive gas removal by supernovae and massive star winds associated with this early period. As a consequence, an important fraction of the stars, especially of low mass, end up moving faster than thescaled-downescapevelocity,andmaybelosttothefield (e.g. Goodwin & Bastian ). This process can dissolve the very young star clusters on time-scale of 4Myr. From the stellar content perspective, the early cluster dissolution depends essentially on the effective starformation efficiency, the total mass converted into stars, and the mass of the more massive stars (e.g. Tutukov 97; Goodwin & Bastian ). However, there is also evidence indicating that the determining factor for cluster survival during gas expulsion is the virial state of the stars just be- Send offprint requests to: C. Bonatto Correspondence to: charles@if.ufrgs.br fore the onset of this phase, so that clusters formed with a dynamically cold stellar component are more likely to survive (Goodwin 9). In any case, the early dissolution of embedded clusters may lead to the formation of OB stellar groups (e.g. Gouliermis et al. ), the subsequent dispersion of which may be an important source of field stars (e.g. Massey, ohnson & Gioia-Eastwood 995). Indeed, Lada & Lada () suggest that only about 5% of the Galactic embedded clusters dynamically evolve into gravitationally bound open clusters (OCs). With such a high dissolution rate, the embedded clusters could be the major contributors of field stars in galaxies for generations. However, recent studies suggest that the early dissolution rate in the Magellanic Clouds is significantly lower (< % - de Grijs & Goodwin ; de Grijs & Goodwin 9) than in the Milky Way (Lada & Lada ), or in other galaxies such as the Antennae (e.g. 7) or M5 (e.g. Bastian et al. 5). Observationally, low-mass star clusters younger than Myr, in general, have Colour-Magnitude Diagrams (CMDs) with an under-populated and developing main sequence (MS) and a more conspicuous population of

2 C. Bonatto and E. Bica: Collinder 97 and vdb 9: Dissolving star clusters? pre-ms (PMS) stars. Typical examples are NGC, NGC4755, NGC9, NGC44, Bochum, Pismis5, NGC9, vdb, and BDBS9 (Bonatto & Bica 9c and references therein). In terms of CMD morphology, bound and un-bound young clusters are expected to present similar evolutionary sequences. On the other hand, the important early changes in the potential that affect the large-scale internal structure of clusters should be reflected on the stellar radial density profile (RDP). Bochum (Bica, Bonatto & Dutra ) and NGC44 (Bonatto & Bica 9b), for instance, appear to be representatives of this scenario (i.e. structures evolving towards dissolution in a few 7 yr), in which an irregularrdp cannot be represented by a cluster-like (i.e. an approximately isothermal sphere) profile. In this context, irregular RDPs in young clusters - when coupled to an abnormally high velocity dispersion - may reflect significant profile erosion or dispersion of stars, and point to important deviations from dynamical equilibrium. The rather complex interplay among environment conditions, effective star-formation efficiency (as defined in Goodwin & Bastian ), and total mass converted into stars, is probably what explains the difference between(dissolving) objects like Bochum and NGC44, and bound young OCs (in which the MS+PMS stars distribute according to a cluster RDP as in NGC and NGC4755). Consistent with the above mass-dependent scenario, the MS+PMS mass of Bochum and NGC44 is a factor - lower than in NGC and NGC4755. In this paper we investigate the nature of the poorlystudied, young (age 5Myr; Sect. 4), large (with radii within pc; Sect. 5), and PMSrich (with stellar masses within 75M ; Sect. ) clusters Cr97 and vdb9. Our main goal is to determine whether such young stellar systems can be characterised as typical OCs or if they are heading towards dissolution. In addition, we will derive their fundamental and structural parameters, most of these for the first time. This paper is organised as follows. In Sect. we recall literature data on both objects. In Sect. we discuss the MASS photometry and build the field-star decontaminated CMDs. In Sect. 4 we derive fundamental cluster parameters. In Sect. 5 we derive structural parameters. In Sect. we estimate cluster mass. In Sect. 7 we build the intrinsic proper motion distribution. In Sect. we compare structural parameters and dynamical state with those of a sample of template OCs. Concluding remarks are given in Sect. 9.. Previous data on Cr97 and vdb9 Collinder 97 (Cr 97) was discovered by Collinder (9) in Vela. It wasalsolisted in the OCcatalogueofAlter et al. (97) as OCl-74, and the ESO/Uppsala Blue Plate Southern (Lauberts 9) as ESO SC. Vogt & Moffat (97) centred the cluster on the B4II star HIP49 at α() = h 44 m 4. s, δ() = 4. This star, which is reported as a double system in SIMBAD, is the brightest in the areawith V = 7.. HIP49 may be the cluster centre in the infrared, since strong absorption is suggested by the B plate image seen to the south/west of HIP49. Vogt & Moffat (97) found the reddening E(B V) =.5 and the distance from the Sun d =.5kpc. Fig.. Top: DSS-II R images centred on Cr97 (left) and vdb 9 (right). Gas emission, dust reflection and/or absorption are present in the fields in varying proportions. Bottom: MASS K s images covering 5 5 (Cr97) and (vdb9). Orientation: North to the top and East to the left. Besides dust absorption and emission in the area, Cr 97 is embedded also in the HII region Gum5 (Gum 955), also known as RCW (Rodgers, Campbell & Whiteoak 9). Assuming thathip 49 is a cluster member, it might responsible for the emission. Indeed, Pettersson & Reipurth (994) have identified lowluminosity emission-line stars in Cr 97 and in the R-association Vela R, which may be part of their lowmass population. They derive d =.kpc and estimate a few yr of age. Battinelli & Capuzzo-Dolcetta (99) and Battinelli, Brandimarti & Capuzzo-Dolcetta (994) estimated.myr of age, d = kpc, a total mass of M, and the integrated absolute magnitude M V = 5.. Found by van den Bergh (9), the embedded cluster vdb9 is related to the reflection nebula vdb-rn9, which is located in Canis Major. The cluster is also listed as FSR (Froebrich, Scholz & Raftery 7) who, based on the H photometry, derived the core and tidal radii R c,h =. and R t,h = 4, and the total number of membersn = 97stars.Magakian() alsoreportsthereflection nebula and related cluster. It is immersed in the large angular size association CMaOB (Clariá 974a; Clariá 974b), and in the reflection nebula association CMaR (van den Bergh 9). The equally extended emission nebula IC77 or Gum (Gum 955) seems to permeate the complex. The relatively bright star at.4 to the northeast of the cluster vdb9 is the Herbig Be star ZCMa, which has a bipolar outflow (Poetzel, Mundt & Ray 99). Undoubtedly, vdb 9 is part of a star-forming complex, and deserves further analyses.

3 Table. Fundamental parameters C. Bonatto and E. Bica: Collinder 97 and vdb 9: Dissolving star clusters? WEBDA This work Cluster α() δ() Age E(B V) d D α() δ() l b Age E(B V) d R SC (hms) ( ) (Myr) (mag) (kpc) ( ) (hms) ( ) ( ) ( ) (Myr) (mag) (kpc) (kpc) () () () (4) (5) () (7) () (9) () () () () () (5) Cr 97 :44:5 4:: :44:4. 4:: ± 4.4 ±..5 ±.. ±.4 vdb 9 7::54 :9: ::5.4 :4: ± 4. ±.. ±..4 ±.9 Table Notes. Col. 7: Optical diameter; Col. : distance from the Sun; col. 5: distance from the Solar circle. WEBDA data for Cr97 is based on photometry from Vogt & Moffat (97). With an approach that differs in various aspects from the present one, Soares & Bica () studied the brighter sequences of vdb9, deriving an age of 5-7Myr, a mean visible absorption = 4.4, and d.5kpc. A few bright stars mixed with nebular gas and/or dust emission are seen in the R images (Fig., taken from LEDAS ). Close-ups of the clusters are shown in the smaller field MASS K s images. Table provides parameters found in the literature and as derived here. The central coordinates were re-computed to match the absolute maxima present in the stellar surface densities (Sect..). We note that the re-computed central coordinates of Cr 97 are almost coincident with HIP 49.. Construction of decontaminated CMDs Since both objects still retain part of the primordial gas and dust (Fig. ), the near-ir provides the adequate depth to study them, especially the faint stellar sequences. For this purpose, we work with the MASS, H, and K s photometry, which provides the spatial and photometric uniformity required for wide extractions that, in turn, result in high star-count statistics. Our group has been developing analytical tools for the MASS photometry that allow us to statistically disentangle cluster evolutionary sequences from field stars in CMDs. Decontaminated CMDs have been used to investigate the nature of star cluster candidates and derive their astrophysical parameters. In summary, we apply (i) field-star decontamination to uncover the intrinsic CMD morphology, essential for a proper derivation of reddening, age, and distance from the Sun, and (ii) colour-magnitude filters, which are required for intrinsic stellar RDPs, as well as luminosity and mass functions (MFs). In particular, the use of field-star decontamination in the construction of CMDs hasprovedto constrainageand distance morethan the raw (observed) photometry, especially for low-latitude and/or bulge-projected OCs (e.g. Bica, Bonatto & Camargo, and references therein). Photometry for both clusters was extracted from VizieR in a wide circular field of radius R ext =, which is adequate to allow determination of the background level (Sect. 5) and to statistically characterise the colour and magnitude distribution of the field stars Leicester Database and Archive Service (LEDAS) DSS/DSS-II service on ALBION; The Two Micron All Sky Survey, All Sky data release (Skrutskie et al. 997) (Sect..). As a photometric quality constraint, only stars with, H, and K s errors lower than.mag were used. Reddening corrections are based on the absorption relations A / =.7, A H / =.7, A KS / =., and A =.7 E( H) given by Dutra, Santiago & Bica (), with R V =., considering the extinction curve of Cardelli, Clayton & Mathis (99). The reddening values are derived from the MASS CMDs (Sect. 4)... Field decontamination Field-star decontamination is usually required for the proper identification and characterisation of star clusters, especially those near the Galactic equator and/or with important fractions of faint stars. Cr97 and vdb9 are located in the rd Galactic quadrant (Table ), which makes field-star contamination a minor issue (e.g. Bica, Bonatto & Camargo ). However, their CMDs are dominated by PMS stars (Sect. 4) and, thus, decontamination is important to avoid confusion with the red dwarfs of the Galactic field. A summary of different decontamination approaches is provided in Bonatto & Bica (9b). In this paper we apply the decontamination algorithm developed in Bonatto & Bica(7), together with an improvement that is described below. Photometric uncertainties are explicitly taken into account, in the sense that stars have a non-negligible probability (assumed to be Gaussian) of having a magnitude and colour significantly different from the average values. The algorithm (i) divides the full range of magnitude and colours of a CMD into a D grid of cells with axes along the magnitudeandthe( H)and( K s )colours,(ii) computes the total (i.e. member+field) number-density of stars within a given cell (η tot ), (iii) estimates the numberdensity of field stars (η fs ) based on the number of comparison field stars with similar magnitude and colours as those in the cell, (iv) computes the expected number-density of member stars, η mem = η tot η fs, (v) converts the numberdensity η fs into the estimated number of field stars and subtracts it from each cell, and (vi) after subtraction of the field stars, the remaining Nclean cell stars in each cell are identified for further use (see below). The wide annulus between 5 < R < (outside the cluster radius - Table ) is used as the comparison field for Cr97; for vdb9 it is located between < R <. The initial cell dimensions are =. and ( H) = ( K s ) =., but cell sizes half and twice the initial values are also used. As a new feature, we also consider shifts in the grid positioning by ±/ofthe respective cell size in the coloursand magnitude axes. When all the different grid/cell size setups are

4 4 C. Bonatto and E. Bica: Collinder 97 and vdb 9: Dissolving star clusters? (-Ks) (-Ks) (-Ks) (-Ks) Fig..Toppanels:number-density(η mem,instarsarcmin )oftheprobablememberstarspresentinthedecontaminated CMDs extracted within R = of Cr97 (left) and R = 5 of vdb9 (right). Bottom: survival frequency of the decontaminated stars. applied, the decontamination takes into account 79 independent combinations. Each setup produces a total number of member stars, N mem = cell Ncell clean,fromwhichwecomputetheexpected total number of member stars N mem by averaging out N mem over all combinations. Stars (identified in step (vi) above) are ranked according to the number of times they survive all runs (survival frequency), and only the N mem highest ranked stars are considered cluster members and transposed to the respective decontaminated CMD (e.g. Figs. 5 and ). The subtraction efficiency, i.e. the difference between the expected number of field stars(which may be fractional) and the number of stars effectively subtracted (which is integer) from each cell, summed over all cells, is 99.4% for vdb9 and 9.7% for Cr97. As a caveat we note that the present decontamination approach implicitly assumes that the field colour-magnitude distribution somewhat matches that of the cluster. Figure shows some gas and dust around both clusters, which might lead to appreciable reddening and differential reddening. While the effect in the foreground contamination may be small, it should be more important in the background. However, the rd Galactic quadrant location of Cr97 and vdb9 is expected to minimise the background contamination. The number-density of the probable member stars (η mem ) are shown in the ( K s ) CMDs in Figure (top panels). In both cases, most of the stars are relatively faint and red, covering a wide range in colour, which is expected of PMS stars somewhat affected by differential reddening. Also shown in Figure (bottom panels) is the survival frequency of the decontaminated stars. Again, the highest membership probabilities occur among the PMS stars... Decontaminated surface density maps Our decontamination approach relies upon differences in the colour and magnitude distribution of stars located in separate spatial regions. For a star cluster, which can be characterised by a single-stellar population projected against a Galactic stellar field, the decontaminated surfacedensity is expected to present a marked excess at the assumed cluster position. Maps of the spatial distribution of the stellar surfacedensity (σ, in units of starsarcmin ), built with fieldstar decontaminated photometry to maximise the cluster/background contrast, are shown in Figs. (Cr 97) and 4 (vdb9). Also shown are the isopleths, in which cluster size and geometry can be observed. The surface density is

5 C. Bonatto and E. Bica: Collinder 97 and vdb 9: Dissolving star clusters? 5 5 Cr97 (R= ) vdb9 (R=5 ) (α cosδ) (α cosδ) Sky ( ) Sky (9.4 ) (α cosδ) Decont Decont. (α cosδ) AV= AV= Fig.. Left: Stellar surface-density σ(stars arcmin ) computed with field-decontaminated photometry to enhance the cluster/background contrast of Cr 97. Right: Average ( Ks ) colour. (α cos(δ)) and in arcmin. AV= (α cosδ) (α cosδ) (α cosδ) (α cosδ) AV=5.5.5 ( KS) ( KS).5 Fig. 5. ( Ks ) CMDs of Cr 97 (left) and vdb 9 (right) showing the observed photometry for representative regions (top) and the equal-area comparison fields (middle). The decontaminated CMDs are shown in the bottom panels. Reddening vectors for AV = 5 are shown in the bottom panels. a blueward dip, due to the few bright and rather blue stars, about the centre (right panel). The stellar distribution of vdb 9, on the other hand, presents a broad and elongated external region (Fig. 4, left panel), with a complex colour distribution (right panel), within. < ( Ks ) <.5. If most of these colour variations (with respect to the average) are due to non-uniform reddening - and not to systematic differences in the stellar content - the upper limits to the differential reddening distribution would be AV <.5 mag for Cr 97 and AV < 5 mag for vdb Fundamental parameters Fig. 4. Same as Fig. for vdb 9. computed in a rectangular mesh with cells.5.5 wide, reaching total offsets of α cos(δ) = with respect to the cluster centre (Table ). Despite the gas and dust associated with the present clusters (Fig. ), the central excesses show up markedly in the decontaminated surface-density distributions. Cr 97 has a well-defined, approximately round, and narrow stellar distribution (Fig., left panel). Its ( Ks ) colour distribution is rather uniform, within. < ( Ks ) <., with ( Ks ) CMDs built with the raw photometry of Cr 97 and vdb 9 are shown in the top panels of Fig. 5. In both cases, the sampled region contains most of the cluster stars (Fig. ). Features typical of very young OCs, such as a relatively vertical and poorly-populated MS, together with a large population of faint and red PMS stars, can be seen when these CMDs are compared to those extracted from the equal-area offset fields4 (middle panels). 4 The equal-area field extractions are used only for qualitative comparisons, since the decontamination uses a wide surrounding ring area (Sect..).

6 C. Bonatto and E. Bica: Collinder 97 and vdb 9: Dissolving star clusters? Cr97 = MS: 5Myr PMS:.Myr PMS: Myr PMS: 5Myr PMS: Myr =5.5.5 ( K S ) = vdb9 = ( K S ) Fig.. Adopted MS+PMS isochrone solutions to the decontaminated CMDs. Reddening vectors as in Fig. 5. Shaded polygons show the MS (dark-gray) and PMS (lightgray) colour-magnitude filters (Sect. 5). As expected, the decontaminated CMDs(bottom panels of Fig. 5) contain essentially the typical stellar sequences of mildly reddened young clusters, with a developing MS and a significant population of PMS stars. Also, in both cases the colour distribution at faint magnitudes ( > ) is wider than the spread predicted purely by PMS models (Fig. ), which implies internal differential reddening. To examine this issue we include in Figs. 5- reddening vectors computed with the MASS ratios (Sect. ) for visual absorptions = to 5. Taking into account the PMS isochrone fit (Fig. ), the differential reddening appears to be lower than = 5, especially for Cr97, which is consistent with the values estimated in Sect... The fundamental parameters are derived by means of the constraints provided by the field-decontaminated CMD morphologies, especially the combined MS and PMS star distribution (Fig. ). We adopt solar metallicity isochrones because both objects are young and located not far from the Solar circle (see below), a region essentially occupied by [Fe/H]. OCs (Friel 995). Padova isochrones (Girardi et al. ) computed with the MASS, H, and K s filters 5 are used to represent the MS. Isochrones 5 These isochrones are very similar to the ohnson-kron-cousins ones (e.g. Bessel & Brett 9), with differences of at most. mag in colour (Bonatto, Bica & Girardi 4). of Siess, Dufour & Forestini () with ages in the range. Myr are used to characterise the PMS sequences. Since the decontaminated CMD morphologies are very similar and typical of young ages (Fig. ), a similar isochrone solution is expected to apply to both objects. As summarised in Naylor & effries (), sophisticated approaches are available for analytical CMD fitting, especially the MS. However, given the poorly-populated MSs, the MASS photometric uncertainties for the fainter stars, the important population of PMS stars, and the differential reddening, we applied the direct comparison of isochrones with the decontaminated CMD morphology. The fits are made by eye, taking the combined MS and PMS stellar distributions as constraint, allowing as well for variations due to photometric uncertainties (which, given the restrictions imposed in Sect., are small in both cases), and differential reddening. Specifically, we start with the MS+PMS isochrones set for zero distance modulus and reddening, and apply shifts in magnitude and colour to them. We take as best-fit the apparent distance modulus and reddening that simultaneously account for the blue border of the MS and PMS stellar distributions. Both clusters present a significant fraction of stars redder than the youngest PMS isochrone. Most of this ( K s ) excess towards the red is probably due to differential reddening. Below we discuss the results, which are shown in Fig.. Cr 97: Given the poorly-populated and nearly vertical MS, acceptable fits to the decontaminated MS are obtained with any isochrone with age in the range Myr. If we take into account the differential reddening, this age spread is consistent with the PMS stars, which are basically contained within the. Myr and Myr isochrones as well. Thus, we take the 5Myr isochrone as a representative solution, and allow for a ±4Myr age-spread in the star formation. With the adopted solution, the fundamental parameters ofcr97arethe near-irreddeninge( H) =.±.5 (E(B V) =.4±. or =.±.5), the observed and absolute distance moduli (m M) =.4±.4 and (m M) O =. ±.4, respectively, and the distance from the Sun d =.5±.kpc. This distanceis compatible with that obtained for the B4II star HIP49 (B = 7.5 and V = 7.), d =.kpc. With R = 7.±.kpc (Bica et al. ) as the Sun s distance to the Galactic centre, the Galactocentric distance of Cr97 is R GC = 7.5±.kpc,whichputsit.kpcoutsidetheSolarcircle. This solution is shown in Fig.. The projected coordinate componentsare(x GC,y GC,z GC ) ( 7.4,.4,.)inkpc, which put Cr 97 at the Orion-Cygnus and Sagittarius- Carina interarm region (e.g. Momany et al. ). We also compute the integrated apparent and absolute magnitudes in the MASS bands for the member stars within R R RDP, together with the reddening corrected colours. The results are m = 5.4 ±.7, M = 4.7 ±., ( H) =.9±.9, ( K s ) =.±.9. vdb 9: Similarly to Cr 97, any isochrone with age within Myr provides an acceptable fit to the nearly vertical MS, while the PMS distribution is essentially contained within the. Myr isochrones. With Derived by means of the Globular Cluster spatial distribution. Recently, Trippe et al. () found R GC =.7±.kpc while Ghez et al. () found R GC =. ±.kpc or R GC =.4 ±.4 kpc, under different assumptions.

7 C. Bonatto and E. Bica: Collinder 97 and vdb 9: Dissolving star clusters? 7 ( H) = Cr97 =5 = vdb9 =5 (Fig. ). Indeed, when transposed to near-ir colour-colour diagrams (Fig. ), the age and reddening solutions derived forcr97andvdb9consistentlymatchmostoftheirfieldstar decontaminated photometry, but a significant fraction of the PMS stars appears to be very reddened. Most of the very red PMS stars occur along the respective reddening vectors. However, a small fraction appears to present an abnormal excess in (H K s ), especially vdb9, in the ( H) (H K s ) diagram, which may come from PMS stars still bearing circumstellar discs and, thus, with excess in H (e.g. Furlan et al. 9; Bonatto et al. ). ( K S ) = = (H K S ) = = (H K S ) Fig. 7. Colour-colour diagrams showing the decontaminated photometry and the 5 Myr PMS track (Siess, Dufour & Forestini ), set with the derived reddening values (Sect. 4). MS stars lie to the blue side of the diagram. Error bars have been omitted for clarity. Reddening vectors as in Figs. 5 and. this solution, the fundamental parameters of vdb9 are E( H) =.7±.4 (E(B V) =.±. or =.7±.4),(m M) =.9±.4,(m M) O =.7±.4, d =. ±.kpc, and R GC =. ±.kpc, thus. kpc outside the Solar circle. Thus, with the coordinate components (x GC,y GC,z GC ) (.,.,.), vdb9 approximately coincides with the Orion-Cygnus arm (e.g. Momany et al. ). SIMBAD provides 5 stars with spectral type and optical photometry to within of the central coordinates of vdb9 (Table ). These are HRW (B7III, B =.7 and V =.9), BD 7 (B.5V, B =.9 and V =.9), BD 7 (BV, B = 9. and V = 9.5), HIP4 (BV, B = 7. and V = 7.4), and NSV4(BV, B = 9.55 and V = 9.4). The average spectroscopicdistanceofthese starsisd =.4±.kpc,in excellent agreement with the photometric value for vdb 9. For vdb9 we derive m =.5±.,M = 4.5±.4, ( H) =.54±.5,( K s ) =.±.4.Thus, Cr97 is intrinsically brighter and somewhat bluer than vdb 9. Besides the similar MS age, Cr97 and vdb9 have in common a significant age spread ( Myr) implied by the PMS stars, which indicates a non-instantaneous star formation. We found a similar scenario in our previous studies of, e.g. NGC4755, NGC and NGC44. Although with a relatively low foreground absorption ( < ), both cases present a significant fraction of PMS stars redder(and fainter) than predicted by the PMS tracks 5. Structural parameters We derive structural parameters by means of the RDPs, which are the projected stellar number density around the cluster centre. Noise in the RDPs is minimised with colourmagnitude filters (Fig. ), which exclude stars with colours not compatiblewith thoseofthe cluster 7. In previousworks we have shown that this filtering procedure enhances the RDP contrast relative to the background, especially in crowded fields (e.g. Bonatto & Bica 7). Rings of increasing width with distance from the cluster centre are used to preserve spatial resolution along the full radial range. The set of ring widths used is R =.5,.5,.,.5, and 5, respectively for R <.5,.5 R <, R < 5, 5 R <, and R. The residual background level of each RDP is the average number-density of field stars. The R coordinate (and uncertainty) of each ring corresponds to the average position and standard deviation of the stars inside the ring. The resulting RDPs are shown in Fig.. As a caveat we note that the projected stellar distribution of vdb 9 is not smooth and radially symmetric, especially at the outskirts (Fig. 4). Thus, the use of circular rings to build the RDP might introduce some noise at large radii, but with little effect in the central parts. We fit the RDPs with the function σ(r) = σ bg +σ /(+ (R/R c ) ), where σ and σ bg are the central and residual background stellar densities, and R c is the core radius. Applied to star counts, this function is similar to that used by King (9) to the surface-brightness profiles in the central parts of globular clusters. To minimise degrees of freedom, only σ and R c are derived from the fit, while σ bg is previously measured in the surrounding field and kept fixed. The best-fit solutions (and uncertainties) are shown in Fig., and the corresponding structural parameters are given in Table. In addition, we estimate the cluster radius (R RDP ), i.e. the distance from the centre where the cluster RDP and field fluctuations are statistically indistinguishable (e.g. Bonatto & Bica 5), and the density contrast parameter δ c = +σ /σ bg (Table ). R RDP may be taken as an observational truncation radius, whose value depends both on the radial distribution of member stars and the field density. 7 They are wide enough to include cluster MS and PMS stars, together with the photometric uncertainties and binaries (and other multiple systems). Because of the relatively low number of stars in the present OCs, fluctuations in surface-brightness profiles are expected to be higher than those in RDPs. Alternative RDP fit functions are discussed in Bonatto & Bica ().

8 C. Bonatto and E. Bica: Collinder 97 and vdb 9: Dissolving star clusters? Table. Derived structural parameters Cluster σ bg σ R c R RDP δ c σ bg σ R c R RDP ( ) ( ) ( ) ( ) (pc) ( pc ) ( pc ) (pc) (pc) () () () (4) (5) () (7) () (9) () () Cr97.7±. 5.±. 4.±. 4.±4..±.5.5.±. 97.9±.9.5±..±. Cr97.5±.4 4.7±. 5.±.9 4.±4..±.5.5.±. 9.±7.9.±..±. vdb9.±..±.7 4.±..±..±..4.5±.4.±4.4.9±.5.±. vdb9.97±..7±. 5.±.9.±..4±..4.±.4 7.±..±.4.±. Table Notes. Col. : cluster/background density contrast parameter (δ c = +σ /σ bg ) computed with the King-like fit parameters. Col. 7: arcmin to parsec scale. ( ): background level kept fixed; ( ): background level allowed to vary. Within uncertainties, the adopted King-like function provides a reasonable description of the RDPs for R > (Fig. ). The dropin the RDP ofcr97for < R( ) < is probably related to enhanced dust absorption. However, the RDP measured in the innermost region (R <. ) in bothcasespresentsasignificant(> σ)excess,orcusp,over the fit. In old star clusters, this feature has been attributed to a post-core collapse structure, as detected in some globular clusters (e.g. Trager, King & Djorgovski 995). Some Gyr-old OCs, e.g. NGC9 (Bonatto & Bica ) and LK (Bonatto & Bica 9a), also display such a dynamical evolution-related feature. Interestingly, this cusp also occurs in the RDP of some very young (< Myr) clusters, e.g. NGC44 (Bonatto & Bica 9b), NGC (Bica, Bonatto & Dutra ), Pismis5 and NGC9 (Bonatto & Bica 9c). Such time-scales are too short for clusters to evolve into a post-core collapse. Thus, the cusp in young clusters is probably related to molecular cloud fragmentation and/or star formation, and may suggest important early deviation from dynamical equilibrium (Sect. ). As a dynamically-related alternative explanation for the RDP cusps, we note that numerical simulations of clusters that form dynamically cool (subvirial) and with a fractal structure (Allison et al. 9), suggest that mass segregation may occur on timescales comparable to their crossing times, typically a few Myr. Also, clusters that form with significant substructure will probably develop an irregular central region, unless such a region collapses and smoothsout the initial substructure. Finally, we also show in Fig. (bottom panels) the RDPs built after isolating the MS and PMS stars by means of the respective colour-magnitude filters (Fig. ). In both clusters, MS stars are more concentrated than the PMS ones and, as expected, the fraction of PMS stars (with respect to the total number) increases with cluster radius, and largely dominates over the MS stars especially for R > 5. The innermost region of vdb 9 is occupied essentially by PMS stars. The different residual background levels betweenthepmsandmsrdpsreflectthedominantpresence of faint and red stars in the field, with respect to the bright and blue ones. As a consequence, the density-contrast parameter of the MS stars in Cr97 is higher than that of the PMS stars. The absence of MS starsin the innermost RDP bin precludes this comparison for vdb 9. Considering the fit parameters (Table ), the densitycontrast parameter is relatively low in both cases,. < δ c <.. However, when the measured innermost RDP value is considered (Fig. ), it increases to δ c = 5.±4. and δ c = 5.7±4.,respectivelyfor Cr97and vdb9. The latter values are more representative of the visual contrast σ (stars arcmin ) σ (stars arcmin ) Cr97 MS+PMS R core R (arcmin) MS PMS RRDP vdb9 MS+PMS R core R (arcmin) Fig.. Top: Stellar RDPs built for the MS and PMS stars combined, together with the best-fit King-like profile (solid line), the σ uncertainty(light-shaded region) and the background level(shaded polygon). Note the pronounced central cusps. Bottom: RDPs built isolately for the MS and PMS stars. produced by the few relatively bright stars in the central regions of both objects. Alternatively, we fitted the RDPs allowing for variations of the background level. Within the uncertainties, the corresponding parameters (Table ) are compatible with the previous (fixed background level) ones. Taken at face values, the core radii - derived from the King-like fits - of Cr97 (R c.5pc) and vdb9 (R c.pc) are somewhat larger than the median value of the R c distribution derived for a sample of relatively nearby OCs by Piskunov et al. (7). R RDP

9 C. Bonatto and E. Bica: Collinder 97 and vdb 9: Dissolving star clusters? 9. Stellar mass estimate By far, the CMDs of Cr97 and vdb9 are dominated in number by PMS stars, followed by the developing, poorlypopulated MS (Fig. ), which implies that most of the cluster mass is stored in the PMS stars. To compute the cluster masses we work with the field-decontaminated photometry within R R RDP and count the number of stars that belong either to the MS or PMS. For the MS, the corresponding mass of each star is taken from the mass-luminosity relation derived from the isochrone solution (Sect. 4). The MS mass range, in both cases, lies within M M. Summing up these values for all stars produces the total number (n MS ) and mass (m MS ) of MS stars. Because of the differential reddening, individual masses cannot be assigned to the PMS stars. Thus, we simply count the number of PMS stars and adopt an average mass value to estimate n PMS and m PMS. To compute the average PMS mass value we use the Kroupa () initial mass function 9 for PMS masses between.m 7M. The result is m PMS =.M. Taking into account the uncertainty in R RDP (Table ), we obtain 4 M, n PMS = for Cr97 n MS = + 7, m MS = , and m PMS = +7 7 M. Thus, Cr97 contains 5 stars and M clu + 59 M of stellar mass. The values for vdb9 are n MS = + 5, m MS = +5 4 M, n PMS = 7 + 5, and m PMS = 4 +7 M. The cluster vdb9 thus contains stars and M clu M of stellar mass. Both objects present similar values of member MS and PMS stars, and total mass. Besides, as anticipated by the CMDs (Fig. ), the MS is poorly-populated in each case, to the point that the mass stored in the PMS stars is the dominant (74% in Cr97 and 5% in vdb9) component of the cluster mass. However, these masses may still be somewhat higher, since because of the differential reddening and the MASS photometric limit, we may not detect the very-low mass PMS stars. Also, given the presence of some dust and gas in Cr97 and vdb9 (Fig. ), their total cluster masses may be a bit higher than the present estimates. Finally, the fact that most of the mass is stored in the faint PMS stars may explain the significant discrepancies with respect to previous mass determinations made in the optical (Sect. ). 7. Proper motions Given the relative proximity of both clusters (Table ), we can use proper motion data to probe kinematical properties of the member stars. For this purpose we use the Third U.S. Naval Observatory CCD Astrograph Catalog (UCAC, Zacharias et al. 9), which also provides the MASS photometry for each star. Proper motions were obtained for the same central coordinates and extraction radius as those used to extract the MASS photometry (Table ). The region analysed is contained within R =, both forcr97and vdb9,with the comparison field located within. Also, we applied the respective colour-magnitude filters(sect. 5) before computing the proper motion distributions. Conversion from 9 Defined as dn/dm m (+χ), it assumes the slopes χ =. ±.5 for. m(m ).5, χ =. ±. for.5 m(m )., and χ =.5 for m(m ) >.. φ(µ) [ stars arcmin /km s ] φ(µ) [ stars arcmin /km s ] Cr97 vdb9 R< Field Subtraction 9 µ α cos(δ) [km/s] Cr97 vdb9 9 9 µ δ [km/s] Fig.9. Proper motion distribution of Cr97 (top) and vdb 9 (bottom). In both cases, the field contribution (light-gray) is subtracted from that of the R < region (heavy-solid line) to produce the intrinsic distribution (heavy-gray). Gaussians are fitted to the resulting distributions (see text for details). For clarity, error bars are only shown in the subtracted profiles. masyr to kms was based on the respective cluster distances (Table ). The proper motion distributions (φ(µ) = dn/dµ, in units of number of stars per arcmin per kms ), for both the R = and field regions, are shown in Fig. 9. The intrinsic distribution corresponds to the subtraction of the field contribution to the R = region. Cr 97 appears to require two components of similar velocity dispersion to describe its intrinsic proper motion distribution, while a single one applies to vdb9. The bestfits, obtained with gaussians, are shown in Fig. 9, and the corresponding parameters are given in Table. In velocity space, the proper motion distributions of Cr97 have peaks separated by kms and 9kms, respectively in right ascension and declination. The spatial velocity dispersion of a cluster that is in virial equilibrium can be expressed as (Spitzer 97) ( )( MD Reff M pc ) (kms ), where G σ = GMD ηr eff.5 is the gravitational constant, M D is the dynamical mass (assumed to be stored only in stars), η 9.75 is a constant (that roughly depends on the density profile), and R eff is the effective, or half-light, radius. Thus, clusters with

10 C. Bonatto and E. Bica: Collinder 97 and vdb 9: Dissolving star clusters? Table. Gaussian fit parameters Cluster µ α cos(δ) σ α µ δ σ δ σ v (kms ) (kms ) (kms ) (kms ) (kms ) Cr97.±4..±.7.±4..4±4. ±5 Cr97 4.5±5..±. +.±. 9.9±. 5± vdb9.±..±. +.5±..±.7 ± Table Notes. σ α and σ δ are the velocity dispersions in right ascension and declination, respectively. Assuming spherical symmetry, we define σ v = (σ α +σ δ). R RDP (pc) (a) R RDP (pc) (c) NGC44 Cr97 vdb9 M D M and R eff pc are expected to present σ kms. In the last column of Table we provide an estimate of the spatial velocity dispersion (σ v ), computed assuming spherical symmetry for Cr97 and vdb9, σv = (σ α +σδ ). Clearly, both clusters present σ v values much higher than the expected dynamical ones. Such a supervirial state has been linked to the impulsive expulsion of the parental molecular cloud residual gas, driven essentially by the strong winds of massive stars and supernovae (Goodwin & Bastian ), which leads to the high rate of dissolution of the young clusters (e.g. Lada & Lada ). A similar scenario has been observed to occur in the dissolving OC NGC 44, which has a spatial velocity dispersion of σ v 5kms (Chen, de Grijs & Zhao 4), which by far exceeds the expected dynamical value (de Grijs, Kouwenhoven & Goodwin ).. Discussion In previous sections we derived a set of astrophysical parameters for Cr97 and vdb9. We now compare them with the same parameters derived (using the same methods and photometry) for a sample of classical OCs, basically gravitationally bound clusters, characterised by a range of properties (e.g. age, mass, size, etc), and located in different environments. The template sample contains some relatively nearby and bright OCs (Bonatto & Bica 5; Bonatto et al. ) together with a few projected towards the central Galaxy (Bonatto & Bica 7). The young OCs NGC (Bonatto, Santos r. & Bica ), NGC (Bica, Bonatto & Dutra ) and NGC9 (Bonatto & Bica 9b) are included for comparison with gravitationally bound objects of similar ages, while the youngngc44(bonatto & Bica9b)mightbedynamically evolving into an OB association. The comparisons are based on diagnostic-diagrams built with template clusters (Fig. ). Panels (a) and (b) deal with the dependence of cluster (R RDP ) and core (R c ) radii on age, respectively. In terms of R RDP, vdb9, and especially Cr 97, present values higher than the other young OCs in the sample. Cr97 is as big as some Gyr-old OCs. With respect to R c, both objects have a core significantly larger than that of the classical (young and old) OCs, comparable to the core of the disrupting NGC 44. Within uncertainties, Cr97 follows the relation R RDP = Because the velocity dispersions considered here are derived from proper motions, the binaries contribution to σ v (Gieles, Sana & Portegies Zwart ) is expected to be minimal. R core (pc). 4 Age (Myr) (b) M cluster (M O ) (d). R core (pc) σ mo =M O pc σ mo =M O pc Fig.. Diagrams dealing with astrophysical parameters of OCs. Gray-shaded circles: template OCs. The young OC NGC 44 is indicated for comparison purposes. The analytical relations in panels (c) and (d) are discussed in the text. (.9±.) R (.±.) C (panel c), derived with the template sample. Similarly to NGC44, vdb9, on the other hand, is an outlier. For gravitationally bound star clusters in which the projected mass density follows a King-like profile (e.g. Bonatto & Bica ) characterised by a central surface mass-density σ M, the mass inside R RDP can be expressed as a function of R c [, σ M, and the ratio R RDP /R c, as M clus = πr C σ Mln +(R RDP /R c ) ]. With the relation between R c and R RDP implied by panel (c), this equation becomes M clus.σ M R C. Thus, the template OCs are constrained within King-like distributions with < σ M (M pc ) < (panel d). Cr97 and vdb9 (together with NGC 44) are outliers with respect to the distribution of bound OCs in the plane R C M clus, in the sense that their locii imply exceedingly low central mass densities (4 < σ M (M pc ) < ). Significantly smaller core radii and/or higher cluster mass would be required for them to be contained within the classical OC boundaries. Together with the kinematical analysis (Sect. 7), a possible conclusion that can be drawn from this comparison with classical bound OCs is that Cr97 and vdb9 may be evolving into OB associations and/or towards dissolution. Given their young ages, the determining factor for this behaviour is certainly related to another mechanism than age-dependent dynamical evolution. Instead,

11 C. Bonatto and E. Bica: Collinder 97 and vdb 9: Dissolving star clusters? it is probably associated with the primordial star formation process, or even earlier, to the molecular cloud fragmentation. Both clusters may have been formed with dynamically hot stellar components, and thus, are unstable against dissolution (e.g. Goodwin 9). A similar conclusion was reached for some equally young objects, but with very different masses: NGC 44 (Bonatto & Bica 9b) and Bochum (Bica, Bonatto & Dutra ) have comparable stellar masses (M clu M ), while NGC9 and Pismis5 (Bonatto & Bica 9c) are at the low-mass end (M clu M, and M clu M, respectively). Thus, early star cluster dissolution appears to leave detectable imprints also on the structure (e.g. RDP) of clusters as massive as several hundreds solar masses. 9. Summary and conclusions Given the complex interplay among environmental conditions, star-formation efficiency, initial dynamical state of the stellar content, and the total mass converted into stars or expelled, the majority of the embedded clusters do not survive the first few tens of Myr, especially the low-mass ones. In this context, it is important to investigate the structural and photometric properties of OCs (with a mass range and in different environments) that are undergoing this early phase. In the present paper we derive astrophysical parameters and investigate the nature of the young OCs Cr97 and vdb9. Their location in the rd Galactic quadrant minimises contamination by field stars. We work with field-star decontaminated MASS photometry (with errors <. mag), which enhances CMD evolutionary sequences and stellar RDPs, yielding more constrained fundamental and structural parameters. The decontaminated CMDs are characterised by similar properties, an under-populated and developing MS, a dominant fraction (> 75%) of PMS stars, and some differential reddening.inbothcases,themsandpmscmdmorphologies consistently imply a time-spread of Myr in the star formation. Thus, we set the age of Cr97 and vdb9 as constrained within 5 ± 4Myr. Their MS+PMS stellar masses are + 59 M (Cr97) and M (vdb9). By means of the proper motion of the member MS and PMS stars, we estimate their velocity dispersions to be in the range 5kms to kms. Compared to a set of classical bound OCs, both Cr97 and vdb9 appear to have core and cluster radii abnormally large, with R c.5,.9pc and R RDP, pc, respectively for Cr 97 and vdb 9. Structurally, the stellar RDPs follow a cluster-like profile for most of the radial range, except in the central region, where they have a pronounced cusp. At less than about Myr, this cusp is probably related to the star formation and/or molecular cloud fragmentation, and not the product of dynamical evolution. A possible conclusion is that Cr 97 and vdb 9 deviate critically from dynamical equilibrium, and similarly to the equally young NGC44 and Bochum (of comparable mass), and Pismis5 and NGC9 (of significantly lower mass), they are heading towards dissolution. This interpretation is also consistent with the supervirial state of both clusters (as well as the dissolving OC NGC 44), which have velocity dispersions much higher than the σ v kms expected for nearly virialised clusters of similar mass and size as Cr97 and vdb9. In this context, Cr97 and vdb9 may be taken as a link between embedded clusters and young stellar associations. We provide evidence that early star cluster dissolution may be detected, for instance, by means of important and systematic irregularities in the stellar radial density profile of clusters as massive as several M. Studies like the present one are important for a better understanding of the crucial early evolution of embedded star clusters - and the dependence on mass and environment - that rarely result in a classical bound OC or, more frequently, lead to their dissolution into the field. Acknowledgements We thank the referee, Simon Goodwin, for interesting comments and suggestions. We acknowledge support from the Brazilian Institution CNPq. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Centre/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This research has made use of the WEBDA database, operated at the Institute for Astronomy of the University of Vienna. References Allison R.., Goodwin S.P., Parker R.., de Grijs R., Portegies Zwart S.F. & Kouwenhoven M.B.N. 9, ApL, 7,99 Alter G., Balazs B., Ruprecht. & Vanysek. 97, in Catalogue of Star Clusters and Associations, Budapest Akademiai Kiado, nd Edition, ed. by Alter G., Balazs B. & Ruprecht. Bastian N., Gieles M., Lamers H..G.L.M., Scheepmaker R.A. & de Grijs R. 5, A&A, 4, 95 Battinelli P. & Capuzzo-Dolcetta R. 99, MNRAS, 49, 7 Battinelli P., Brandimarti A. & Capuzzo-Dolcetta R. 994, A&AS, 4, 79 van den Bergh S. 9, A, 7, 99 Bessel M.S. & Brett.M. 9, PASP,, 4 Bica E., Bonatto C. & Dutra C., A&A, 49, 9 Bica E., Bonatto C., Barbuy B. & Ortolani S., A&A,, 5 Bica E., Bonatto C. & Camargo D., MNRAS, 5, 49 Bonatto C., Bica E. & Girardi L. 4, A&A,, 57 Bonatto C. & Bica E. 5, A&A, 47, 4 Bonatto C., Santos r..f.c. & Bica E., A&A, 4, 57 Bonatto C., Bica E., Ortolani S. & Barbuy B., A&A,, Bonatto C. & Bica E., A&A, 5, 9 Bonatto C. & Bica E. 7, MNRAS, 77, Bonatto C. & Bica E., A&A, 477, 9 Bonatto C. & Bica E. 9a, MNRAS, 9, 4 Bonatto C. & Bica E. 9b, MNRAS, 94, 7 Bonatto C. & Bica E. 9c, MNRAS, 97, 95 Cardelli.A., Clayton G.C. & Mathis,.S. 99, Ap,, Chen L., de Grijs R. & Zhao.L. 7, A, 4, Clariá.. 974a, Ap, 79, Clariá.. 974b, A&A, 7, 9 Collinder P. 9, AnLun,, Dutra C.M., Santiago B.X. & Bica E., A&A,, 9 Friel E.D. 995, ARA&A 995,, Froebrich D., Scholz A. & Raftery C.L. 7, MNRAS, 74, 99 Furlan E., Watson D.M., McClure M.K., Manoj, P., Espaillat C., D Alessio P., Calvet N., Kim K.H. et al. 9, Ap, 7, 94 Ghez A.M., Salim S., Weinberg N.N., Lu.R., Do T., Dunn.K., Matthews K., Morris M.R. et al., Ap, 9, 44 Gieles M., Sana H. & Portegies Zwart S.F., MNRAS, 4, 75 Girardi L., Bertelli G., Bressan A., Chiosi C., Groenewegen M.A.T., Marigo P., Salasnich B. & Weiss A., A&A, 9, 95 Goodwin S.P. 9, Ap&SS, 4, 59 Goodwin S.P. & Bastian N., MNRAS, 7, 75 Gouliermis D., Kontizas M., Korakitis R., Morgan D.H., Kontizas E. & Dapergolas A., A, 9, 77