Current data on the globular cluster Palomar 14 are consistent with MOND
(Research Note)
G. Gentile^{1,2}  B. Famaey^{3,4}  G. Angus^{5,6}  P. Kroupa^{4}
1  Institut d'Astronomie et d'Astrophysique, Université Libre de Bruxelles, CP 226, Bvd du Triomphe, 1050, Bruxelles, Belgium
2 
Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281, 9000 Gent, Belgium
3 
Observatoire Astronomique, Université de Strasbourg, CNRS UMR 7550, 67000 Strasbourg, France
4 
AIfA, Universität Bonn, 53121 Bonn, Germany
5 
Dipartimento di Fisica Generale ``Amedeo Avogadro'', Università degli studi di Torino, via P. Giuria 1, 10125, Torino, Italy
6 
Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Torino, Torino, Italy
Received 15 September 2009 / Accepted 15 November 2009
Abstract
Context. Certain types of globular clusters have the very
important property that the predictions for their kinematics in the
Newtonian and modified Newtonian dynamics (MOND) contexts are
divergent.
Aims. Here, we advice regarding the recent claim that the
stellar kinematics data (using 17 stars) of the globular cluster
Palomar 14 are inconsistent with MOND.
Methods. We compare the observations to the theoretical
predictions using a KolmogorovSmirnov test, which is appropriate for
small samples.
Results. We find that with the currently available data the MOND
prediction for the velocity distribution can only be excluded with a
very low confidence level, clearly insufficient to claim that MOND is
falsified.
Key words: globular clusters: individual: Palomar 14  gravitation  methods: statistical
1 Introduction
A plethora of observational data on various astronomical scalesseem to support the idea that the amount of visible matter in the Universe is several times smaller than the total amount of matter (e.g. Hinshaw et al. 2009). The current paradigm of structure formation and evolution in the Universe is known as the Cold Dark Matter (CDM) model. However, as long as the dark matter particle has not been discovered (and its cosmological abundance confirmed), it is worth considering alternative theories to explain the current data.
For instance, MOND (modified Newtonian dynamics) was proposed by Milgrom (1983) as an alternative to galactic dark matter. MOND stipulates that below a certain gravitational acceleration cm s^{2} the actual gravitational acceleration g is stronger than expected in Newtonian gravity . Asymptotically, in MOND it reaches the value for . This allows it to naturally explain various galaxy scaling relations (e.g., Faber & Jackson 1976; Tully & Fisher 1977; McGaugh et al. 2000; McGaugh 2004; Gentile 2008; Donato et al. 2009; Gentile et al. 2009).
The effects of MOND and dark matter are however often rather degenerate and modeldependent since the gravitational potential predicted by MOND can almost always be attributed by a Newtonist to an ad hoc dark matter distribution. Objects for which the predictions of the two theories are unambiguously different are unfortunately rare. We note that galaxy clusters are not good discriminant tests: indeed, in galaxy clusters, the acceleration predicted by MOND is not large enough (e.g., Sanders 1999), but some form of hot dark matter can be added within the MOND context to make the data consistent with the predictions (e.g., Angus et al. 2009). If there were systems where MOND predicted more gravity than observed, they would make a strong case against MOND.
One example of such objects that should be (almost) devoid of cold dark matter in the CDM cosmological model, and for which MOND predicts a significantly stronger gravity than that attributable to visible matter, are tidal dwarf galaxies (TDGs). However, from the observations of three young TDGs around NGC 5291, Gentile et al. (2007) found that MOND fits the data remarkably well with zero free parameters, whereas CDM fails to explain them.
Apart from TDGs, another type of object that has recently been put forward as a distinguishing test for CDM and MOND are the globular clusters of our own Galaxy (Baumgardt et al. 2005). Of particular interest are those which are diffuse and distant from the Milky Way, to ensure that the internal acceleration probes the deep MOND regime (which is not always necessarily the case: see NGC 2419 in Baumgardt et al. 2009) and that, simultaneously, the gravitational acceleration due to the Milky Way (external field) is weak enough. In MOND the deviation from a Newtonian behaviour should start appearing around a_{0}, whereas in the ``Newtonian gravity plus CDM'' picture no discrepancy is expected since globular clusters should contain (almost) no cold dark matter.
In a recent paper, Jordi et al. (2009) analysed the data of 17 stars from the globular cluster Palomar 14 and claimed that (within the assumption of Palomar 14 being on a circular orbit) MOND is inconsistent with the observed velocity dispersion: the MOND prediction is too high, whereas the data are consistent with the Newtonian prediction (with no dark matter). They discuss two separate cases, depending on the inclusion or exclusion of a star (Star 15), which based on its lineofsight velocity is not a definite member of the cluster. However, given such a small sample size, an appropriate test should be used to discern between the hypotheses. An example of such a test is the KolmogorovSmirnov (KS) test (e.g. Soong 2004). An illustrative example of the uncertain value of the velocity dispersion measured with a small number of stars is NGC 2419: the original analysis by Olszewski et al. (1993) yielded a velocity dispersion for 12 stars of km , which was reevaluated by Baumgardt et al. (2009) with 40 stars to be km , i.e. about two sigma greater than the original. This gives a perfectly reasonable M/L in MOND, although the velocity dispersion profile, if confirmed, could be a problem for MOND (Sollima & Nipoti 2009). This will allow for a very interesting test of MOND when data points at larger radii will be obtained in NGC 2419. In any case, this is an excellent example of an underestimation of the true velocity dispersion due to a sample size originally too limited.
Therefore, the question we ask in the present research note is the following: do the Palomar 14 data really exclude MOND, given the small number of stars used in the analysis of Jordi et al. (2009)? We then also investigate the minimum number of stars needed to exclude MOND in a globular cluster similar to Palomar 14. Or, equivalently, we ask how many globular clusters like Palomar 14 would be needed to exclude MOND.
2 Method
Given the small sample size, the formal error on the velocity dispersion is not sufficient to distinguish between various models, so we use the KStest to redo the analysis of Jordi et al. (2009). The KStest compares two cumulative distribution functions (cdfs), then the maximum difference D between these two cdfs yields a Pvalue.
Here we compare the cdf of the data (separately for the samples with and without Star 15) to the cdf of a Gaussian with a dispersion equal to 1.27 km s^{1} (the MOND prediction with M/L=2, Baumgardt et al. 2005; Jordi et al. 2009). We note that there is no evidence, apart from the dynamics itself, about the centre of the Gaussian, therefore we chose the centres that minimise D (in other words, the centres that maximise the Pvalue).
The MOND prediction for the velocity dispersion was derived by Baumgardt et al. (2005). First, they calculated the internal and external field gravitational accelerations ( and , respectively) of a number of globular clusters, including Pal 14. Then, based on results obtained by Milgrom (1986), they found that the MOND prediction is simply given by the Newtonian one multiplied by , because Pal 14 is in the deep MOND regime and is larger than . Baumgardt et al. (2005) assume M/L=2 from the observed M/L ratios of globular clusters: Mandushev et al. (1991) find M/L=1.21 from 32 clusters, whereas Pryor & Meylan (1993) find M/L=2.3 from 56 clusters. Baumgardt et al. (2005) state that the latter value is more plausible because the modelling used in Pryor & Meylan (1993) takes mass segregation into account. To get an estimate of the uncertainty on the M/L, the data compilation by Dabringhausen et al. (2008) suggest that M/L ratios of globular clusters have a spread of roughly 0.5 dex.
Table 1: Velocities and cumulative distribution functions (observed and Gaussian, with and without Star 15).
Figure 1: Cumulative distribution function (cdf) vs. radial velocity for the sample with Star 15. The solid steplike line represents the data, whereas the dotted line is the cdf of a Gaussian with a standard deviation of 1.27 km s^{1} (the MOND prediction) and the centre at 72.30 km s^{1}. The red segment shows D, the maximum difference between the two cdfs. 

Open with DEXTER 
3 Results
In Table 1 we list the observed and predicted cdfs (Gaussians with a standard deviation of 1.27 km s^{1}) separately for the sample with and without Star 15. The two predicted cdfs are different because they have different centres (72.30 km s^{1} for the sample with Star 15 and 72.41 km s^{1} for the sample without Star 15). We also show them in Figs. 1 and 2. In both cases, there are two regions where the difference almost reaches the maximum: around 71.75 km s^{1} and around 72.65 km s^{1}, respectively. And in both samples, the difference between the two observed and predicted cdf are very similar at these two velocities. With Star 15, the maximum difference is 0.215, and without Star 15 the maximum difference is 0.239. These values of D correspond to Pvalues of 0.360 and 0.273, respectively. This means that using the KStest, the data presented in Jordi et al. (2009) can exclude the MOND with M/L=2 hypothesis only with 64% and 73% confidence, depending on the inclusion of Star 15. These confidence levels are clearly not sufficient to exclude MOND. We note that if the M/L is not 2, but rather 1.1 (a possibility considered by Jordi et al. 2009), then MOND is perfectly consistent with a velocity dispersion of 0.85 km s^{1}.
Figure 2: Cumulative distribution function (cdf) vs. radial velocity for the sample without Star 15. The lines are the same as in Fig. 1, but the Gaussian's centre is at 72.41 km s^{1}. 

Open with DEXTER 
Figure 3: Mean Pvalue vs. number of stars in the sample for a number of MonteCarlo realisations (10 per sample size) of distributions of stars following a Gaussian with a standard deviation of 0.85 km s^{1} (equivalent to the sample with Star 15). 

Open with DEXTER 
Now we can ask ourselves: if one wanted a Pvalue of 0.1 (i.e. an exclusion confidence of 90%), how many stars would one need (with the hypothesis that they follow a Gaussian with a standard deviation equal to the measured one)? We measure a dispersion of 0.85 km s^{1} with Star 15 and 0.63 km s^{1} without Star 15. Hence, we create a series of MonteCarlo realisations, using the software package R (Ihaka & Gentleman 1996), of an increasing sample size. We create 10 realisations per sample size. Then, we compare these mock data sets to the MOND prediction using the same method as above, take the mean Pvalue for each sample size, and look for the minimum number of stars necessary to obtain . Our results (shown in Figs. 3 and 4) are that for the mock data sets equivalent to the sample with Star 15 (i.e., the mock data sets that follow a Gaussian with a standard deviation = 0.85 km s^{1}) a minimum of about 80 stars (or alternatively 80/17 5 clusters) would have been needed to exclude the MOND hypothesis with a confidence level of 90%, whereas without Star 15 this minimum number of stars decreases to about 30.
4 Discussion
In this research note, we showed that current observational data on Pal 14 are not significantly discrepant with the theoretical prediction of MOND. Let us however also note that this theoretical prediction is not really unique and depends on many factors, such as stellar M/L, rotation and anisotropy. For instance, in Angus (2008, thesis), the velocity dispersion of Pal 14 was discussed before any velocity dispersion was published: the line of sight velocity dispersion for various M/Ls and velocity anisotropies were computed (radial, isotropic and centrally isotropic with increasing tangentially biased orbits). It was found that lower velocity dispersions (consistent with 0.85 km ) were attainable with very radial orbits, and in that case the line of sight velocity dispersion profile would moreover not be flat in the outskirts. Rotation could also play a role although the round appearance of Pal 14 justifies the norotation assumption. Of course, the most important uncertainty comes from the M/L ratio itself.
Figure 4: Same as in Fig. 3, but for a Gaussian with a standard deviation of 0.63 km s^{1} (equivalent to the sample without Star 15). 

Open with DEXTER 
Then, another possible oversight may be the mass segregation of stars within the globular cluster: if significant energy partitioning has occured due to the short MOND relaxation time (Ciotti & Binney 2004), this would mean that the lowmass, unobserved stars could have a larger spread in configuration space and a higher velocity dispersion. Also, as already noted in Jordi et al. (2009), the MOND theoretical prediction was based on assuming a purely circular orbit for Pal 14 around the Milky Way, while if it is on an eccentric orbit, it could (i) have lost many low mass stars at perigalacticon leading to a decrease of the theoretical stellar M/L, and (ii) be ``frozen" in the Newtonian regime due to the period of recovery while transiting from the large external field endured at perigalacticon. Finally, we note that the fact that MOND cannot be excluded also implies, in the Newtonian plus dark matter cosmological framework, that the current data of Pal 14 cannot exclude a certain amount of unseen mass in the cluster. With the assumption of a stellar M/L of 2, a velocity dispersion of 1.27 km s^{1} would imply a totaltoluminous mass ratio of 2 to 4.
5 Conclusion
Even assuming an isotropic MOND model with no rotation for a globular cluster on a circular orbit with M/L=2, we showed that, based on a KS test, which is the relevant statistical test for small samples, the currently available data are insufficient to distinguish between Newtonian gravity and MOND in Palomar 14, contrary to the claim of Jordi et al. (2009). While the objects proposed by Baumgardt et al. (2005) provide one of the best distinguishing tests between MOND and cold dark matter plus Newtonian dynamics, more observations would be needed to exclude the aforementioned MOND model if the observed velocity dispersion is representative of the true one: about 80 stars in Pal 14, or about five similarly problematic globular clusters for this model. In this respect, even though it might still not be conclusive, velocity data of 21 stars in Pal 3 and 24 stars in Pal 4 will be of prime interest for testing fundamental physics in our Galactic backyard.
AcknowledgementsG.G. is a postdoctoral researcher of the FWOVlaanderen (Belgium). B.F. acknowledges support from the Hulmboldt foundation and the CNRS. G.W.A. is supported by the Univsersity of Torino and Regione Piemonte, and by the INFN grant PD51.
References
 Angus, G., Famaey, B., & Diaferio, A. 2009, MNRAS, in press[arXiv:0906.3322] [Google Scholar]
 Baumgardt, H., Grebel, E. K., & Kroupa, P. 2005, MNRAS, 359, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Baumgardt, H., Côté, P., & Hilker, M. 2009, MNRAS, 396, 2051 [NASA ADS] [CrossRef] [Google Scholar]
 Ciotti, L., & Binney, J. 2004, MNRAS, 351, 285 [NASA ADS] [CrossRef] [Google Scholar]
 Dabringhausen, J., Hilker, M., & Kroupa, P. 2008, MNRAS, 386, 864 [NASA ADS] [CrossRef] [Google Scholar]
 Donato, F., Gentile, G., & Salucci, P. 2009, MNRAS, 397, 1169 [NASA ADS] [CrossRef] [Google Scholar]
 Faber, S., & Jackson, R. 1976, ApJ, 204, 668 [NASA ADS] [CrossRef] [Google Scholar]
 Gentile, G. 2008, ApJ, 684, 1018 [NASA ADS] [CrossRef] [Google Scholar]
 Gentile, G., Famaey, B., & Combes, F. 2007, A&A, 472, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gentile, G., Famaey, B., Zhao, H. S., & Salucci, P. 2009, Nature, 461, 627 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Harris, W. E., & van den Bergh, S. 1984, AJ, 89, 1816 [NASA ADS] [CrossRef] [Google Scholar]
 Hilker, M. 2006, A&A, 448, 171 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hinshaw, G., Weiland, J. L., & Hill, R. S. 2009, ApJS, 180, 225 [NASA ADS] [CrossRef] [Google Scholar]
 Holland, S., & Harris, W. E. 1992, AJ, 103, 131 [NASA ADS] [CrossRef] [Google Scholar]
 Ihaka, R., & Gentleman, R. 1996, Journal of Computational and Graphical Statistics, 5(3), 299, URL: http://www.rproject.org/ [CrossRef] [Google Scholar]
 Jordi, K., Grebel, E. K., & Hilker, M. 2009, AJ, 137, 4586 [NASA ADS] [CrossRef] [Google Scholar]
 Mandushev, G., Staneva, A., & Spasova, N. 1991, A&A, 252, 94 [NASA ADS] [Google Scholar]
 McGaugh, S. S. 2004, ApJ, 609, 652 [NASA ADS] [CrossRef] [Google Scholar]
 McGaugh, S. S., Schombert, J. M., Bothun, G. D., & de Blok, W. J. G. 2000, ApJ, L99 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Milgrom, M. 1983, ApJ, 270, 365 [NASA ADS] [CrossRef] [Google Scholar]
 Milgrom, M. 1986, ApJ, 302, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Olszewski, E. W., Pryor, C., & Shommer, R. B., ASP Conf. Ser., 48, 99 [Google Scholar]
 Pryor, C., & Meylan, G. 1993, in Structure and Dynamics of Globular Clusters, ASP Conf., 50, 357 [Google Scholar]
 Sanders, R. H. 1999, ApJ, 512, L23 [NASA ADS] [CrossRef] [Google Scholar]
 Sollima, A., & Nipoti, C. 2009, MNRAS, 401, 131 [NASA ADS] [CrossRef] [Google Scholar]
 Soong, T. T. 2004, Fundamentals of Probability and Statistics for Engineers (Publisher: John Wiley & Sons) [Google Scholar]
 Tully, R. B. & Fisher, J. R. 1977, A&A, 54, 661 [NASA ADS] [Google Scholar]
All Tables
Table 1: Velocities and cumulative distribution functions (observed and Gaussian, with and without Star 15).
All Figures
Figure 1: Cumulative distribution function (cdf) vs. radial velocity for the sample with Star 15. The solid steplike line represents the data, whereas the dotted line is the cdf of a Gaussian with a standard deviation of 1.27 km s^{1} (the MOND prediction) and the centre at 72.30 km s^{1}. The red segment shows D, the maximum difference between the two cdfs. 

Open with DEXTER  
In the text 
Figure 2: Cumulative distribution function (cdf) vs. radial velocity for the sample without Star 15. The lines are the same as in Fig. 1, but the Gaussian's centre is at 72.41 km s^{1}. 

Open with DEXTER  
In the text 
Figure 3: Mean Pvalue vs. number of stars in the sample for a number of MonteCarlo realisations (10 per sample size) of distributions of stars following a Gaussian with a standard deviation of 0.85 km s^{1} (equivalent to the sample with Star 15). 

Open with DEXTER  
In the text 
Figure 4: Same as in Fig. 3, but for a Gaussian with a standard deviation of 0.63 km s^{1} (equivalent to the sample without Star 15). 

Open with DEXTER  
In the text 
Copyright ESO 2010