Movatterモバイル変換


[0]ホーム

URL:


EDP Sciences logo
0
Subscriber Authentication Point
EDPS Account
Astronomy & Astrophysics (A&A)
Search
Menu
HomeAll issuesVolume 558 (October 2013)A&A, 558 (2013) A53Full HTML
Free Access
Issue
A&A
Volume558, October 2013
Article NumberA53
Number of page(s)8
SectionCatalogs and data
DOIhttps://doi.org/10.1051/0004-6361/201322302
Published online03 October 2013
A&A 558, A53 (2013)

Global survey of star clusters in the Milky Way

II. The catalogue of basic parameters

N. V. Kharchenko1,2,3,A. E. Piskunov1,2,4,E. Schilbach1,S. Röser1 andR.-D. Scholz2

1 Astronomisches Rechen-Institut, Zentrum für Astronomie der Universität Heidelberg, Mönchhofstraße 12-14, 69120 Heidelberg, Germany
2 Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany
e-mail: rdscholz@aip.de
3 Main Astronomical Observatory, 27 Academica Zabolotnogo Str., 03680 Kiev, Ukraine
4 Institute of Astronomy of the Russian Acad. Sci., 48 Pyatnitskaya Str., 109017 Moscow, Russia

Received:17 July 2013
Accepted:9 August 2013

Abstract

Context.Although they are the main constituents of the Galactic disk population, for half of the open clusters in the Milky Way reported in the literature nothing is known except the raw position and an approximate size.

Aims.The main goal of this study is to determine a full set of uniform spatial, structural, kinematic, and astrophysical parameters for as many known open clusters as possible.

Methods.On the basis of stellar data from PPMXL and 2MASS, we used a dedicated data-processing pipeline to determine kinematic and photometric membership probabilities for stars in a cluster region.

Results.For an input list of 3784 targets from the literature, we confirm that 3006 are real objects, the vast majority of them are open clusters, but associations and globular clusters are also present. For each confirmed object we determined the exact position of the cluster centre, the apparent size, proper motion, distance, colour excess, and age. For about 1500 clusters, these basic astrophysical parameters have been determined for the first time. For the bulk of the clusters we also derived the tidal radius. We estimated additionally average radial velocities for more than 30% of the confirmed clusters. The present sample (called MWSC) reaches both the central parts of the Milky Way and its outer regions. It is almost complete up to 1.8 kpc from the Sun and also covers neighbouring spiral arms. However, for a small subset of the oldest open clusters (log t ≳ 9) we found some evidence of incompleteness within about 1 kpc from the Sun.

Key words:globular clusters: general / open clusters and associations: general / Galaxy: stellar content / galaxies: photometry / galaxies: fundamental parameters / galaxies: star clusters: general


The MWSC survey is only available at the CDS via anonymous ftp tocdsarc.u-strasbg.fr (130.79.128.5) or viahttp://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/558/A53

© ESO, 2013

1. Introduction

Star clusters are the main building blocks of the stellar populations in our Galaxy. They are found in high numbers among different Galactic populations, and their astrophysical parameters can be determined with relatively high precision. In the literature about 4000 objects are currently known that can be regarded as Galactic star clusters of various types or candidates. However, about half of them are cited only with a name, approximate coordinates (RA and Dec), and a rough angular size. We call all clusters mentioned in the literature, or more precisely, the clusters of our input sample (defined below)namedclusters henceforth to distinguish this sample from our final sample of clusters for which we determined astrophysical parameters within our Milky Way Star Clusters (MWSC) project.

We started the MWSC survey a few years ago and aimed to build a comprehensive sample of Galactic star clusters with well-determined parameters, spatially complete enough to enable an unbiased study of the content and evolution of the star clusters of our Galaxy. For a reliable cluster membership construction and a reliable cluster parameter determination we used a combination of uniform kinematic and accurate near-infrared (NIR) photometric data gathered in the all-sky catalogue PPMXL (Röser et al. 2010). In the first paper (Kharchenko et al. 2012), hereafter called Paper I, we introduced the survey, explained the underlying motivation, provided a short review of similar studies, described the observational basis of the survey and our data processing pipeline, and presented preliminary results obtained in the second Galactic quadrant. The present paper summarises the results of the full survey carried out for a compiled list of all named clusters, covering the whole sky. The MWSC catalogue of 3006 clusters, the corresponding stellar data with membership probabilities, as well as supplementary material on the full input list of clusters will be available from the CDS. We do not plan, however, to restrict the survey to named clusters only. Our second aim is extending the cluster sample by detecting hitherto unknown clusters on the basis of the PPMXL catalogue. This work is still in progress, and we will present the results in a following paper.

In Sect.2 we briefly describe the basic input data set and the pipeline of the cluster parameter determination and give the current general statistics of the outcome. In Sect.3 we characterise the derived cluster parameters. In Sect.4 we discuss the general properties of the MWSC sample of star clusters. Section5 summarises the results achieved so far.

2. Survey constituents and output statistics

2.1. Observational basis

The basic stellar data for our study were taken from the all-sky catalogues PPMXL (Röser et al. 2010) and 2MASS (Skrutskie et al. 2006). PPMXL gives coordinatesα, δ, and proper motion componentsμα andμδ in the ICRS, and low-accuracy photometry from USNO-B1.0 (Monet et al. 2003) for about 900 million objects down toV ≈ 20. For some 400 million entries the catalogue contains accurateJ,H,Ks magnitudes from 2MASS (Skrutskie et al. 2006). For the current survey we used an improved merged version of these catalogues called hereafter 2MAst, see Paper I for a description of its construction. We used 2MAst to confirm clusters from our input list and to determine cluster parameters in the astrometric and photometric systems that are homogeneous over the whole sky.

Table 1

Comparison of the COCD and MWSC surveys.

In Table1 we compare the basics of the current survey with our previous work on open clusters (Kharchenko et al. 2005a,b), here referred to as COCD (catalogue of open cluster data), which was based on the shallower ASCC-2.5 catalogue. The comparative depth of the new survey has increased by at least 3 mag by considering the transitions from the optical to the infra-red. This provides the respective increase of the limiting and completeness distances. Including the NIR photometric data has, as a consequence, led to an increased relative weight of the clusters of extreme age (both young, normally heavily reddened, and old, which are relatively faint in the optical).

The target list was compiled from sources available in the literature. As the primary source we used the data from COCD. For additional optical clusters and associations the data were taken from the catalogue ofDias et al. (2002), called hereafter as DAML02 (version 3.1, 24/nov/2010), and fromMelnik & Dambis (2009). For clusters detected in the NIR the information came fromBica et al. (2003a,b),Dutra et al. (2003),Froebrich et al. (2007,2010),Bukowiecki et al. (2011), and other sources. We described the input list in detail in Paper I. Furthermore, we included embedded and globular clusters taken from the catalogues byLada & Lada (2003) andHarris (1996)(edition 2010)1 in our target list.

2.2. Pipeline

The pipeline has been described in detail in Paper I and here we only refer to its basic features. The main purpose of the pipeline is i) to reveal clusters from the fore-/background using kinematic, photometric, and spatial criteria; ii) to construct a list of cluster members; and iii) to determine basic cluster parameters. For each cluster, the direct outcome of the pipeline were the improved coordinates of the cluster centre and the apparent sizes of different morphological parts of the cluster, the average proper motion, the distance, the reddening, and the age. As a theoretical basis, we used recent Padova stellar models with isochrones computed with the CMD2.2 on-line server2, whereas the pre-MS isochrones were based on the models ofSiess et al. (2000) transformed by us to theJHKs photometric system. In the pipeline, we used several diagrams of kinematic (vector point diagram, magnitude-proper motion relation) and photometric (colour–magnitude, two-colour andQJHKs-colour diagrams) data. The membership probabilities of stars in the diagrams were determined from their location with respect to the reference sequences (represented by isochrones in photometric diagrams, or the average cluster proper motion in kinematic diagrams), which themselves depend on the cluster parameters we intended to find. Hence, this required an iterative approach, allowing us to successively improve both cluster membership and cluster parameters. As an initial approximation we either used data from the input list or made our own estimates by eye if the input parameters were obviously incorrect or not available. As a rule, the process converged after a few iterations. The pipeline provided spatial, kinematic, and photometric membership probabilities for each star in the sky area around a cluster. In total, about 64 million stars were retrieved. For each confirmed cluster (see below), a coherent set of basic cluster parameters was determined.

The final information on membership was then also used to determine of other characteristics of a cluster, such as parameters of theKing (1962) density profile and the radial velocities (RVs). King parameters (core and tidal radii, and normalisation factors) were computed for almost all clusters. The RV of a cluster was determined by averaging RVs of individual cluster members, for which these data were available in the literature. Altogether, RVs were obtained for about 30% of the MWSC objects.

Since open clusters belong to a relatively young Milky Way population, their metallicities are high (of the order of the solar value), and, as the stellar model computations show (see for example the Padova isochrones of the same age and differentZ), their variations only mildly influence the derived cluster parameters such as age and distance. Therefore we adopted the isochrones with solar metallicity (Z = 0.019) to determine the parameters of open clusters. When a metallicity measurement was found in the literature for a given cluster, its value was copied to our catalogue. Since we used the data from the bibliographic data collection of DAML02 as the main source of open cluster metallicities, the MWSC metallicities, unlike other cluster parameters, should be regarded as highly heterogeneous. For globular clusters, the distances, reddenings and ages determined from colour–magnitude diagrams (CMDs) strongly depend on the adopted metallicity, which means that the approach with solar metallicity is not applicable. Moreover, for the majority of globular clusters we observed only the brightest red giant members in 2MAst. The steepness of the giant branch makes it difficult to estimate their distances reliably. Therefore, we adopted the published distances and reddenings from theHarris (2010) catalogue, and determined ages and metallicities from the theoretical isochrones that provided the best fit to the observed CMD.

Table 2

Statistics of MWSC objects.

2.3. Object statistics

In total the input list of named clusters comprised 3784 entries. In our analysis, however, we did not find real counterparts for 21% of the targets in 2MAst. The corresponding statistics is given in Table2. About one half of the uncomfirmed targets seem to be dubious objects that show no reasonable sequences either in the kinematic or photometric diagrams, therefore we concluded that they are not physical clusters. Therefore, no membership probability was provided for the stars in the corresponding sky areas. Another 4% of entries from the input list were not seen in 2MAst since they are too faint for this catalogue. Moreover, very many targets in the input list turned out to be double or even multiple entries. In different sources in the literature, the same object was referred to by different names and/or with slighly different positions. We kept unconfirmed entries in the output list and marked them with corresponding flags and notes.

About 79% of the objects considered were confirmed as real clusters. In total, the 3006 confirmed clusters contain about 400 000 stars that are most probably cluster members, that is, stars with kinematic and photometric membership probabilities higher than 60%. On average, this corresponds to about 130 stars per cluster (to be compared with the about 20 highly probable members per cluster in COCD).

The MWSC survey includes clusters of different types, which are listed in Table3. There are 142 globular clusters fromHarris (2010), which make up about 90% of entries of this catalogue. Five new globular cluster candidates are star clusters observed in the direction of the Galactic centre. We classified them as globular clusters because of their CMDs. We will discuss globular clusters from the MWSC survey in more detail in a dedicated paper (Kharchenko et al. 2013, in prep.).

Associations in the MWSC are defined as plain assemblies of early-type stars that appear like star clusters, but do not show a pronounced central concentration in the density profile of the member distribution on the sky. We stress here that both our pipeline and the basic input data are highly unsuitable for a study of extended star-forming regions with multi-centre structures, which are usually referred to by the term association. Therefore, many classical associations cannot be processed with our pipeline, and were flagged as dubious objects in the MWSC. For 21 compact cluster-like assosiations fromMelnik & Dambis (2009) the pipeline provided acceptable results however. We refer to targets as candidates in Table3 if they were classified as clusters in the original input list, but did not show a clear central concentration although they contained early-type stars. Finally, individual clusters observed within large and complex associations that successfully passed our pipeline were classified as confirmed open clusters.

Furthermore, we considered several sub-groups of open clusters (see Table3). The descriptors “cluster with nebulosity”, “moving group”, and “remnant” were adopted from the original sources. We also marked old scarce clusters with irregular density profiles as “remnant candidates”.

Table 3

Cluster types in the MWCS survey.

Table 4

Statistics of cluster parameters in the MWSC, the first seven parameters are called basic parameters, the other three additional parameters.

3. Cluster parameter scope

One of the main goals of our survey is to determine a uniform and homogeneous set of astrophysical cluster parameters for all named clusters.

Our set of cluster parameters includes positions (coordinates of the centre, distance) that can be converted into 3D coordinates, structure data (apparent sizes of distinct cluster parts, tidal parameters), kinematic information (proper motions, RVs) related to 3D velocity components, and astrophysical parameters (age, reddening, and sometimes cluster metallicity). In Table4 we show the number of clusters for which a given parameter was determined in the MWSC survey (Col. 2), and the number of clusters for which the parameter was estimated for the first time within the MWSC project (Col. 3). Column 4 gives the corresponding percentage of these clusters. Except for most of the metallicities, all parameters listed (including coordinates and apparent radii) were re-determined or newly determined in our study.

We distinguish in Table4 basic parameters determined via the MWSC pipeline for all confirmed objects (first seven rows), and additional ones (last three rows) that required a different approach and were obtained for a fraction of clusters only. To determine the King tidal parameters we use information on the projected density distribution of cluster members, which is described in Paper I and the references therein. However, not all of the confirmed clusters follow King profiles. Especially for apparently large loose objects with poor radial density profiles one fails to fit a King model.

Originally, RVs were known for only 670 clusters of our sample (Kharchenko et al. 2007b;Dias et al. 2002). To update the RV data, we used information on the membership and looked for the relevant RV measurements in CRVAD-2 (Kharchenko et al. 2007a) and in SIMBAD. This led to 75 and 14 additional cluster RVs, respectively. For another 40 clusters, new RVs were determined byConrad et al. (2013) from dedicated RAVE observations proposed earlier on the basis of COCD data. Additionally, we cross-matched about 63.6 million entries in 3006 MWSC cluster areas with the about 2.7 million spectroscopic measurements in the SDSS DR9 (Ahn et al. 2012) and found only a small overlap of about 20 000 measurements in 93 cluster areas. Only in 62 of these areas at least one cluster member had an RV measured in SDSS. Finally, we identified about one hundred cluster members as IRAS point sources with measured RVs (Bronfman et al. 1996) and determined RVs for another 92 star clusters.

Metallicities for open clusters were taken from the DAML02 andConrad et al. (2013). For globular clusters we had to vary the cluster metallicity parameter to achieve the best agreement between the observed colour–magnitude diagrams and isochrones. The corresponding pipeline modifications for globular clusters will be described in more detail in the forthcoming paper (Kharchenko et al. 2013, in prep.). From Table4 one notes that the cluster metallicity is the most rarely known parameter in the MWSC survey.

In Table5 we give the results of a comparison of the basic MWSC parameters with the corresponding data found in the literature. For the listed samples of 3σ-clipped deviants, we applied the method of a least-squares bisector (Isobe et al. 1990), which provides the coefficients of a regressiony = a + b x, and the standard deviationσ around the bisector. The parametersa andb indicate the magnitude of a possible bias between the compared data sets, andσ is the measure of the spread of the data around the regeression line. We found (Table5) that the scale coefficientb is close to unity for all basic parameters, whereas the zero-point differencesa are formally significant.The bias in zero-points can be caused by uncertainities of our parameters as well as of the literature values. Taking into account the accuracy of the input data and the spreadσ, the differences are not dramatic however.

Table 5

Comparison of the bacic cluster parameters in the MWSC survey and data in the literature

Table 6

Estimated accuracy of the derived MWSC cluster parameters

thumbnailFig. 1

Distribution functions of MWSC clusters with respect to the different astrophysical parameters determined in this paper.

thumbnailFig. 2

Distribution of clusters of the MWSC survey projected onto the GalacticXY-plane. Blue dots are open clusters and associations, red triangles mark globular clusters. The dotted circle shows the border of the Galactic disk (diameter ~20 kpc) as indicated by distant open clusters. The thick solid cyan and magenta sections of the spirals indicate the position of local spiral arms as defined by the COCD clusters (Piskunov et al. 2006). They are formally extended to the edge of the disk with thin solid curves. The light (yellow) thick circle around the Sun with radius of 1.8 kpc marks the completeness limit of the survey. A cross at (X,Y) = (8.5,0) kpc indicates the position of the Galactic centre.

Our estimates on the accuracy of the cluster parameters are given in Table6. There are two ways to estimate the achieved accuracy. For parameters that are average or fitted values (such as proper motions, RVs, or tidal radii, but also ages, if at least two stars were used in fitting to the evolved portion of the isochrone), one can estimate the spread of cluster members around the fit. This does not take into account the accuracy of the pipeline itself and gives us an internal estimate (lower limit) of the error only. For the overall characterisation of the internal errors of the full cluster sample, we used the most frequent value (the mode) of the error distribution. These values are shown in the second column of Table6. In the third column we show the external estimates of the derived parameter errors, which come from the comparison of our values with those published in the literature. For this comparison we used the dispersions given in the fourth column of Table5, and assumed that our data have at least the same accuracy as those from the literature. This comparison includes inaccuracies originating from our pipeline as well as from the technique/data used to determine the literature parameters. Therefore, these external errors represent an upper limit of our error estimates.

thumbnailFig. 3

Observed (derived) kinematic data of open clusters within | b| ≤ 20° versus Galactic longitude for four ranges of projected distancesdXY (as indicated at the top). Theupper panelsshow proper motions, themiddle panelstangential velocities, and thebottom panelsRVs. Individual clusters are marked by thin horizontal lines with vertical error bars (in the RV panel, error bars are given only for clusters whose rms-errors could be estimated). For a given distance range, the curves show the systematic velocity components due to solar motion and local Galactic rotation computed with COCD-based values (Piskunov et al. 2006). The dashed curves correspond to the contributions at the smallest distance of a given distance range, whereas solid curves show this contribution at the largest distance. For example, in theupper left panelthe dashed curve shows systematic variations versus longitude in proper motions expected for a cluster at 0.1 kpc from the Sun, whereas the solid curve shows these variations for a cluster at 1.0 kpc.

As one can see from Table6, the mean cluster proper motions are accurate within 1 mas/y, and proper motions of northern clusters are somewhat more accurate than those of the southern clusters. Since distances and reddenings were determined from fitting the isochrones by eye to the observed colour–magnitude distributions, we were not able to determine the internal accuracy of the procedure. Therefore, we relied on the comparison with previously known values, which indicates a good quality of our derived distances and reddenings, even better than our preliminary quality assessment in Paper I. For clusters with well-observed turn-off points, we were able to estimate the spread of cluster members around the selected isochrone and consequently the internal uncertainty of the derived ages caused by this spread. We found that this uncertainty depends on age and is smaller for older clusters with more populated and clear turn-offs. The external age errors are considerably larger than the internal ones. One possible reason for this discrepancy are the different theoretical isochrones applied. But a particularly strong impact on the age determination can be produced by wrongly including non-members (or vice versa by neglecting true members) above the turn-off. This occurs especially when only photometric membership is applied. Because we combined photometric and proper motion membership, our results are probably in general more reliable than those from pure photometric membership studies (as discussed in Paper I). Tidal radii are accurate to within 25%. The distribution of the RV-errors peaks near 1 km s-1 and most of the clusters haveεRV better than ±5 km s-1. However, there is a caveat on the mean error estimation of the RV of a cluster. The uncertainty may be larger if RVs of only a few stars are available, and if these stars also have low membership probability. For 446 out of 953 clusters with RVs, more than one member was measured, whereas the cluster RVs of 279 and 228 clusters were based on only one or an unknown number of members (if not given in the literature), respectively. The RVs of 183 clusters have no reported errors.

In Fig.1 we show histograms of cluster parameters determined in the MWSC survey. Note that these histograms characterise the distribution of all 3006 clusters in our catalogue, not the actual cluster population in the solar neighbourhood. We found that the typical angular radiusr2 of the clusters in our sample isr2 ≈ 3...12 arcmin, or that the typical tidal radiirt are between 5...10 pc. Most proper motions lie within ±10 mas/y, and the RVs within ±100 km s-1. Though the typical distance is about 2 kpc, a long tail is observed in the distribution of distant clusters up to more than 10 kpc from the Sun. The reddening distribution is almost flat forE(J − K) < 0.3 (E(B − V) < 0.7) mag with a slowly decreasing tail to higher values. The MWSC survey is dominated by older clusters. More than 45% of the clusters have ages between 400 Myr and 2 Gyr. This can be attributed to the NIR basis of the survey since the older clusters with red giants are sufficiently prominent even at large distances from the Sun. On the other hand, the local maxima in the age distribution both for extremely old and extremely young clusters are a consequence of including globular clusters on the one hand and dedicated samples of embedded NIR clusters on the other.

4. General description of the cluster sample

The spatial distribution of clusters from our MWSC survey in the GalacticXY −plane is shown in Fig.2. The plot includes both open and globular clusters (note that the five most distant globular clusters do not appear in Fig.2 since they are outside the plotted frame). The open clusters cover a wide range of galactocentric distances and reach the central region of the Galaxy as well as its fringe at about 20 kpc from the Galactic centre.

Figure3 illustrates the kinematic properties of open clusters selected in the MWSC survey within 4 kpc from the Sun. We show the distributions of proper motions and tangential and RVs versus Galactic longitude for different ranges of distances. As a reference, we also plot the systematic components of the velocities due to solar motion and Galactic rotation for heliocentric distances given at the top of each panel. The curves were computed with the corresponding parameters from COCD (seePiskunov et al. 2006) based on the Hipparcosdata. The top row of Fig.3 shows that the cluster proper motions coincide well with the predicted proper motion distribution. The systematic contribution of solar motion is relevant only for nearby clusters (dXY = d·cosb ≪ 1.0 kpc). At larger distances the main systematic contributor is the Galactic rotation, and, on average, the observed proper motions of clusters follow this closely. This behaviour can be regarded as independent proof that cluster proper motions from the MWSC survey are in the inertial system represented by Hipparcos. The tangential velocities (the middle row of Fig.3) and RVs (the bottom row of Fig.3) also show an expected behaviour, which confirms that the distances and RVs of clusters in the MWSC survey are sufficiently reliable and can be used for studies of the Galactic disk kinematics. However, with increasing distances (dXY > 2.0 kpc) we found stronger deviations of the observed velocities from the simple rotation model, especially for proper motions and tangential velocities. This effect is probably also present in the RVs, but the number of data points is too low for a conclusive discussion. We cannot exclude the possibility that the accuracy of the cluster parameter determination may decrease at larger cluster distances, or that neglecting the radial metallicity gradient in the disk can bias the cluster distances. On the other hand, this effect may indicate that the adopted aproximation of the rotation law is valid for the solar neighbourhood, but is not sufficient to describe the Galactic rotation at larger distances from the Sun.

Completeness is one of the most important characterictics of a sample used in a statistical study. In Fig.4 we show the surface density of open clusters in the MWSC survey as a function of their distancedXY in the Galactic plane. The upper panel shows young (log t < 7.9, blue curve) and moderately young clusters (log t = 7.9...8.3, green curve), the bottom panel moderately old (log t = 8.3...9.0, brown curve) and old clusters (log t > 9.0, red curve). Additionally, we plot the total distribution of all open clusters in both panels (black curve). The total distribution is almost flat up to 1.8 kpc and steadily declines towards larger distances. This shows the increasing incompleteness of the MWSC survey at projected distances larger than 1.8 kpc. The average surface density of clusters withdXY < 1.8 kpc is 121 kpc-2, which is close to the 114 kpc-2 determined previously from the sample of COCD clusters within the smaller completenes area of 0.85 kpc (Piskunov et al. 2006).

The distributions of clusters of different ages follow the general trend, although with a few exceptions. A significant excess of the youngest clusters at about 400 pc from the Sun is related to the Orion star formation complex. For the youngest and oldest clusters the distributions decrease more slowly atdXY > 1.8 kpc because these clusters may contain absolutely bright stars, and therefore they can be observed in the NIR at larger distances than clusters of moderate ages. Another prominent feature appears in the distribution of the oldest (red curve in the bottom panel of Fig.4) clusters: the surface density of nearby clusters increases with increasingdXY up to distances of about 1.1 kpc from the Sun. This effect hints at an incompleteness of the input list at high Galactic latitudes. Indeed, one of the most important sources of input data for the MWSC survey was the cluster list byFroebrich et al. (2007) obtained from a systematic search for clusters in the 2MASS catalogue at low Galactic latitudes | b| < 20°. Due to the strong concentration of younger clusters to the Galactic plane, this limitation has no serious consequences for their completeness in the MWSC survey. The oldest clusters show a larger scale height and a larger scattering in theZ-direction however (see below). Therefore, they could be incompletely represented in the MWSC survey, especially in the solar neighbourhood. Assuming that the real surface density of the oldest clusters is equal to the average density observed in the rangedXY = 1.1...2 kpc, we estimate that ≈40 old clusters are missing within\hbox{$d_{XY}\leqslant 1$} kpc. This effect must be taken into account in statistical studies based on the MWSC survey.

thumbnailFig. 4

Distribution of the surface density ΣXY of star clusters versus distancedXY from the Sun projected onto the Galactic plane. The density distribution of all clusters is given in black in both panels. Inthe upper and lowerpanels we show the distributions of younger (blue: log t < 7.9, green: log t = 7.9...8.3) and older clusters (brown: log t = 8.3...9.0, red: log t > 9.0), respectively. The dotted vertical line indicates the adopted completeness limit, the dashed horizontal lines correspond to the average surface density within the completeness limit of clusters in different age groups.

thumbnailFig. 5

Distribution of clusters along the GalacticZ-axis as a function of age. Open and globular clusters are marked with small blue crosses and open red circles, respectively. Theright panelshows a zoomed-in segment of the left panel.

The cluster distribution from the MWSC survey along the GalacticZ-axis is shown in Fig.5 as a function of cluster age. Here we consider both open and globular clusters with the parameters determined in this paper. The most striking feature of the distribution is the smooth transition from the oldest open clusters to the youngest globulars at ≈6 Gyr. Within the age ranget ≈ 3...7 Gyr, open and globular clusters show a similar vertical distribution. For older globular clusters, the vertical scattering increases steadily with age.

5. Summary and conclusions

Our MWSC project aimed at completing a full survey of star clusters in the wider neighbourhood of the Sun. We provided a catalogue containing the basic astrophysical data for all clusters in this survey, that is, exact positions of the cluster centre, proper motions, apparent radii, distances, reddenings, and ages. All these quantities were determined from data in the stellar all-sky surveys PPMXL and 2MASS. Because these latter surveys give homogeneous data sets all over the sky, and because a uniform pipeline has been applied to all objects, the astrophysical quantities derived in this paper have a uniform and homogeneous nature. A reliable membership determination is the basic issue; all the derived astrophysical parameters rely on this.

The starting point of the MWSC survey was an input list of targets compiled from the literature. The 3784 objects in the input list were called named clusters throughout the paper. From our analysis based on PPMXL and 2MASS data, we found that 3006 (79%) of the named clusters are related to real objects. The vast majority of confirmed objects are open clusters, but due to the content of the input list, stellar associations and globular clusters were also found. In addition to the basic astrophysical parameters mentioned above, we also determined tidal parameters for the bulk of the clusters (98%). From RV measurements of cluster members available in the literature, we estimated additional mean RVs for more than 30% of the MWSC clusters.

For about 1500, or 50% of the confirmed MWSC clusters, we presented basic astrophysical parameters for the first time. We also compared our results with literature data for subsets of clusters where possible. No severe systematic differences were found. Considering the size of the cluster sample and the uniform and homogeneous nature of the cluster parameters, the MWSC survey is unprecedented.

Our sample of MWSC clusters covers a large section of the Galactic disk and reaches the very centre of the Milky Way as well as its outer regions. The area of data completeness now reaches the neighbouring spiral arms, which allows comparative studies of the cluster population in the inter- and intra-arm regions. Our sample of open clusters is almost complete up to a distance of about 1.8 kpc from the Sun, except for the subset of the oldest open clusters (log t ≳ 9), where we found evidence of incompleteness within 1 kpc from the Sun. We attribute this effect to the incompleteness of our input list, that is, the data from the literature. This stimulated a search for clusters missing in the literature. This work is in progress at the moment of writing, but we will soon be able to publish a list that completes the current data.


Acknowledgments

This study was supported by DFG grant RO 528/10-1, and by Sonderforschungsbereich SFB 881 “The Milky Way System” (subproject B5) of the German Research Foundation (DFG). We acknowledge the use of the Simbad database, the VizieR Catalogue Service and other services operated at the CDS, France, and the WEBDA facility, operated at the Institute for Astronomy of the University of Vienna. We thank the referee, Bruce A. Twarog, for his useful comments.

References

  1. Ahn, C. P., Alexandroff, R., Allen de Prieto, C., et al. 2012, ApJS, 203, 21 [NASA ADS] [CrossRef][Google Scholar]
  2. Bica, E., Dutra, C., Soares, J., & Barbuy, B. 2003a, A&A, 404, 223 [NASA ADS] [CrossRef] [EDP Sciences][Google Scholar]
  3. Bica, E., Dutra, C. M., & Barbuy, B. 2003b, A&A, 397, 177 [NASA ADS] [CrossRef] [EDP Sciences][Google Scholar]
  4. Bronfman, L., Nyman, L.-A., & May, J. 1996, A&AS, 115, 81 [NASA ADS][Google Scholar]
  5. Bukowiecki, Ł., Maciejewski, G., Konorski, P., & Strobel, A. 2011, Acta Astron., 61, 231 [NASA ADS][Google Scholar]
  6. Conrad, C., Scholz, R.-D., Kharchenko, N. V., et al. 2013, A&A, accepted, DOI:10.1051/0004-6361/201322070[Google Scholar]
  7. Dias, W. S., Alessi, B. S., Moitinho, A., & Lépine, J. R. D. 2002, A&A, 389, 871 [NASA ADS] [CrossRef] [EDP Sciences][Google Scholar]
  8. Dutra, C., Bica, E., Soares, J., & Barbuy, B. 2003, A&A, 400, 533 [NASA ADS] [CrossRef] [EDP Sciences][Google Scholar]
  9. Froebrich, D., Scholz, A., & Raftery, C. L. 2007, MNRAS, 374, 399[Google Scholar]
  10. Froebrich, D., Schmeja, S., Samuel, D., & Lucas, P. W. 2010, MNRAS, 409, 1281 [NASA ADS] [CrossRef][Google Scholar]
  11. Harris, W. E. 1996, AJ, 112, 1487 [NASA ADS] [CrossRef][Google Scholar]
  12. Harris, W. E. 2010 [arXiv:1012.3224][Google Scholar]
  13. Isobe, T., Feigelson, E. D., Akritas, M. G., & Babu, G. J. 1990, ApJ, 364, 104 [NASA ADS] [CrossRef] [EDP Sciences][Google Scholar]
  14. Kharchenko, N. V., Piskunov, A. E., Röser, S., Schilbach, E., & Scholz, R.-D. 2005a, A&A, 438, 1163 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet][Google Scholar]
  15. Kharchenko, N. V., Piskunov, A. E., Röser, S., Schilbach, E., & Scholz, R.-D. 2005b, A&A, 440, 403 [NASA ADS] [CrossRef] [EDP Sciences][Google Scholar]
  16. Kharchenko, N. V., Scholz, R.-D., Piskunov, A. E., Roeser, S., & Schilbach, E. 2007a, VizieR Online Data Catalog: III/254[Google Scholar]
  17. Kharchenko, N. V., Scholz, R.-D., Piskunov, A. E., Roeser, S., & Schilbach, E. 2007b, VizieR Online Data Catalog: III/2801[Google Scholar]
  18. Kharchenko, N. V., Piskunov, A. E., Schilbach, E., Röser, S., & Scholz, R.-D. 2012, A&A, 543, A156 [NASA ADS] [CrossRef] [EDP Sciences][Google Scholar]
  19. King, I. 1962, AJ, 67, 471 [NASA ADS] [CrossRef][Google Scholar]
  20. Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57 [NASA ADS] [CrossRef][Google Scholar]
  21. Melnik, A. M., & Dambis, A. K. 2009, MNRAS, 400, 518 [NASA ADS] [CrossRef][Google Scholar]
  22. Monet, D. G., Levine, S. E., Canzian, B., et al. 2003, AJ, 125, 984 [NASA ADS] [CrossRef][Google Scholar]
  23. Piskunov, A. E., Kharchenko, N. V., Röser, S., Schilbach, E., & Scholz, R.-D. 2006, A&A, 445, 545 [NASA ADS] [CrossRef] [EDP Sciences][Google Scholar]
  24. Röser, S., Demleitner, M., & Schilbach, E. 2010, AJ, 139, 2440 [NASA ADS] [CrossRef][Google Scholar]
  25. Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 [NASA ADS][Google Scholar]
  26. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef][Google Scholar]

All Tables

Table 1

Comparison of the COCD and MWSC surveys.

Table 2

Statistics of MWSC objects.

Table 3

Cluster types in the MWCS survey.

Table 4

Statistics of cluster parameters in the MWSC, the first seven parameters are called basic parameters, the other three additional parameters.

Table 5

Comparison of the bacic cluster parameters in the MWSC survey and data in the literature

Table 6

Estimated accuracy of the derived MWSC cluster parameters

All Figures

thumbnailFig. 1

Distribution functions of MWSC clusters with respect to the different astrophysical parameters determined in this paper.

In the text
thumbnailFig. 2

Distribution of clusters of the MWSC survey projected onto the GalacticXY-plane. Blue dots are open clusters and associations, red triangles mark globular clusters. The dotted circle shows the border of the Galactic disk (diameter ~20 kpc) as indicated by distant open clusters. The thick solid cyan and magenta sections of the spirals indicate the position of local spiral arms as defined by the COCD clusters (Piskunov et al. 2006). They are formally extended to the edge of the disk with thin solid curves. The light (yellow) thick circle around the Sun with radius of 1.8 kpc marks the completeness limit of the survey. A cross at (X,Y) = (8.5,0) kpc indicates the position of the Galactic centre.

In the text
thumbnailFig. 3

Observed (derived) kinematic data of open clusters within | b| ≤ 20° versus Galactic longitude for four ranges of projected distancesdXY (as indicated at the top). Theupper panelsshow proper motions, themiddle panelstangential velocities, and thebottom panelsRVs. Individual clusters are marked by thin horizontal lines with vertical error bars (in the RV panel, error bars are given only for clusters whose rms-errors could be estimated). For a given distance range, the curves show the systematic velocity components due to solar motion and local Galactic rotation computed with COCD-based values (Piskunov et al. 2006). The dashed curves correspond to the contributions at the smallest distance of a given distance range, whereas solid curves show this contribution at the largest distance. For example, in theupper left panelthe dashed curve shows systematic variations versus longitude in proper motions expected for a cluster at 0.1 kpc from the Sun, whereas the solid curve shows these variations for a cluster at 1.0 kpc.

In the text
thumbnailFig. 4

Distribution of the surface density ΣXY of star clusters versus distancedXY from the Sun projected onto the Galactic plane. The density distribution of all clusters is given in black in both panels. Inthe upper and lowerpanels we show the distributions of younger (blue: log t < 7.9, green: log t = 7.9...8.3) and older clusters (brown: log t = 8.3...9.0, red: log t > 9.0), respectively. The dotted vertical line indicates the adopted completeness limit, the dashed horizontal lines correspond to the average surface density within the completeness limit of clusters in different age groups.

In the text
thumbnailFig. 5

Distribution of clusters along the GalacticZ-axis as a function of age. Open and globular clusters are marked with small blue crosses and open red circles, respectively. Theright panelshows a zoomed-in segment of the left panel.

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.


[8]ページ先頭

©2009-2025 Movatter.jp