Issue 
A&A
Volume 598, February 2017



Article Number  A48  
Number of page(s)  22  
Section  Galactic structure, stellar clusters and populations  
DOI  https://doi.org/10.1051/00046361/201629239  
Published online  31 January 2017 
A revised moving cluster distance to the Pleiades open cluster^{⋆}
^{1} Univ. Grenoble Alpes, IPAG, 38000 Grenoble, France
email: galli@astro.iag.usp.br
^{2} CNRS, IPAG, 38000 Grenoble, France
^{3} Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Universidade de São Paulo, Rua do Matão, 1226, Cidade Universitária, 05508900, São Paulo – SP, Brazil
^{4} Centro de Astrobiología, INTACSIC, PO Box 78, ESAC Campus, 28691 Villanueva de la Cañada, Madrid, Spain
^{5} Dpto. de Inteligencia Artificial, UNED, Juan del Rosal, 16, 28040 Madrid, Spain
Received: 4 July 2016
Accepted: 28 September 2016
Context. The distance to the Pleiades open cluster has been extensively debated in the literature over several decades. Although different methods point to a discrepancy in the trigonometric parallaxes produced by the Hipparcos mission, the number of individual stars with known distances is still small compared to the number of cluster members to help solve this problem.
Aims. We provide a new distance estimate for the Pleiades based on the moving cluster method, which will be useful to further discuss the socalled Pleiades distance controversy and compare it with the very precise parallaxes from the Gaia space mission.
Methods. We apply a refurbished implementation of the convergent point search method to an updated census of Pleiades stars to calculate the convergent point position of the cluster from stellar proper motions. Then, we derive individual parallaxes for 64 cluster members using radial velocities compiled from the literature, and approximate parallaxes for another 1146 stars based on the spatial velocity of the cluster. This represents the largest sample of Pleiades stars with individual distances to date.
Results. The parallaxes derived in this work are in good agreement with previous results obtained in different studies (excluding Hipparcos) for individual stars in the cluster. We report a mean parallax of 7.44 ± 0.08 mas and distance of 134.4^{+2.9}_{2.8} pc that is consistent with the weighted mean of 135.0 ± 0.6 pc obtained from the nonHipparcos results in the literature.
Conclusions. Our result for the distance to the Pleiades open cluster is not consistent with the Hipparcos catalog, but favors the recent and more precise distance determination of 136.2 ± 1.2 pc obtained from Very Long Baseline Interferometry observations. It is also in good agreement with the mean distance of 133 ± 5 pc obtained from the first trigonometric parallaxes delivered by the Gaia satellite for the brightest cluster members in common with our sample.
Key words: stars: distances / open clusters and associations: individual: Pleiades
Full Table B.2 is only available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/598/A48
© ESO, 2017
1. Introduction
The problem of distance determination has always played a central role in astronomy. For example, distances are necessary to determine the dimensions of celestial objects, the physical properties of stars, and their true motion (i.e., spatial velocity). In this context, determining the distance to open clusters is also vitally important when calibrating the astronomical distance scale. The stars associated with open clusters are assumed to have similar properties (age, distance, kinematics, and chemical composition), which makes them ideal targets to refine stellar evolution models.
As the closest open cluster to the Sun in terms of age and richness, the Pleiades is an important cornerstone to many studies related to star formation from physical models to observational properties of young stars. Given the crucial role of the Pleiades for calibration purposes, one would expect its distance to be well established. However, there is still a current debate in the literature regarding the distance to the cluster.
The Pleiades distance controversy began when the first studies in the preHipparcos era using the main sequence (MS) and/or isochrone fitting methods (Nicolet 1981; Giannuzzi 1995) delivered a distance estimate for the cluster that roughly exceeds by 10−15% the distance determination of 118.3 ± 3.5 pc (van Leeuwen 1999) that was obtained from the trigonometric parallaxes of the Hipparcos catalog (ESA 1997). Later studies using similar methods but different samples of stars (Pinsonneault et al. 1998; Stello & Nissen 2001; Percival et al. 2005; An et al. 2007) confirmed the previous results for the distance of the cluster, and the new distance estimate of 120.2 ± 1.9 pc (van Leeuwen 2009) that was based on the recalibrated parallaxes of 53 stars from the new version of the Hipparcos catalog (van Leeuwen 2007) did not change the situation.
In the meantime, two independent studies returned trigonometric parallaxes for a few Pleiades stars. Gatewood et al. (2000) measured the parallax of seven stars in the cluster yielding a weighted mean distance of 130.9 ± 7.4 pc, and Soderblom et al. (2005) find a distance of 134.6 ± 3.1 pc based on trigonometric parallaxes for another three stars.
Distance determination obtained from orbital modeling of eclipsing binaries is another possibility for estimating the distance to the Pleiades, and the results vary from 132.0 ± 5.0 pc (David et al. 2016) to 139.1 ± 3.5 pc (Southworth et al. 2005). Unfortunately, this method is restricted to a few systems that have been discovered to date and that can be used in this regard (Zwahlen et al. 2004; Munari et al. 2004; Groenewegen et al. 2007).
Recently, Melis et al. (2014) measured the trigonometric parallaxes of four stars (or stellar systems) in the Pleiades based on multiepoch Very Long Baseline Interferometry (VLBI) observations. The weighted mean distance of 136.2 ± 1.2 pc obtained in their work represents the most precise distance determination for the cluster to date, and the parallaxes for individual stars have a precision that is better than 1%. Indeed, the VLBI trigonometric parallax method is the most precise and accurate technique for measuring distances nowadays, and it has already proved to deliver good results in nearby associations of young stars (see e.g., Loinard et al. 2007, 2008).
In a recent paper, Mädler et al. (2016) use a different approach to derive the distance of 15 stars in the cluster based on stellar twins. This method assumes that the difference in brightness between two stars with different locations in the sky, but identical physical properties (i.e., stellar twins), results directly from the difference in their distances. Thus, by knowing the distance of one star from other methods, they calculate the distance to its stellar twin. Doing so, they found an average distance of 134.8 ± 1.7 pc for the Pleiades open cluster that clearly supports the nonHipparcos results.
The moving cluster method is another alternative to compute distances of comoving groups of stars. The method takes the observables proper motions and radial velocities into account to compute the distance to individual cluster members (see e.g., Galli et al. 2012). Although a variant of this methodology has already been used in previous studies to investigate the distance of the Pleiades, the admittedly large error bars in the distances results of 130.6 ± 10.8 pc (Narayanan & Gould 1999) and 125.9 ± 7.5 pc (Röser & Schilbach 2013) make it difficult to resolve the current dispute between the Hipparcos and VLBI distance determinations of the cluster.
In this context, we decided to revisit the Pleiades based on the moving cluster method by using a refurbished implementation of the method, an updated census of cluster members and more precise data to provide a revised distance determination to the cluster with this technique. The results obtained in this work are independent of previous distance determinations of the cluster in the literature and they will be useful for doublechecking the upcoming stateoftheart parallaxes delivered by the Gaia space mission in the near future.
This paper is organized as follows. In Sect. 2 we describe the moving cluster method and present a new implementation of the method that is used in this work. In Sect. 3 we present the sample of Pleiades stars and the dataset (proper motions and radial velocities) used in this analysis. Section 4 describes our convergent point analysis using different samples of cluster members to investigate our solution. In Sect. 5, we present the individual parallaxes (i.e., distances) for Pleiades stars obtained in this work, and discuss our results in Sect. 6. We summarize our results and conclusions in Sect. 7.
2. The moving cluster method
2.1. Kinematic parallaxes
Moving groups are kinematic aggregates of stars whose members share the same space motion. The low velocity dispersion, typically of a few km s^{1} (Mathieu 1986), allows them to be easily identified as overdensities in the velocity space (see e.g., Antoja et al. 2008). As discussed in the seminal work of Blaauw (1964), the proper motions of comoving stars in such kinematic groups converge to a vertex referred to as the convergent point (CP) of the moving group, because of their parallel space motion or the state of linear expansion.
In this context, individual parallaxes for members of young moving groups can be derived from their proper motions, radial velocities and the CP position. They are given by, (1)where A = 4.74047 km yr/s is the ratio of one astronomical unit in km to the number of seconds in one Julian year, μ_{∥} is the proper motion component directed parallel to the great circle that joins the star and the CP (see Galli et al. 2012, for more details), V_{r} is the radial velocity of the star, and λ is the angular distance between the star and the CP position.
Alternatively, for binaries (or multiple systems) and group members with unknown (or poor) radial velocities one can derive the kinematic parallax of the star from the spatial velocity of the moving group under the assumption that all members share the same space motion. In this case, it is given by (2)where V_{space} is the average spatial velocity of the moving group. In both cases, the parallax uncertainty is calculated from error propagation of the equations above, and takes the error budgets owing to proper motions, radial velocities, (or the spatial velocity) and the CP into account. The parallaxes derived in this way are obviously not as precise as those obtained from VLBI observations, but errors in the same order of the Hipparcos catalog and recent trigonometric parallax results from the ground (Weinberger et al. 2013; Ducourant et al. 2014) can be achieved (see e.g., Bertout & Genova 2006; Galli et al. 2013).
2.2. A new implementation of the convergent point search method based on Markov chain Monte Carlo methods
The computation of kinematic parallaxes using the moving cluster method requires prior knowledge of the CP position. The technique that we use to calculate the CP of a moving group is based on methods developed by de Bruijne (1999) and Galli et al. (2012; see these papers for a more detailed description about the basic concepts, algorithms, and implementation of the method). As discussed in these papers, the convergent point search method (CPSM) simultaneously determines the most likely CP position and performs a membership analysis to select the moving group members.
The CPSM takes stellar proper motions, the velocity dispersion, and a distance estimate of the cluster as input parameters. However, in general, the results are not sensitive to the assumed distance of the cluster (see also Mamajek 2005). We anticipate that the results and conclusions for the Pleiades cluster presented in the upcoming sections are rather insensitive to the assumed distance of the cluster ranging from the Hipparcos distance (~120 pc) to the VLBI distance (~136 pc). On the other hand, the velocity dispersion is one important input parameter in the CPSM to identify the moving group members. While a low velocity dispersion value would not be sufficient to recover all cluster members, a high velocity dispersion would enable the method to include field stars (interlopers) in the solution.
Fig. 1 CP coordinates (upper panel) and associated errors (lower panels) for the TW Hydrae association as a function of the number of walkers (with 200 iteration steps) used in the MCMCbased version of the CPSM. Each point is an average of 1000 MCMC realizations of the method and the error bars indicate the scatter of the solution given for each ensemble of walkers. 
The original implementation of the CPSM (de Bruijne 1999; Galli et al. 2012) uses an analytical minimization routine to solve the leastsquare equations and return the CP position with its uncertainty. Although this procedure has already proved to deliver good results for nearby associations with tens of stars (Galli et al. 2013; Ducourant et al. 2014), it was clearly less attractive when dealing with rich clusters, such as the Pleiades, which require more computing time to converge to a final solution. Thus, we developed a new version of the CPSM that exploits parallelism and exhibits good performance at low computational cost. This was achieved using the emcee package (ForemanMackey et al. 2013) written in Python programming language, which implements the Markov chain Monte Carlo (MCMC) method proposed by Goodman & Weare (2010). The algorithm was adapted to our purposes and applied to the general problem of computing the CP coordinates (α_{cp},δ_{cp}) of a moving group. Briefly, the new implementation of the CPSM that we use in this work exploits the parameter space using a number of socalled walkers and iteration steps, which are defined by the user to search for the CP position via Bayesian inference. The walkers move around the parameter space (α,δ) and take tentative steps towards the lowest valley of X^{2} that roughly defines the CP position (see de Bruijne 1999; Galli et al. 2012). Our final CP solution obtained after one MCMC realization is computed from the distribution of individual solutions given by the ensemble of walkers.
Fig. 2 Distribution of the CP coordinates after 1000 MCMC realizations of the CPSM for the TW Hydrae association. The red solid line indicates the kernel density estimator. 
In the following, we apply the CPSM to the TW Hydrae association whose CP position is well established (Mamajek 2005; Ducourant et al. 2014) to calibrate our new implementation of the CPSM and illustrate its application. In this context, we use the sample of 30 stars identified as kinematic members of the association by Ducourant et al. (2014). We vary the number of walkers from 10 to 200 and run the MCMC version of the CPSM using 50, 100, and 200 iteration steps for each walker to compare the results. Figure 1 shows the CP coordinates with their errors obtained for 200 iteration steps. We conclude that the formal errors on the CP coordinates are significantly larger than the observed scatter in the individual CP solutions obtained with the ensemble of walkers used in each case. This makes our choice of the minimum number of walkers to be used in our analysis rather arbitrary. In practice, we verified that our final solution is not sensitive to the number of walkers in the range of 100 to 200 if we run a significant number of MCMC realizations. We also confirmed that convergence of the Markovchains of the ensemble of walkers was attained after 50 iterations where the mean of both CP coordinates is clearly bounded by the variance of the sample. So, we decided to work with 100 walkers and 200 iteration steps as a more conservative approach. In addition, we emphasize that all CP solutions presented in this paper are calculated from a total of 1000 MCMC realizations of the CPSM, and the results listed are averaged values obtained from the distribution of the CP coordinates.
Fig. 3 Number of positions (left panel) and timebase (right panel) used to compute the stellar proper motions for Pleiades stars in the DANCe project. The red solid lines indicate the kernel density estimator. The blue dashed lines indicate the selection criteria applied to our sample as described in Sect. 4.1. 
Doing so, we run the CPSM for the TW Hydrae association using a velocity dispersion of 0.8 km s^{1} and distance estimate of 50 pc (as described in Sect. 4.1 of Ducourant et al. 2014). The resulting CP solution is located at (see also Fig. 2) We conclude that our result is in very good agreement with the more recent CP solution obtained by Ducourant et al. (2014), (α_{cp},δ_{cp}) = (100.1°,−27.1°) ± (3.0°,1.3°), but the formal errors of our solution are smaller. Thus, we have tested and calibrated our new version of the CPSM based on the MCMC method, and it will be used in the upcoming sections to investigate the CP position of the Pleiades cluster.
Fig. 4 Distribution of the relative errors on proper motions (i.e., μ/σ_{μ}) in right ascension (left panel) and in declination (right panel). The solid lines indicate the normal components of our mixture model (with k = 4), and the dashed line indicates the observed distribution of proper motions errors. 
3. Sample of Pleiades stars
3.1. An updated census of the cluster (Bouy et al. 2015)
In a recent work, Bouy et al. (2015) apply a probabilistic method based on multivariate data analysis (see Sarro et al. 2014) to select high probability members of the Pleiades open cluster in the DANCe (Bouy et al. 2013) and Tycho2 (Høg et al. 2000) catalogs. They identify a total of 2107 high probability members of the Pleiades open cluster making it the most complete census of the cluster to date. This list includes most cluster members identified in previous surveys (e.g., Stauffer et al. 2007; Lodieu et al. 2012) and 812 previously unknown members. As discussed in Sect. 2 of their paper, the absence of some cluster members in this sample is justified by their different properties implying low membership probabilities, while a few of them may have been missed because they fall in the domain of incompleteness between the DANCe survey (i ≥ 14 mag) and the Tycho2 catalog (V ≤ 12 mag).
The catalog of Pleiades stars given by Bouy et al. (2015) contains stellar positions, proper motions, multiband photometry (filters u,g,r,i,z,Y,J,H,K_{s}) and membership probabilities. This is the initial sample of stars that we use in our forthcoming analysis to investigate the distance of the Pleiades.
3.2. Proper motions
Our initial sample of Pleiades stars consists of 2010 stars (prob. ≥75%) from the DANCe survey and 207 stars (prob. ≥48%) from the Tycho2 catalog with 110 stars in common (see Bouy et al. 2015). Thus, the proper motions that we use in this work come mostly from the DANCe survey (1900 stars), and we use the proper motions from Tycho2 for the remaining stars (i.e., the brightest stars in our sample). The proper motions from the DANCe catalog provided by the Vizier/CDS tables are not anchored to an absolute reference system. We applied the procedure described in Sect. 7.10 of Bouy et al. (2013) to tie these proper motions to the International Celestial Reference System (ICRS). The corrected measurements are presented in the upcoming tables, together with other results from our analysis.
The DANCe project combines multiepoch panchromatic images collected with different instruments at various observatories (see Bouy et al. 2013, for more details) to derive accurate proper motions over a wide field from the ground. The proper motions were computed using, on average, about 40 stellar positions collected at different epochs and a timebase longer than 10 yr for most stars in our sample (see Fig. 3). The stars have a median proper motion of (μ_{α}cosδ,μ_{δ}) = ( + 19,−44) mas/yr and a median precision better than 1 mas/yr in each component.
Figure 4 shows the relative errors on proper motions in each component and reveals the existence of different subgroups in our sample regarding the quality of their proper motions measurements. One reason for the observed distribution of proper motion errors is the different observing campaign for each star in the sample, which is based on archival data collected in the DANCe project.
To identify the stars with best proper motion measurements in our sample (to be used in Sect. 4.2) we model the observed distribution of proper motion errors as a mixture of kcomponents. In this case, each component follows a normal distribution with mean μ and variance σ^{2}. The mixture model probability distribution is given by (3)where φ_{i} is the probability that a star belongs to the kcomponent of the model. We use the mixtools package in R programming language to calculate the mixture model parameters (φ_{i},μ_{i},σ_{i}) based on the expectationmaximization algorithm implemented in the normalmixEM procedure.
We perform a KolmogorovSmirnov (KS) test between the observed distribution of proper motion errors (see Fig. 4) and the empirical probabilitydensity function that is given by our mixture model to investigate the number of components (i.e., the value of k) that best represent our data. To assess the robustness of this finding, we construct a total of 1000 synthetic samples from Eq. (3) using k = 2,3, and 4. Then, we perform a KS test between each synthetic dataset and the observed distribution of proper motion errors. We compute the fraction f_{0.05} of synthetic realizations of our data sample with a pvalue higher than the adopted significance level of α = 0.05 (see e.g., Feigelson & Jogesh Babu 2012). Our results from this analysis are summarized in Table 1 and they indicate that the mixture model with four components (k = 4) is the one that best describes the observed distribution of proper motions errors. These results will be used in Sect. 4 to construct a control sample to investigate the CP position of the cluster.
Results of the mixture model applied to the distribution of proper motions errors in right ascension and declination.
3.3. Radial velocities
The radial velocity is one important parameter that is needed to derive the kinematic parallax of cluster members in our analysis. In the following, we describe our search in the literature for radial velocity information for the stars in our sample.
First, we searched the Apache Point Observatory Galactic Evolution Experiment (APOGEE) catalog in the Sloan Digital Sky Survey III (SDSSIII). The APOGEE program completed a systematic and homogeneous spectroscopic survey to build a highresolution (R ~ 22 500), nearinfrared spectra database for more than 10^{5} stars in our Galaxy (Majewski et al. 2016). We identified a total of 174 stars from our sample in the APOGEE survey. The zeropoint of the APOGEE radial velocities is estimated to be −0.355 ± 0.033 km s^{1} from comparison with other studies (see Nidever et al. 2015). As discussed in Appendix B of Galli et al. (2013), a small shift in the radial velocities could lead to a more significant offset in the parallaxes derived from the moving cluster method. We verified that a small variation of 0.355 km s^{1} in the radial velocities would produce an offset of 0.45 mas in parallaxes using typical values of proper motions, radial velocities, and the CP position of the Pleiades cluster that will be presented in the forthcoming analysis. Thus, to avoid a systematic bias in the resulting parallaxes we decided to correct the APOGEE radial velocities from their absolute zeropoint before using them in this work.
Then, we searched the CDS databases to access more radial velocity data for the remaining stars in the sample. The search made use of a query for radial velocity information in script mode using the webbased data mining tools available in the SIMBAD database (Wenger et al. 2000). Our search results for radial velocities is based on: Wilson (1953), Gontcharov (2006), Kharchenko et al. (2007), White et al. (2007), Mermilliod et al. (2009) and Kordopatis et al. (2013). Doing so, we retrieved radial velocity information for 241 stars. We note that 21 stars in this list are in common with the APOGEE sample and we calculated the weighted mean of the multiple radial velocity values.
Thus, our search for radial velocity data returned a sample of 394 stars with at least one radial velocity measurement published in the literature. However, this sample also includes stars with poor radial velocity measurements and obvious outliers. After removing obvious outliers from the radial velocity distribution by a 3σ elimination, we end up with a sample of 340 stars. Figure 5 shows the distribution of radial velocity data for the stars in our sample. The average radial velocity of this sample is 5.6 ± 0.1 km s^{1} with a median value of 5.4 km s^{1}. The radial velocities collected in this work will be used in Sect. 5 to derive individual kinematic parallaxes of Pleiades stars from the moving cluster method.
Fig. 5 Distribution of radial velocities for Pleiades stars collected from the literature. 
3.4. Binaries and multiple systems
Binaries and multiple systems are common features during star formation and various surveys in different regions and clusters have been performed to measure the frequency of such systems (see e.g., Ghez et al. 1993; Duchêne et al. 1999; Daemgen et al. 2015). Previous studies on the Pleiades cluster indicate that the binary fraction varies from 28% to 44%, depending on the mass range, orbital period, and angular separation of the binary systems (Bouvier et al. 1997; Lodieu et al. 2007).
Despite previous efforts, the binarity/multiplicity census of the Pleiades cluster is still far from complete and would require a large observing programme, which is clearly beyond the scope of this paper. However, unresolved binaries and highorder multiple systems (n ≥ 3) appear brighter than single stars and can be identified in a colormagnitude diagram (CMD). To do so, we plot different CMDs using the i,J,H,K_{s} filters that are available for most stars in our sample, and identified those stars that define a binary sequence as potential binaries and multiple system candidates. Figure 6 illustrates this procedure for one CMD used in our analysis. We consider the star to be a binary or multiple system if it is classified as such in more than one CMD. Thus, we found 510 binary or multiple systems that amounts to 24% of our initial sample. Although this procedure is more likely to remove only the equalmass binaries, it serves as a first estimate to characterize the binarity/multiplicity in our sample, and it will be useful in our forthcoming analysis to investigate the distance of the cluster. We note that we are not rejecting these stars from the sample. As explained in Sect. 2.2, the parallaxes of binaries and multiple systems derived in this work will be inferred from the average spatial velocity of the cluster (see Sect. 5).
Fig. 6 CMD to illustrate the locus of potential binaries and multiple systems that might be present in our sample. 
4. Convergent point analysis
In the following, we investigate the CP position of the Pleiades open cluster obtained from different subsets of our initial sample of stars. First, we build two control samples as a quality check of our results and to discuss the velocity dispersion of the cluster. Then, we apply the CPSM to the full sample of stars and investigate our solution based on Monte Carlo simulations.
4.1. Preliminary inspection of proper motion data
The accuracy and precision of proper motions are of ultimate importance in the CP analysis. To ensure the quality of our results presented in the forthcoming sections, we proceed as follows.
First, we select from our initial sample of 2107 stars only those stars whose proper motion measurements have been computed with a minimum of n ≥ 10 positions and a time base Δt ≥ 10 yr (see Fig. 3). These selection criteria reduce the sample to 1657 stars.
Then, we remove the outliers from the proper motion distributions based on a 3σ elimination in each component. Doing so, we end up with a sample of 1551 stars that will be used in the remaining of this paper to investigate the distance of the cluster.
4.2. Control sample 1: the high probability cluster members
The DANCe catalog is not a homogeneous dataset because the stars in our sample have different astrometric observing campaigns and the multiband photometry (u,g,r,i,z,Y,J,H,K_{s}) is not complete for all stars. The membership probabilities computed from the DANCe project take both astrometric and photometric data into account, and they are more accurate for stars with a complete dataset.
With this in mind, we apply the following selection criteria (in addition to those described in Sect. 4.1) to identify the most likely cluster members and the stars with more precise proper motion measurements in our sample:

1.
stars with complete data (astrometry + photometry);

2.
proper motion errors (μ_{α}/σ_{μα}) ≥ 47 and (μ_{δ}/σ_{μδ}) ≥ 128;

3.
membership probability p ≥ 0.9975 (3σ).
Our selection criterium #2 refers to the stars that roughly define the last component of the mixture model applied to the distribution of proper motion errors in Fig. 4. Our choice of using the membership probability threshold of 99.75% (see criterium #3) is justified to minimize, as much as possible, the contamination rate in our control sample (see Table 4 of Sarro et al. 2014, for more details). Doing so, we end up with a sample of 296 stars.
Then, we applied the CPSM with σ_{v} = 0 to identify the stars that show strict convergence to the CP of the moving group. We find a sample of 118 stars that defines a core moving group of the cluster (thus, our control sample 1). The corresponding CP is located at with chisquared statistics (i.e., χ^{2}/ν = 133.0/116) and correlation coefficient of ρ = −0.99. We note that our first CP estimate is already consistent with the solution of (α_{cp},δ_{cp}) = (92.5°,−47.9°) ± (5.4°,5.3°) derived by Makarov & Robichon (2001). However, our result is more precise, which is a result of the more precise proper motion data available nowadays and the improved methodology (described in Sect. 2.2) to calculate the CP position.
We repeated this analysis after removing the binary (and multiple system) candidates from our sample (as described in Sect. 3.4), and confirmed that the resulting CP position is in perfect agreement with the solution given above. For clarity of presentation, we provide only the results, including binaries and multiple systems, because they contain more stars in the solution. In Appendix A we provide an alternative approach to estimate the CP coordinates of the Pleiades that also supports the results presented in this section.
4.3. Control sample 2: the nuclear members of the cluster
The DANCe catalog for the Pleiades open cluster covers an area of ~ 80deg around the cluster center. Ideally, one would expect the contamination rate by field stars to increase with increasing distance to the cluster center. To minimize the existence of field stars (interlopers), we decided to build a second control sample by retaining only the nuclear members of the cluster within 2deg from the cluster center (α ≃ 57°, δ ≃ + 24°) out of the sample of 1551 stars, which was selected as described in Sect. 4.1. Doing so, we end up with a sample of 914 stars that will be used to investigate the typical velocity dispersion in the center of the cluster. In the outer regions, the velocity dispersion required to identify all cluster members might be larger owing to a combination of several effects (e.g. mass segregation, stellar encounters, tidal disruption, and cluster evaporation). Thus, we argue that the cluster members identified by the CPSM in this work represent a minimum moving group of the Pleiades cluster. As illustrated in Fig. 7, this sample contains most of the stars with known distances in the literature that can be used to compare with our results.
Fig. 7 Distribution of stellar positions (left panel) and proper motions (right panel) for the Pleiades stars in the DANCe catalog. The various symbols and colors indicate the nuclear members of the cluster that define our control sample 2, and the stars with individual distances reported in the literature. 
Fig. 8 Results obtained from the CPSM (CP coordinates and number of moving group members) as a function of the velocity dispersion of the cluster. The solid and dashed lines indicate the CP solution derived from our control sample 1 at 1σ and 3σ (see Sect. 4.1). 
Fig. 9 CP solutions (black dots) for 1000 Monte Carlo realizations of the Pleiades open cluster overlaid on the X^{2} contour levels (solid lines) indicated by the color code. The red triangle indicates the final CP solution given in Sect. 4.4. 
Then, we run the CPSM using different values for the velocity dispersion of the cluster from 0 to 2 km s^{1} in steps of 0.1 km s^{1}. In each step we calculate the CP position and count the number of stars identified by the CPSM as kinematic moving group members. The results of this analysis are shown in Fig. 8. We observe that the solution is rather unstable at low values for the velocity dispersion, although the resulting CP positions are still statistically compatible with the results obtained from control sample 1 (see Sect. 4.1). The scatter reduces significantly after σ_{v} = 0.5 km s^{1} with the inclusion of more cluster members in the solution. We note that with σ_{v} = 0.8 km s^{1}, our solution retrieves the maximum number of cluster members (732 stars) that can be identified by the CPSM from this sample of stars, and that the CP coordinates converge to the solution obtained with the control sample 1. Thus, our methodology confirmed 80% of the selected sample as moving group members, as expected for the CPSM (see Fig. 8 of Galli et al. 2012). Increasing the velocity dispersion in the CP analysis for values higher than σ_{v} = 0.8 km s^{1} only increases the noise (i.e., errors) in the CP coordinates with no benefit of including more cluster members (see Fig. 8). Thus, based on this analysis, we infer the velocity dispersion of the cluster to be σ_{v} = 0.8 ± 0.1 km s^{1} (see also discussion Sect. 4.5). The corresponding CP solution is located at with chisquared statistics (i.e., χ^{2}/ν = 753.0/730) and correlation coefficient of ρ = −0.99. Figure 7 illustrates the distribution of the stellar proper motion for the 732 stars selected by the CPSM. These stars define our control sample 2, and the results obtained in this section will be used as a quality check of our final CP solution for the Pleiades cluster (see below).
4.4. The CPSM applied to the full sample of Pleiades stars
Our final analysis consists in running the CPSM on the selected sample of 1551 stars used in this work (see Sect. 4.1). Using a velocity dispersion of σ_{v} = 0.8 km s^{1}, we identify a moving group with 1210 stars and CP located at with chisquared statistics (i.e., χ^{2}/ν = 1142.8/1208) and correlation coefficient of ρ = −0.99. The resulting CP solution is in good agreement with the results presented in Sects. 4.2 and 4.3. Based on this investigation, we consider the 1210 stars identified directly by the CPSM to be confirmed members of the cluster.
To gain more confidence in the results presented above, we constructed a number of 1000 synthetic samples of the Pleiades cluster based on Monte Carlo simulations. These samples are generated by resampling the stellar proper motions of individual stars from a Gaussian distribution with mean and variance that corresponds to the proper motion measurements and their uncertainties. Then, we run the CPSM for each synthetic realization of the Pleiades cluster and compute the CP position. The results of this analysis are shown in Fig. 9. The centroid of the Monte Carlo simulations is located at with and ρ = −0.99. This confirms our CP solution showing that it is indeed representative of the Pleiades open cluster.
All the investigations reported in this section, using different samples and techniques to derive the CP position of the Pleiades cluster, make us confident that our solution is well constrained. It will be used in Sect. 5 to calculate individual distances of cluster members.
4.5. A posteriori assessment of the velocity dispersion
An alternative approach for investigating the velocity dispersion of the cluster consists of analyzing the distribution of the proper motion component that results from the peculiar motions of the stars. To do so, we calculate the proper motion component μ_{⊥} that is directed perpendicular to the great circle that joins the star and the CP of the moving group following the procedure outlined in Sect. 2.1 of Galli et al. (2012). The distribution of μ_{⊥} values for the 1210 cluster members identified in Sect. 4.4 is shown in Fig. 10. The average and standard deviation of the distribution are, respectively, 0.03 mas/yr and 1.45 mas/yr. As expected, the average value of μ_{⊥} is close to zero, which reflects the good convergence of the stellar proper motions to our CP solution.
The onedimensional velocity dispersion of the cluster can be roughly estimated from σ_{v} ≃ 4.74 dσ_{μ}, where we consider σ_{μ} to be the observed scatter that comes from the μ_{⊥} statistics. The resulting velocity dispersion depends on the assumed distance d of the cluster that is needed to convert the dispersion of proper motions (in units of mas/yr) to velocity (in units of km s^{1}). Using the Hipparcos and VLBI distance estimates of ~ 120 pc and ~ 136 pc for the Pleiades we find a velocity dispersion of 0.82 km s^{1} and 0.93 km s^{1}, respectively. These numbers confirm the value of σ_{v} = 0.8 ± 0.1 km s^{1} inferred from the analysis presented in Sect. 4.3.
Fig. 10 Distribution of the μ_{⊥} values calculated from the CP solution given in Sect. 4.4 for the 1210 cluster members. The red solid line indicates the kernel density estimator. 
5. Movingcluster distance to the Pleiades
In the following we discuss the distance of the Pleiades based on the CP solution given in Sect. 4.4 and the sample of 1210 cluster members identified in that analysis. First, we use the sample of stars with known radial velocity to calculate individual parallaxes and the spatial velocity of each star. Then, we calculate approximate parallaxes for other group members from the spatial velocity of the cluster.
5.1. Kinematic parallax and spatial velocity for cluster members with known radial velocity
We calculate the individual kinematic parallax for the stars with known radial velocity from Eq. (1), and use them to calculate the threedimensional Galactic UVW velocities from the procedure described by Johnson & Soderblom (1987). In a recent paper, BailerJones (2015) presents the common approach of estimating distances from parallaxes as an inference problem that requires the use of prior assumptions. Consequently, we first apply the methodology described in Sect. 7 of his paper, using an exponentially decreasing volume density prior (with L = 10^{4}) to avoid biases in the computed UVW velocities. One direct consequence of this approach is that the UVW error bars of individual stars are not symmetric. Then, we perform an iterative clip on the distribution of spatial velocities using Chauvenet’s criterion (see e.g. Bevington & Robinson 2003) to remove obvious outliers and spot possible errors on the parallaxes derived from radial velocities. This leaves us with a sample of 64 stars with known radial velocity, which is used to accurately constrain the spatial velocity of the cluster. Binaries and multiple systems (see Sect. 3.4) are excluded from this count, because it is not possible to derive their kinematic parallaxes based on a single radial velocity measurement.
Table B.1 lists the parallaxes derived in this work for cluster members with known radial velocity. The mean parallax is 7.44 ± 0.08 mas (median of 7.38 mas and standard deviation of 0.66 mas). This is consistent with a distance estimate of pc given a confidence interval of 95%. This result is in good agreement with the VLBI distance determination of 136.2 ± 1.2 pc reported by Melis et al. (2014), but it is not consistent with the Hipparcos distance of 120.2 ± 1.9 pc (van Leeuwen 2009).
Individual values for the space velocity of each star are also listed in Table B.1. Figure 11 illustrates the distribution of the Galactic space velocity in each component and Table 2 lists their average values. In a previous paper, Robichon et al. (1999) report the space velocity of (U,V,W) = (−6.4,−24.4,−13.0) ± (0.5,0.7,0.4) km s^{1} (not corrected for the Solar motion) based on Hipparcos parallaxes. The difference of (ΔU,ΔV,ΔW) = (0.2,−4.3,−1.7) ± (0.5,0.8,0.4) km s^{1} between these two estimates is significant and comes from the different parallaxes used for each star in both studies.
Fig. 11 Histogram of the Galactic space velocity components for the 64 cluster members with known radial velocities. The red solid line indicates the kernel density estimator. 
Spatial velocity of the Pleiades cluster derived in this work from the sample of 64 stars with known radial velocities.
5.2. Kinematic parallax for other cluster members
Once the spatial velocity of the cluster is defined we can derive approximate parallaxes to the remaining group members under the assumption that all stars in a moving group share the same space motion. As a consistency check, we first compare the parallaxes derived from Eqs. (1) and (2) for the sample of 64 stars with known radial velocities. Figure 12 illustrates this comparison. We find a mean difference and root mean square (rms) of 0.01 mas and 0.58 mas, respectively, between the two methods to compute kinematic parallaxes. These numbers are smaller than the typical error (about 1.0 mas) of the parallaxes listed in Table B.1 and confirm that both techniques yield consistent results. We performed a consistency check of the parallaxes derived from Eq. (2) using the mean, median, and mode of the cluster spatial velocity. We verified that the agreement between both methods to compute kinematic parallaxes is better when we use the mode instead of the mean (or median) of the cluster spatial velocity owing to the asymmetry observed in the distribution for the total spatial velocity V_{space} (see Fig. 11).
Fig. 12 Comparison of the parallaxes derived from individual radial velocities and the mode of the cluster spatial velocity for the sample of 64 stars. The red solid line indicates perfect correlation. 
Fig. 13 Probability density function for the parallaxes obtained from the sample of members with known radial velocities (64 stars) and all cluster members (1210 stars). 
Then, we calculate from Eq. (2) the kinematic parallax for the remaining 1146 cluster members with unknown (or poor) radial velocities, binaries and multiple systems in our sample of stars selected by the CPSM. Table B.2 lists the resulting kinematic parallaxes for these stars. As explained in Sect. 2.1, the associated uncertainties take into account the errors on proper motions, the CP position, and the average error of the individual spatial velocities derived for the 64 stars with known radial velocities (see Sect. 5.1). The rms of 0.58 mas that results from the analysis shown in Fig. 12 is added quadratically to the parallax uncertainties to consider for possible systematic errors that could arise from using the spatial velocity of the cluster to derive individual kinematic parallaxes.
The mean parallax of the complete sample of cluster members (1210 stars) is 7.44 ± 0.01 mas (median of 7.46 mas and standard deviation of 0.42 mas). This result is consistent with our first estimate using radial velocities (see Sect. 5.1), but it is more precise, which comes naturally from the more significant number of stars used to derive the standard error of the mean. Figure 13 illustrates the density distribution of parallaxes obtained for these two samples. We note that the observed dispersion of parallaxes is also smaller for the complete sample. This mostly reflects the spread of the stellar proper motions since radial velocities are available for only a small fraction of the sample. However, we retain the results given in Sect. 5.1 as our final solution in this work, because the individual parallaxes computed from radial velocities, when available, are better constrained.
6. Discussion
6.1. Parallaxes for different samples of the Pleiades cluster
In Sect. 4 we built two control samples based on the membership probability of cluster members and their distance to the cluster center to guide our CP analysis. Here, we perform a posteriori assessment of our results to investigate the distance of the Pleiades given by different subsets of our sample of cluster members selected by the CPSM (see Sect. 4.4). We construct these samples by selecting the stars within a given radius from the cluster center and membership probability threshold. Table 3 lists the basic properties and Fig. 14 illustrates the parallax distribution of each sample.
Fig. 14 Probability density function for the parallaxes of different samples of cluster members selected from their membership probability (left panel) and distance to the cluster center (right panel). 
It is apparent from this analysis that the parallax results obtained with different samples of cluster members are consistent between themselves. In particular, we note that the average parallaxes listed in Table 3 are all compatible at the 3σ level. This gives us more confidence in the distance determination of the Pleiades open cluster derived in this paper (see Sect. 5.1) and confirms that it is not dependent on the adopted selection criteria.
We also observe from Table 3 that the standard deviation of parallaxes increases with the distance to the cluster center. Figure 15 shows the parallax dispersion as expected from the cluster depth along the line of sight as a function of the apparent magnitudes (from the ifilter) and radius from cluster centre. This diagram provides information on the spatial distribution and segregation of stars within the cluster. We note a clear decrease of the stellar magnitudes with increasing radius and parallax dispersion. This confirms the effect of mass segregation in the cluster (using the stellar magnitudes as a proxy for spectral types) as anticipated in previous studies (see e.g. Raboud & Mermilliod 1998; Palmer et al. 2014).
Kinematic parallax of the Pleiades derived from different subgroups in our sample of cluster members selected by the CPSM.
Fig. 15 Apparent magnitudes of cluster members as a function of the distance to the cluster center and parallaxes for individual stars. 
6.2. Comparison with previous parallaxes of individual stars
Here, we compare the parallaxes derived in this work with previous results published in the literature for individual stars. Fortunately, our sample of 1210 cluster members contains 64 stars from various sources in the literature that applied different techniques to measure the stellar parallaxes. Their parallaxes are listed in Table B.3. This comparison includes results from Gatewood et al. (2000), Melis et al. (2014), Mädler et al. (2016), Hipparcos catalog (van Leeuwen 2007) and the first results delivered by the Gaia satellite. In this context, it is important to mention that the binary system HII 2147 NW+HII 2147 SE from the VLBI study of Melis et al. (2014) is not resolved in our sample. To make a better comparison with our results we decided to work with the weighted average of the proper motions and parallaxes of the two components.
Fig. 16 Comparison between the kinematic parallaxes derived in this work and parallaxes from the literature (van Leeuwen 2007; Melis et al. 2014; Mädler et al. 2016; Gatewood et al. 2000; Gaia Collaboration 2016). The black dashed line indicates the expected distribution for equal results. 
Figure 16 compares our kinematic parallaxes with published results in the literature, and Table 4 presents the mean difference and rms that result from the comparison with each source. This comparison confirms the good accuracy of our results when compared to the individual parallaxes of twin stars obtained by Mädler et al. (2016) and the very precise multiepoch VLBI trigonometric parallaxes from Melis et al. (2014) given that these works are independent. The comparison with the trigonometric parallaxes measured by Gatewood et al. (2000) returned that the parallax for one star is not consistent with our results, although we observe a good agreement with the other two stars in common. In general our results do not reproduce the parallaxes from the Hipparcos catalog. But, interestingly, a few stars in our sample with parallaxes between 7.0 mas and 8.0 mas show a reasonable agreement with the Hipparcos parallaxes. The question arises whether the discrepancy of Hipparcos parallaxes indeed applies to all cluster members as commonly reported in the literature.
Comparison of the parallaxes derived in this work with other studies in the literature.
During the submission process of the current manuscript, the first Gaia data release (Gaia DR1) became available providing astrometry and photometry for 1 billion sources brighter than magnitude 20.7 in the satellite’s whitelight photometric band (Gaia Collaboration 2016). Gaia DR1 also includes trigonometric parallaxes and mean proper motions for about 2 million stars from the TychoGaia astrometric solution (TGAS). We use this first dataset of trigonometric parallaxes, which has just been delivered by the Gaia collaboration, to make a comparison with the results obtained from the moving cluster method in this work. As expected, this is the source of trigonometric parallaxes with the largest number of stars in common with our sample of Pleiades stars (see Table 4). This gives us the possibility to perform a more detailed comparison of the results delivered by these studies.
As illustrated in Fig. 16, the TGAS trigonometric parallaxes are in good agreement with the kinematic parallaxes derived from the movingcluster method and the other independent studies discussed before. Indeed, this confirms the longterm suspicion that the Hipparcos parallaxes for Pleiades stars are biased. However, there are still a few stars in our sample that exhibit discrepant proper motions and/or parallaxes with respect to the TGAS results (at the 3σ level). We summarize their properties in Table 5.
Comparison of the proper motions and parallaxes with TGAS results from Gaia DR1 for discrepant stars.
Fig. 17 Comparison between the proper motions measurements in right ascension (left panel) and declination (right panel) from the DANCe project and the TGAS catalog. The gray squares indicate the discrepant stars listed in Table 5. 
We note from Fig. 16 that the kinematic parallaxes derived in this work for three stars (HIP 17999, TYC2 12468081, and TYC2 18079551) are not consistent with TGAS. It is clear from Table 5 that the proper motions obtained for HIP 17999 in both surveys are consistent between each other and also in good agreement with the observed proper motions for other cluster members (see, for example, Fig. 7). Interestingly, the TGAS trigonometric parallax of π = 3.40 ± 0.73 mas is also not consistent with the results reported by the Hipparcos satellite: π = 9.83 ± 1.00 mas (ESA 1997) and π = 9.93 ± 0.75 mas (van Leeuwen 2007). So, if we assume that HIP 17999 is indeed a cluster member, as suggested by its proper motion, then the TGAS trigonometric parallax for this target needs further clarification. On the other hand, TYC2 12468081 and TYC2 18079551 exhibit different proper motion values in the right ascension component that justifies the different parallax results in both studies. We searched for proper motion measurements for these two stars in the UCAC4 (Zacharias et al. 2013) and PPMXL (Roeser et al. 2010) catalogs, and verified that they are more consistent with the TGAS results. As discussed in previous papers of the DANCe project (see e.g. Bouy et al. 2013, 2015) the proper motion measurements close to the saturation limit of the instruments (that is considered to be the limit of sensitivity of the Tycho2 catalog) suffer from biases and incompleteness, such that they should be considered with caution. Replacing the DANCe proper motions for these two stars by the TGAS results will force them to be rejected in our membership analysis. We verified that removing TYC2 12468081 and TYC2 18079551 from our sample has negligible impact on the results presented in the previous sections.
Fig. 18 Summary of distance measurements in the literature for the Pleiades open cluster obtained from different methods. The results presented here come from: (i) moving cluster – this work (see Sect. 5.1), Röser & Schilbach (2013), Narayanan & Gould (1999); (ii) trigonometric parallax – van Leeuwen (1999), van Leeuwen (2009), Melis et al. (2014), Soderblom et al. (2005), Gatewood et al. (2000), TGAS (Gaia Collaboration 2016) – using result reported in this work (see Sect. 6.2); (iii) orbital modeling – Groenewegen et al. (2007), Southworth et al. (2005), Zwahlen et al. (2004), Munari et al. (2004), Pan et al. (2004), David et al. (2016); (iv) twin stars Mädler et al. (2016), (v) MS and isochrone fitting – An et al. (2007), Percival et al. (2005), Stello & Nissen (2001), Pinsonneault et al. (1998), Giannuzzi (1995) and Nicolet (1981). 
Figure 17 compares the DANCe and TGAS proper motions for the stars in common. In general there is a good agreement between the proper motions in both surveys, but another four stars (TYC2 180016211, TYC2 18001181, HIP 17921, and HIP 18559) exhibit proper motion measurements that are not consistent at the 3σ level (see also Table 5). However, the parallaxes derived in this work are still compatible with the TGAS results. The shift (or absolute difference) in the observed proper motions is smaller compared, for example, to TYC2 12468081 and TYC2 18079551 to produce a more significant bias in the parallaxes. In addition, we note that the proper motions for TYC2 12466171 are still compatible within 3σ, because of the large error bars of the DANCe measurement, but the difference in the right ascension component is 10.0 mas/yr and it clearly stands out in the comparison shown in Fig. 17. The systematic shift (or mean difference) between the two proper motion datasets (after removing the discrepant stars) is −0.2 mas/yr in right ascension and 0.3 mas/yr in declination. The rms of the comparison is 1.2 mas/yr in right ascension and 1.1 mas/yr in declination. This represents a good agreement between the DANCe and TGAS proper motions.
After removing the discrepant stars listed in Table 5 from our sample of stars in common with TGAS, the mean difference and rms given in Table 4 become 0.12 mas and 0.49 mas, respectively. The mean parallax derived from the TGAS results is 7.53 ± 0.19 mas with median of 7.58 mas and standard deviation of 0.38 mas. However, as reported in the paper describing the summary of the Gaia DR1 properties, there is a systematic uncertainty of about 0.3 mas, which is correlated over small spatial scales and must be added to the parallax uncertainties. Thus, we report the mean parallax of 7.53 ± 0.30 mas as a more conservative solution obtained from the TGAS results for the stars in common with our sample. This is consistent with a distance estimate of 133 ± 5 pc for the cluster. However, this result is somewhat different from the value of 134 ± 6 pc reported in Sect. 5.5 of the Gaia DR1 paper (see Gaia Collaboration 2016), because of the different sample of stars used in each case. So, we conclude that the distance to the Pleiades cluster derived from the TGAS trigonometric parallaxes is consistent with the distance determination obtained in this work and previous studies.
To summarize, the investigations reported above using parallaxes derived from independent studies and different samples of Pleiades stars show that there is a good agreement with the results obtained in this work, and this make us confident that our distance estimate (see Sect. 5.1) is indeed representative of the cluster.
6.3. Revisiting the distance of the Pleiades open cluster
Figure 18 summarizes the distances obtained for the Pleiades cluster in the literature by each method. This is an updated version of Fig. 1 from Melis et al. (2014) to include the more recent results.
We note that all methods employed so far yield a distance determination that is larger than the result obtained by Hipparcos. However, in many cases the error bars do not allow us to unambiguously reject the Hipparcos result. In this context, it is important to distinguish between the methods that yield distances to individual cluster members and the ones that can only be used to estimate the average distance of the cluster. Among the methods included in this discussion, the MS and isochrone fitting are the only methods that cannot be used to deliver individual distances to cluster members. Although these methodologies yield mostly smaller errors in the distance determination, the resulting distance modulus is representative of the cluster photometry, but not individual stars. The error on the average distance obtained in this work is sometimes greater compared to the other distance determinations of the cluster. This reflects the dispersion of individual parallaxes (see Sect. 5.1) in our solution, which is a direct result of the spread in proper motions and radial velocities.
As illustrated in Fig. 18, the average distance for the Pleiades cluster obtained in this work clearly supports the recent VLBI distance determination for the cluster, based on the very precise and accurate astrometry by Melis et al. (2014). Our results are also very consistent with other nonHipparcos distances reported by different studies over the last two decades. The weighted mean distance of the nonHipparcos results listed in Fig. 18 (including this work) yields a mean distance of 135.0 ± 0.6 pc (median of 133.8 pc and standard deviation of 2.7 pc). This confirms the historical discrepancy between the Hipparcos distance determination and the results given by other methods (see also Sect. 1).
Despite our efforts to investigate the distance of the cluster, we emphasize that individual distances for Pleiades stars are available for only a few cluster members in the literature (as discussed in Sect. 1). The results obtained in this work with the moving cluster method significantly increase the number of Pleiades stars with known individual distances, thus yielding a sample that is the largest one to date (even if we only consider the 64 stars with parallaxes derived from radial velocities).
7. Conclusions
We have implemented a new version of the CPSM based on the MCMC method that exploits parallelism and exhibits good performance when dealing with rich clusters, such as the Pleiades, with several hundreds of stars. This new implementation of the CPSM was calibrated on the TW Hydrae association and applied to the Pleiades open cluster.
We performed extensive investigations to search for the CP position of the Pleiades cluster using (i) different samples of stars, (ii) an alternative approach based on the intersection of great circles defined by the stellar proper motions in the sky, and (iii) Monte Carlo simulations. Our analysis allowed us to accurately constrain the CP position, infer the velocity dispersion, and confirm 1210 stars as kinematic members of the cluster.
The moving cluster method applied to our sample returned individual parallaxes for 64 stars based on their radial velocities, and approximate parallaxes for the remaining 1146 group members from the spatial velocity of the cluster. We emphasize that this is the largest sample of Pleiades stars with individual distances to date. Our results are in good agreement with previous parallaxes for individual stars published in the literature (excluding Hipparcos), and demonstrate the good accuracy of the distance determination obtained in this work. We report in this paper a distance estimate of pc that is consistent with previous distance determinations of the cluster given by different methods. In particular, our results support the recent distance determination of 136.2 ± 1.2 pc, obtained from the trigonometric parallax method that is based on VLBI multiepoch observations and that represents the most precise and direct alternative to measure distances nowadays. Our distance determination is also in good agreement with the distance of 133 ± 5 pc obtained from the TGAS trigonometric parallaxes for the stars in common with our sample.
Our results are not consistent with the Hipparcos distance determination of 120.2 ± 1.9 pc for the cluster, and they confirm the discrepancy of these parallaxes as reported extensively in the literature. The weighted mean distance of the nonHipparcos methods, including the results of this work, is 135.0 ± 0.6 pc, which exceeds the Hipparcos distance by roughly 12%.
This study represents one more step towards better constraining the distance of the Pleiades cluster before the upcoming and very precise trigonometric parallaxes from the Gaia space mission. Our results are already consistent with the first trigonometric parallaxes included in the TGAS catalog for the brightest stars in the cluster and they will be useful in the next future for a comparative analysis of other cluster members that were not included in the Gaia DR1.
Acknowledgments
We thank the referee for constructive comments that helped us to improve the manuscript. This work was funded by the São Paulo State Science Foundation (FAPESP). This research has made use of the computing facilities of the Laboratory of Astroinformatics at IAG/USP (São Paulo, Brazil). It has also made use of the SIMBAD database operated at the CDS (Strasbourg, France).
References
 Abad, C., Vieira, K., Bongiovanni, A., Romero, L., & Vicente, B. 2003, A&A, 397, 345 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 An, D., Terndrup, D. M., Pinsonneault, M. H., et al. 2007, ApJ, 655, 233 [NASA ADS] [CrossRef] [Google Scholar]
 Antoja, T., Figueras, F., Fernández, D., & Torra, J. 2008, A&A, 490, 135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 BailerJones, C. A. L. 2015, PASP, 127, 994 [NASA ADS] [CrossRef] [Google Scholar]
 Bertout, C., & Genova, F. 2006, A&A, 460, 499 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bevington, P. R., & Robinson, D. K. 2003, in Data reduction and error analysis for the physical sciences (Boston, MA: Mc GrawHill) [Google Scholar]
 Blaauw, A. 1964, ARA&A, 2, 213 [NASA ADS] [CrossRef] [Google Scholar]
 Bouvier, J., Rigaut, F., & Nadeau, D. 1997, A&A, 323, 139 [NASA ADS] [Google Scholar]
 Bouy, H., Bertin, E., Moraux, E., et al. 2013, A&A, 554, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bouy, H., Bertin, E., Sarro, L. M., et al. 2015, A&A, 577, A148 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Daemgen, S., Bonavita, M., Jayawardhana, R., Lafrenière, D., & Janson, M. 2015, ApJ, 799, 155 [NASA ADS] [CrossRef] [Google Scholar]
 David, T. J., Conroy, K. E., Hillenbrand, L. A., et al. 2016, AJ, 151, 112 [NASA ADS] [CrossRef] [Google Scholar]
 de Bruijne, J. H. J. 1999, MNRAS, 306, 381 [NASA ADS] [CrossRef] [Google Scholar]
 Duchêne, G., Bouvier, J., & Simon, T. 1999, A&A, 343, 831 [NASA ADS] [Google Scholar]
 Ducourant, C., Teixeira, R., Galli, P. A. B., et al. 2014, A&A, 563, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 ESA 1997, VizieR Online Data Catalog: I/239 [Google Scholar]
 Feigelson, E. D., & Jogesh Babu, G. 2012, in Modern Statistical Methods for Astronomy (UK: Cambridge University Press) [Google Scholar]
 Fix, E., & Hodges, J. 1951, Discriminatory analysis, Nonparametric discrimination: Consistency properties, Technical Report 4, Project Number 2149004, USAF School of Aviation Medicine, Randolph Field, Texas [Google Scholar]
 ForemanMackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [CrossRef] [Google Scholar]
 Gaia Collaboration (Brown, A. G. A., et al.) 2016, A&A, 595, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Galli, P. A. B., Teixeira, R., Ducourant, C., Bertout, C., & BenevidesSoares, P. 2012, A&A, 538, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Galli, P. A. B., Bertout, C., Teixeira, R., & Ducourant, C. 2013, A&A, 558, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gatewood, G., de Jonge, J. K., & Han, I. 2000, ApJ, 533, 938 [NASA ADS] [CrossRef] [Google Scholar]
 Ghez, A. M., Neugebauer, G., & Matthews, K. 1993, AJ, 106, 2005 [NASA ADS] [CrossRef] [Google Scholar]
 Giannuzzi, M. A. 1995, A&A, 293 [Google Scholar]
 Gontcharov, G. A. 2006, Astron. Lett., 32, 759 [NASA ADS] [CrossRef] [Google Scholar]
 Goodman, J., & Weare, J. 2010, Comm. Appl. Math. Comput. Sci., 5, 65 [CrossRef] [MathSciNet] [Google Scholar]
 Groenewegen, M. A. T., Decin, L., Salaris, M., & De Cat, P. 2007, A&A, 463, 579 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27 [Google Scholar]
 Johnson, D. R. H., & Soderblom, D. R. 1987, AJ, 93, 864 [NASA ADS] [CrossRef] [Google Scholar]
 Kharchenko, N. V., Scholz, R.D., Piskunov, A. E., Röser, S., & Schilbach, E. 2007, Astron. Nachr., 328, 889 [NASA ADS] [CrossRef] [Google Scholar]
 Kordopatis, G., Gilmore, G., Steinmetz, M., et al. 2013, AJ, 146, 134 [NASA ADS] [CrossRef] [Google Scholar]
 Lodieu, N., Dobbie, P. D., Deacon, N. R., et al. 2007, MNRAS, 380, 712 [NASA ADS] [CrossRef] [Google Scholar]
 Lodieu, N., Deacon, N. R., & Hambly, N. C. 2012, MNRAS, 422, 1495 [NASA ADS] [CrossRef] [Google Scholar]
 Loinard, L., Torres, R. M., Mioduszewski, A. J., et al. 2007, ApJ, 671, 546 [NASA ADS] [CrossRef] [Google Scholar]
 Loinard, L., Torres, R. M., Mioduszewski, A. J., & Rodríguez, L. F. 2008, ApJ, 675, L29 [NASA ADS] [CrossRef] [Google Scholar]
 Mädler, T., Jofré, P., Gilmore, G., et al. 2016, A&A, 595, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2016, Astron. Nachr., 337, 8 [Google Scholar]
 Makarov, V. V., & Robichon, N. 2001, A&A, 368, 873 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mamajek, E. E. 2005, ApJ, 634, 1385 [NASA ADS] [CrossRef] [Google Scholar]
 Mathieu, R. D. 1986, Highlights of Astronomy, 7, 481 [CrossRef] [Google Scholar]
 Melis, C., Reid, M. J., Mioduszewski, A. J., Stauffer, J. R., & Bower, G. C. 2014, Science, 345, 1029 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Mermilliod, J.C., Mayor, M., & Udry, S. 2009, A&A, 498, 949 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Munari, U., Dallaporta, S., Siviero, A., et al. 2004, A&A, 418, L31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Narayanan, V. K., & Gould, A. 1999, ApJ, 523, 328 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Nicolet, B. 1981, A&A, 104, 185 [NASA ADS] [Google Scholar]
 Nidever, D. L., Holtzman, J. A., Allen de Prieto, C., et al. 2015, AJ, 150, 173 [NASA ADS] [CrossRef] [Google Scholar]
 Palmer, M., Arenou, F., Luri, X., & Masana, E. 2014, A&A, 564, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pan, X., Shao, M., & Kulkarni, S. R. 2004, Nature, 427, 326 [NASA ADS] [CrossRef] [Google Scholar]
 Percival, S. M., Salaris, M., & Groenewegen, M. A. T. 2005, A&A, 429, 887 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pinsonneault, M. H., Stauffer, J., Soderblom, D. R., King, J. R., & Hanson, R. B. 1998, ApJ, 504, 170 [NASA ADS] [CrossRef] [Google Scholar]
 Raboud, D., & Mermilliod, J.C. 1998, A&A, 329, 101 [NASA ADS] [Google Scholar]
 Robichon, N., Arenou, F., Mermilliod, J.C., & Turon, C. 1999, A&A, 345, 471 [NASA ADS] [Google Scholar]
 Röser, S., & Schilbach, E. 2013, in Advancing the Physics of Cosmic Distances, ed. R. de Grijs, IAU Symp., 289, 66 [Google Scholar]
 Roeser, S., Demleitner, M., & Schilbach, E. 2010, AJ, 139, 2440 [NASA ADS] [CrossRef] [Google Scholar]
 Sarro, L. M., Bouy, H., Berihuete, A., et al. 2014, A&A, 563, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Soderblom, D. R., Nelan, E., Benedict, G. F., et al. 2005, AJ, 129, 1616 [NASA ADS] [CrossRef] [Google Scholar]
 Southworth, J., Maxted, P. F. L., & Smalley, B. 2005, A&A, 429, 645 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Stauffer, J. R., Hartmann, L. W., Fazio, G. G., et al. 2007, ApJS, 172, 663 [NASA ADS] [CrossRef] [Google Scholar]
 Stello, D., & Nissen, P. E. 2001, A&A, 374, 105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 van Leeuwen, F. 1999, A&A, 341, L71 [NASA ADS] [Google Scholar]
 van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 van Leeuwen, F. 2009, A&A, 497, 209 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Venables, W. N., & Ripley, B. D. 2002, in Modern Applied Statistics with S, fourth edition (Springer) [Google Scholar]
 Weinberger, A. J., AngladaEscudé, G., & Boss, A. P. 2013, ApJ, 762, 118 [NASA ADS] [CrossRef] [Google Scholar]
 Wenger, M., Ochsenbein, F., Egret, D., et al. 2000, A&AS, 143, 9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 White, R. J., Gabor, J. M., & Hillenbrand, L. A. 2007, AJ, 133, 2524 [NASA ADS] [CrossRef] [Google Scholar]
 Wilson, R. E. 1953, in General Catalogue of stellar radial velocity, Carnegie Institute Washington D.C. Publication, 601 [Google Scholar]
 Zacharias, N., Finch, C. T., Girard, T. M., et al. 2013, AJ, 145, 44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zwahlen, N., North, P., Debernardi, Y., et al. 2004, A&A, 425, L45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Great circle intersections
An alternative approach for investigating the CP position of a moving group is based on the idea of representing the stellar proper motions of individual cluster members by great circles on the celestial sphere and calculating their intersections. One would expect the density of great circle intersections to be higher close to the CP position and its mirrorpoint (see e.g. Abad et al. 2003; Galli et al. 2012).
As described in Sect. 3 of Galli et al. (2012), the motion of each cluster member over a great circle on the celestial sphere can be described by the polar vector p that is given by (A.1)where r = (x,y,z) = (cosαcosδ,sinαcosδ,sinδ) is the position of a star with coordinates (α,δ) and proper motion (μ_{α},μ_{δ}). The intersection of great circles is defined by p_{i} × p_{j}, where i and j denote the polar vectors of two cluster members. Thus, the number of great circle intersection points is N_{int} = n(n−1) for a sample of n stars.
We calculate the great circle intersections for the stars in our control sample 1 (see Sect. 4.1) and display the results in Fig. A.1. We note that the great circle intersections roughly define a line connecting the CP with its mirror point. The two solutions correspond to the regions of highest density of great circle intersections in the diagram. They appear as an extended region (instead of a point), because of the proper motion errors and the velocity dispersion of the cluster, which prevent perfect parallelism of the stellar motions from occurring. It is clear from this analysis that a preferred direction of motion exists and the search for the CP position can be limited to a small area of the sky to optimize the methodology outlined in Sect. 2.1. Moreover, the diagram itself confirms the existence of a moving group structure in this sky region, and it can also be used to provide an independent estimate of the CP coordinates.
In this context, one simple approach consists of running a kNN regression. In this method, the kNN algorithm (Fix & Hodges 1951; Venables & Ripley 2002) is used to predict the CP position, which is taken to be the weighted mean of the k nearest neighbors (i.e., great circle intersections). We perform a weighted regression to weight the contribution of the neighbors by the inverse of their distance. The resulting CP position obviously depends on the number of nearest neighbors k that is taken into account. To overcome this problem, we varied the value of k in the range of 100 ≤ k ≤ N_{int}, and run the kNN regression to compute the CP position as a function of k. The lower limit for k is chosen to avoid a noisy CP solution that is derived with only a few neighbors. Figure A.2 shows the distribution of the CP coordinates obtained from the kNN regression. The average CP position of (α_{cp},δ_{cp}) = (92.2°,−48.6°) ± (0.8°,0.9°) is in good agreement with the solution derived from the CPSM.
The methodology described in this section, which is based on a kNN regression analysis of the great circle intersections provides an interesting tool to detect moving group structures in stellar catalogs and to estimate their CP coordinates. It is used to support and confirm our first CP solution obtained for a control sample of the Pleiades cluster. However, we emphasize that the CP position derived directly from the CPSM is likely to be more accurate, because it takes into account the proper motion errors and the velocity dispersion of the cluster.
Fig. A.1 Diagram of great circle intersections for the stars in control sample 1. The colors indicate the regions with highest density of great circle intersections that are close to the CP position (and its mirror point). Each point (black dots) represents the intersection of two great circles in the sample. 
Fig. A.2 Histogram of the CP coordinates derived from the kNN analysis of the great circle intersections using different number of neighbors (i.e., values of k). The red solid line indicates the kernel density estimator. 
Appendix B: Additional tables
Parallaxes and spatial velocities for the sample of 64 stars with known radial velocities.
Approximate parallaxes for the sample of 1146 stars.
Parallax results for cluster members in common with other studies.
All Tables
Results of the mixture model applied to the distribution of proper motions errors in right ascension and declination.
Spatial velocity of the Pleiades cluster derived in this work from the sample of 64 stars with known radial velocities.
Kinematic parallax of the Pleiades derived from different subgroups in our sample of cluster members selected by the CPSM.
Comparison of the parallaxes derived in this work with other studies in the literature.
Comparison of the proper motions and parallaxes with TGAS results from Gaia DR1 for discrepant stars.
Parallaxes and spatial velocities for the sample of 64 stars with known radial velocities.
All Figures
Fig. 1 CP coordinates (upper panel) and associated errors (lower panels) for the TW Hydrae association as a function of the number of walkers (with 200 iteration steps) used in the MCMCbased version of the CPSM. Each point is an average of 1000 MCMC realizations of the method and the error bars indicate the scatter of the solution given for each ensemble of walkers. 

In the text 
Fig. 2 Distribution of the CP coordinates after 1000 MCMC realizations of the CPSM for the TW Hydrae association. The red solid line indicates the kernel density estimator. 

In the text 
Fig. 3 Number of positions (left panel) and timebase (right panel) used to compute the stellar proper motions for Pleiades stars in the DANCe project. The red solid lines indicate the kernel density estimator. The blue dashed lines indicate the selection criteria applied to our sample as described in Sect. 4.1. 

In the text 
Fig. 4 Distribution of the relative errors on proper motions (i.e., μ/σ_{μ}) in right ascension (left panel) and in declination (right panel). The solid lines indicate the normal components of our mixture model (with k = 4), and the dashed line indicates the observed distribution of proper motions errors. 

In the text 
Fig. 5 Distribution of radial velocities for Pleiades stars collected from the literature. 

In the text 
Fig. 6 CMD to illustrate the locus of potential binaries and multiple systems that might be present in our sample. 

In the text 
Fig. 7 Distribution of stellar positions (left panel) and proper motions (right panel) for the Pleiades stars in the DANCe catalog. The various symbols and colors indicate the nuclear members of the cluster that define our control sample 2, and the stars with individual distances reported in the literature. 

In the text 
Fig. 8 Results obtained from the CPSM (CP coordinates and number of moving group members) as a function of the velocity dispersion of the cluster. The solid and dashed lines indicate the CP solution derived from our control sample 1 at 1σ and 3σ (see Sect. 4.1). 

In the text 
Fig. 9 CP solutions (black dots) for 1000 Monte Carlo realizations of the Pleiades open cluster overlaid on the X^{2} contour levels (solid lines) indicated by the color code. The red triangle indicates the final CP solution given in Sect. 4.4. 

In the text 
Fig. 10 Distribution of the μ_{⊥} values calculated from the CP solution given in Sect. 4.4 for the 1210 cluster members. The red solid line indicates the kernel density estimator. 

In the text 
Fig. 11 Histogram of the Galactic space velocity components for the 64 cluster members with known radial velocities. The red solid line indicates the kernel density estimator. 

In the text 
Fig. 12 Comparison of the parallaxes derived from individual radial velocities and the mode of the cluster spatial velocity for the sample of 64 stars. The red solid line indicates perfect correlation. 

In the text 
Fig. 13 Probability density function for the parallaxes obtained from the sample of members with known radial velocities (64 stars) and all cluster members (1210 stars). 

In the text 
Fig. 14 Probability density function for the parallaxes of different samples of cluster members selected from their membership probability (left panel) and distance to the cluster center (right panel). 

In the text 
Fig. 15 Apparent magnitudes of cluster members as a function of the distance to the cluster center and parallaxes for individual stars. 

In the text 
Fig. 16 Comparison between the kinematic parallaxes derived in this work and parallaxes from the literature (van Leeuwen 2007; Melis et al. 2014; Mädler et al. 2016; Gatewood et al. 2000; Gaia Collaboration 2016). The black dashed line indicates the expected distribution for equal results. 

In the text 
Fig. 17 Comparison between the proper motions measurements in right ascension (left panel) and declination (right panel) from the DANCe project and the TGAS catalog. The gray squares indicate the discrepant stars listed in Table 5. 

In the text 
Fig. 18 Summary of distance measurements in the literature for the Pleiades open cluster obtained from different methods. The results presented here come from: (i) moving cluster – this work (see Sect. 5.1), Röser & Schilbach (2013), Narayanan & Gould (1999); (ii) trigonometric parallax – van Leeuwen (1999), van Leeuwen (2009), Melis et al. (2014), Soderblom et al. (2005), Gatewood et al. (2000), TGAS (Gaia Collaboration 2016) – using result reported in this work (see Sect. 6.2); (iii) orbital modeling – Groenewegen et al. (2007), Southworth et al. (2005), Zwahlen et al. (2004), Munari et al. (2004), Pan et al. (2004), David et al. (2016); (iv) twin stars Mädler et al. (2016), (v) MS and isochrone fitting – An et al. (2007), Percival et al. (2005), Stello & Nissen (2001), Pinsonneault et al. (1998), Giannuzzi (1995) and Nicolet (1981). 

In the text 
Fig. A.1 Diagram of great circle intersections for the stars in control sample 1. The colors indicate the regions with highest density of great circle intersections that are close to the CP position (and its mirror point). Each point (black dots) represents the intersection of two great circles in the sample. 

In the text 
Fig. A.2 Histogram of the CP coordinates derived from the kNN analysis of the great circle intersections using different number of neighbors (i.e., values of k). The red solid line indicates the kernel density estimator. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.