A&A 416, 125-136 (2004)
DOI: 10.1051/0004-6361:20034238

Proper motion surveys of the young open clusters Alpha Persei and the Pleiades[*]

N. R. Deacon - N. C. Hambly

Institute for Astronomy, University of Edinburgh, Blackford Hill, Edinburgh EH9 3HJ, UK

Received 28 August 2003 / Accepted 28 October 2003

Abstract
In this paper we present surveys of two open clusters using photometry and accurate astrometry from the SuperCOSMOS microdensitometer. These use plates taken by the Palomar Oschin Schmidt Telescope giving a wide field ($5^{\circ}$ from the cluster centre in both cases), accurate positions and a long time baseline for the proper motions. Distribution functions are fitted to proper motion vector point diagrams yielding formal membership probabilities. Luminosity and mass functions are then produced along with a catalogue of high probability members. Background star contamination limited the depth of the study of Alpha Per to R=18. Due to this the mass function found for this cluster could only be fitted with a power law ( $\xi(m) = m^{-\alpha}$) with $\alpha=0.86^{+0.14}_{-0.19}$. However with the better seperation of the Pleiades' cluster proper motion from the field population results were obtained down to R=21. As the mass function produced for this cluster extends to lower masses it is possible to see the gradient becoming increasingly shallow. This mass function is well fitted by a log normal distribution.

Key words: astrometry - stars: low mass, brown dwarfs - Galaxy: open clusters and associations: individual: Alpha Per, Pleiades

1 Introduction

The young open cluster Alpha Per lies at a distance of 183 pc ( $(M-m)_{\rm o}=6.31$, van Leeuwen 1999) with an age of 90 Myr (Stauffer et al. 1999). We assume in this paper that AR=0.23 (O'Dell et al. 1994). It has been reasonably well studied in the past using various techniques. Prosser's catalogue (Prosser 1993) contains a comprehensive list of candidate member stars identified by photometry, spectroscopy and radial velocity measurements. The faintest of these has $V \approx 21.8$( $M_{V} \approx 15.5$). These candidate stars have been supplemented by those identified in similar studies by Zapatero Osorio et al. (1996), Stauffer et al. (1999) and Rebolo et al. (1992). In addition Prosser et al. (1998) identified candidate stars which were optical counterparts to X-ray sources within the cluster.

An early proper motion study by Heckmann et al. (1956) identified stars with proper motions consistent with cluster membership. However this study only went down to $V \approx 12$ and did not calculate formal membership probabilities. Prosser (1992) also used stellar proper motions as an initial selection technique. Although this study was deeper than the earlier one (V=18.8) it also did not calculate formal membership probabilities.

Most recently Barrado y Navascués et al. (2001) produced a deep CCD survey of the cluster. This identified several candidate low mass stars and Brown Dwarfs down to $I_{\rm c}=22$. They then used this data to produce a mass function. The mass function was fitted with a power law whose slope was found to be $\alpha = 0.59$.

The Pleiades is the best studied open cluster, here we assume $(M-m)_{\rm o}=5.4$ (van Leeuwen 1999) and AR=0.15 (O'Dell et al. 1994). In recent years several deep CCD surveys of the cluster such as Dobbie et al. (2002) and Moraux et al. (2003) have produced candidates down to $0.03~ M_{\odot}$. These have been complimented by proper motion studies such as those of Hambly et al. (1999) and Moraux et al. (2002). However the only wide field, high precision proper motion survey is that of Hambly et al. (1991). Adams et al. (2001) used 2MASS data along with proper motions to produce a wide field study. This contained many candidate stars outside the tidal radius of the cluster showing background star contamination is a problem for many techniques.

  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{applot.ps}
\end{figure} Figure 1: A plot showing the sky coverage of the Alpha Per study. The solid line represent the later POSSII fields while the dotted lines are the earlier POSSI fields. The dashed circles marks out a five degree radius from the cluster centre (marked by an $\times $).
Open with DEXTER

2 Observational data and reduction

2.1 Observations

The SuperCOSMOS facility at the Royal Observatory Edinburgh has (using plates from the United Kingdom Schmidt Telescope) produced complete southern sky surveys in BJ, R and I with an additional second epoch R survey. These surveys are now publicly available (Hambly et al. 2001). The scanning program has now moved on to the northern hemisphere, using film and glass copies of plates taken by the Oschin Schmidt Telescope on Mount Palomar, California. These data will soon be publicly available. The declination of Alpha Per ( $\delta \approx 49^{\circ}$) and the Pleiades ( $\delta \approx 24^{\circ}$) means this prerelease northern hemisphere data must be used. Details of these plates are given in Table 1 for Alpha Per and Table 2 for the Pleiades. Unfortunately due to saturation caused by the bright core stars of the cluster and a satellite track, 3.1% (2.3 sq. deg) of the Pleiades survey area was unusable. The region of the CO cloud near Merope was used in the proper motion survey, it is assumed that the increased reddening in this region will not affect the results. The areas covered by both surveys are shown in Figs. 1 and 2.

  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{plplot.ps}
\end{figure} Figure 2: A plot showing the sky coverage of the Pleiades study. The solid line represent the later POSSII fields while the dotted lines are the earlier POSSI fields. The dashed circles marks out a five degree radius from the cluster centre (marked by an $\times $).
Open with DEXTER

2.2 Sample selection

The sample of stars used was selected to maximise completeness and reliability (e.g. by minimising astrometric errors). Highly elliptical objects, non-stellar objects and poor quality images near bright stars were all removed. Deblended images were included. This is because Alpha Per (lying at a galactic latitude of  $-7^{\circ}$) is in a crowded field meaning many stellar images will be merged with others on the plate. Unfortunately the deblending algorithm is rather crude and this can lead to increased astrometric and photometric errors. It is necessary to calculate the astrometric errors caused by the deblending alogorithm before deciding whether to include deblended objects.

Table 1: Schmidt plates used in the Alpha Per study.

Table 2: Schmidt plates used in the Pleiades study.

2.3 Errors in the proper motions

In the Alpha Per data set there were several overlap regions between the plates used in this study. Hence objects in this region would have two different measurements of their position. These were used to calculate the rms error on the proper motions. The errors were calculated in both axes ($\xi$ and $\eta $) and for objects that had been deblended and those that had not. Figure 3 shows how these errors vary with magnitude. Notice the higher errors at the bright end and the faint end. The overlap regions in the Pleiades were too small to be used for a similar calculation. Hence a small selection of faint blue stars was selected. These are expected to have negligible proper motions, so any measured proper motion would be the result of random errors. This yielded an error estimates of $\sigma_{x}=6.2$ mas/yr and $\sigma_{y}=6.3$ mas/yr. It should be noted that these errors are high because the stars used are faint, Fig. 3 shows that for fainter stars errors are higher. When a study was done using candidate stars from other studies for Alpha Per it was found that 172 out of 323 were merged with other images. As it has been shown that there is no significant increase in the errors for deblended stars it was decided to include these to maximise completeness.

2.4 Photometric calibration

The photometric calibration used by the SuperCOSMOS survey is subject to large systematic errors especially at the bright end. This situation can be improved when large numbers of photoelectric measurements are available, as is typically the case in open clusters, using a process of photometric recalibration. To do this the photographic magnitudes from the sample were compared with their photoelectric magnitudes. The photoelectric data for recalibrating the Alpha Per data set was provided by John Stauffer and with others being taken from Barrado y Navascués et al. (2003). Figure 4 shows a plot of photoelectric magnitude vs. photographic, the line shown is a fourth degree polynomial fit to the data. The best fit line shown in the graph was then used to recalibrate the photographic magnitudes. Unfortunately the data from Barrado y Navascués' only includes stars which fall in one field (200). Hence the recalibration was carried out on this plate and propagated through to the other plates. The scatter from the fitted polynomial was used to estimate the photometric errors to be, $\sigma_{R} = 0.20$ and $\sigma_{I} = 0.20$. A similar process was carried out for the Pleiades using data from Jameson & Skillen (1989) and Stauffer (1984). Figure 5 shows a plot of photoelectric magnitude vs. photographic, again the line shown is a fourth degree polynomial fit to the data. The photometric errors were found to be, $\sigma_{R} = 0.17$ and $\sigma_{I} = 0.21$.
  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{errplot.ps}
\end{figure} Figure 3: A plot showing how the errors in the proper motion vary with magnitude. The larger errors for bright stars in the $\eta $ axis are an artifact of the the plate measuring process due to saturated images.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{APphot.ps}
\end{figure} Figure 4: Comparing the photographic plate magnitudes with the measured photoelectric magnitudes for the Alpha Per data set. The line shown is a quartic fit to the data. The photoelectric data has been naturalised to the photographic passband system.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{PLphot.ps}
\end{figure} Figure 5: Comparing the photographic plate magnitudes with the measured photoelectric magnitudes for the Pleiades data set. The line shown is a quartic fit to the data. Again the photoelectric data has been naturalised to the photographic passband system.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{aa0238f6.eps}
\end{figure} Figure 6: A colour magnitude diagram for stars in the initial sample for Alpha Per. Only stars above the line shown were included in the proper motion survey. For symbols see text.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{aa0238f7.eps}
\end{figure} Figure 7: A colour magnitude diagram for stars in the initial sample for the Pleiades. Only stars above the line shown were included in the proper motion survey. For symbols see text.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{prob.ps}
\end{figure} Figure 8: Proper motion vector point diagrams and probability histograms for each magnitude interval in the Pleiades study. In each VPD the cluster is the group in the bottom right hand corner, separate from the field stars around the origin.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{APprob.ps}
\end{figure} Figure 9: Proper motion vector point diagrams and probability histograms for each magnitude interval in the Alpha Per study. In each VPD the cluster is the group in the bottom right hand corner, separate from the field stars around the origin.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{comp1.ps}
\end{figure} Figure 10: An estimate of the incompleteness of the study of Alpha Per. The line shown is a least squares fit.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{compPL.ps}
\end{figure} Figure 11: An estimate of the incompleteness of the study of the Pleiades. The line shown is a least squares fit.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8cm,clip]{cLF.ps}
\end{figure} Figure 12: The Luminosity Function derived for Alpha Per.
Open with DEXTER

2.5 Colour selection

After photometric recalibration, data cuts could be made on the basis of colour and brightness. Firstly it was decided that as stars brighter than R=10 produce high centroiding errors these would be excluded. In the Alpha Per study the low separation between the cluster and the field in the Proper Motion Vector Point Diagram produced a large amount of background star contamination at the faint end. Hence stars fainter than R=18 were excluded. With it's better separation between the cluster and field such a cut was not necessary for the Pleiades data set. To reduce the amount of background star contamination colour selections were made. Figures 6 and 7 show colour magnitude diagrams for each sample. The dots represent stars in the sample while the star symbols are records from the sample which have been identified as candidates from both our study and previous studies, the triangles are stars thought to be candidiates by previous studies but not by our study. It can clearly be seen that the candidate stars from other surveys form a main sequence-like line distinct from the large mass of background stars. Only stars above the line shown were used in this study. In the case of the Pleiades Fig. 7 shows the colour magnitude diagram. The large bulk of background stars are easily identified and removed. The line shown is again the colour cut used. Both lines used in the colour cuts were chosen by inspection As astrometric errors vary with magnitude the data was divided into the following magnitude ranges, for Alpha Per 10<R<12, 12<R<14, 14<R<16 and 16<R<18 and for the Pleiades, 10<R<12, 12<R<14, 14<R<16, 16<R<17.5 and R>17.5. Each of these was treated independently in the fitting process.

2.6 The fitting process

The membership probabilities were calculated by a process similar to that outlined by Sanders (1971). The full mathematical detail is outlined in Appendix A. To ensure that the data fitting process is robust a series of twenty sets of simulated data were produced and run through the fitting program. The data sets were generated by taking uniform random number distributions. These were converted into Gaussian and exponential distributions with the correct scale lengths, means and variances to model the expected distributions for the field and cluster stars. Each data set included 350 cluster stars and 350 field stars, both typical numbers for the real data. Table 3 shows the results obtained. It is apparent that there is no significant offset in any of the calculated parameter values.
  \begin{figure}
\par\includegraphics[width=8cm,clip]{PLcLF.ps}
\end{figure} Figure 13: The Luminosity Function derived for the Pleiades. The striking feature is the peak at R=16-17.
Open with DEXTER

Table 3: The results of running 20 sets of simulated data through the parameter fitting program. With the exception of f all parameters have units of milliarcseconds per year.

3 Results

The proper motion vector point diagrams for each of the magnitude ranges for both clusters are shown in Fig. 8 (the Pleiades) and Fig. 9 (Alpha Per) along with probability histograms. In the case of Alpha Per the cluster is especially apparent in the vector point diagram 16<R<18 and consequently it also has the highest number of high probability cluster members. As the Pleiades study goes far deeper we can see the cluster becoming more apparent with increasing magnitude before becoming engulfed by the mass of field stars in the lowest magnitude interval. Accordingly the number of high probability members increases with increasing magnitude, the high number of field stars in the lowest magnitude interval means there are no stars with membership probabilities >$90\%$. A full list of stars with calculated membership probabilities greater than $60\%$ is available for both clusters from CDS in Strasbourg. Example tables are shown in Appendix B.

4 Analysis

4.1 Luminosity functions

As Schmidt plates have poor detection rates close to the plate limit it was necessary to carry out a completeness estimate. This was done by measuring how the number of stars in the full stellar sample increased with Rmagnitude. For a uniform distribution of stars the number in each sample should rise exponentially with increasing magnitude with any dropoff being caused by incompleteness. An estimate of this incompleteness can be found by taking the logarithm of the number of stars in each magnitude interval, fitting a best fit line up to the point where the dropoff begins and then using the deficit to estimate the incompleteness. Figure 10 shows such a fit for the Alpha Per data set and Fig. 11 that for the Pleiades sample.

The luminosity functions were found by summing the membership probabilities of stars in each one magnitude wide interval. Each interval was then corrected for incompleteness. Figure 12 shows the Luminosity Function produced for Alpha Per. Binarity was not taken into account when producing this. The luminosity function is smooth and rises with increasing magnitude. The Pleiades luminostiy function (Fig. 13) by contrast has distinct features, the most obvious of which is the peak at 16<R<17. The errors shown are Poisson errors which also take into account errors in the number of backround stars.

  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{mlplot.ps}
\end{figure} Figure 14: The mass-luminosity relation used for Alpha Per. The line is a sixth degree polynomial fit.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=6.9cm,clip]{APmlplot.ps}
\end{figure} Figure 15: The mass-luminosity relation used for the Pleiades. The line is a sixth degree polynomial fit.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=6.9cm,clip]{mfall.ps}
\end{figure} Figure 16: The derived cluster Mass Function for the Pleiades. Note the increasingly shallow gradient towards lower masses. The line shown is a log normal fit to the data. The solid circles are data points from this survey while the crosses are from Hambly et al. (1999) and references therein.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{APmfall.ps}
\end{figure} Figure 17: The derived cluster Mass Function for Alpha Per. The line show is a power law fit.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{mlerr.ps}
\end{figure} Figure 18: A plot showing possible systematic errors in the mass luminosity relation used for the Pleiades.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{APcrosscat.ps}
\end{figure} Figure 19: A histogram of the membership probabilities calculated in this study of Alpha Per for stars previously identified as candidate members.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{crosscat.ps}
\end{figure} Figure 20: A histogram of the membership probabilities calculated in this study of the Pleiades for stars previously identified as candidate members.
Open with DEXTER

4.2 Mass function

To convert the luminosty function to a mass function a mass luminosity relationship is required. This is provided by the models of Barraffe et al. (1998). Figure 14 shows the mass luminosity relation taken from these models for a 90 Myr old population such as Alpha Per. The line shown is a sixth degree polynomial least squares fit to the data which is used to find the mass of a star of a particular luminosity. The mass function is found from the luminosity function, $\Phi(M)$. In some particular magnitude interval of width ${\rm d}M$ there will be a number of stars d $N=\Phi(M){\rm d}M$. In the corresponding mass interval the same number dN is given by d $N=\xi(m){\rm d}m$, where dm is the width of the interval and $\xi(m)$ is the mass function. This can be calculated by taking the number of stars in a particular luminosity interval and dividing by the width of that interval in mass. Again binarity was not taken into account. The mass function for Alpha Per is shown in Fig. 17. The mass function was then fitted by a power law fit of the form,

\begin{displaymath}\xi(m) \propto m^{-\alpha}.
\end{displaymath} (1)

It was found that $\alpha=0.86^{+0.14}_{-0.19}$. As mentioned earlier Barrado y Navascués' et al. (2002) derived a mass function for Alpha Per at lower masses (down to $M=0.035~M_{\odot}$). They also fitted a power law function of the form, finding that $\alpha = 0.59$. It therefore appears that the mass function flattens off significantly towards lower masses.

A mass function was also produced for the Pleiades. This time the mass luminosity relation was that for a 120 Myr old population and is shown in Fig. 15. Again this was taken from Barraffe et al. (1998). The mass function is shown in Fig. 16, the data points from this survey are supplemented by those from Hambly et al. (1999) and references therein. It has been suggested by Adams & Fatuzzo (1996) that the initial stellar mass function can be approximated by a log normal function as shown below,

\begin{displaymath}\log_{10}\xi(m) = a_{0} + a_{1}\log_{10}m + a_{2}(\log_{10}m)^{2}.
\end{displaymath} (2)

It is clear that a power law fit is not appropriate for the mass function derived here; a log normal was fitted to the Mass function. The parameters are a0=2.213, a1=-2.069 and a2=-0.745. Hambly et al. (1999) produced a similar log normal fit to the mass function with parameters, a0=2.341, a1=-2.313 and a2=-1.191. These are in generally good agreement with the discrepancy in a2 probably due to the poorer constaint of the Hambly study at the faint end. The referee has pointed out the possibilty of systematic errors in the mass-luminosity relationship from the Baraffe et al. (1998) models. We checked for such errors by taking an empirical BCR vs. $T_{\rm eff}$relation derived using the same data and techniques as described in Hambly et al. (1999) and using it to estimate the R magnitude from  $M_{\rm bol}$ for each of the points taken from Baraffe's model. The mass of these points was then compared to that given by the mass-luminosity relation for the recalculated R. The results of this are shown in Fig. 18.

4.3 Comparison with other studies

As mention earlier there are several other methods to determine cluster membership. Examining the membership probabilities of stars from other studies can help both to confirm individual star's cluster membership and to examine the relevance of both study's results. It would be worrying if candidates from previous studies did not have, by and large, high membership probabilities. The candidate stars were paired with the appropriate records in the catalogue. The maximum pairing radius was set to 6 arcseconds. A table of those candidate stars for both clusters with calculated membership probabilities is availible from CDS in Strasbourg, examples are provided in Appendix B. Figure 19 shows a histogram of the membership probabilities for these candidate stars for Alpha Per while Fig. 20 shows the same graph for the Pleiades study. It is encouraging that both show that most stars thought to be members in previous studies have high membership probabilities and that the large number of probable nonmembers shown in the membership probability histograms (see Figs. 8 and 9) do not appear here. Both Figs. 19 and 20 show a sharp increase in the number of candidate stars from previous studies at $p=60\%$. For this reason it was decided to include only stars with a membership probability greater than this in the catalogues of probable members for both clusters (available at the CDS). It should be noted that several HHJ objects (see Hambly et al. 1991) do not have calculated membership probabilities. This is because these objects appeared highly elliptical on the first epoch plate SuperCOSMOS measurements and hence were not included in the present proper motion survey.

5 Conclusion

We have presented the first proper motion survey of the cluster Alpha Persei to calculate formal membership probabilities. A Luminosity Function for the range R=10 to R=18 has been produced and a Mass Fuction derived from this. The Mass Function was fitted with a power law distribution with $\alpha=0.86^{+0.14}_{-0.19}$. A catalogue of 339 high probability members has been created. This has been cross-correlated with previous studies of the cluster showing that a large number of stars previously thought to be cluster members have high calculated membership probabilities. A similar study has also been undertaken on the Pleiades and while not the first of it's type it has produced several low mass star candidates and yielded a better constrained Mass Function at the faint end.

Acknowledgements
The authors would like to thank John Stauffer for supplying data on candidate member stars from other surveys and for a prompt and thorough referee's report.

   
Appendix A: Mathematics

In order to calculate membership probabilities, distributions were fitted to both the cluster stars and the field stars in the vector point diagram. The cluster stars were fitted with a circularly symmetric Gaussian distribution as in Sanders (1971) however the field stars were fitted with a decaying exponential in the direction of cluster proper motion and a Gaussian perpendicular to this direction as in Hambly et al. (1995). In order to fit a distribution for the field stars the vector point diagram was rotated so that the cluster lay on the y axis. The total distribution function $\Phi$ is given by Eq. (A.1) where  $\Phi_{\rm f}$ is the field star distribution, $\Phi_{\rm c}$ is the cluster star distribution and f is the fraction of stars which are field stars. $\Phi_{\rm f}$ is given in (A.2) and  $\Phi_{\rm c}$ is given in (A.3). $c_{\rm o}$ is the normalisation for the exponential between the limits $\mu_{1}$ and $\mu_{2}$and is given by (A.5).

 \begin{displaymath}\Phi = f \Phi_{\rm f} + (1-f) \Phi_{\rm c}
\end{displaymath} (A.1)


 \begin{displaymath}\Phi_{\rm f} = \frac{c_{\rm o}}{\sqrt{2 \pi} \Sigma_{x}} \exp...
...{x\rm f})^{2}}{2
\Sigma_{x}^{2}} - \frac{\mu_{y}}{\tau}\right)
\end{displaymath} (A.2)


 \begin{displaymath}\Phi_{\rm c} = \frac{1}{2 \pi \sigma^{2}} \exp \left(-\frac{(...
...\rm c})^{2}+(\mu_{y} - \mu_{y\rm c})^{2}}{2 \sigma^{2}}\right)
\end{displaymath} (A.3)


 \begin{displaymath}c_{\rm o} = \frac{1}{\tau\left({\rm e}^{-\frac{\mu_{2}}{\tau}} - {\rm e}^{-\frac{\mu_{1}}{\tau}}\right)}
\end{displaymath} (A.4)


 \begin{displaymath}c_{\rm o} \int_{\mu_{1}}^{\mu_{2}}{\rm e}^{-\mu_{y}/\tau} {\rm d}\mu_{y}=1.
\end{displaymath} (A.5)

Using the maximum likelihood method

 \begin{displaymath}\sum_{i} \frac{\delta \ln \Phi_{i}}{\delta \Theta} =0
\end{displaymath} (A.6)

(where $\Theta$ is some parameter) the following set of nonlinear equations are found.

 \begin{displaymath}f: \sum_{i} \frac{\Phi_{\rm f} - \Phi_{\rm c}}{\Phi}=0
\end{displaymath} (A.7)


 \begin{displaymath}\sigma: \sum_{i} \frac{\Phi_{\rm c}}{\Phi}
\left (\frac{(\mu_...
...m c})^{2} + (\mu_{y} - \mu_{y\rm c})^{2}}{\sigma^2}-2\right)=0
\end{displaymath} (A.8)


 \begin{displaymath}\Sigma_{x}: \sum_{i} \frac{\Phi_{\rm f}}{\Phi} \left( \frac{(\mu_{x} -
\mu_{x\rm f})^{2}}{\Sigma_{x}^{2}} -1\right)=0
\end{displaymath} (A.9)


 \begin{displaymath}\mu_{x\rm f} : \sum_{i} \frac{\Phi_{\rm f}}{\Phi} (\mu_{x} - \mu_{x\rm f})=0
\end{displaymath} (A.10)


 \begin{displaymath}\mu_{x\rm c} : \sum_{i} \frac{\Phi_{\rm c}}{\Phi} (\mu_{x} - \mu_{x\rm c})=0
\end{displaymath} (A.11)


 \begin{displaymath}\mu_{y\rm c} : \sum_{i} \frac{\Phi_{\rm c}}{\Phi} (\mu_{y} - \mu_{y\rm c})=0
\end{displaymath} (A.12)


 \begin{displaymath}\tau : \sum_{i} \frac{\Phi_{\rm f}}{\Phi}\left( \frac{\mu_{y}...
...mu_1}{\tau}} - \mu_{2}{\rm e}^{-\frac{\mu_2}{\tau}})\right)=0.
\end{displaymath} (A.13)

These are solved by a simple bisection algorithm. Applying the calculated values of the parameters, the probability of the ith star being a cluster member is,

\begin{displaymath}p_{i}=\frac{(1-f)\Phi_{c i}}{\Phi_{i}}\cdot
\end{displaymath} (A.14)

The fitted values for the parameters for each magnitude interval are shown in Table A.1 for Alpha Per and Table A.2 for the Pleiades.

Table A.1: The calculated parameters for each magnitude interval for Alpha Per. With the exception of f all parameters have units of milliarcsenonds per year.

Table A.2: The calculated parameters for each magnitude interval for the Pleiades. With the exception of f all parameters have units of milliarcseconds per year.

   
Appendix B: Example tables

Example tables of the catalogue produced. All coordinates are J2000 and all magnitudes are recalibrated photographic magnitudes in the natural photographic (R59F, IIV-N) system. Some membership probabilities may be underestimated at the bright end due to the low separation of the main sequence from the background stars at such magnitudes in the colour magnitude diagram.

Table B.1: An example table of the catalogue of high probability members of Alpha Per.

Table B.2: An example table of the catalogue of candidate members of Alpha Per from other studies with calculated membership probabilities. Several of these stars have alternative names, for comprehensive cross-identifications, see the Stauffer and Prosser Catalogue of open cluster data (Stauffer, private communication).

Table B.3: An example table of the catalogue of high probability members of the Pleiades.

Table B.4: An example table of the catalogue of candidate members of the Pleiades from other studies with calculated membership probabilities. Several of these stars have alternative names, see the Stauffer and Prosser Catalogue of open cluster data (Stauffer, private communication).

References



Copyright ESO 2004