A&A 470, 585-596 (2007)
DOI: 10.1051/0004-6361:20067003
L. Balaguer-Núñez1,2,3 - D. Galadí-Enríquez4 - C. Jordi1,5
1 - Departament d'Astronomia i Meteorologia, Universitat de
Barcelona, Avda. Diagonal 647, 08028 Barcelona, Spain
2 -
Shanghai Astronomical Observatory, CAS Shanghai 200030, PR China
3 -
Institute of Astronomy, Madingley Road, CB3 OHA Cambridge, UK
4 -
Instituto de Astrofísica de Andalucía (CSIC).
Camino Bajo de Huétor 50, 18008 Granada, Spain
5 -
Institut d'Estudis Espacials de Catalunya - IEEC, Edif. Nexus
Gran Capità 2-4, 08034 Barcelona, Spain
Received 21 December 2006 / Accepted 13 April 2007
Abstract
Following deep astrometric and photometric study of the cluster
NGC 2682 (M 67), we are able to accurately determine its fundamental
parameters. Being an old and rich cluster, M 67 is relevant
for the analysis of the Galactic disk evolution.
M 67 is well studied but lacks a wide and deep
Strömgren photometric study. The brightest stars of the open cluster M 67 were used as
standard stars in our studies of NGC 1817
and NGC 2548, and the extension of the field covered, as well as the
number of observations, allowed us to obtain the best set of
Strömgren data ever published for this cluster.
We discuss the results of our CCD
intermediate-band
photometry, covering an area of about 50
50
down to
19. Moreover, a complete membership segregation
based on astrometric and photometric criteria is obtained.
The photometric analysis of a selected sample of stars yields a reddening value of E(b-y) = 0.03
0.03, a distance modulus of V0-MV = 9.7
0.2 and [Fe/H] = 0.01
0.14.
Through isochrone fitting we found an age of
= 9.6
0.1
(4.2
0.2 Gyr). A clump of approximately 60 stars around
V = 16, (b-y) = 0.4 could be interpreted as a population
of pre-cataclysmic variable stars (if members), or as a stream of
field G-type stars placed at twice
the distance of the cluster (if non-members).
Key words: Galaxy: open clusters and associations: individual: NGC 2682 - techniques: photometric - astrometry - methods: observational - methods: data analysis
NGC 2682 (C0847+120), also known as M 67, in Cancer
[
= 8
,
;
,
]
is probably the most
thoroughly studied old open cluster
in the Galaxy, thanks to its small distance from us
(estimated to be
900 pc). Typically quoted
values for the age of the cluster (
4 Gyr) place it among the
oldest open clusters. In our astrometric and photometric
study of several open clusters (Balaguer-Núñez et al. 2004a,2005),
the photometric standard stars were taken from this cluster (Nissen et al. 1987),
thus high-quality Strömgren wide-field CCD photometry of it also
resulted from our observations. The existence of a proper motions study (Zhao et al. 1993) of similar quality and characteristics to those used in our study of other clusters
(based on the same plate material and methods, by the Shanghai
Astronomical Observatory), allowed us to apply the same techniques and analysis tools
to this cluster. This paper on M 67 closes the series at this stage of our open
cluster programme.
Photometric studies of M 67 have been performed by
several authors. Montgomery et al. (1993)
have presented a deep (
20)
CCD photometry of 1468 stars
within 15
of the centre of the cluster, and Fan et al. (1996) have studied
spectrophotometry of similar depth in nine BATC intermediate-band filters
for stars in a 1
92
1
92 area centred on the cluster.
In addition, variability studies have been performed in this cluster, most
notably by Gilliland et al. (1993), who conducted a very sensitive, highly
temporally sampled study of stars in the central few arcmin of the cluster.
Stassun et al. (2002), using differential CCD photometry to search for
variability, obtained sensitive photometry of 990 stars in a roughly
square region one-third of a degree on a side centred 5
north of the cluster centre. Sandquist (2004), also as a by-product of variability studies, conducted a high precision
colour-magnitude diagram analysis.
Nissen et al. (1987) obtained accurate
photoelectric
photometry of a sample of 79 stars within a radius of 10
from the centre, from the subgiant branch down to the unevolved main sequence. For this reason
stars in M 67 were used by us as photometric standard stars for the
observations of NGC 1817 and NGC 2548
(Balaguer-Núñez et al. 2004a,2005). The number of observations taken in our programme allowed us to obtain
high quality photometry of this cluster.
In spite of being a well studied cluster,
our Strömgren photometry deserves an analysis due both to its
deepness and to the extension of the covered area.
Table 1: Log of the observations used in this study.
There are also a number of proper motion and radial velocity studies.
Sanders (1977) calculated probabilities of membership based on
relative proper motions of 1866 stars in the cluster field with a limiting photographic magnitude of 17. From 10 plate pairs of different origin and with a maximum epoch difference of 68 years, he found 649 probable members. Girard et al. (1989) gave relative proper motions for 663 stars in an area
of 42
34
from 44 plates with a maximum separation of 66 years. Although the plates used were of different depths, the deepest plates provided a limiting visual
magnitude down to 16. In 1993, Zhao et al. derived relative proper motions for 1046 stars
within a 1
5
1
5 area in the region from PDS measurements
of 9 plates with a maximum epoch difference of 80 years and a magnitude
limit
.
On the other hand, the radial velocity studies by Mathieu et al. (1986, 1990) gave precise radial-velocity measurements for 170 stars, including all main-sequence stars brighter than V = 12.8.
In this paper we discuss the results of our CCD photometric study,
covering an area of
about 50
50
(Fig. 1)
down to V
19. Section 2 contains the details of the CCD observations and their
reduction and transformation to the standard system. In Sect. 3 we discuss
a new membership segregation based on the combination of
parametric and non-parametric methods applied to the above-mentioned
proper motions from Zhao et al. (1993).
This study was selected because it covers
the largest area and epoch difference.
Also, this material is fully homogeneous with that used by us for the other
clusters in this series, from the point of view of telescope and plates
and proper motion treatment and reduction.
This section includes a discussion of the limits on proper motion detection
of binarity.
In Sect. 4 we discuss the colour-magnitude diagram and identify the sample of probable cluster members using astrometric as well as photometric criteria. Section 5 contains the derivation of the fundamental cluster parameters of reddening, distance, metallicity and age. Section 6 studies the multiple star systems, blue stragglers and a new feature on the colour-magnitude diagram. Section 7 summarizes our conclusions.
The Strömgren CCD photometry of the area was obtained at Calar
Alto Observatory (Almería, Spain) in January 2000
using the 1.23 m telescope of Centro Astronómico Hispano-Alemán (CAHA).
Further data were obtained at Observatorio del Roque de los Muchachos
(ORM, La Palma, Canary Islands, Spain) in February 2000 using the 2.5 m
Isaac Newton Telescope (INT) of ING (equipped with the Wide-Field Camera, WFC),
and in December 1998 and February 2000 using the 1 m Jakobus Kapteyn
Telescope (JKT) of ING, with the H filters. A log of the
observations, the total number of frames,
exposure times and seeing conditions is given in Table 1.
We obtained photometry for a total of 1843 stars in an area of
50
50
in M 67 region, down to a limiting magnitude V
19.
The area covered is shown in the finding chart
of the cluster (Fig. 1).
Due to the lack of H
filter at the WFC-INT, it was only possible to
measure it at the JKT and CAHA telescopes, thus limiting the spatial coverage
with this filter.
Only 288 stars (in the central region) have H
values.
![]() |
Figure 1:
Finding chart of the area under study. The covered area
is marked in black on an image of a plate (POSSI.E -optical R-.DSS1.486-LOW) plotted with
Aladin (Bonnarel et al. 2000). The H![]() |
Open with DEXTER |
The reduction of the photometry is explained in length at Balaguer-Núñez et al. (2004a). Having the standard stars in the same images, the process is, in this case, equivalent to performing differential photometry.
Table 2:
Number of stars observed (N) and mean internal errors ()
as a function of apparent visual magnitude.
Our general procedure has been to routinely obtain twilight
sky flats for all the filters
and a sizeable sample of bias frames (around 10) before and/or after every run.
Flat fields are typically fewer in number, from five to ten per filter.
Two or three dark frames of 2000 s were also taken. IRAF routines were used for the reduction process.
The bias level was evaluated individually for each frame by averaging the counts
of the most stable pixels in the overscan areas. The 2D structure of the bias
current was evaluated from the average of a number of dark frames with
zero exposure time. Dark current was found to be negligible in all the cases.
Flatfielding was performed using sigma clipped, median stacked, dithered twilight flats.
Ten short exposures in each filter were taken every night with a magnitude limit of V
19.
At least one exposure per filter was pointed to make the centre of the
cluster fall on each of the WFC chips.
Our fields are not crowded. Thus, the synthetic aperture technique provides the most efficient measurements of relative fluxes within the frames and from frame to frame. We used the appropiate IRAF packages, and DAOPHOT and DAOGROW algorithms (Stetson 1990,1987). We analyzed the magnitude growth curves and determined the aperture correction with the IRAF routine MKAPFILE.
For WFC images from the INT, we employ the pipeline specifically developed by the Cambridge Astronomical Survey Unit. The process bias subtracts, gain corrects and flatfields the images. Catalogues are generated using algorithms described in Irwin (1985). The pipeline gives accurate positions in right ascension and declination linked to the USNO2 Catalogue (Monet et al. 1998), and instrumental magnitudes with their corresponding errors. A complete description can be found in Irwin & Lewis (2001) and in http://www.ast.cam.ac.uk/~wfcsur/index.php.
The coefficients of the transformation equations were computed
by a least squares method using the instrumental magnitudes of the
standard stars and the standard magnitudes and colours in the
system. Up to 68 standard stars in the field of the cluster
(Nissen et al. 1987) were used depending on the size of the frame.
Those standard stars with residuals greater than 2
were rejected.
The reduction was performed for each night independently and in two steps.
The first step is to determine the extinction coefficients for
each passband from the standard stars. With the extinction coefficients
fixed, the transformation coefficients to the standard system
were fitted. The final errors as a function of apparent visual magnitude
are given in Table 2 and plotted in Fig. 2.
![]() |
Figure 2: Mean internal errors of magnitude and colours (including the effects both from the instrumental measurements and the transformation uncertainties) as a function of the apparent visual magnitude, V, for all observed stars in the cluster region. The structure in the magnitude dependence is mainly due to different numbers of measurements for individual stars since the centre of the cluster was shifted to be observed by each of the four CCD chips. |
Open with DEXTER |
Table 3 lists the
data for all 1843 stars observed in a region of 50
50
around the open cluster M 67 (Fig. 1). Star positions
are given as frame (x,y) and equatorial (
,
) coordinates.
An identification number was assigned to each star following
the order of increasing right ascension. Column 1 is the ordinal star number; Cols. 2 and 3
are
and
;
Cols. 4 and 5 are the respective x, y coordinates in arcmin (arbitrary origin, close to cluster centre);
Cols. 6 and 7 are the (b-y) and its error, Cols. 8 and 9 the V magnitude and its error,
Cols. 10 and 11 the m1 and its error, Cols. 12 and 13 the c1 and its error,
and Cols. 14 and 15 the H
and its error.
In Col. 16, stars considered candidate members (see Sect. 4.1) are labelled "M'', while those classified as non-members have the label "NM''.
The cross-identification of stars in common with the astrometry (Zhao et al. 1993), BDA (http://obswww.unige.ch/WEBDA), Sanders (1977), Hipparcos (ESA 1997), Tycho-2 (Høg et al. 2000) and USNO-2 (Monet et al. 1998) catalogues is provided in Table 4.
After Nissen et al. (1987), few studies of this cluster have been carried out
in Strömgren photometry. Anthony-Twarog (1987) used a CCD to cover a larger
area with up
to 44 stars in common with our sample. Comparing Anthony-Twarog's
data with ours, the mean differences in the sense ours minus others,
we get -0.01( = 0.04) in V,
0.00(0.02) in b-y, 0.00(0.04) in m1 and 0.00(0.04) in c1.
The data published by Joner & Taylor (1997), with 14 stars in common with our catalogue,
up to a limit of
13, give differences of
0.01(0.03) in V, 0.01(0.02) in b-y, -0.02(0.03) in m1 and
-0.01(0.04) in c1.
Before Nissen et al. (1987), the study of Strom et al. (1971) with 19 common stars
up to a limit of
13.5 gives:
-0.02(0.03) in V, -0.03(0.02) in b-y, 0.03(0.05) in m1,
-0.02(0.06) in c1 and -0.02(0.07) in H
,
but a large scatter
and a probable colour term in their results was already noted by Nissen et al. (1987).
The V magnitude derived from the y filter can be compared to the published broadband data. Several studies of M 67 give us mean differences in V, in the sense ours minus others (see Table 5). Looking at the differences, one concludes that there are no apparent systematic trends in our photometry.
Table 5: Comparison of the V magnitude to the published broadband data.
We have studied the membership segregation using parametric and
non-parametric approaches in the same way as in Balaguer-Núñez et al. (2005,2004b).
Thanks to the proper motions obtained by Zhao et al. (1993) from high-quality plates
taken with the double astrograph at the Zo-Sè station of the
Shanghai Observatory, as in the case of NGC 1817 and NGC 2548,
we have been able to use the same procedure in a very homogeneous way.
The focal length of the Gautier 40 cm double astrograph at the
Zo-Sè station is 6.9 m (hence a plate scale
of 30
mm-1). All the plates were measured using the Photometric
Data Systems model 1010 automatic measuring machine at the Purple Mountain
Observatory in Nanjing. Zhao et al. give relative proper motions of
1046 stars within a 1
5
1
5 area,
from measurements of 9 plates. The plates have a maximum epoch difference of 80 years.
Zhao et al. (1993) applied the plate-pair technique adopted many times
at Shanghai Observatory to derive the relative proper motions.
All linear and quadratic coordinate-dependent terms and the coma term were included
in the plate-pair solution. The quoted mean errors of the relative proper motions vary
from 0.4 mas yr-1 for bright stars in the inner part of the cluster
field to some 1.5 mas yr-1 for faint stars in the outer part of the cluster.
The comparison with the proper motions of Sanders (1977) and Girard et al. (1989)
shows a satisfactory agreement.
The magnitude limit of Zhao et al.'s study is
.
There are 50 stars in common with the Tycho-2 Catalogue (Høg et al. 2000).
From their absolute proper motions we can calculate the transformation of
Zhao et al.'s relative proper motions (
)
to the absolute reference frame:
Membership determination in Zhao et al. (1993) was calculated
with an 8-parametric Gaussian model, and a list of stars with
probability higher than 0.8 and a distance to the centre less than 45
contained 282 cluster members. For coherence with our analysis of NGC 1817 and NGC 2548
(Balaguer-Núñez et al. 2005,2004b),
we apply a 9-parametric Gaussian model and a non-parametric method
to the proper motions obtained by Zhao et al.
In contrast to their approach to the segregation of cluster members,
we do not use any spatial information. Selecting members on kinematic and
photometric criteria leaves the spatial information untouched, ready for
a clean interpretation of radial and spatial trends.
Zhao et al. give x,y coordinates, relative proper motions with their errors, number of plates used for proper motion determination, and cross-identifications with Sanders (1977). The authors state that their x,y coordinates are taken from one of the plates, but a detailed analysis of the data shows that this is the case for the majority of the stars, but there is a small subset (61 stars) whose x,y coordinates seem to have been taken from at least two other plates, and this introduces shifts in their quoted positions. Fortunately, Zhao et al.'s cross-identifications with Sanders (1977) were correct, which allowed us a complete cross-identification of our photometric catalogue with the proper motion data, not losing any stars. This inhomogeneity in the Zhao et al. (1993) x,y data does not affect our results, since we do not use these coordinates in our study. All tests performed show that Zhao et al.'s proper motions are not affected. The mistake in the x,y coordinates of 61 stars seems to have been introduced in the final preparation of the tables by Zhao et al. (1993).
We consider any observational plane (a, b), where a and b stand for any physical magnitudes.
The distribution of any stellar sample on this plane can be described
by its frequency function, ,
measured in units of number
of stars per unit area. The integral of the frequency function all over
the plane, i.e., its volume, gives the number of members of the
sample, N.
If the frequency fuction is normalised to unit volume,
then we get the probability density function (PDF),
,
measured in units of fractions of the sample per unit area.
If two populations are superposed on the same region of the
observational plane, but following different frequency functions
and
,
then, for any point (a,b) on that plane, it is possible to compute the probability P for a star placed at that point of belonging to population
:
Table 6:
Distribution parameters and their uncertainties
from a 9-parametric Gaussian model applied to M 67 cluster and the field.
measures the volume of the cluster frequency function
(
),
and
are the mean relative proper motions,
,
and
and
are the dispersions of the
Gaussians fitted to the cluster and field distributions,
and
is the correlation for the field function (it fixes the
orientation of the field elliptical function). The units of
and
are mas yr-1.
The parametric method assumes the existence of two populations
in the VPD: cluster and field. The corresponding frequency functions
are modelled as parametric Gaussian functions: a circular Gaussian model
is adopted for the cluster distribution, while a bivariate (elliptical)
Gaussian is admitted to describe the field. We apply
a 9-parametric Gaussian model following Balaguer-Núñez et al.
(2004b, 2005)
and the results are shown in Table 6.
Our implementation takes into account the individual errors of
the proper motion measurements (Zhao & He 1990). If we transform the obtained
mean relative proper motion of the cluster to the absolute
system by the equations above, we get
(
) = (-7.1
0.8, -7.6
0.4) mas yr-1.
![]() |
Figure 3:
Empirical PDFs in the vector-point diagram. Top:
![]() ![]() ![]() |
Open with DEXTER |
The cluster/field segregation from astrometry has also been analysed
with a non-parametric approach, as explained in length in Balaguer-Núñez et al. (2005,2004b).
We perform an empirical determination of the frequency functions
from the vector-point diagram (VPD),
without relying on any previous assumption about their profiles.
In the area occupied by the cluster, the frequency function
is
made up of two contributions: cluster,
,
and field,
.
To disentangle the two populations, we studied the VPD for the plate area
outside a circle centred on the cluster. To this end, the centre of the cluster was chosen as the point of highest spatial density. The only two assumptions that we need to apply are (1) that it is possible to determine the frequency function of the field stars from some spatial
area free of cluster members (outside the circle in our case); and, (2) that
this frequency function is representative of the field frequency function
in the area occupied by the cluster.
We did tests with circles of very different radii searching for a reasonable tradeoff between cleanness (absence of a significant amount of cluster members)
and signal-to-noise ratio (working area not too small, to avoid small number statistics).
The kernel density estimator technique (Hand 1982) was
applied in the VPD to these data. Details of the procedure
can be found in Galadí-Enríquez et al. (1998). Circular Gaussian kernel functions were used, with
Gaussian dispersion (also called smoothing parameter) h elected according to
Silverman's rule (Silverman 1986).
This way, the empirical frequency functions were computed
for a grid with cell size of 0.2 mas yr-1, well below the proper
motion errors. The procedure was tested for several subsamples applying different
proper motion cutoffs and the adopted one, as in the cases of NGC 1817 and
NGC 2548, is of
15 mas yr-1.
We finally find that the area outside a circle with a radius of 35
centred on the cluster yields a clean frequency function with low cluster contamination and low noise.
We can scale this field frequency function to represent
the field frequency function in the inner circle,
,
by simply applying a factor linked to the area. The cluster
empirical frequency function can then be determined as
=
.
Normalizing the frequency functions to volume unity we get the
probability density functions for the mixed
population (inside the circle), for the field
(outside the circle) and for the cluster (non-field)
population. These PDFs are displayed in Fig. 3.
We estimate the typical noise level,
,
present in the result
and we restricted the probability calculations to the stars with cluster PDF
3
.
The maximum of the cluster PDF is located at (
) = (-0.6
0.2, 0.4
0.2) mas yr-1, which coincides well with the values obtained in the parametric method.
If we transform this value to the absolute system, we get an absolute
proper motion for the cluster of (
) = (-7.1
0.7, -7.7
0.5) mas yr-1. Thirty six out of the 50 Tycho-2 stars are classified as members giving a mean value of (
) = (-8.4
2.1, -6.1
2.2) mas yr-1. This is fully compatible with our result, but we consider
our figures more reliable, since they are derived from the whole sample
of stars.
The non-parametric approach does not take
into account the errors of the individual proper motions.
But the full width at half maximum (FWHM) of the PDF of the
cluster gives an estimation of the width of the distribution.
We obtained a FWHM of 4.1
0.2 mas yr-1.
Taking into account the Gaussian dispersion owed to the smoothing
parameter h = 1.34 mas yr-1, this would correspond to a rms error on proper motions of 1.55 mas yr-1. But from Zhao et al. (1993) we know that the mean proper motion
precision is 1.24 mas yr-1,
what gives us an intrinsic dispersion component of 0.93 mas yr-1,
of the same order that the value obtained by the
parametric membership determination (
= 0.89
0.05 mas yr-1,
4 km s-1 at the distance of 900 pc).
![]() |
Figure 4:
The histogram of cluster membership probability
of M 67. The solid line gives the results for classical parametric
method ![]() ![]() |
Open with DEXTER |
Membership probabilities are computed according to Eq. (1), both from the classical, parametric frequency functions ()
and from the non-parametric, empirical
frequency funcions (
). The cluster membership probability histogram (Fig. 4)
shows a clear separation between cluster members and field stars
in both approaches (the solid line being the classical parametric method,
dotted line the non-parametric approach). The non-parametric technique yields
an expected number of cluster members from the integrated volume of the
cluster frequency function in the VPD areas of high cluster density,
where PDF
3
.
The expected number of cluster members in this VPD area is
of 393. Sorting the sample in order of decreasing
non-parametric membership probability,
,
the first 393 stars
are the most probable cluster members. The minimum value of the non-parametric
probability (for the 393-rd star) is
= 0.87.
To decide where to set the limit among members and non-members in the list sorted in order of
decreasing parametric membership probability, ,
we accept
the size of the cluster predicted by the non-parametric method, 393 stars. Thus we consider that the 393 stars of highest
are the most probable members, according to the results of the parametric technique. The minimum value of the parametric probability
(for the 393-rd star) is
.
Contrary to what happened with the previous two clusters studied by us, and in other cases
found in the literature (Galadí-Enríquez et al. 1998), in this case the original parametric
segregation (Zhao et al. 1993) at
was underestimating, and not
overestimating, the number of probable members.
This is because M 67 is much more populated, and displays
a higher constrast against the field than the other clusters.
The need to reach a
cutoff much lower than the
cutoff
illustrates the effects of an imperfect modelling of the wings of the real
distribution in the parametric method.
Also, the importance of having some objective criterion to
decide the
cutoff is evident from our results.
With these limiting probabilities (
;
),
we get a 97
(1014 stars) agreement
in the segregation yielded by the two methods.
Table 7 lists the cross-identifications with Zhao et al. (1993)
and Sanders (1977), the
and the
for the 1046 stars.
As in Balaguer-Núñez et al. (2004b), we accept as probable members those stars classified as member by at least one of the two methods. This way we get a list of 412 astrometric probable member stars. Figure 5 shows the proper motion VPD and the sky distribution for all the measured stars, where empty circles denote an astrometric probable member of M 67, and all other stars are considered field stars and they are indicated by dots.
![]() |
Figure 5: The proper motion vector-point diagram (VPD) and spatial distribution of stars in M 67 region (empty circles for the astrometric probable members of M 67 and dots for field stars). Units of relative proper motions are mas yr-1. |
Open with DEXTER |
A comparison of the cluster/field segregation for the 143 stars
in common with the radial velocity studies
by Mathieu et al. (1990,1986) is given
in Tables 8 (candidate members) and 9 (non-members).
The radial velocities have typical standard deviations
for a set of measurements of any given star that range
between 0.5 and 0.8 km s-1.
To quantify the differences in the segregation, we do as in
Balaguer-Núñez et al. (2005,2004b) and set an agreement
index
to 1 if the parametric probability,
,
agrees with the radial velocity
segregation, 2 if the non-parametric probability,
,
agrees,
3 if both probabilities,
and
,
agree and 0 if none does.
We find 125 out of 143 stars with
,
that is 87
agreement
with the radial velocities segregation. This 13
disagreement consists of
15 stars out of 41 being considered non-members on the basis
of proper motions, while only 3 out of 102 were found to be astrometric
members but considered non-members on the basis of radial velocities.
Those three have only one measurement of their radial velocities
and two of them are known to be spectroscopic binaries (SB),
therefore the measurement may be different to the systemic velocity and thus
they cannot be completely ruled out as members.
If we compare the parametric and non-parametric methods, the behaviour
is similar. For the
parametric method we find a total of 121 stars (85)
whose membership
assignation coincides with the radial velocity criterion, while
for the non-parametric method this amounts to 124 stars (87
).
![]() |
Figure 6:
UCAC2 proper motion vs. magnitude for M 67 area. Top row:
![]() ![]() ![]() |
Open with DEXTER |
The internal dynamics of non-resolved multiple systems could induce a motion of the measured photocentre different from the motion of the barycentre. This effect has to be null for equal-mass binaries and almost undetectable for binaries with components of extremely different masses, and should reach a maximum somewhere in between. When looking at the histogram of proper motion modules, this effect would skew the distribution inducing a longer tail towards high values.
Bica & Bonatto (2005) have made an attempt to
identify high-velocity stars as unresolved binary cluster members and, as a result, to map and quantify the binary component in colour-magnitude diagrams.
They analysed 9 open clusters using UCAC2 proper motions (Zacharias et al. 2004)
and 2MASS photometry (Skrutskie et al. 2006). As a test case of the method, they use M 67.
Their study of the modulus of the proper motion in the cluster area yields a double-peaked
distribution instead of the skewed Gaussian that, in principle,
could be expected. The authors associate the peaks with high-velocity and low-velocity populations in the cluster, and the high-velocity group is interpreted as produced by
unresolved binary systems. This behaviour is found in the core and in the halo
of the cluster. The authors claim internal dispersions of the order of
6 km s-1 for single stars and around
11 km s-1for unresolved binaries, and a difference between the two peaks
of more than 20-30 km s-1. These strikingly unrealistic numbers and the
lack of any error estimation adds suspicion to the fact that
all the stars from the high-velocity peak are
dimmer than
12 (Fig. 6 in their paper), when the abundance
of binaries among the bright members of M 67 is well known.
A similar study using the proper motions by Zhao et al. (1993) gives as a result a bell-shaped distribution with a tail towards higher velocities. We were not able to find a double-peaked distribution in the core nor in the halo. Different tries of dividing the distribution did not show any relation to binarity when checked on the corresponding colour-magnitude diagram.
As is well depicted by the study of Platais et al. (2003), the presence of
residual colour/magnitude terms in any proper motion analysis is almost
unavoidable. This is especially true when one deals with an all-sky catalogue,
where the proper motions are derived from different first epoch sources. Aware
of these facts, we have checked in UCAC2 (Zacharias et al. 2004) proper motions for magnitude terms.
We plot proper motion vs. magnitude in Fig. 6.
The UCmag magnitude in the top row is in the UCAC bandpass
(579-642 nm, between V and R) and should be considered approximate.
UCAC2 proper motions are derived by using over 140 ground- and space-based
catalogues, including Tycho-2 (Høg et al. 2000) and the AC2000.2 (Urban et al. 1998).
In the bottom row, we have plotted the 2MASS J magnitude as a double check of the results.
Cluster member candidates from our selection
are marked as empty circles. The UCAC2 proper motions of these stars
show a trend with magnitude for magnitudes fainter than from ,
highlighted by the regression fit marked by a thick line.
The behaviour detected by Bica & Bonatto (2005) is not related to binarity, but due
to systematics in the proper motion data.
The use of proper motion data from UCAC2 Catalogue should be carefully
performed to avoid systematics.
![]() |
Figure 7: Colour-magnitude diagrams of the M 67 area. Figure on the right shows the astrometric probable members as empty circles. Thick line is a shifted ZAMS, with the chosen margin for candidate members (V+0.5, V-1) in thin lines. See text for details. |
Open with DEXTER |
We use the V vs. (v-y) colour-magnitude diagram for our study (Fig. 7 right and centre) because it defines the main-sequence of a cluster significantly better than the traditional V vs. (b-y) diagram (Meibom 2000) in the presence of extinction. The colour-magnitude diagram of all the stars in the area displays a fairly well defined main sequence. The advanced age of this cluster is obvious at first sight. Moreover, the sequence of binaries can be easily followed.
Unfortunately, proper motions are only available for the
brightest stars in the area. Photometric
measurements help to reduce the possible field contamination in the
proper motion membership - among bright stars -, as well as to enlarge the selection of members towards faint magnitudes. Our astrometric segregation of member stars has a limiting magnitude
of
15.5. From this magnitude down to
we construct a ridge line
following a fitting of the observational zero-age main sequence (ZAMS)
(Hilditch et al. 1983; Crawford 1979; Olsen 1984; Crawford 1978,1975), on the V vs. (v-y) diagram.
A selection of stars based on the distance to this ridge line is then
performed. The width of the sequence defined by astrometric candidate members allows us to establish the margins for the selection of fainter stars. The chosen margin includes all the stars
between V+0.5 and V-1 from the ridge line, as shown in the
right panel of Fig. 7. This magnitude range allows for photometric errors in the two axes,
intrinsic dispersion around the ZAMS and multiplicity effects.
The interval of 1 mag at the bright side comprises the width of the sequence and is large enough to allow for equal-mass binaries (
)
and for most triples (
for equal-mass systems). The interval of 0.5 mag on the faint side
agrees with the observed width of the main sequence and is larger
than the expected photometric errors.
Luckily, the high density of M 67 and its outstanding contrast against the
population of field stars (shown by the comparison of the density functions
in Fig. 3) makes the fine-tuning of these margins
non-critical: the number of non-members introduced by our choice
is low, and the loss of true members is minimised.
![]() |
Figure 8: The colour-colour diagrams of M 67. Empty circles denote candidate members of M 67, chosen with astrometric and non-astrometric criteria as explained in Sect. 4.1. The thick line is the solar ZAMS standard relation shifted by E(b-y) = 0.03 when necessary. Evolutionary and binarity effects are noticeable in these diagrams. Some of the candidate members, known to be multiple stars, have very discordant values. |
Open with DEXTER |
Extreme outliers have been rejected making use of the colour-colour
diagrams (Fig. 8) with the help of the standard relations
(Hilditch et al. 1983; Crawford 1979; Olsen 1984; Crawford 1978,1975). The astrometric and photometric selection gives
a final set of 776 candidate members. They are plotted
in Fig. 8 as filled circles in the
[m1] - [c1],
m1 - (b-y),
c1 - (b-y) and
diagrams. This photometric selection is known to contain a certain amount of field contamination, so we prefer to talk about "candidate members'',
rather than "bona fide members''. The selection done leads to a set of
stars overwhelmingly dominated by true members:
the methods applied in Sect. 5 to derive the physical parameters of the cluster are readily capable of extracting reliable mean values from such a sample.
The stars of the area selected as candidate cluster members were
classified into photometric regions and their physical parameters determined,
following the algorithm described in Masana et al. (2007) and
Jordi et al. (1997). The algorithm uses
photometry and standard relations among colour indices
for each of the photometric regions of the HR diagram.
![]() |
Figure 9:
The histograms of the distance modulus, reddening and metallicity of the selected member stars of M 67 with H![]() |
Open with DEXTER |
Only 251 stars among the 776 candidate members have H measurements. Thus, the computation of physical parameters is only possible for that subset. The results are shown in Fig. 9.
Excluding peculiar stars and those with inconsistence among their
photometric indices and applying an average with a 2
clipping to that subset,
we found a mean reddening value of E(b-y) = 0.03
0.03
(corresponding to E(B-V) = 0.04) and a mean distance modulus of V0-MV = 9.4
0.4.
The distance modulus may be biased towards short values
due to the presence of multiple stars treated as if they
were single ones. A substantial fraction of binaries tend to have equal mass companions.
These binaries would be well-separated from the cluster main-sequence, so
they can be easily detected and removed. However, low mass ratio binaries would
remain, and might systematically reduce the estimated distance since they are
brighter and redder than single cluster members.
If we set a conservative limit for stars with lower
value than 9.4 in the distance modulus, we can see from the V vs. b-y diagram that these lower values mainly correspond to stars above the main sequence of the cluster
which are most probably multiple stars. Excluding those we can get
a value of V0-MV = 9.7
0.2 from the study of 66 stars that
can be confidently considered singles.
Metallicity is better calculated studying only the 81 F and G type stars in
our sample following Masana (1994).
We find a mean value of [Fe/H] = 0.01 0.14.
The publication by Clem et al. (2004) of empirically constrained colour-temperature relations in the Strömgren system makes possible an isochrone fitting to our results. The least
unsatisfactory fitting is found for the Pietrinferni et al. (2004,2006) tracks. Figure 10
shows the isochrones from Pietrinferni et al. (2004) for scaled solar models
of solar metallicity and ages of 4.0, 4.2, 4.4 and 4.6 Gyr for models with overshooting
and of 3.4, 3.6, 3.8 Gyr for canonical models. The adopted reddening and distance modulus are E(b-y) = 0.03 and V0 - MV = 9.7; from Pietrinferni et al. (2006) for -enhanced models
of solar metallicity and ages of 3.6, 3.8, 4.0 and 4.2 Gyr for models with overshooting
and of 3.0, 3.2, 3.4 Gyr for canonical models.
For this case, the adopted reddening and distance modulus are E(b-y) = 0.05
and V0 - MV = 9.7. Empty circles are candidate members.
Canonical models give lower ages.
-enhanced models reproduce better the behaviour of the red giants clump and red giant branch, but seem to behave worse in the lower giants branch and
subgiants branch where models with overshooting seem to give a better fit.
None of the models provides a good fit over all the areas of the colour-magnitude diagram.
We decided to adopt a compromise between these isochrones and give an estimation of the age of t = 4.2
0.2 Gyr (
= 9.62
0.02)
in agreement with previous estimates by other authors.
![]() |
Figure 10: Isochrones from Pietrinferni et al. (2004) for scaled solar models ( top row) of solar metallicity and ages of 4.0, 4.2, 4.4 and 4.6 Gyr for models with overshooting ( left panel) and ages of 3.4, 3.6, 3.8 Gyr for canonical models ( right panel). Isochrones from Pietrinferni et al. (2006) for alpha-enhanced models ( bottom row) of solar metallicity and ages of 3.6, 3.8, 4.0 and 4.2 Gyr for models with overshooting ( left panel) and ages of 3.0, 3.2, 3.4 Gyr for canonical models ( right panel). |
Open with DEXTER |
The high binary content of M 67 was already noted by Racine (1971):
"more than half of the M 67 main-sequence stars between G2 and K5 appear
to be unresolved multiple stars''.
Montgomery et al. (1993) studied the distribution of binaries and calculated a 22
of
equal-mass component binaries. Trying to account for the binaries with
low mass ratios when studying the distribution of multiple stars from a fiducial
main sequence, they found a ratio of 38
.
Due to the possible presence
of very low mass ratio binaries, this percentage is a lower limit.
Fan et al. (1996) reconsidered the mass ratios
and number of binaries in M 67, and using models leaving the mass ratio
to randomly vary from zero to one, they found that the
true binary fraction in this cluster depends critically on how to account for
the contribution of low mass-ratio binaries. From the models, around 50
of binaries seems a plausible scenario, with a binary mass ratio distribution
more consistent with being random than double-peaked.
![]() |
Figure 11: Multiple star systems (marked as crosses) and blue stragglers (diamonds) in a colour magnitude diagram in b-y ( left panel) and u-y ( right panel) of M 67. Empty circles are candidate member stars. Dots are non-member stars. A ZAMS is shown as reference. |
Open with DEXTER |
A list of the multiple star systems compiled by Sandquist (2004) and present in our photometry is shown in Table 10 with our membership segregation, including in the last two rows the recently discovered eclipsing systems by Sandquist (2006). Figure 11 shows the distribution of multiple star systems and blue stragglers in the colour-magnitude diagram.
The stars S757, S1036 (EV Cnc), S1282 (AH Cnc) and ES379 (ET Cnc), are W UMa contact binaries (where "S'' refersto Sanders's 1977 study and "ES'' refers to Eggen & Sandage's 1964 study). W UMa variables are one class of binary star that can produce blue stragglers after angular momentum loss causes the two stars to coalesce.
Many studies have been done in M 67 and its rich and varied population
of binaries and blue stragglers (BS). It has a super-blue straggler (S977), a single star
with a mass of
,
more than twice the cluster turn-off mass.
Gilliland et al. (1993) discovered two BSs with low amplitude
Scuti pulsations (S1280 and S1282) and evidence of longer period variations in other stars. Sandquist & Shetrone (2003, hereafter SSh03) discuss that S968, S1066 and S1263 could be long period variables.
Sandquist et al. (2003) and van den Berg et al. (2001) point out that S1082 is part of a triple system with two components BSs; one in a close binary and the other one, the brightest component, showing some
evidence of being a
Scuti star. The close binary is an RS CVn system
and also an X-ray source. Other BSs detected in X-rays, S997 and S1072, have wide eccentric
orbits whose nature is still not understood (van den Berg et al. 2004). Also with an eccentric orbit, S1284 is an X-ray source, too.
We consider BS candidates those stars that, according to our photometry, are located above and to the blue of the main sequence turnoff (Bailyn 1995; Stryker 1993), and, at the same time, have not been classified as non-members in our astrometric analysis. Those member stars that are above the main sequence turnoff but slightly to the red are usually also included in the BSs list. Our selection of BSs is then enlarged by the so-called yellow stragglers (YS), astrometric member stars located between the turnoff and the giant branch, excluding those known to be unresolved binaries with feasible standard components (e.g. van den Berg 2001; Portegies Zwart et al. 1997). This allows a re-evaluation of the membership for not all, but a good part of the stars formerly proposed as BS members of M 67. This helps to clean out the list of bona fide BSs in this cluster. We can perform a further check on the reliability of the selection, drawing the V vs. (u-y) diagram (see right panel of Fig. 11). Since red giants are faint in UV, the photometric blends, which mimic BSs in visible colour-magnitude diagrams are less problematic (see for example Sabbi et al. 2004). We find a total of 23 BS candidates, listed in Table 11. Two of them, S489 and S1466, had not been proposed as BSs before. S1466 was considered non-member by some authors (Sanders 1977; Girard et al. 1989), but it shows X-ray emision as detected by Chandra (van den Berg et al. 2004).
The work by Deng et al. (1999, hereafter DCLC99), based on the extensive photometry by Fan et al. (1996), proposes a list of 24 BSs. Of them, 17 are astrometric members contained in Table 11, and 2 others (S977 and S2226) appear in our BS selection but they are not covered by our astrometric study. BSs S1031 and S1440 are not included in DCLC99's list, but they had been also considered as BSs by Ahumada & Lapasset (1995) and SSh03.
The list by DCLC99 includes five stars not in Table 11. The reason is that one of them (S975) has an astrometry incompatible with being a member, two other (S1267 and S1434) are astrometric members but are not covered by our photometry, and the other two (S145 and S277) are neither in our astrometric nor photometric studies.
Other stars (S740, S821, S856 S1165, S1183) proposed as BSs by Ahumada & Lapasset (1995) are already considered normal turnoff member stars by SSh03 and DCLC99, a conclusion supported by our results. S1947, also proposed by Ahumada & Lapasset (1995), is not covered by our studies.
The BS population of M 67 is believed to be abnormally large
with respect to other open clusters. Ahumada & Lapasset (1995)
give a ratio of the number of BSs to that of main-sequence
stars within 2 mag below the turnoff,
,
of 30/200 for
this cluster, with a mean of 0.025 for clusters with 6.5
8.6.
However, we find in our sample a ratio of
= 23/289, similar to
the 24/286 found by DCLC99,
in closer agreement with the average found among other open clusters.
![]() |
Figure 12: Isochrones from Pietrinferni et al. (2004) from left to right: and 2.5, 1.8, 1.2, 1.0, 0.7, 0.5, 0.4 Gyr for M 67. |
Open with DEXTER |
We can derive photometric estimates of BS masses by comparing their
position in the colour-magnitude diagram with
the appropriate isochrones (see Fig. 12) from
the set of Pietrinferni et al. (2004) with the corresponding metallicities, reddening and
distance modulus. We have plotted
the standard isochrones from 2.5 Gyr (
= 1.46
), 1.8, 1.2,
1.0, 0.7, 0.5 and 0.4 Gyr (
= 2.70
). The reported
isochrones encompass the whole distribution of BSs, thus constraining the range
of masses covered
to 1.5
2.7
for M 67.
BSs with
1.5
,
i.e. formed by the merging of
low mass main-sequence stars, are still hidden in the main sequence.
There are some BSs that cannot lie on the main sequence of any reasonable
isochrone, but which can be well fit by the sub-giant branch sequences.
These stars, mentioned above,
have evolved to the thinning hydrogen burning shell phase and
are rapidly moving to the base of the RGB, i.e., they are YS,
and also some of these are candidates to be evolved
blue stragglers (E-BS), i.e, a BS in its helium-burning phase
(see Bellazzini et al. 2002, and references therein).
The presence of a striking clump of stars around
16,
0.4,
1.1, can be noticed in Fig. 7.
This accumulation of stars in the same area cannot be a coincidence and deserves
further consideration. We have dubbed this group Lola's Bunch, after
one of the authors of this paper. There are
approximately 60 stars in this clump on the V vs. b-y diagram,
where this accumulation reaches a density more than twice that
of the surrounding areas of the photometric space.
Furthermore, the spatial distribution of the stars
(Fig. 13) shows a lack of them
in the centre of the cluster. The
average distance from the centre of these stars is r = 19.94
7.85
.
![]() |
Figure 13:
Spatial distribution of clump stars (crosses) and members from our
astrometric study (empty circles) in M 67 area. Dots represent the stars
present in our CCD photometric sample: their spatial coverage does not reach
the lower right nor the upper left corners. Clump stars, selected exclusively
on a photometric basis, are disposed around the cluster centre, at an average
distance of 20
![]() |
Open with DEXTER |
We have looked for photometry of different sources to check the reality of the bunch, discarding any bizarre photometric effect in our data. We have seen that the only existing study covering an area wide enough to include the stars of this feature is the one by Fan et al. (1996) using the BATC system, that covers almost two square degrees.
In Fig. 14 we can see a colour-magnitude diagram
of the central 50
50
area of M 67 in two different colours:
mag3890 vs. 3890-6075 and 3890-9190. Although the contrast is lower than the factor two found in the Strömgren passbands, the bunch also appears as a significant accumulation in the diagrams,
with a local density equal to 1.5 times that of the surronding area in the
colour-magnitude diagram mag3890 vs. 3890-9190, and 1.75 times in the colour-magnitude
diagram mag3890 vs. 3890-6075. The bunch
is more noticeable when using the BATC filter 3890, quite similar to
Strömgren v. The lack of clump members at the central region of the cluster
is also observed in BATC data.
Stars placed in this photometric region of the colour-magnitude diagram
are found all over the two square degrees field covered by BATC data,
but we have noticed that the bunch gives a better constrast
when selecting the central 50
50
.
This is also approximately
the area covered by our Strömgren photometry, and coincides with the spatial area
over which most astrometric members are concentrated (Fig. 5).
![]() |
Figure 14:
The BATC colour-colour diagrams in the central 50 ![]() ![]() |
Open with DEXTER |
Considering the main observational facts, we see that either of the two possible options (clump stars being members or non-members) are interesting. If those stars were mainly cluster members, then their photometry could be explained if they were detached binaries composed of a white dwarf and a red dwarf, i.e: pre-cataclysmic variable stars (pre-CVs). Current knowledge of white dwarfs and pre-CV stars indicates that the combinations required to obtain the observed colours are physically feasible (Bergeron et al. 1995). Stassun et al. (2002) have found that five stars in this photometric area of the M 67 colour-magnitude diagram are variables of kinds compatible with the pre-CV status, but they failed to perceive the relevance of this concentration of stars in their data. Confirming the existence of a significant population of pre-CV stars in such a well studied open cluster (and thus, of known distance, age and metallicity) would have very deep implications for the understanding of the nature and evolution of these binary systems. First, their clumping in a region of the colour-magnitude diagram, and not in a locus following the trend that would be expected from the white dwarf cooling sequence, would impose tight restrictions on the evolutionary history of the precursor binary systems. Second, their distribution in the outer region of the cluster could only be explained in terms of mass segregation, implying that the pairs have had enough time to experience this dynamical process since the mass-loss that made them as light as they would be now.
The other option, having a clump made up from field stars, is no less challenging. The relatively narrow intervals of apparent magnitude and colour means that all these stars have similar spectral type and, thus, similar absolute magnitude and distance. In this case, we would have to explain how an accumulation of field G stars could be formed, at a definite distance: we would have a curtain or stream of G stars placed twice as far as M 67. Furthermore, their spatial distribution (that seemingly avoids the centre of M 67) would require some explanation.
The scarcity of the data available make any further consideration peculation: we have already started spectroscopic observations with the aim of clarifying the nature of this peculiar group of stars.
In this paper we give a catalogue of accurate Strömgren
and J2000 coordinates for 1843 stars in an area of 50
50
around M 67.
We give a selection of probable members of M 67 combining this
photometric study with astrometric analysis using parametric as well
as non-parametric approaches. A better determination of the cluster
physical parameters based on our accurate photometry gives:
E(b-y) = 0.03
0.03, [Fe/H] = 0.01
0.14,
a distance modulus of V0-MV = 9.7
0.2 and an age of
4.2
0.2 Gyr (
= 9.6
0.1). The values are coherent with
previous studies (see Chen et al. 2003, and references therein).
An anomalous bunch of more than 60 stars in the colour-magnitude diagram around V = 16, (b-y) = 0.4 is tentatively interpreted as composed of pre-CVs belonging to the cluster, or as a stream of G stars placed twice as far as the cluster, two alternatives to be tested with further observations.
Acknowledgements
We thank Simon Hodgkin and Mike Irwin for their inestimable help with the reduction of the images taken at the WFC-INT. L.B-N, also wants to thank Gerry Gilmore and Floor van Leeuwen for their continuous help and valuable comments, as well as all the people at the IoA (Cambridge) for a very pleasant stay. L.B-N. gratefully acknowledges financial support from EARA Marie Curie Training Site (EASTARGAL) during her stay at IoA. Based on observations made with the INT and JKT telescopes operated on the island of La Palma by the RGO in the Spanish Observatorio del Roque de Los Muchachos of the Instituto de Astrofísica de Canarias, and with the 1.52 m telescope of the Observatorio Astronómico Nacional (OAN) and the 1.23 m telescope at the German-Spanish Astronomical Center, Calar Alto, operated jointly by Max-Planck Institut für Astronomie and Instituto de Astrofísica de Andalucia (CSIC). This study was also partially supported by the contracts No. AYA2003-07736, AYA2006-15623-C02-02, with MCYT. This research has made use of Aladin, developed by CDS, Strasbourg, France. This research has made use of the WEBDA database, operated at the Institute for Astronomy of the University of Vienna.
Table 8:
The cross-identification of stars considered as candidate members by us,
with the radial velocities sudies by Mathieu et al. (1990,1986) and
the comparison of those with the parametric ()
and
non-parametric (
)
probabilities.
See text for explanation of the agreement index
.
Table 9: The same as Table 8, for stars considered as non-members by us.
Table 10: The cross-identification of known multiple star systems and their astrometric segregation from this work (last column). First column is our identification number (Table 3), second column that from BDA and third that from Sanders (1977).
Table 11: Our BS candidates in M 67, and their photometry. First column is our identification number (Table 3), second column that from BDA, and third that from Sanders (1977). Last column gives the astrometric membership (void values for stars without astrometric membership determination).