A&A 370, 87-99 (2001)
DOI: 10.1051/0004-6361:20010230
J. Sanner - M. Geffert
Sternwarte der Universität Bonn, Auf dem Hügel 71, 53121 Bonn, Germany
Received 14 September 2000 / Accepted 7 February 2001
Abstract
We studied the fields of nine nearby open star clusters based on
the Tycho-2 catalogue. We determined membership probabilities for the stars
in the cluster fields from the stellar proper motions and used the Tycho-2
photometry to compute the initial mass function (IMF) of the clusters from
the main sequence turn-off point down to approx.
.
We found IMF
slopes ranging from
down to
(when the Salpeter
1955 value would be
). We also studied the
membership of individual stars of special astrophysical interest. In some
cases previous results had to be revised. As a by-product, we
investigated some general properties of the Tycho-2 catalogue; we
confirmed that the Tycho-2 proper motions show only marginal deviations
from the Hipparcos catalogue. On the other hand, in some regions the
completeness of the catalogue seems to decrease at magnitudes brighter than
claimed by Høg et al. (2000).
Key words: open clusters and associations: general - catalogs - astrometry - stars: kinematics - Hertzsprung-Russell (HR) and C-M diagrams - stars: luminosity function, mass function
The Tycho-2 catalogue (Høg et al. 2000) provides positions and stellar proper motions in conjunction with two colour photometric data. These were computed from the original Tycho (star mapper) data obtained during the Hipparcos mission (ESA 1997) and ground based first epoch data, including the Astrographic Catalogue (Wood 1970). This work was triggered by the findings of Makarov et al. (2000) who recently found that the precision of the Tycho-2 proper motions are of the same order as the Hipparcos data. With its significantly larger number of stars, the Tycho-2 catalogue therefore is a suitable source of data for studies of stellar kinematics. In earlier investigations Hipparcos data were used for a search for open clusters (Platais et al. 1998) and checks of open cluster candidates (Baumgardt 1998).
In this investigation we used fields of nine nearby open star clusters. We separated the field stars from the cluster members with the help of the proper motions provided by Tycho-2 and determined colour magnitude diagrams (CMDs) from its photometry. From the CMDs we derived the age, distance modulus, and reddening by fitting isochrones and compared the results with, e.g., distances derived from Hipparcos parallaxes (van Leeuwen 1999a, hereafter referred to as vL99a). In a final step, we computed the initial mass function (IMF) of the clusters with a maximum likelihood method. The results of this part of the study can give us clues about the completeness of the Tycho-2 data.
Especially for open clusters, where the contamination of the data with field stars poses a major challenge, knowledge of membership is essential to derive the IMF. The most precise method to determine membership is by studies of stellar kinematics, as one can assume that all cluster stars move in the same way, whereas the distribution of the motions for the field stars is much wider and differently centred. From this point of view, a catalogue like Tycho-2, including photometry and proper motions, provides an ideal opportunity for an IMF study.
The shape of the IMF is an important parameter to understand the fragmentation of molecular clouds and therefore the formation and development of stellar systems. Besides studies of the Solar neighbourhood (Salpeter 1955; Tsujimoto et al. 1997), work on star clusters plays an important role (Scalo 1986,1998), as age, metallicity, and distance of all cluster stars can be assumed to be the same.
Restricted to certain mass intervals, most studies agree that the IMF can
be described by a power law in the form
![]() |
(1) |
Up to now, thorough IMF studies are available only for a few of the nearby star clusters, e.g. the Pleiades (van Leeuwen 1980; Meusinger et al. 1996; Hambly et al. 1999) or Praesepe (Williams et al. 1995; Pinfield et al. 1999). Some of the clusters were also studied by Tarrab (1982), but her work was based on an inhomogeneous data sample with sometimes weak information on membership. Discrepancies of her results with later work have been found (see, e.g., Sanner et al. 2000 for NGC1960 or Credner 1996 for NGC2323), so that revisiting the clusters studied by Tarrab (1982) is promising.
Nowadays, although CCD cameras provide a good tool for IMF studies down to low
masses (
), the limited field size of typical imagers of the
order of arc minutes makes comprehensive studies a difficult task: Nearby
clusters have a large angular diameter in the sky, so that mosaics of lots of
different fields would have to be observed to obtain a general overview of the
stellar content of a cluster. For the determination of a cluster IMF a
snapshot of one or only few fields is not enough, since spatial
inhomogeneities (mass segregation; see, e.g., Raboud & Mermilliod
1998) are likely to falsify the results. With a uniform photometric
and astrometric all sky survey, this goal can be achieved at least down to a
reasonable mass limit.
Section 2 contains a description of the data and methods used for this work as well as some general results, followed in Sect. 3 by the results sorted by cluster. The final section contains a summary and conclusions.
Since the Tycho-2 data are claimed to be complete down to
(see Høg et al. 2000 and the discussion below), we were
restricted to a small magnitude (and hence mass) interval. To reach stars with
masses as low as possible we chose
to study nearby (
or
)
open star clusters only.
Robichon et al. (1999) published a list of astrometric data of
nearby open clusters based on the Hipparcos catalogue. From this information
the nine targets listed in Table 1 were chosen. A main
criterion for the selection was a clear separation between the field and
cluster proper motion centres to make the membership determination an easy
task and to avoid the presence of numerous field stars with the same proper
motion as the cluster. A maximum of 10 to 15 stars can be expected to belong
to the field but to have (within the errors) the same proper motion as the
corresponding cluster. The Hipparcos parallaxes of seven of them were recently
studied by vL99a to derive the distances of the clusters. Similar computations
can also be found in Robichon et al. (1999). Due to the very high
quality of the Hipparcos data and the reliability of the results, we preferred
these values for the distance moduli and reddenings of these clusters compared
with those of previous studies. The results of vL99a and Robichon et al. (1999) agree well within the errors for our clusters, except for
Blanco1 with a deviation of
.
We
will discuss this point later (Sect. 3.5).
cluster |
![]() |
![]() |
[h m] | [![]() ![]() |
|
Blanco1 | 00 04.3 | -29 56 |
Stock2 | 02 15.0 | +59 16 |
![]() |
03 22.0 | +48 37 |
Pleiades | 03 47.0 | +24 07 |
NGC2451A | 07 45.4 | -37 58 |
IC2391 | 08 40.2 | -53 04 |
Praesepe | 08 40.4 | +19 40 |
IC2602 | 10 43.2 | -64 24 |
NGC7092 | 21 32.2 | +48 26 |
cluster | field | ||||||||
cluster |
![]() |
![]() |
![]() |
![]() |
![]() |
a | b | ![]() |
r |
[mas yr-1] | [mas yr-1] | [mas yr-1] | [mas yr-1] | [mas yr-1] | [mas yr-1] | [mas yr-1] | [![]() |
[![]() |
|
Blanco1 | +19.27 | +2.81 | 2.02 | +10.98 | -7.71 | 16.48 | 11.93 | +24.5 | 11.1 |
Stock2 | +17.15 | -12.94 | 4.61 | +0.36 | -1.58 | 6.76 | 6.40 | -31.8 | 14.0 |
![]() |
+22.63 | -25.74 | 1.52 | -0.74 | -3.53 | 13.11 | 8.14 | -37.0 | 3.7 |
Pleiades | +19.96 | -45.09 | 1.78 | +5.43 | -11.42 | 17.62 | 14.31 | -53.1 | 7.6 |
NGC2451A | -21.96 | +14.59 | 1.17 | -5.40 | +4.59 | 11.21 | 7.33 | -57.8 | 5.4 |
IC2391 | -24.64 | +23.25 | 1.13 | -6.82 | +6.29 | 12.34 | 7.28 | -48.3 | 6.0 |
Praesepe | -35.88 | -13.18 | 1.93 | -6.40 | -9.20 | 11.47 | 11.32 | +6.5 | 10.6 |
IC2602 | -17.63 | +10.30 | 1.50 | -8.55 | +2.43 | 9.62 | 5.68 | -26.6 | 4.5 |
NGC7092 | -8.34 | -19.16 | 1.19 | -2.00 | -2.26 | 10.93 | 6.52 | +41.1 | 2.7 |
The Tycho project made use of its own photometric system, referred to as BTand VT. These two filters are close to Johnson B and V, respectively,
and ESA (1997) provides transformation equations to the Johnson system:
After extracting sufficiently large areas around the centres of the targets - taking at least twice the diameter given in the Lyngå (1987) catalogue to avoid a biasing due to mass segregation - we derived the membership probabilities P of the individual stars from their proper motions. We fitted a combination of a sharp (for the members) and a wider (for the field stars) two-dimensional Gaussian distribution to the vector point plot diagram. The optimal parameters (location of the field and cluster proper motion centres in the diagram and the standard deviations of the two Gaussians) were determined with a maximum likelihood procedure as described in Sanders (1971).
Due to the large fields, the vector point plot diagrams were dominated by the
field stars. For some of the clusters, this caused the software to
misinterpret the field stars as the cluster members and vice versa. Since the
stars determined as members are located within a circle or an ellipse, and the
non-members consist of all the rest, the correct membership probabilities
could not be computed as the complement of the determined ones: Any star at a
sufficient distance from the field centre would obtain a high probability with
this method, even if that star is located at the opposite end of the
vector point plot diagram, compared with the actual cluster. This problem
could be solved by removing all stars with a very large deviation from the
cluster proper motion (>
). Since the clusters were
easily detectable in the vector point plot diagrams and the Robichon et al. (1999) data gave additional clues about their location, it can be
excluded that the region covering the target was accidentally eliminated. Also,
cross checks of the positions of the stars selected as members in the vector
point plot diagram showed that the membership determination was successful.
The vector point plot diagrams (see the diagram of the
Per cluster
derived from the Tycho-2 data in Fig. 7 as an example)
showed a standard deviation of the cluster stars in the region of 1.5-2 mas yr-1.
Since these values correspond to proper motion errors
rather than to internal motions, a good coincidence between our results and
the values given in Høg et al. (2000, their Table 2) for the
Tycho-2 accuracy was found. The only exception is the cluster Stock2 with a
width of the distribution of approx.
.
We will return
to this outlier in Sect. 3.2.
From the stars with a sufficiently high membership probability (P>0.5 to
P>0.8, depending on the object) we derived CMDs of the clusters. We fitted
isochrones based on the stellar models of Bono et al. (1997) and
provided by Cassisi (private communication) to the CMDs to derive estimations
for the distance moduli, reddenings, and ages of the clusters. Our initial
values were those of vL99a and Robichon et al. (1999) which in
most cases led to a sufficiently good fit. For some clusters, however, the
parameters had to be slightly adjusted, since the isochrones based on these
results clearly did not represent the CMDs well. In those cases, which will
individually be discussed in the corresponding sections, we used our own
results. The metallicity was in all cases assumed to be Solar, which is,
according to the Lyngå (1987) catalogue, appropriate for all
clusters. The only exception may be NGC2451A which is
metal-deficient. Values for its metallicity as given in the literature range
from Z=0.006 to Z=0.013 (Strobel 1991), which, however, in terms
of masses, leads to an error of a few per cent at most. The CMDs, with
isochrones overplotted, are displayed in Figs. 1,
3, and 5. The parameters which are listed in
Table 3 were necessary for the transformation of the VTmagnitudes to the initial stellar masses. From
plots of the
members we determined lower limits of the cluster diameters: We selected
diameters outside of which almost no more members are present, assuming that
these outliers are field stars with the same proper motion as the
corresponding cluster. These results are added to Table 3. The
"true'' sizes might be larger due to mass segregation and the limited
magnitude range of our study.
![]() |
Figure 1:
CMDs of the members of six of the eight open star
clusters from Tycho-2 data. We overplotted the best-fitting isochrones. The
CMD of Stock2 is discussed in more detail in Sect. 3.2, the CMD
of the Pleiades is plotted in Fig. 3. Note the lack of stars
fainter than
![]() |
Open with DEXTER |
cluster | (m-M)0 | EB-V | t |
![]() |
[mag] | [mag] | Myr | [
![]() |
|
Blanco1* | 6.8 | 0.03 | 50 | 105 |
Stock2 | 7.5 | ![]() |
100 | 260 |
![]() |
6.3 | 0.09 | 20 | 255 |
Pleiades* | 5.6 | 0.05 | 75 | 300 |
NGC2451A* | 6.4 | 0.00 | 20 | 140 |
IC2391* | 5.8 | 0.00 | 20 | 110 |
Praesepe* | 6.0 | 0.00 | 650 | 195 |
IC2602* | 5.8 | 0.03 | 10 | 185 |
NGC7092 | 7.6 | 0.12 | 70 | 170 |
We found that the members of Stock2 are differentially reddened, whereas the main sequences of the other objects either show a uniform reddening or are too sparsely populated for a distinct proposition. The variations of the reddening in the field of Stock2 will be discussed in detail in Sect. 3.2.
We computed the masses from the VT magnitudes with the help of the
mass-luminosity relation provided together with the isochrones. We used
6th order polynomials as an approximation. Attempts with polynomials
of lower orders failed since the low masses were not sufficiently well
reproduced. A comparison of our isochrones basing on the Bono et al. (1997) models with the Geneva isochrones (Schaller et al. 1992) revealed that the difference between the masses as functions
of the V magnitudes are less than .
We derived the IMF slopes of
the clusters from the masses of the members with a maximum likelihood
method, i.e. we maximised the likelihood function
![]() |
(4) |
A critical point of this work was
the completeness of the Tycho-2 data which would have contaminated the
results towards too high
values, i.e. too shallow slopes. Høg
et al. (2000) state in their Table 1 that the Tycho-2 data are 99% complete down to
and still 90% down to
.
However, it cannot be excluded that these figures vary for
different fields, depending on crowding or the first epoch material used to
build Tycho-2. We tried to get an idea about the completeness of the data by
the optical impression of the CMDs and by computing the IMF for different mass
ranges: We assumed that the IMF indeed has the shape of a single power
law within the range of the Tycho-2 photometry. If the lower mass limit is
below the completeness limit, the resulting IMF would be too shallow, since
stars at the low mass end of the function would be missed. When we vary the
lower mass limit up towards the completeness limit,
would approach its
"correct'' value. Above this point,
would remain on a stable
level. This phenomenon is illustrated with the Stock2 data in
Fig. 2. An alternate interpretation of this behaviour of the
slope as a function of the lower mass limit is that the shape of the IMF
deviates from a power law with constant slope, e.g. due to a gap in the main
sequence - a phenomenon which is not uncommon for open star cluster CMDs
(see, e.g., the recent work of Rachford & Canterna 2000 and
references therein). Therefore, the assumption that the IMF indeed follows a
"perfect'' power law may lead to underestimations of the completeness limit:
Read the figures given in Table 4 as upper limits (in
terms of magnitudes) for the completeness limit in the fields of the
clusters. Both interpretations recommend us to stop the determination of the
IMF at the point where the slope becomes variable: We would either compute
the IMF in a mass interval in which the data are incomplete, or we would fit a
single power law to an IMF which does not follow such a function.
![]() |
Figure 2:
IMF slope ![]() ![]() ![]() ![]() |
Open with DEXTER |
We found varying completeness limits by this investigation. For two of
the clusters, the star counts drop significantly at magnitudes as bright as
(see Table 4).
Most of the clusters are sparsely populated in the mass range under consideration. This made the choice of the correct completeness limit a challenging task. Figure 10 shows an overview of the IMFs of six of the clusters; the Pleiades IMF is sketched separately in Fig. 4. The results are discussed in the following section and summarised in Table 4.
An alternative way of determining the completeness limit would be the comparison of the Tycho-2 data with a second catalogue with fainter limiting magnitude, like the GSC 1.2 (Röser et al. 1998b) or USNO A2.0 (Monet et al. 1998): As long as no stars of the other catalogue are missing in Tycho-2, we can assume Tycho-2 to be complete. However, for this procedure it is required that the brightest stars missing in Tycho-2 are present in the second catalogue. This in general cannot be expected, since stars might be missing in both catalogues due to the presence of close neighbours, crowding, or similar effects. With this method we would only be able to obtain a limit at which Tycho-2 definitely is incomplete, which is not enough for a thorough IMF study.
Adding the supplement 1 of the Tycho-2 catalogue (see Høg et al. 2000 for details on the two supplements) did not affect the
completeness in general: Although we found from 7 (Blanco1) to 141 (Per) stars in supplement 1, only less than 10 stars per cluster (except for
Per, with 13 objects) could be found for which proper motions were
given. No more than four of them in each field also contained both BT and
VT photometry. For the stars in Tycho-2 and its supplement 1 which are
listed without proper motions, we also attempted to use the proper motions
given in other Hipparcos-based catalogues, namely ACT (Urban et al. 1998) and TRC (Høg et al. 1998). For the Pleiades, this
resulted in 52 additional stars. 21 of them can be considered as members
according to their proper motions, but only 7 objects could be added to the
IMF computation. The others either were located well below our completeness
limit (4 stars), or BT was not given (4 stars), or the data points were so
far away from the main sequence that we can assume that these data suffer from
bad photometry or the stars are non-members with the same proper motions as
the Pleiades (6 stars). An IMF computation with the 7 stars added led to a
difference in
of 0.02 compared with the IMF based on the Tycho-2
proper motions only. This difference is so much smaller than the errors that
we did not follow this procedure for the other clusters.
Comparing the proper motion centres of the clusters with the values
determined by Robichon et al. (1999) using Hipparcos data, we
found that the results coincide very well in both
and
within
.
This indicates that the
absolute proper motions of Tycho-2 correspond with the Hipparcos
catalogue and that at least in the fields of the open clusters studied there
are no systematic errors caused, e.g., by the first epoch material.
The first object under investigation is the Pleiades cluster. The Pleiades
show a proper motion centre which is located more than
away from the average field proper motion. In addition,
previous IMF studies are available, so that we can externally check the
quality of our results.
In the case of the Pleiades, the standard deviation of the member proper motions might not only reflect the inaccuracy of the Tycho-2 catalogue, but also intrinsic properties of the Pleiades, like internal velocity dispersion or perspective effects. However, since the members under consideration are located in a comparably small area (see Table 3), these effects are of minor importance.
Fitting isochrones to the CMD of the members, we found that a distance
modulus of
as given by vL99a is not the best
parameter in terms of quality of the fit: The isochrone passes right in the
middle of the main sequence, whereas it is our strategy to fit the isochrone
towards the blue side of the main sequence to account for binary stars or Be
stars, which broaden the "true'' main sequence towards the red. It seems that
an isochrone with
provides a better fit. For both
isochrones, a reddening of
and an age of the
Pleiades of 75 Myr were assumed. The CMD with both isochrones overplotted is
shown in Fig. 3. There has been a recent discussion about the
reliability of the Hipparcos parallaxes on which vL99a based his distance
determinations: Pinsonneault et al. (1998) as well as
Narayanan & Gould (1999) claim that small systematic errors in
the Hipparcos parallaxes occur, especially in the field of the Pleiades, so
that they derived a distance which corresponds to our value, while vL99a
excludes these errors as a consequence of extensive tests of the Hipparcos
data (see also van Leeuwen 1999b).
We cannot resolve this problem with our results, therefore we computed the
IMF based on both parameter sets leading to results of
and
,
with our and vL99a's distances, respectively. This
means that within the errors the IMF determination does not depend on the
uncertainties of the distance determination. The two mass functions are
sketched in Fig. 4. In Tables 3 and 4 we retain our values.
In a study of the IMF of the Pleiades, almost in the same mass range as ours,
van Leeuwen (1980) arrived at a result of (transformed to our
notation)
based on the members known at that
time. Within the errors, this is in good agreement with our values.
Meusinger et al. (1996) studied the IMF down to approx.
and stated that globally the Pleiades IMF cannot be represented by a
single power law. Therefore they did not provide a specific slope for
,
but they noted their IMF as being "somewhat steeper'' than
van Leeuwen's (1980) result in the corresponding mass interval,
so that their findings seem to support our result very well.
![]() |
Figure 3:
Colour magnitude diagram of the Pleiades
members. We overplotted the Bono et al. (1997) isochrones for
![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
Initial mass functions of the Pleiades
cluster. The upper diagram is based on a distance modulus of
![]() ![]() ![]() ![]() |
Open with DEXTER |
A region with a radius of
around the centre of the Pleiades contains
397 stars which are common to both the Hipparcos and Tycho-2
catalogues. This gave us a good opportunity to compare the proper motions of
the two catalogues. We found a good agreement with
,
,
so that the absolute reference systems of both
catalogues appear to be almost identical.
![]() |
Figure 5:
Colour magnitude diagrams of the members of
Stock2. The leftmost diagram shows all members with their apparent
magnitudes. Note the broad main sequence which inspired us to study the
differential reddening of the cluster. The middle diagram shows the subfield
of Stock2 which is accentuated in Fig. 6. This field
provided the evidence for the distance modulus of Stock2, determined
as
(m-M)0=7.5 mag. Overplotted is the isochrone which best fit
these data, the reddening is
EB-V=0.55 mag, the adopted
age t=100 Myr. The right panel is the de-reddened CMD
of Stock2 with absolute magnitudes. We overplotted the 100 Myr
isochrone. The red (![]() |
Open with DEXTER |
One of the three stars mentioned in the list of Ahumada & Lapasset
(1995) as blue stragglers of the Pleiades can be found in the
Tycho-2 main catalogue; the other two are too bright, so that no proper
motions were determined for them. The object included in Tycho-2, HD23338,
has a proper motion of
,
,
corresponding to a membership
probability of P=0.98, so that there is clear evidence that this object is a
member of the Pleiades. We marked this star in the CMD in
Fig. 3 with a cross. The star is located to the left of the
best fitting isochrone, however it is still well below the turn-off point of
the main sequence so that we do not see evidence for the blue straggler nature
of this object. For the other two objects, we can refer to the Hipparcos
proper motions instead: As stated before, the proper motions of Hipparcos and
Tycho-2 are almost identical, so that the pure Hipparcos figures can be taken
as an indication about membership without any transformation or other special
treatment, although we do not give explicit membership probabilities. In the
case of the Pleiades, the two stars not in Tycho-2 are HD23302 (HIP17499)
and HD23630 (HIP17702). Their proper motions are close to the centre of
proper motions of the Pleiades so that we can assume that these objects indeed
belong to the Pleiades. Note that since these two objects do not appear in the
Tycho-2 catalogue with complete sets of data, we did not include them in our
membership analysis and therefore they are not present in the CMD of the
Pleiades members (Fig. 3).
Slettebak (1985) provides data on Be and shell stars in open star
clusters. He reports four members of the Pleiades. The proper motions of
two of these stars, HD23302 (17Tau) and HD23630 (
Tau) are not
listed in Tycho-2, but are in the Hipparcos catalogue. Both objects are likely
members of the Pleiades according to their Hipparcos proper motions. The two
remaining stars of the Slettebak (1985) list, HD23480 (23Tau)
and HD23862 (28Tau) are cluster members with membership probabilities of
P=0.97, so that Tycho-2 supports the proposition that these objects are
members of the Pleiades, too. It is interesting that HD23302 appears in both
lists. However, it should be remarked that its category in Ahumada & Lapasset
(1995) is 3, i.e. only weak evidence for a blue straggler nature,
whereas Slettebak's (1985) discussion seems to be more sure.
Stock2 is by far the most populous cluster (down to the limiting magnitude
of Tycho-2) of our sample. There has been a controversial discussion about the
distance modulus and reddening of the cluster (see the overview in the
recently published work of Foster et al. 2000). Here we just
summarise that distance moduli from
(Piskunov
1980) up to
(Pandey et al. 1989) are
proposed; the reddenings found are of the order of
.
Our CMD of the members (Fig. 5, left diagram) shows a broad
(0.3 mag) main sequence over the entire magnitude interval,
so that we decided to study the CMDs of smaller areas to find out whether
Stock2 is differentially reddened. We found that indeed for small portions
of the data the main sequences are much narrower, and we found different
reddening values: Large portions of the field of Stock2 are reddened by
,
however, we found areas with reddenings from
up to
.
The results of this
investigation are presented in the reddening map of Stock2
(Fig. 6). Since the single fields do not cover many objects, due
to the bright limiting magnitude of Tycho-2, a much more detailed reddening
map could be achieved with a deeper photometric study of Stock2.
![]() |
Figure 6:
Reddening map of Stock2. The largest dots
represent
![]() ![]() ![]() |
Open with DEXTER |
While for most areas isochrones with a large range of distance moduli (
)
led to good fits, one field
(
,
,
marked in Fig. 6) could not be
sufficiently well fitted after applying a distance modulus of more than
.
The CMD of this field together with the
corresponding isochrone is sketched in Fig. 5, middle
diagram. On the other hand, all the other fields do correspond to a
distance modulus of
.
Therefore we adopt this
value to represent the distance modulus of Stock2. The right diagram of
Fig. 5 shows the composition of the de-reddened CMDs of the
individual sub-fields on the basis of absolute magnitudes. It is evident that
the main sequence is much more pronounced than the original CMD (left
diagram in Fig. 5), at least down to
.
With
a more detailed reddening map the CMD could still be improved.
The CMD shows three bright very blue (
,
marked in
the right diagram of Fig. 5) as well as numerous red (
)
stars that do not fit to the selected isochrone. These objects
can be explained as randomly having the same proper motions as Stock2, but
in fact being field stars. The objects are uniformly distributed over the
field in an
plot; for example, two of them can be found in
the CMD of the selected region of Stock2 mentioned before
(Fig. 5, middle diagram). Another explanation for the blue
objects would be that they are less reddened than proposed for the field in
which they are located.
The IMF was computed from 262 main sequence stars within a field of view of
,
after some 70 too red and blue stars discussed
before had been removed. 204 objects are located above the completeness limit
of
,
so that the slope of
finally is based on these stars.
We had mentioned in Sect. 2 that the standard deviation of the
members of Stock2 was larger than for all the other clusters and higher than
expected from the general accuracy of the Tycho-2 catalogue. Stock2 is
located within the Vatican declination zone (
to
)
of
the Astrographic Catalogue. A remark of Eichhorn (1974, p. 297)
that the accuracy of the Vatican measurements are comparably low could explain
this substandard behaviour.
Five stars in the field of Stock2 are mentioned by Ahumada & Lapasset
(1995) as blue stragglers belonging to the cluster. Three of
them - their stars 6, 31 (HD13402), and 45 - show membership
probabilities of P=0. The other two - stars 110 (HD13735) and 138
(HD13909) - have high probabilities of
,
however, after the
de-reddening process, they are located well in the main sequence of the
CMD. Therefore, we propose that all five objects are not blue straggler
members of Stock2.
Astrometric and photometric aspects of Praesepe have been well studied before (e.g. Jones & Cudworth 1983; Jones & Stauffer 1991), including IMF studies (Williams et al. 1995; Pinfield et al. 1999).
The Tycho-2 stars which were determined as members of Praesepe mostly are
located within a distance of
of the centre of the cluster. The CMD
shows a scattered main sequence, so that a differential reddening seems
likely. Since the total number of detected member stars brighter than
is only 79, a meaningful reddening map as determined for Stock2
(Sect. 3.2) could not be obtained. The isochrone shifted according to
the vL99a results fits the main sequence quite well, however, the isochrone
lies slightly too red, so that we propose a reddening of
.
Due to the age of Praesepe, which we determined to be 650 Myr, only a small range of main sequence stars could be covered. Therefore, it was not possible to derive a trustworthy IMF from the Tycho-2 data.
Ahumada & Lapasset's (1995) catalogue of blue stragglers contains five stars in the region of Praesepe. Three of them, HD73618, HD73666 (40Cnc), and HD73819 have proper motions leading to high membership probabilities (P>0.9). HD73711 shows an intermediate probability of P=0.79, but it is still likely that this object belongs to Praesepe. The final object, HD73120, is a clear non-member with P=0.
Although numerous photometric and kinematic studies are available for the
Per cluster (e.g. Mitchell 1960; Fresneau
1980; Prosser 1992,1994), the IMF of this object has
not been investigated since the work of Tarrab (1982).
The vector point plot diagram of the
Per cluster is presented in
Fig. 7. The good separation between field and cluster proper
motions is clearly visible: Only very few field stars might erroneously be
identified as cluster members. The width of the distribution of the cluster
stars is
which coincides with the
proposed standard deviation for individual objects included in the Tycho-2
catalogue. A circle indicates which stars were selected as the members of the
Per cluster.
![]() |
Figure 7:
Vector point plot diagram of the stars within
![]() ![]() ![]() |
Open with DEXTER |
vL99a found values of
and
for the
Per cluster. We could confirm these values with
the Tycho-2 data. On the other hand, together with the proposed age of
,
the upper part of the main sequence is not well represented by
the fitted isochrone. We prefer a younger isochrone with
,
which is plotted in Fig. 1. A higher age could only be explained
by the brightest members detected being blue stragglers.
Three stars are listed as blue stragglers by Ahumada & Lapasset
(1995). Two of them, HD22928 and HD24760, are too bright to be
included in the Tycho-2 catalogue with proper motions. One of them,
HD21428, has a membership probability of P=0.34 according to its Tycho-2
proper motion, so that it can be assumed that this object is a field star. The
proper motions of the other two, HD22928 (HIP17358) and HD24760
(HIP18532), can again be taken from the Hipparcos catalogue: Both have a
proper motion which indicates they are not members of
Per, since for
HD22928 the deviation of
from the cluster centre is approx.
,
and for HD24760,
is too far
away from the value for the members.
In Slettebak (1985), three objects are listed as Be or shell stars
in the region of the
Per cluster. HD21551 has a membership
probability of P=0.39 and HD22192 (
Per) of P=0.27 which means
that both objects are likely non-members. HD25940 (48Per) is not included
in the Tycho-2 catalogue, but it is a Hipparcos star (HIP19343) with a
proper motion of
,
.
Comparing this with the proper motion
centre of the cluster (Table 1), this means that this star,
too, shows a proper motion which does not coincide with the motion of the
cluster. Only the fourth star, HD21455, mentioned in the text of Slettebak
(1985) as a non-member, has a membership probability of P=0.76which makes it a likely cluster member from the proper motion point of view,
since the probability at which we distinguished members and non-members in the
case of
Per was P=0.6. Since Slettebak (1985)
mentions rotational velocities from Kraft (1967) and radial
velocities from Petrie & Heard (1969) as further membership criteria,
it is still possible that this object in fact is not a member of the cluster.
Blanco1 (also referred to as the
Scl cluster) turned out to be the
cluster with the second faintest completeness limit of our sample with
,
following the rough method presented earlier in
Sect. 2. Nevertheless, the IMF was computed on the basis of
34 stars only. The distance moduli of Blanco1 according to vL99a -
- and Robichon et al. (1999) -
- differ by almost
.
We found that an isochrone on the
basis of the vL99a value fits the CMD of the members slightly better, but
the Robichon et al. (1999) distance modulus cannot be excluded,
either. Instead of a reddening of
,
as proposed by
vL99a, we suggest a lower value of
.
This value is
in good agreement with other studies of Blanco1, e.g. Westerlund
et al. (1988). The cluster age as given in vL99a could be well
reproduced. The luminosity function of Blanco1 (see Fig. 8)
shows a gap in the range from approx.
to
(corresponding to masses from
to
). When experimenting with the lower mass limit during the IMF
determination, this led to a temporary increase of
from -2.27 to as
much as -1.12, which normalised after passing the gap. It must be remarked
that only 11 members have masses higher than
,
so that the IMF can
most reliably be described with the help of the lower mass stars. The gap
seems not to be a consequence of incomplete Tycho-2 data, since it has been
observed before (Eggen 1983). In addition, these gaps are not
uncommon for open clusters (Rachford & Canterna 2000). As a
consequence, Blanco1 is an example of a star cluster for which a power law
does not well describe the shape of the entire IMF. This is also illustrated
by the comparably high error of the maximum likelihood fit. It would be
interesting to obtain a deeper photometric study of this object to find out
whether the IMF shows a "normal'' power law shaped behaviour for lower mass
stars or if further gaps appear. An IMF computed on the basis of the
vL99a distance modulus leads to a slope of
.
The star
Scl itself is, according to Ahumada & Lapasset
(1995), a blue straggler cluster member. For this object we find a
membership probability of P=0 as a consequence of its proper motion of
,
.
![]() |
Figure 8:
VT luminosity function of Blanco1. Note the
dip around
![]() |
Open with DEXTER |
NGC2451 has been the subject of a discussion regarding whether it really is
an open star cluster or not (Eggen 1983,1986; Röser & Bastian
1994). Recently, Carrier et al. (1999) found
evidence that this object actually consists of two clusters in the same line
of sight, as proposed earlier by Röser & Bastian
(1994). Baumgardt (1998) was only able to confirm
existence of one of the two clusters on the basis of the Hipparcos and ACT
catalogues. Carrier et al. (1999) named the two objects
NGC2451A and B. We will stick to this nomenclature in the
following. Cluster A, which is the target of our study, is located at a
distance of approx. 190 pc, whereas component B is 360 pc (
)
away from the Sun. NGC2451B will not be studied here,
as its proper motion is too close to the centre of the field proper motion.
For NGC2451A, we found that adopting the vL99a reddening of
leads to an isochrone which is too far redwards to properly
represent the main sequence. vL99a does not distinguish between components A
and B. However, the values given for the proper motions (his Table 1) indicate
that he - as well as Robichon et al. (1999) - really does
investigate component A. Carrier et al. (1999) found the reddening
to be
,
which better corresponds to our isochrone
fit. Due to the lack of evolved stars, the age determination is a critical
task: A
isochrone fits the data as well as a 10 Myr
isochrone. Carrier et al. (1999) seem to overestimate the age of
NGC2451A with their value of 50 Myr.
With
,
NGC2451A shows the shallowest IMF of our
sample. Since the usable portion of the main sequence consists of 27 stars
only, this result is quite uncertain. Still, we can conclude that the IMF of
this object is not as steep as Scalo's (1998) average slope of
would suggest.
Ahumada & Lapasset (1995) mention two stars in the field of
NGC2451 as blue stragglers. One of them, HD61831, is a likely
cluster member with a probability of P=0.9. However, since it is located
slightly right of the fitted isochrone in the CMD, its nature as a blue
straggler is doubtful. On the other hand, the second object, HD63465, has a
proper motion of
,
which leads to a membership probability
of P=0. Following Carrier et al. (1999) who indirectly give
evidence for the proper motion of cluster B with the statement that
HR2968 (
,
)
belongs to NGC2451 B, HD63465 does not seem to be a
member of either of the clusters.
As for NGC2451A, IC2391 seems to be less reddened than proposed by
vL99a. Our result that IC2391 is unreddened can also be found, e.g.,
in Perry & Hill (1969). The main sequence of IC2391 (see
Fig. 1) is not very evenly populated: Several gaps of
can be found around
,
,
and
.
Below
only
very few stars with the proper motion of the cluster are present in the CMD,
and they are widely spread in colour. It seems that the data are complete down
to
.
IC2391 in general is sparsely populated, so
that we had to compute its IMF on the basis of only 29 stars. The result,
,
again is at the shallow end of the IMF slope
interval.
IC2602 shows a clear and narrow main sequence down to
.
Below this point, almost no stars are found, and the few detected members
are widely scattered in the CMD. In a positional plot of all stars with
of the Tycho-2 data (see
Fig. 9) one finds that in a
wide
area around the centre of IC2602, almost no stars are found. An equivalent
plot of, e.g., the GSC, Version 1.2, (Röser et al. 1998b) does not show
this decrease in stellar density. This might mean that for some reason the
stars fainter than
in this area - and hence members of
the cluster - are not included in the Tycho-2 catalogue. The total stellar
density, as one can see in images of this region, is low enough so that neither
crowding nor neighbouring bright stars might have disturbed the measurements
of either the Tycho star mapper or the first epoch data. We also checked the
stars in Tycho-2 that are listed without proper motions and the two
supplements of Tycho-2. Supplement 1 contains 31 stars, but only for one of
them is a proper motion given. There is an above-average number of these stars
in the region of IC2602, but when we added all these objects, the total
stellar numbers are still too low. The original Tycho catalogue does not show
this lack of stars.
corresponds to a mass of
,
and the
number of members which could be used for the IMF computation was no more
than ten, too small a number to compute the mass function of IC2602.
![]() |
Figure 9:
Plots of the stars in the region of IC2602 in
the magnitude range
![]() ![]() |
Open with DEXTER |
Ahumada & Lapasset (1995) list a "well known blue straggler''
for IC2602:
Car. It indeed has a Hipparcos proper motion close
to the centre of the cluster proper motion.
For the sparse open cluster NGC7092, reddenings of a few hundredths of a
magnitude can be found in the literature (e.g., Mohan & Sagar
1985; Manteiga et al. 1991; Nicolet 1981),
however, we found that no stars in the entire field studied (members and
non-members) justify an isochrone reddened by less than
.
Concerning the distance modulus, our result of
is in good agreement with other studies. The determination of the
IMF slope of NGC7092 led to a high error of
.
This cannot
only be explained with the low number of 25 members above the estimated
completeness level, which, with
,
is the brightest of the entire study besides IC2602. Another point is the shape of the IMF
which cannot be well represented by a power law: The CMD of NGC7092 shows a
magnitude region with (almost) no stars, in this case for
.
In terms of masses, this interval corresponds
to
.
This again underlines the advantages
of a maximum likelihood estimation compared with a least square fit to an IMF
histogram, which would have led to an error of
only,
so that the unusual shape of the IMF would not have been striking.
![]() |
Figure 10: Initial mass functions of six of the open star clusters using Tycho-2 data. The Pleiades IMF can be found in Fig. 4. For Praesepe and IC2602, no reliable IMF could be determined and their histograms are not shown. For more details about the plots, see the explanations in Fig. 4 |
Open with DEXTER |
cluster | # | ![]() |
mass range | VT range |
stars | [![]() |
[mag] | ||
Blanco1 | 34 |
![]() |
[1.1; 4.8] | [6.1; 11.4] |
Stock2 | 204 |
![]() |
[1.5; 4.1] | [7.6; 11.0] |
![]() |
70 |
![]() |
[1.1; 6.8] | [5.0; 10.5] |
Pleiades | 127 |
![]() |
[1.0; 4.1] | [5.0; 10.9] |
NGC2451A | 27 |
![]() |
[1.3; 6.8] | [4.8; 10.0] |
IC2391 | 29 |
![]() |
[1.1; 8.1] | [3.5; 10.7] |
NGC7092 | 25 |
![]() |
[1.4; 3.4] | [6.5; 9.9] |
We used the information of the Tycho-2 catalogue for the determination of
the IMF of nearby open star clusters in the
mass range and to
demonstrate the capabilities of Tycho-2.
The IMF of seven open star clusters could be represented by power laws with
exponents ranging from
to
.
Within the errors,
this is in good agreement with the range of IMF slopes as given, e.g., by
Scalo (1998). An average slope from our results could be computed as
.
The standard deviation of the sample of slopes equals
,
whereas the standard deviation as derived from the individual
errors of the slopes is only
.
Two of the clusters, Blanco1 and NGC7092, show significant gaps over
almost
in luminosity or a few tenths of a Solar mass,
influencing the slope determined and their errors. For these clusters, further
studies with fainter limiting magnitudes will be needed to find a more precise
IMF, since not only would the mass range be larger, but also the numbers of
members to be detected.
With a Monte Carlo simulation, we tested whether our assumption that the IMF of
the clusters can be represented with a (single) power law within the mass
range covered by our study: We repeatedly simulated the mass functions of a
cluster with the same number of members as found from Tycho-2 and with the
slope derived from the data. From these samples, we determined the slopes with
the same method as for the actual objects. A probable difference between the
errors as found from the maximum likelihood method and the standard deviation
of the
values found from the simulation might be an indication
that the IMF of that object does not follow a power law. All clusters
show an error smaller by 0.15 in the simulated clusters. This could not only
be explained by the fact that the "true'' IMF of the clusters does not follow
a power law, but also by the slight contamination of the data with field
stars (See the remarks in Sect. 2). The influence of these
objects at arbitrary distances cannot be well simulated in the Monte Carlo
experiments. The only exception is NGC7092, with a difference of 0.4. The
gaps in the main sequence mentioned in Sect. 3.9 and the low
number of total stars might be the cause for this behaviour. However, we do
not find evidence in any of the clusters that the choice of a power law as
the representative of the IMF was clearly incorrect.
Comparing our results with the study of Tarrab (1982), we find good
agreement for four of the five clusters which are discussed in both
publications. Only the slope of NGC7092, for which Tarrab found - in our
notation -
,
dramatically differs from our value of
.
Tycho-2 allows us to discuss the membership of individual objects
mentioned previously in the literature. We inspected the catalogue of blue
stragglers of Ahumada & Lapasset (1995). Many of their objects
belong to their "class 3'', which means that there is only weak evidence (like
the location in a CMD) for stars being blue stragglers and members of the
clusters. For these objects, if they are bright enough, proper motions are
now available, so that the list can be improved by cleaning it of field
stars. The second paper under consideration was Slettebak's (1985)
list of Be and shell stars in open clusters. While we could confirm the list
for the Pleiades, the nature of the stars mentioned in the context of the
Per cluster has to be revised.
Our analysis showed that the completeness given for the Tycho-2 catalogue by
Høg et al. (2000) is too optimistic, at least in several fields of
the clusters studied: A field with a diameter of approx.
around the
centre of IC2602 contains just a few stars with
,
in
contrast to other star catalogues. It should be remarked that this result is
not in contradiction with the statement about the global completeness of
Tycho-2 in Høg et al. (2000). However, for studies in small
excerpts of the Tycho-2 data, one has to be aware of these regional variations
in the completeness. On the other hand, the quality of the Tycho-2 proper
motions appears to be as reliable as the Hipparcos data.
Regarding upcoming missions of astrometry satellites - FAME (Horner et al. 1998), DIVA (Röser et al. 1998a), or GAIA (Lindegren & Perryman 1998) - studies based on uniform all-sky data samples like the IMF investigations demonstrated in this work will play a more important role in the future: The limiting magnitudes of these instruments will be fainter, proper motions will be available for all detected objects (in contrast to, e.g., the Hipparcos catalogue, the proper motions of which were determined for the stars of an input catalogue only) and will become independent of ground based observations, so that their data quality will be higher. This means that - concerning the determination of mass functions - the IMF of nearby clusters will be covered down to lower masses. Due to the fainter limiting magnitude and the higher accuracy of the proper motions (and hence a more precise membership determination) a larger number of open star clusters can then be studied.
Acknowledgements
The authors acknowledge Santi Cassisi for providing the isochrones necessary for our studies. For valuable discussions, we thank Andrea Dieball (especially for the hints concerning the differential reddening and the study of the stellar density) and Ram Sagar. Thanks to Klaas S. de Boer and Andrea Dieball for carefully reading the manuscript of this publication. This research has made use of NASA's Astrophysics Data System Bibliographic Services, the CDS data archive in Strasbourg, France, and J.-C. Mermilliod's WEBDA database of open star clusters.