Astronomy Reports, Vol. 47, No. 4, 2003, pp. 263–275. Translated from Astronomicheski˘ı Zhurnal, Vol. 80, No. 4, 2003, pp. 291–303. Original Russian Text Copyright c 2003 by Kharchenko, Pakulyak, Piskunov.
The Subsystem of Open Clusters in the Post-Hipparcos Era: Cluster Structural Parameters and Proper Motions
N. V. Kharchenko1,L.K.Pakulyak1,andA.E.Piskunov2 1Main Astronomical Observatory, ul. Akademika Zabolotnogo 27, Kiev, 03680 Ukraine 2Institute of Astronomy, Russian Academy of Sciences, ul. Pyatnitskaya 48, Moscow, 109017 Russia Received August 23, 2002; in final form, October 10, 2002
Abstract—The wide neighborhoods of 401 open clusters are analyzed using the modern, high-precision, homogeneous ASCC-2.5 all-skycatalog. More than 28 000 possible cluster members (including about 12 500 most probable members) are identified using kinematic and photometric criteria. Star counts with the ASCC-2.5 and USNO-A2.0 catalogs are used to determine the angular and linear radii of the cluster cores and coronas, which exceed the previouslypublished values byfactors of two to three. The segregation (differing central concentration) of member stars bymagnitude is observed. The mean proper motions are determined directlyin the Hipparcos system for 401 clusters, for 183 of them for the first time. The heliocentric distances of 118 clusters are determined for the first time based on color–magnitude diagrams for their identified members. c 2003 MAIK “Nauka/Interperiodica”.
1. INTRODUCTION velocity, color–magnitude diagram, stellar content, and age. Open star clusters are typical Galactic disk- The last decade has witnessed sweeping qualita- population objects. At the same time, the spatial tive and quantitative developments in the creation of coordinates, velocities, and ages of open clusters can all-skycatalogs of positional, kinematic, and pho- be determined much more accuratelythan for other tometric stellar data. The results of the Hipparcos Population-I objects. The reason is that, usually, up space mission formed the basis for several catalogs to several hundred stars can be observed in each (Hipparcos, Tycho-1 [1], and Tycho-2 [2]), which cluster, and the spatial–kinematic parameters of all contain high-precision stellar data and establish an the cluster stars are either virtuallyidentical or obey inertial reference frame and two-color photometric certain relations (spectroscopic and photometric da- system over the entire sky. Many ground-based cata- ta). This makes it possible to select cluster members logs have alreadybeen reduced to these systems,in- among all stars seen in the direction of a cluster, cluding the currentlymost massive and deepest cat- making these objects an especiallyvaluable tool alog, containing more than 500 million stars down to for analyzing the structure and evolution of stellar a limiting magnitude of B 22m—the USNO-A2.0 systems such as individual clusters or the Galactic catalog [3] (hereafter, the USNO catalog). Recently, disk as a whole. Kharchenko [4] combined all the data from modern, The principal drawback of clusters as objects for high-precision, astrometric catalogs containing stars such studies is the inhomogeneity, insufficient ac- down to V =12− 14m into the “All-SkyCombined curacy, and incompleteness of stellar data. This is Catalog of Astrometric Data for 2.5 Million Stars” due to the fact that the available stellar data—proper (hereafter, the ASCC). motions and photometric parameters being the most About 1500 open clusters have been discovered in abundant among them—have been obtained in clus- the Galaxyso far; however, onlyafter the Hipparcos ter neighborhoods byvarious authors using di fferent mission had been successfullycompleted and its re- instruments and different instrumental systems with sults implemented in practice did it become possible different degrees of completeness and accuracy. Fur- to create catalogs of open-cluster parameters deter- thermore, the apparent surface distribution of cluster mined in systems that are uniform over the entire sky members can be distorted byinterstellar clouds along in terms of the observational material and/or methods the line of sight and/or inside young open clusters. used. All this makes the correct identification of cluster The most complete source of geometrical and as- member stars rather difficult, distorting the principal trophysical parameters of open clusters is the Lund parameters of the cluster: its spatial structure, space Catalog of Open Cluster data (LCOC), compiled by
1063-7729/03/4704-0263$24.00 c 2003 MAIK “Nauka/Interperiodica” 264 KHARCHENKO et al. Lynga˚ [5] and containing a total of 1150 objects. precision stellar data obtained in all-sky, homoge- These data are the result of compiling an extensive neous, astrometric and photometric systems, having bibliographyand contain, among other things, the in view the subsequent use of these parameters for distances and ages of 420 open clusters. The dis- analyzing individual open clusters and the Galactic tances, color excesses, and ages of 340 open clusters, disk as a system of such clusters. We address some recentlyexpanded to include a total of 425 open clus- of these tasks in this paper. Section 2 describes the ters, can be found in [6–8]. sources of the observational data, and Section 3, The studyof the nearest clusters was among the the open-cluster sample used and the technique em- primarytasks of the Hipparcos [1] mission. Robichon ployed. Section 4 summarizes the underlying princi- et al. [10], Baumgardt et al. [11], and Loktin and ples of identifying probable cluster members based on Beshenov [12] determined the trigonometric paral- proper-motion data and BV spaceborne photometry laxes of 50, 205, and 45 open clusters based on pre- and the results of applying these principles. Sec- selected lists of 550, 715, and 638 members, respec- tion 5 analyzes the structural parameters of the open tively. Robichon et al. [10] and Baumgardt et al. [11] clusters obtained, and Section 6, specific features obtained the proper motions of the open clusters in of the distributions of stars of various masses in- the Hipparcos reference frame. Dias et al. [13, 14] side these open clusters. Section 7 reports the open- determined these parameters for 205 clusters us- cluster proper motions in the Hipparcos frame deter- ing data for about 5900 Tycho-2 [2] stars identi- mined in this paper. fied as cluster members using the Sanders statistical method. As a result, open-cluster proper motions 2. CATALOGS OF STELLAR DATA were first determined directlyin an inertial reference frame. Glushkova et al. [15] determined the proper Our principal tools for determining the parameters motions of 331 open clusters using 2980 photomet- of the open clusters were the ASCC and USNO star ricallyidenti fied cluster members selected from the catalogs. “Catalog of About Four Million Stars” [16] created The ASCC, which is based on large, modern, at the Sternberg Astronomical Institute and then re- high-precision catalogs and covers the magnitude duced the resulting proper motions to the Hipparcos interval V =12− 14m, is most complete in terms reference frame. In all these studies, the authors have of its list of stars and the composition of the stel- adopted previouslypublished structural parameters lar data reduced to homogeneous, all-sky, standard for the open clusters (the central coordinates and systems. The input data for the ASCC include both sizes) without redetermining them. catalogs based directlyon data from the Hippar- The size of a cluster is one of its main parame- cos and Tycho [1] projects (the files hip_main.dat ters, which is related to the dynamical states of the and h_dm_com.dat, with multiple-star components, cluster and the Galactic disk and plays a key role in and tyc_main.dat), and catalogs based on earlier da- determining whether particular stars belong to the ta from the Astrographic Catalog: ACT RC [19], cluster or not. Often, however, onlythe size of the TRC [20], and Tycho-2 [2]. The list of Hipparcos central region of the cluster has been considered, stars is incomplete, especiallyin crowded fields, where although manystudies (see Kholopov’s [17] mono- some of the clusters are located. The same draw- graph for references) have found that clusters pos- back is also characteristic of catalogs based on much sess extended coronas whose identification is made more complete observations performed in the Tycho difficult bysuperposed rich stellar fields. It is the project, due to the specifics of the TUC 1/2 input size of the core—even if the cluster was known to star lists, which were based primarilyon the Guide possess a corona—that Lynga˚ tended to include in Star Catalog [21]. Therefore, an appreciable fraction the LCOC, and manyauthors mistook these sizes for of stars in such areas of the skyis provided bythe those of the entire clusters. Researchers have rarely ground-based CMC11 [22] and PPM [23–25] cata- redetermined cluster sizes, and veryfew of them have logs, which we also used in compiling the ASCC. done this on a massive scale using a homogeneous The ASCC includes 2 501 968 stars with J2000 technique. Danilov and Seleznev [18] (hereafter DS) equatorial coordinates for the observing epoch of created the most homogeneous catalog of structural 1991.25 and compiled proper motions in the Hip- characteristics of 103 Northern clusters based on star parcos system. The trigonometric parallaxes are counts down to B 16m performed on wide-angle computed as weighted averages over the catalogs [1]. plates obtained with a single instrument—the SBG With few exceptions, the magnitudes are adopted camera of the Kourovka Observatory. from the Tycho catalogs and reduced to the Johnson Our aim is to determine the structural and kine- B,V system using a photometric reference system matic parameters of all open clusters whose mem- defined by20 000 stars never suspected to be variable. bers can be found in large modern catalogs of high- The final catalog gives the standard errors for all these
ASTRONOMY REPORTS Vol. 47 No. 4 2003 SUBSYSTEM OF OPEN CLUSTERS 265 quantities. The MK and/or HD spectral types; mul- and µy ≡ µδ), list of cluster members, distance d,and tiplicityand variability flags; and the Hipparcos, HD, color excess E(B − V ). We carried out several iter- and DM numbers were taken from the corresponding ations for each cluster and in each performed a joint source catalogs. Its completeness and the accuracyof analysis of the corresponding sky chart, vector dia- its data makes the ASCC a unique tool for selecting gram of the proper motions, magnitude dependence cluster members and analyzing star clusters. of the proper-motion components, V vs. (B − V ) The USNO catalog is currentlythe deepest all- color–magnitude diagrams, the location of the zero- skycatalog, containing J2000 equatorial coordinates age main sequence (ZAMS) in the color–magnitude and magnitudes in systems close to the B and R diagram, and the radial distribution F (r) of the stellar bands for a total of 526.3 million stars down to a density. limiting magnitude of B 22m. These data are based on Schmidt plate scans taken from 1950 to the 1980s As a zeroth approximation, we used the area with at the Palomar and European Southern Observato- the LCOC coordinates and a diameter equal to three ries. The USNO positions were computed using the times the LCOC value (but not greater than 1◦). We ACT RC [19] catalog, whose positions and proper then searched this area for cluster members using motions are, in turn, tied to the Hipparcos system. the technique described in Section 4. We adopted the The magnitudes of bright stars are tied to the Tycho-1 point of maximum surface densityof cluster members [1] data, while those of faint stars were calibrated us- as the cluster center, which we subsequentlyused to ing CCD observations. The use of the USNO catalog compute the function F (r) to determine the radius for analyzing star clusters is hindered by the fact that of the region of skyonto which the members are it contains no objects in the vicinityof bright stars and projected. We then used the resulting spatial param- in crowded areas of the sky. However, for sufficiently eters as input data for the next approximation. As a compact open clusters, the USNO catalog can be rule, two approximations were sufficient to identify used to determine the structural cluster parameters the cluster members and determine the structural and unbiased bydetector-size limitations. kinematic parameters; however, some clusters with In addition, we made extensive use of the “Mil- sparse structure or erroneous initial central coordi- lennium Star Atlas” [26], which contains more than nates required three and more approximations. one million stars of the Tycho-1 catalog [1] and gives some idea about how clusters appear against field Asarule,thedifferences between the refined open- stars and about the appearance of specific features cluster coordinates and those listed in [5] are small due to the presence and distribution of neighboring and can be explained bythe consideration of di fferent clusters, field stars, interstellar gas, and dustymate- distributions of stars with different magnitudes. These ◦ rial. discrepancies exceed 0.3 for onlyten clusters, due to errors in the open-cluster coordinates in the old catalogs or to problems with identifying the centers of 3. LIST OF OPEN CLUSTERS large open clusters, such as Mel 20 (α Per), Mel 111 AND TECHNIQUE EMPLOYED (Coma), and NGC 2632 (Praesepe). We compiled a preliminarylist of about 500 open clusters byidentifying all the LCOC clusters in the We adopted the functional form of the ZAMS ASCC and the “Millennium Star Atlas” [26], pro- specified bySchmidt-Kaler [27] and determined its position on the color–magnitude diagram based on vided that at least ten ASCC stars were located within − an area of skywith a diameter three times the value the distance d, color excess E(B V ), and extinction − − given in the LCOC. We did not consider the Hyades AV =3.1 E(B V ).Weadoptedthed and E(B and similar open clusters, which span verylarge ar- V ) values for 231, 17, 22, and 13 open clusters from eas, in order to preserve the homogeneityof the meth- Loktin et al. [7, 8], with allowance for the correc- ods used to determine the cluster parameters. We tions therein; Robichon et al. [10]; the LCOC; and discarded some clusters in the process of our work, the WEBDA open cluster database [28], in order of and our final list contains a total of 401 open clusters, decreasing priority. For some clusters, primarily those each represented in the ASCC byat least five probable with variable extinction, we corrected the E(B − V ) members. values in accordance with the shape of the color– We used a multifaceted approach to determine magnitude diagram. For 118 clusters with unknown and/or refine the structural and kinematic parameters d and E(B − V ) values, we determined these pa- of the open clusters—the central equatorial coordi- rameters byZAMS fitting, taking into account the nates (α, δ)J2000,radiir1 of the core and r2 of the appearance of the cluster color–magnitude diagram corona (which corresponds to the radius of the entire and the absolute magnitudes of stars with known cluster), proper-motion components (µx ≡ µα cos δ spectral types.
ASTRONOMY REPORTS Vol. 47 No. 4 2003 266 KHARCHENKO et al.
δ 2000 V –65.5° 4 (‡) (b) 5 6 –66.0 7 8 9 –66.5 10 11 12 –67.0 13 0 0.5 1.0 1.5 2.0 B – V
h α 14.9 14.8 14.7 2000 µ x 20 (d) µ y 10
20 0 (c) –10 10 –20 µ y 20 0 (e) 10
–10 0
–10 –20 –20 –20 –10 0 10 20 5678910111213 µ x V
Fig. 1. The cluster С1439-661 (vdB-Hagen 164) according to the ASCC data. The circles on the equatorial-coordinates map (a) show the most probable cluster members and the triangles all remaining stars. The size of each symbol corresponds to the magnitude of the corresponding star. The pluses and crosses indicate the positions of the cluster center according to the LCOC and our analysis, respectively. The dashed circles with radii of 0◦.35 and 0◦.95 bound the core and corona regions, respectively. The large and small circles in plots (b)–(e) indicate probable cluster members located in the core and corona, respectively, while the points show all remaining stars. The proper motions are in units of 0.001 /yr. The position of the ZAMS on the color–magnitude diagram (b) corresponds to d = 541 pc and E(B − V )=0m. 15. The dashed lines in (d) and (e) indicate the average motion of the cluster: (µx,µy)=(−7.1 ± 0.4, −10.9 ± 0.4) × 0.001 /yr.
Figure 1 shows the result of our iterative procedure 4. IDENTIFICATION OF CLUSTER for identifying the members of the cluster С1439- MEMBERS 661 (vdB-Hagen 164) and determining its previously Unbiased determinations of cluster parameters are unknown distance and E(B − V ). possible onlyif the cluster members are correctly
ASTRONOMY REPORTS Vol. 47 No. 4 2003 SUBSYSTEM OF OPEN CLUSTERS 267
log F(r) 2 Stock 2
1
0
–1 2 NGC 2632
1
0
–1 2 NGC 7092
1
0
–1 0 1°2°3°4°5° r
Fig. 2. Distribution of the surface density F (r) of ASCC stars with V ≤ 11m. 5 as a function of the angular distance from the cluster center in areas of skywith three selected open clusters. The vertical axis has a logarithmic scale and F (r) is in units of square degree. The thin, medium, and bold lines show F (r) computed for all stars, possible cluster members (P ≥ 1%), and most probable cluster members (P ≥ 61%), respectively. The dashed lines show the radii of the cluster cores (r1)and coronas (r2).
identified. The most objective methods for identify- and µx,y are the components of the mean proper ing members are based on kinematic parameters, motion of the cluster. The factor Ci is equal to 1 or first and foremost, proper motions. However, because 0 depending on whether the star satisfies the photo- the proper motions of clusters and of the disk stars metric criterion for cluster membership. We adopted against which the open clusters are observed are usu- allycomparable to each other and to their standard errors, we used photometric data as a supplementary Table 1. Statistical data for ASCC stars projected onto criterion. 401 clusters We used the following iterative procedure to iden- tifycluster members in each of the approximations Cluster areas described above. For all stars located within a circular Sample region of a given radius, we computed the probability 0 ...r r ...r 0 ...r P of cluster membership for the ith star: 1 1 2 2 0.5 µi − µ 2 All stars 9497 53372 62869 P i = Ci × exp − x x (1) 2 εi µx P ≥ 61% 4513 7910 12423 i − 2 µy µy 14% ≤ P<61% 1861 8464 10325 + × 100%, εi µy 1% ≤ P<14% 586 4775 5361 i i where µx,y and εµ are the components of the proper x,y P<1% 6960 21149 28109 motion of the ith star and their standard errors (1σ)
ASTRONOMY REPORTS Vol. 47 No. 4 2003 268 KHARCHENKO et al.
0.5 0 (‡) (c) 0 –0.5 (DS) (ASCC) 2 1 r r –0.5 log log –1.0
–1.0 –1.5 –1.0 –0.5 0 log r1(USNO)
(b) (d) 0.5
0 (ASCC) 2 r –0.5 log
–1.0 –1.5 –1.0 –0.5 0 0.5 1.0 –1.0 –0.5 0 log r2(USNO) log r2(DS)
Fig. 3. Comparison of angular radii of open-cluster coronas r2 (a–c) and cores r1 (d) determined from star counts using the ASCC, USNO catalog, and DS. The scales are logarithmic and radii are in degrees. The straight lines correspond to loci of equal radii.
Ci =1if the star’s position in the color–magnitude P ≥ 61%.Glushkovaet al. [15] identified, on aver- i age, nine stars from catalog [16] per cluster. Thus, our diagram is no more than 1σ of the color index ε(B−V ) to the left of the ZAMS. multifaceted approach has enabled us to most com- pletelyidentifyopen-cluster members using available i Stars with C =1and proper motions deviating high-precision stellar data. from the cluster proper motion byless than one, We did not applythe identi fication procedure to two, or three σ have cluster membership probabilities the USNO stars, due to the lack of proper motions, greater than 61%, 14%,and1%, and we therefore uncertaintyof the photometric system,and the large considered them to be most probable, probable, and errors of the R magnitudes, whose systematic errors possible cluster members, respectively. We computed were 1m and standard errors were 0m. 5, stronglydis- ≥ the components µx,y based on stars with P 61% torting the B vs. (B − R) color–magnitude diagram. within the core, then computed refined probabilities in accordance with (1). We stopped this iterative pro- 5. STRUCTURAL PARAMETERS cedure when two consecutive steps yielded identical OF THE CLUSTERS lists of cluster members. The measured size of a cluster depends on many Table 1 lists the total numbers of ASCC stars factors: the method and qualityof member identi fica- projected onto the core and corona regions of 401 tion, the limiting magnitude and densityof the data clusters, bounded bycircles with radii r1 and r2.We used, and the presence and distribution of absorbing identified a total of 12 500 most probable members matter. with P ≥ 61%, which, on average, corresponds to 31 We determined the structural parameters of the stars per cluster. Robichon et al. [10] and Baumgardt open clusters from star counts using the ASCC and et al. [11] identified 11 and 3.5 stars per cluster, re- USNO catalog, assuming that the distribution of the spectively. However, they used the obviously incom- cluster members is centrallysymmetricand consists plete Hipparcos [1] catalog. Dias et al. [13, 14] con- of two components with radii r1 and r2—the core sidered, on average, 28.7 Tycho-2 [2] stars per clus- and corona. These radii correspond to the cluster- ter with membership probabilities greater than 51%; centric distances where the absolute value of the stel- this corresponds to 25 such stars for the criterion lar surface densitygradient decreases abruptly( r1)
ASTRONOMY REPORTS Vol. 47 No. 4 2003 SUBSYSTEM OF OPEN CLUSTERS 269
log r(LCOC) 0.5
0
–0.5
–1.0
–1.5
–2.0 –1.0 –0.5 0 0.5 –1.5 –1.0 –0.5 0 0.5 –1.0 –0.5 0 log r2(ASCC) log r2(USNO) log r2(DS)
Fig. 4. Comparison of the corona angular radii r2 determined from star counts using the ASCC, USNO catalog, and DS with the values given byLyng a˚ [5] in the LCOC. The scales are logarithmic and the radii are in degrees. The straight lines show loci of equal radii.
and ceases to change (r2). Our method is less formal We tried to overcome these problems in two ways. than determining the structural parameters in the First, in the region of each cluster, we computed context of a specific model (e.g., a King model [29]); F (r) for stars brighter than the limiting completeness however, it can be applied to poorer clusters, thereby magnitude Bc of the catalog, which we determined by considerablyexpanding the sample studied. This is analyzing the integrated magnitude function N(B). whywe preferred this method for our analysis. For clusters with dense cores, we estimated the sizes We constructed the functions F (r) corresponding of the possible drop in the surface distribution of stars to the radial densitydistributions of stars brighter and analyzed the surface distribution taking into ac- than V =11m. 5, where the catalog reaches its com- count the presence or absence of such a drop. Second, pleteness limit, out to angular distances from the we constructed F (r) both for the entire ring and sep- cluster center of 10◦ in annular zones with steps of aratelyfor the four quadrants. This enabled us to take 0◦.05. We used three stellar samples: (i) all stars, and into account possible errors in the coordinates of the cluster members with probabilities (ii) P ≥ 1% and cluster center or asymmetry of the stellar distribution (iii) P ≥ 61%. Figure 2 illustrates the degree of vari- due to the origins above. We determined the corona ation of F (r) with radial distance for various groups sizes for 363 and the core radii for 140 open clusters of stars using three clusters (Stock 2, NGC 2632 by carrying out a joint analysis of all five F (r) func- (Praesepe), and NGC 7092) as examples. tions, which we computed over the distance interval ≤ In addition to the ASCC, we used the USNO r ρ in annular zones with steps of 0.5. catalog when possible. We determined the structural Figure 3 compares the angular radii of cluster parameters of the selected open clusters based on the cores and coronas determined directlyfrom star USNO catalog data, using the central coordinates counts in this paper and byDS. It is clear that the ASCC ◦ and corona sizes determined from the ASCC data. radii of small cores (r1 < 0.1) are overestimated To this end, we selected from the USNO catalog all due to the insufficient angular resolution of the star ASCC stars within a radius ρ =3r2 . The USNO catalog cannot be used to determine the sizes of some nearby open clusters with large angular diameters and low Table 2. Characteristic linear radii (in pc) of the cluster surface densities for their member stars, which do not cores and coronas in the Galaxy substantiallyexceed the overall stellar densityin a min min max max crowded field. Certain problems also arise in regions Source R1 R2 R1 R2 R1 R2 with high stellar surface densities located in the vicin- ACSS 0.6 1.5 7.1 19.2 1.7 5.0 ityof bright stars and nebulae where there are no stars in the USNO catalog and in some areas where the USNO 0.5 0.9 4.8 15.9 1.5 3.5 limiting magnitude of the USNO catalog is several DS 0.5 3.0 4.5 17.2 2.7 7.6 magnitudes brighter than the average value for the catalog. LCOC 0.2 – 10.0 – 2.0 –
ASTRONOMY REPORTS Vol. 47 No. 4 2003 270 KHARCHENKO et al.
2.0
1.5 , pc 2
R 1.0
log 0.5
0 ASCC USNO DS
1.0
, pc 0.5 1 R 0 log –0.5 ASCC USNO LCOC
012345012345012345 d, pc
Fig. 5. Linear sizes of the cluster cores (R1) and coronas (R2) determined from star counts using the ASCC, USNO catalog, and DS and adopted from the LCOC as a function of the distances d of the clusters.
ASCC counts using the ASCC. The corona sizes r2 catalog, and DS, respectively. To prevent further se- are larger than those determined from the USNO lection effects, we used onlythose LCOC and DS ASCC ± × USNO clusters that are included in our list. data (r2 =(1.4 0.1) r2 ) throughout the entire radius interval, despite the significantlydeeper ThefactthatourclusterlistoverlapsthatofDS limiting magnitude of the USNO star counts. The and that there is a satisfactoryagreement of results excess over the results of DS is not statistically makes it possible to establish a relation between the ASCC ± × DS significant, r2 =(1.1 0.1) r2 ,andisdue “empirical” sizes inferred from star counts and the to the fact that identification of cluster members sizes obtained byDS by fitting the results of star is not carried out in either the USNO catalog or counts to model King curves. We used the core radius DS, and the outer regions of the clusters were lost DS rc computed byDS using the King formula as r1 against the rich surrounding stellar fields. Figure 2, (we transformed these values into angular units using which shows the F (r) curves based on the ASCC distances taken from DS). data computed for all stars, stars with P ≥ 1%,and stars with P ≥ 61%, demonstrates how the degree The relation between the sizes of the cluster of selection affects the form of F (r) and therebythe coronas and cores is more or less the same whether theyare determined using the ASCC or USNO cluster sizes. Thus, the selection of cluster members ASCC ± × ASCC USNO in the ASCC enabled us to derive the most reliable catalog: r2 =(2.4 0.2) r1 and r2 = ± × USNO cluster sizes, which exceed those determined using (2.6 0.1) r1 . Because the King radius refers the deepest catalogs but without selection of cluster DS to the most central region of the cluster, whereas r2 members byseveral tens of percent. approximatelycorresponds to our de finition of the Figure 4 compares our radii with the extensive cluster radius, the correlation between the data of DS compilation of cluster sizes given in the LCOC. For DS ± × DS hasadifferent form: r2 =(4.5 0.4) r1 .Acom- comparison, we also give data adopted from the ho- parison of these relations yields a relation between mogeneous sample of DS. Figure 4 shows how biased the King parameter and the empirical core radius (mostlyunderestimated) are the cluster sizes used by rDS =(0.5 ± 0.1) × rASCC =(0.6 ± 0.1) × rUSNO. manyauthors, who adopt them from the LCOC and 1 1 1 forget that Lynga˚ [5] himself aimed to include in his DS also computed the tidal radii rt of open clus- catalog onlythe sizes of the cluster cores. ters, which we transformed into angular units in the DS These data are, on average, underestimated by same wayas for r1 andcomparedwiththeempirical factors of 3.0 ± 0.1, 2.2 ± 0.1,and2.9 ± 0.2 com- cluster radii r2. As expected, the tidal radii are upper pared to the sizes of the cluster coronas (and, con- limits for r2 and, on average, exceed the latter bya sequently, the sizes of the clusters themselves) de- factor of 1.5 ± 0.1, 2.1 ± 0.1,and1.6 ± 0.1 for the termined from star counts using the ASCC, USNO ASCC, USNO catalog, and DS, respectively.
ASTRONOMY REPORTS Vol. 47 No. 4 2003 SUBSYSTEM OF OPEN CLUSTERS 271
MV –6 –4 logA 7.5 logA 7.5 –2 0 2 4 –6 –4 7.5 < logA 8.0 7.5 < logA 8.0 –2 0 2 4 –6 –4 8.0 < logA 8.5 8.0 < logA 8.5 –2 0 2 4 –6 –4 logA > 8.5 logA > 8.5 –2 0 2 4 012345012345 R/R1
Fig. 6. Absolute magnitude MV of the brightest cluster member (P ≥ 61%) within a ring of a given radius as a function of the relative radius R/R1 for clusters in four age intervals. The left and right plots show stars that are still on the main sequence and those that have left it.
As pointed out above, distance estimates d are For the USNO sample, on the contrary, it becomes available for all the open clusters, enabling us to de- systematically more difficult to identifynearbyopen termine the linear radii of the cores (R1) and coronas clusters against the rich background stellar field. The (R2). Figure 5 shows the dependence of these radii trends in the upper boundaries of the linear radii with on d. These distributions are characterized bystrong increasing d are much weaker for the other samples. selection effects near the lower boundaryfor the radii, We computed the lower and upper boundaries due to the fact that, as the distance increases, small Rmin and Rmax for the linear sizes of the structural open clusters become increasinglyinconspicuous components byaveraging over the ten clusters with against the background stellar field. This selection the smallest and largest radii in each sample (ex- effect is evident for both R and R . All the samples 1 2 h χ RASCC ≈ 80 suffer from this effect, even the homogeneous DS cept and Per, which have 2 pc; they sample, whose objects are located at characteristic correspond to the two uppermost points in the top distances d =1–2kpc. left graph in Fig. 5). The mean linear radii R are determined from clusters located at distances smaller The upper boundaryof the distribution of linear than 1 kpc, where selection effects are comparatively radii is less subject to selection effects. The system- ASCC weak. These data are summarized in Table 2, and, atic effect in R1 maybe due to the overestimation although theysu ffer from various selection effects, ASCC of small values of r1 discussed above, i.e., to theygive some idea of the characteristic sizes of problems in resolving the cores of distant clusters. Galactic open clusters.
ASTRONOMY REPORTS Vol. 47 No. 4 2003 272 KHARCHENKO et al.
µ x(ASCC) 20 (‡) (b) (c) (d)
0
–20
–40 –40 –20 0 20 –40 –20 0 20 –40 –20 0 20 –40 –20 0 20 µ µ y(ASCC) x 20 (‡) (b) (c) (d)
0
–20
–40 –40 –20 0 20 –40 –20 0 20 –40 –20 0 20 –40 –20 0 20 µ y
Fig. 7. Comparison of cluster proper motions (in 0.001 /yr) derived here with those of (a) Robichon et al. [10], (b) Baumgardt et al. [11], (c) Dias et al. [13, 14], and (d) Glushkova et al. [15]. The straight lines are the loci of equal mu components in the two catalogs compared. The circles show the data for clusters included in all five catalogs.
6. SPATIAL SEGREGATION main-sequence stars and of the stars that have left the OF STARS IN CLUSTERS main sequence. We can see that the area occupied by Our data can be used to analyze differences be- stars expands substantiallywith decreasing bright- tween the spatial distributions of stars of different ness: bright stars are observed exclusivelyin central masses in clusters, i.e., the mass segregation, for cluster regions, whereas fainter, i.e., less massive, which we have used the absolute magnitude M as stars are distributed over the entire cluster volume. V The observed distributions were constructed using an analog. We selected two of the brightest most manyclusters with reliable and mutuallyindependent probable cluster members (P ≥ 61%) in each annu- data. This suggests that the observed trends are lar zone around the cluster centers used to compute statisticallysigni ficant. We eliminated the potential F (r): one located in the main-sequence band in the danger of systematic effects arising as a result of the color–magnitude diagram, which deviates to the left mutual normalization of the cluster sizes bycompil- ε − and right of the ZAMS byno more than (B V ) and ing a homogeneous sample of clusters and choosing 0m. 2+ε − , and the second located to the right (B V ) thecoreradiusR1 as the normalizing coefficient, of this band, in the red part of the color–magnitude which is less subject to selection effects than R2.The diagram corresponding to late evolutionarystages. realityof the observed e ffect is also supported bythe We determined the absolute magnitudes MV of these absence of segregation in the group of old clusters. stars from their known V and d values. To ensure The data in Fig. 6 can also be interpreted in a the homogeneityof the sample, we selected clusters different way: stars from the entire mass interval, from ≈ with d<1 kpc and mean ratios R2/R1 4.5 whose lowest to highest, are observed at the center of the members span a magnitude interval of no less than cluster, whereas onlythose with masses not exceed- m 4 . We subdivided the sample into four groups in age ing a certain threshold can be seen at the periphery. ≤ ≤ ≤ A log A 7.5, 7.5 < log A 8.0, 8.0 < log A 8.5, It is this effect that we call mass segregation. The log A>8.5 consisting of 10, 27, 14, and 9 clusters, younger the cluster, the stronger the segregation. A respectively. comparison of results for clusters of different ages Figure 6 shows the dependence of MV for these shows that the form of the initial distribution of stars brightest stars on the relative distance R/R1 for with radial distance does not change with time, and a number of selected open clusters. The left-hand the observed flattening of this distribution is due to and right-hand graphs show the distributions of the the burning out of massive stars at the cluster center.
ASTRONOMY REPORTS Vol. 47 No. 4 2003 SUBSYSTEM OF OPEN CLUSTERS 273
Table 3. Parameters of 401 open clusters
Parameter Number of Identifier Unit of Comment clusters measure IAU number 401 LCOC Name 401 LCOC J2000 right ascension 401 α h ASCC J2000 declination 401 δ deg ASCC ASCC Core radius (ASCC) 401 r1 deg ASCC ASCC Corona radius (ASCC) 401 r2 deg ASCC 1 USNO Core radius (USNO) 140 r1 deg USNO 1 USNO Corona radius (USNO) 363 r2 deg USNO 1 Completeness B mag. limit of the catalog 363 Bc mag. USNO Number of most probable members 401 N(P ≥ 61%) ASCC
Number of stars used to determine the 401 Nµ ASCC mean proper motion