Issue 
A&A
Volume 527, March 2011



Article Number  A33  
Number of page(s)  7  
Section  Galactic structure, stellar clusters and populations  
DOI  https://doi.org/10.1051/00046361/201015573  
Published online  21 January 2011 
Distant star clusters of the Milky Way in MOND
^{1}
Institute for Advanced Studies in Basic Sciences (IASBS),
PO Box 113659161,
Zanjan,
Iran
email: haghi@iasbs.ac.ir
^{2}
Argelander Institute for Astronomy (AIfA),
Auf dem Hügel 71,
53121
Bonn,
Germany
email: pavel@astro.unibonn.de
^{3}
School of Mathematics and Physics, University of Queensland,
Brisbane,
QLD
4072,
Australia
email: h.baumgardt@uq.edu.au
Received:
12
August
2010
Accepted:
1
December
2010
We determine the mean velocity dispersion of six Galactic outer halo globular clusters, AM 1, Eridanus, Pal 3, Pal 4, Pal 15, and Arp 2 in the weak acceleration regime to test classical vs. modified Newtonian dynamics (MOND). Owing to the nonlinearity of MOND’s Poisson equation, beyond tidal effects, the internal dynamics of clusters is affected by the external field in which they are immersed. For the studied clusters, particle accelerations are much lower than the critical acceleration a_{0} of MOND, but the motion of stars is neither dominated by internal accelerations (a_{i} ≫ a_{e}) nor external accelerations (a_{e} ≫ a_{i}). We use the Nbody code NMODY in our analysis, which is a particlemeshbased code with a numerical MOND potential solver developed by Ciotti et al. (2006, ApJ, 640, 741) to derive the lineofsight velocity dispersion by adding the external field effect. We show that Newtonian dynamics predicts a lowvelocity dispersion for each cluster, while in modified Newtonian dynamics the velocity dispersion is much higher. We calculate the minimum number of measured stars necessary to distinguish between Newtonian gravity and MOND with the KolmogorovSmirnov test. We also show that for most clusters it is necessary to measure the velocities of between 30 to 80 stars to distinguish between both cases. Therefore the observational measurement of the lineofsight velocity dispersion of these clusters will provide a test for MOND.
Key words: gravitation / methods: numerical / globular clusters: general
© ESO, 2011
1. Introduction
Observable matter in galaxies and in clusters of galaxies cannot produce sufficient gravity to explain their dynamics. Cold dark matter (CDM) scenarios or alternative theories of gravitation are therefore invoked to resolve the problem. Nowadays, the CDM hypothesis is the dominant paradigm. This hypothetical matter does not interact with electromagnetic radiation and only shows its presence through its gravitational interaction. Even though the dark matter hypothesis has successfully explained the internal dynamics of galaxy clusters, gravitational lensing, and the standard model of cosmology within the framework of general relativity (GR) (Spergel 2003), much experimental effort has failed to yield a detection of dark matter particles. Moreover the results of highresolution simulations of structure formation do not reproduce some observations on galactic scales, such as the central structures of rotation curves, the prevalence of low bulgetodisc ratios, and the numbers and spatial distribution of the subhalos (Klypin et al. 1999; Moore et al. 1999; Metz et al. 2008; Kroupa et al. 2010). Even the ability of the dark matter theory to account for the TullyFisher and Freeman relations is controversial (Bosch & Dalcanton 2000; Governato et al. 2010). These shortcomings have not led to the rejection of the theory only because on galactic scales baryons are at least nonnegligible contributors to the mass density, consequently simulations that include the complex physics of star formation are essential for reliable predictions. Currently such simulations are still at an experimental stage and are usually substituted by “semianalytic” arguments that have weak theoretical underpinnings.
One of the alternative theories to CDM is the socalled modified Newtonian dynamics (MOND) theory, which was originally proposed by Milgrom (1983) to explain the flat rotation curves of spiral galaxies at large distances by a modification of Newton’s second law of acceleration below a characteristic scale of a_{0} ≃ 1.2 × 10^{10} m s^{2} = 3.6 pc/(Myr)^{2} without invoking dark matter (Bekenstein & Milgrom 1984). In MOND, Newton’s second law is modified to μ(a/a_{0})a = a_{N} + ∇ × H, where ρ is the massdensity distribution, a_{N} is the Newtonian acceleration vector, a is the MONDian acceleration vector, a = a is the absolute value of MONDian acceleration, μ is an interpolating function for the transition from the Newtonian to the MONDian regime, which runs smoothly from μ(x) = x at x ≪ 1 to μ(x) = 1 at x ≫ 1 (Bekenstein & Milgrom 1984). Different interpolating functions have been suggested, such as the simple function, μ(x) = x/(1 + x) (Famaey & Binney 2005) and the standard interpolation function, (Milgrom 1983). Because the simple function fits galactic rotation curves better than the standard function (Gentile et al. 2011), in this paper we use the simple function. The value of the curl field H depends on the boundary conditions and the spatial mass distribution and vanishes only for some special symmetries (Bekenstein & Milgrom 1984). The nonlinearity of the MOND field equation leads to difficulties for standard Nbody codes and makes the use of the usual Newtonian Nbody simulation codes impossible in the MOND regime.
It has been shown that on galactic scales MOND can explain many phenomena better than CDM (Begmann 1989; Begmann et al. 1991; Sanders & McGaugh 2002; SanchezSalcedo & Hernandez 2007; Haghi et al. 2006; Malekjani et al. 2009; Gentile et al. 2007; Milgrom 1994; Brada & Milgrom 2000; Milgrom 1995; Wu et al. 2008; Zhao & Famaey 2006; Tiret & Combes 2007; Hasani et al. 2010). MOND has been generalized to a generalrelativistic version (Bekenstein 2004; Sanders 2005; Zlosnik et al. 2007; Milgrom 2009), making it possible to test its predictions for gravitational lensing. Dynamics of galaxies in clusters (2002) and the merging of galaxy clusters, where the baryonic mass is clearly separated from the gravitational mass (Clowe et al. 2006) cannot be completely explained by MOND without invoking some kind of hot dark matter, perhaps in the form of massive (active or sterile, 2 to 11 eV) neutrinos (Angus et al. 2006, 2010).
In order to decide whether MOND is a comprehensive theory to explain the dynamics of the universe, it is desirable to study MOND for objects in which no dark matter is supposed to exist and where the characteristic acceleration of the stars is less than the MOND critical acceleration parameter a_{0}. Globular clusters (GCs) are a perfect candidate since they are the largest virialized structure that do not contain dark matter (Moore 1996).
In the distant halo of our Milky Way there exist several lowmass GCs where both internal and external accelerations of stars are significantly below the critical acceleration parameter a_{0} of MOND. Because GCs are assumed to be darkmatterfree, if MOND is true, the motions of stars must deviate from the standard Newtonian dynamics. It has been proposed by Baumgardt (2005) that some of these distant Galactic GCs are perfect tools to test gravitational theories in the regime of very weak accelerations. For MOND, the internal velocity dispersion among the stars in these clusters would be significantly higher than in Newtonian dynamics.
The mean velocity dispersion of stellar systems for the two extreme cases of internal (a_{i} ≫ a_{e}) or external (a_{e} ≫ a_{i}) field domination have been derived analytically by Milgrom (Milgrom 1986, 1994), assuming that the systems are everywhere in the deepMOND regime (a_{e},a_{i} ≪ a_{0}). Many systems that can be used to test MOND are not completely either internally or externally dominated. Globular clusters or dwarf galaxies of the Milky Way for example have internal and external accelerations that are of the same order (Baumgardt et al. 2005), consequently one has to determine the velocity dispersion numerically for intermediate cases. Sollima & Nipoti (2010) constructed selfconsistent, spherical models for stellar systems in MOND, neglecting the external field effect and presented a dynamical model for six galactic GCs. The presence of the external field effect breaks the spherical symmetry and validity of their model. Recently Haghi et al. (2009, hereafter HBK09) investigated the dynamics of star clusters by numerically modeling them in MOND, assuming circular orbits. They performed Nbody simulations and presented analytical formulae for the velocity dispersion of stellar systems in the intermediate MOND regime, which are useful for a comparison with observational data of several GCs and dSph galaxies (for details on the numerical calculations see HBK09). In a followup paper, Jordi et al. (2009) determined the velocity dispersion (using 17 stars) and massfunction slope of Pal 14 and showed that MOND can hardly explain the lowvelocity dispersion of this system. However, Gentile et al. (2010) showed that with the currently available data, the KolmogorovSmirnov (KS) test is still unable to exclude MOND with a sufficiently high confidence level. Moreover, the low density of Pal 14 suggests that binary stars may be an important issue for interpreting its measured velocity dispersion (Küpper & Kroupa 2010), and the true velocity dispersion of Pal 14 could be much lower than the value reported by Jordi et al. (2009), thereby possibly posing an even larger challenge for MOND, but also for Newtonian gravity and for any understanding of the dynamics of this object as being in equilibrium.
In this paper we calculate the prediction of MOND and Newtonian dynamics on the velocity dispersion of six other distant clusters of the MW (Table 1). In order to see the pure MONDian effects, we concentrate on systems in which the tidal radius is much larger than the gravitational radius^{1} and therefore tidal effects are unimportant. In other words, this paper provides the basis for further observational efforts. The measurements of a low (Newtonian) velocity dispersion would mean that MOND in its present form is in severe trouble and that globular clusters do not possess dark matter. In contrast, a highvelocity dispersion would either favor MOND or could be a hint to the existence of dark matter in GCs (Baumgardt et al. 2009).
The paper is organized as follows: in Sect. 2 we give a brief review of the external field effect (EFE) in MOND. The simulation setup is explained in Sect. 3. The numerical results for six clusters are discussed in Sect. 4. We present our conclusions in Sect. 5.
2. External field effect in MOND
In classical Newtonian dynamics, a uniform external field does not affect the internal dynamics of a stellar system. In other words, in the frame of the system relative motions of objects are the same as in an isolated system.
In MOND, the situation is entirely different. Owing to the nonlinearity of Poisson’s equation, the strong equivalence principle (SEP) is violated (Bekenstein & Milgrom 1984), and consequently the internal properties and the morphology of a stellar system are affected by both the internal and external field. This socalled external field effect (EFE) significantly affects nonisolated systems and can provide a strict test for MOND.
The EFE is indeed a phenomenological requirement of MOND and was postulated by Milgrom (1983) to explain the dynamical properties of open clusters in the MW that do not show MONDian effects. The EFE allows highvelocity stars to escape from the potential of the Milky Way (Famaey et al. 2007; Wu et al. 2007) and implies that rotation curves of spiral galaxies should fall where the internal acceleration becomes equal to the external acceleration (Gentile et al. 2007; Wu et al. 2008).
For a star cluster with a density distribution ρ_{c}, which is embedded in a host galaxy with a density distribution ρ_{ext}, the acceleration of stars in the cluster satisfies the modified Poisson equation (1)where ∇Φ is the MONDian potential generated by the total matter density. For star clusters or dwarf galaxies far out in the halo of the Milky Way the local density of the Milky Way is negligible (i.e. ρ_{ext} ≪ ρ_{c}). One way to solve Eq. (1) is then to assume that the total acceleration is the sum of the internal a_{i} and the external a_{e} acceleration, which both satisfy the modified Poisson equation as (Bekenstein & Milgrom 1984; Wu et al. 2007, 2008; Angus 2008) (2)In Eq. (2), the direction of the external acceleration is important for the dynamics of stars in the clusters. For example, for two stars the acceleration would be different if a_{e} is parallel or antiparallel to their internal acceleration.
Several attempts have been made considering the EFE. Wu et al. (2007) simulated an isolated object with a static MONDian potential solver, but changed the boundary condition on the outermost grid point to be nonzero. Famaey et al. (2007) estimated the escape velocity of galaxies in MOND and assumed μ(a_{e} + a_{i}/a_{0})a_{i} = a_{N}. As a first order test they replaced a_{e} + a_{i} with (a_{e} + a_{i}), or . This is an approximation, since they neglected the direction of a_{e} (i.e. a possible angular difference between a_{i} and a_{e}).
In order to study the evolution of the systems we need to calculate the acceleration at each step by adding the constant external gravitational field to the internal acceleration inside the μfunction, considering the different angle between external and internal acceleration for all stars throughout the evolution.
3. Simulation setup
We performed a large set of Nbody simulations of star clusters with the NMODY code, which has been developed by Ciotti et al. (2006). NMODY is a parallel, threedimensional particlemesh code for the timeintegration of collisionless Nbody systems (Londrillo & Nipoti 2009). The code numerically solves the nonlinear MOND field equations, which can be used to perform numerical experiments in either MONDian or Newtonian dynamics. The potential solver of NMODY is based on a grid in spherical coordinates and is best suited for modeling isolated systems. NMODY uses the leapfrog method to advance the particles. The code and the potential solver have been presented and tested by Ciotti et al. (2006) and Nipoti et al. (2007).
In the present study we used a spherical grid (r,θ,ϕ) made of N_{r} × N_{θ} × N_{ϕ} = 128 × 128 × 128 grid cells for the integration. The total number of particles is N_{p} = 10^{5}. The details of the scaling of the numerical MOND models and code units are discussed in Nipoti et al. (2007).
Our treatment of the external field effect is based on the idea that the star clusters we study are much smaller than they are distant from the Milky Way, so the external field is nearly constant over the cluster area. We carry out the simulations in a rotating reference frame centered on the cluster. We also assume that the clusters are on a circular orbit, so the external field of the Milky Way is also constant with time in the clusters frame.
We assume that tidal forces arising from a gradient of the external field and the Coriolis acceleration caused by the rotating reference frame can be neglected. This is justified as long as the size of the cluster is much smaller than the tidal radius. We then solve Eq. (2) numerically with NMODY, using a constant acceleration a_{e} as the boundary condition.
In order to include the EFE for nonisolated systems, we applied several changes to the NMODY code. The changes were encoded in the source file mondlib.f90, which contains all relevant routines implementing the MOND potential solver. The Mondian potentials and accelerations are assigned to stars in subroutines: “gmond” and “mondacc”, and we placed the vector of a_{e} in the interpolating function for all stars, i.e., a_{e} is added to a_{i} via vector summation and the code solves Eq. (2) by iteration. At each step of potential solving, we added the constant external field for all grid points according to Eq. (2). This is a first approximation that allows us to focus on the effects of a constant external field, which in turn allows us to deal with the intermediate MOND regime. We would like to stress that our numerical solution agrees very well with analytic ones for the extreme cases a_{i} ≪ a_{e} ≪ a_{0}, a_{e} ≪ a_{i} ≪ a_{0} and a_{e},a_{i} ≫ a_{0} (HBK09).
Our method is a first attempt to take into account the constant EFE in NMODY, in order to avoid solving the modified Poisson equation for the large area including both galaxy and cluster, which is impossible with the current version of NMODY. This can only be studied with a realistic, highresolution MOND simulation that includes the density distribution causing the external field.
In this work we start from Newtonian equilibrium Plummer models and produce MONDian equilibrium initial systems. In order to have a MONDian equilibrium initial system, we increased the initial velocity of the particles. This method is useful when the external field is important. The details of the method are described in HBK09.
Fig. 1
Lineofsight velocity dispersion for clusters listed in Table 1 for various masses as found by NMODY. In order to compare with the real cluster, the halfmass radii of all models are fixed. For low masses, which mean low internal accelerations, the prediction of Newtonian dynamics differs from the numerical solution. As the mass increases (at about 10^{6} M_{⊙}), the internal acceleration grows and the system enters the Newtonian regime and the numerical solutions get close to the Newtonian prediction. In order to focus on a reasonable range of mass for these clusters, we show the results in the range of 500−2 × 10^{5}. Dashed vertical lines show the expected cluster mass, calculated from the absolute luminosity by assuming a stellar masstolight ratio of M/L_{V} = 2. This value agrees with observed masstolight ratios of galactic GCs (Mieske et al. 2008). 
Globular clusters modeled in this paper.
4. Modeling the distant globular clusters
We calculated the velocity dispersion for the Galactic GCs mentioned in Table 1. These clusters are best suited for testing MOND. A large set of dissipationless NMODY computations with MONDian equilibrium initial conditions were performed for stellar systems that are embedded in the outer MW halo and affected by different values of the external field. Clusters in Table 1 are generally far out in the Galactic halo so that the external acceleration caused by the MW is small. In order to produce different internal acceleration regimes, we changed the cluster mass from 500 M_{⊙} to 10^{7} M_{⊙} and assume the halfmass radius to be constant. The models are evolved for several crossing times to reach the equilibrium state. The EFE on the predicted lineofsight velocity dispersion for nonisolated stellar systems with different internal accelerations have been carried out in HBK09.
For the GCs listed in Table 1, observational efforts are underway to determine their velocity dispersion and to constrain their mass (Grebel et al. 2009). Because the halfmass relaxation times of most clusters in Table 1 are on the order of a Hubble time or even larger, dynamical evolution does not play an important role for these clusters, therefore we can assume that the halflight radius is equal to the halfmass radius. The projected halfmass radii, R_{hp}, and galactocentric distances, R_{G}, are taken from Harris (1996). The projected halfmass radii relate to the threedimensional halfmass radii as R_{hp} = γR_{h} with γ ≈ 0.74. For all clusters in Table 1, the tidal radii are larger than the halfmass radii by a factor of 5 to 10, which means that tidal effects play no significant role for the internal dynamics of the clusters in our sample. We perform numerical modeling to obtain the mean velocity dispersion.
Fig. 2
Top panel: the ratio of the velocity dispersion in the direction of the external field compared to the perpendicular direction for Pal 3 as an example plotted against the radius. The system is isotropic throughout the cluster. Bottom panel: the anisotropy profile, β, for Pal 3 as an example plotted against the radius. The system is isotropic throughout the cluster. 
Fig. 3
Upper panel: radial velocity distribution of stars for a Newtonian (red triangles) and MONDian (black dots) system in which both have the same internal velocity dispersion. The smoothed line shows a Gaussian fitted to the distributions. Both distributions follow the Gaussian fit very closely. Lower panel: the relative difference between Newtonian and Mondian velocity distributions. There is a difference between the two models, but it is too small (less than 1%) to be detected observationally. 
Because the cluster masses are not known from observations, it is useful to calculate lineofsight velocity dispersions for different values of cluster mass. If the velocity dispersion is determined observationally, one can constrain the cluster mass. In Fig. 1 the resulting global lineofsight velocity dispersion as a function of mass is plotted for the six stellar clusters and is compared with the Newtonian results. The black line shows the Newtonian prediction for the velocity dispersion, which we calculate using (Haghi et al. 2009). The red line shows the numerical calculation of the velocity dispersion in the MONDian regime. For each cluster the external acceleration of the Galaxy is given in Table 1. All clusters are in the intermediate regime, for which there is no analytical prediction in MOND. In the lowacceleration region, the numerical solutions show a considerable relative difference with Newtonian results. As the mass increases, the internal acceleration grows and then gets close to the Newtonian results (at about M = 10^{6} M_{⊙}). The vertical dashed lines show the expected mass for each cluster calculated from the absolute luminosities by assuming a stellar masstolight ratio of M/L = 2 (Mieske et al. 2008). It should be noted that in the simulations we do not assume any M/L, but we vary the mass and therefore the M/L ratio of the clusters. The velocity dispersions corresponding to this mass for both cases are indicated in Table 1. The M/L ratio of 2 is only used in Table 1 to illustrate there is a MOND effect, which makes a difference for the studied clusters. Our predictions for the MONDian case are lower than those by Baumgardt et al. (2005). The absolute differences between MONDian and Newtonian results are largest for the clusters Pal 4, Pal 3, AM 1, and Eridanus. Therefore these clusters would be the best cases to test MOND. The same qualitative results were also obtained by Sollima & Nipoti (2009), who have studied some of the clusters of our list (i.e., Eridanus, AM1, Pal 3, Pal 4, and Pal 14, which we have studied in HBK09). They have obtained the cluster’s velocity dispersion profile. Their predictions for the velocity dispersion are generally higher than ours. For AM1 for example, the mean value of the velocity dispersion is ~1.7 km s^{1} since the predicted MOND velocity dispersion profile is flat throughout the cluster area (Fig. 5 of Sollima & Nipoti 2009). This is higher than the value of 1.5 km s^{1} we find, which is to be expected because they did not include the EFE. If we take the EFE into account, the systems tend to be in the quasiNewotonian regime, where the gravity is weaker than in the deepMondian regime (HBK09).
In MOND the external field introduces an anisotropy so that clusters could in principle become elongated. In order to test for a possible elongation of the velocity dispersion in the direction of the external field (σ_{x}), we define the parameter . Figure 2 shows that the lineofsight velocity dispersion in the direction of the external field differs only randomly compared to that perpendicular to the external field. This small a difference would be unobservable.
Fig. 4 Mean Pvalue vs. number of stars in the sample for a number of random realizations of distributions of stars following a Gaussian with a standard deviation of σ_{N} for clusters listed in Table 1. The minimum number of stars necessary to reject MOND is defined where the mean Pvalue drops below 0.05. 
Figure 3 shows the distribution of radial velocities of the N_{total} = 10^{5} stars in the velocity range –10 to +10 km s^{1} for both Newtonian (red triangles) and Mondian (black dots) cases together with the best Gaussian fit to this distribution using the least square method. We created the Newtonian model with a mass of 177 000 M_{⊙} and the Mondian one with a mass of 320 00 M_{⊙}, so that the velocity dispersions are the same. The bottom panel shows the relative deviation, (N_{Newtonian} − N_{Mond})/N_{total}, where N is the number of stars in a velocity bin. Evidently the observed radial velocity distribution is well approximated by this Gaussian. The Newtonian one is slightly different from the Mondian distribution, but the relative difference is very small (less than 1%), therefore it is impossible to find a detectable deviation observationally.
5. Distinguishing between MOND and Newtonian gravity models
Usually, it is not feasible to measure the velocity dispersion of a large sample of stars in distant and lowmass GCs. If only a few stars are measured, the observed velocity dispersion has a large error, which makes it difficult to distinguish between MOND and Newtonian dynamics. Recently, Gentile et al. (2010) have shown that the small sample size for Pal 14 (i.e., radial velocity of 17 stars were used to determine the observed velocity dispersion) can only reject MOND with a low confidence level.
In order to see how many stars are necessary to reject MOND with high statistical confidence if the underlying velocity dispersion is Newtonian, we created artificial data sets of increasing sample size from the Newtonian simulations for each of the six clusters. For each sample size we produced 20 random realizations of a mock sample, and each time applied a KS test with the null hypothesis being that MOND is correct (thus comparing the mock data set from the Newtonian simulation with the velocity dispersion of the MONDian simulation, σ_{MOND}). We approximate the cluster mass from the Vband luminosity by assuming a stellar masstolight ratio of M/L = 2 (vertical dashed line in Fig. 1). Figure 4 shows the mean Pvalue from the KS test as a function of the number of stars in the sample. The number of stars changes in the range of [10–160]. For any given cluster and a fixed number of stars a 1sigma range of Pvalues are shown as error bars. The minimum number of stars, N_{min}, which is necessary to achieve P < 0.05 (i.e. exclusion confidence of 95%), are obtained for each cluster and are given in Table 1. Measurement of N_{min} stars should therefore be enough to rule out MOND if the cluster mass is simultaneously determined by independent measurements, possibly through e.g. star counts.
6. Conclusions
We compared global lineofsight velocity dispersions of six distant lowdensity GCs of the Milky Way in MONDian and Newtonian dynamics and showed that they have a significantly higher velocity dispersion in MOND than the prediction of Newtonian dynamics.
Using the NMODY code, we obtained a large set of dissipationless numerical solutions for GCs with the MONDian initial conditions as endproducts of the Nbody computations by adding the external field. In order to produce different internal acceleration regimes, we changed the mass of the system over a wide range for each cluster. In addition, our results show that the clusters Pal 4, Pal 3, AM 1, and Eridanus are the best cases to test MOND owing to their larger absolute difference between MONDian and Newtonian velocity dispersions. These results will allow us to test MOND more rigorously than was possible so far, and this will enable us to compare MOND with observational data in the future.
Recently Sollima & Nipoti (2010) have performed a test for MOND by constructing selfconsistent dynamical models for outer galactic clusters. These authors neglected the external field effect of the Milky Way, and hence obtained higher estimates for the velocity dispersion of clusters compared with our results. This difference is reasonable. Indeed, the external field in MOND is leading the system to move from the deepMONDian regime toward the quasiNewtonian regime, and thus to reduce the velocity dispersion. In other words, the larger the external field, the smaller the internal acceleration, which implies a lower velocity dispersion for the cluster (Haghi et al. 2009).
It should be noted that the clusters could be on eccentric orbits, which means that the MOND predictions would be different because of the variation of the external field (in direction and amplitude) along the orbit.
Using a KS test, we calculated the minimum number of stars that are sufficient to exclude MOND (under the hypothesis that these GCs are on circular orbits) at the 95% confidence level. We found that between 30 to 80 stars are necessary for most clusters to distinguish between both cases. This number of stars can be observed with current 8 m class telescopes. Additional observational efforts to determine the velocity dispersions of these clusters and constraining the mass of the clusters by star counts would be highly important and provide a strict test of MOND. On the other hand, if MOND is the correct theory, these observations could be used to constrain the external field and consequently to put constraints on the potential in which the systems are embedded. According to the anisotropy profile, the simulated systems are isotropic throughout.
The gravitational radius is a measure of the size of the system and is related to the mass and potential energy as given in Eqs. (2)–(132) in Binney & Tremaine (1987). In many stellar systems the gravitational radius can be approximated by the threedimensional halfmass radius r_{h} as r_{g} = 1.25r_{h} if the assumption of virial equilibrium is valid.
Acknowledgments
We would like to thank C. Nipoti for providing us with the NMODY code and his help in using it. H.H. thanks the stellar dynamics group of the Argelander Institute for Astronomy for giving financial support for this research. H.B. acknowledges support from the German Science foundation through a Heisenberg Fellowship and from the Australian Research Council through Future Fellowship grant FT0991052.
References
 Angus, G. W. 2008, MNRAS, 387, 1481 [NASA ADS] [CrossRef] [Google Scholar]
 Angus, G. W., Famaey, B., & Zhao, H. S. 2006, MNRAS, 371, 138 [NASA ADS] [CrossRef] [Google Scholar]
 Angus, G.W., Famaey, B., & Diaferio, A. 2010, MNRAS, 402, 395 [NASA ADS] [CrossRef] [Google Scholar]
 Baumgardt, H., Grebel, E. K., & Kroupa, P. 2005, MNRAS, 359, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Baumgardt, H., Cote, P., Hilker, M., et al. 2009, MNRAS, 396, 205 [Google Scholar]
 Begmann, K. G. 1989, A&A, 233, 47 [NASA ADS] [Google Scholar]
 Begmann, K. G., Broeils, A. H., & Sanders, R. H. 1991, MNRAS, 249, 523 [NASA ADS] [CrossRef] [Google Scholar]
 Bekenstein, J. 2004, PRD, 70, 083509 [Google Scholar]
 Bekenstein, J. D., & Milgrom, M. 1984, ApJ, 286, 7 [Google Scholar]
 Binney, S., & Tremaine, S. 1987, Galactic Dynamics (Princeton, NJ: Princeton Univ. Press) [Google Scholar]
 Bosch, V. D., & Dalcanton, J. J. 2000, ApJ, 534, 146 [NASA ADS] [CrossRef] [Google Scholar]
 Brada, R., & Milgrom, M. 2000, ApJ, 541, 556 [NASA ADS] [CrossRef] [Google Scholar]
 Clowe, D., Brada, M., Gonzalez, A. H., et al. 2006, ApJ, 648, L109 [Google Scholar]
 Ciotti, L., Londrillo, P., & Nipoti, C. 2006, ApJ, 640, 741 [NASA ADS] [CrossRef] [Google Scholar]
 Famaey, B., & Binney, J. 2005, MNRAS, 361, 633 [NASA ADS] [CrossRef] [Google Scholar]
 Famaey, B., Bruneton, J., & Zhao, H. S. 2007, MNRAS, 377, L79 [NASA ADS] [CrossRef] [Google Scholar]
 Gentile, G., Famaey, B., Combes, F., et al. 2007, A&A, 472, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gentile, G., Famaey, B., Angus, G., & Kroupa, P. 2010, A&A, 509, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gentile, G., Famaey, B., & de Blok, W. J. G. 2011, A&A, in press [arXiv:1011.4148] [Google Scholar]
 Governato, F., Brook, C., & Mayer, L., et al. 2010, Nature, 463, 203 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Haghi, H., Rahvar, S., & Hasani Zonoozi, A. 2006, ApJ, 652, 354 [NASA ADS] [CrossRef] [Google Scholar]
 Haghi, H., Baumgardt, H., Kroupa, P., et al. 2009, MNRAS, 395, 1549 (HBK09) [NASA ADS] [CrossRef] [Google Scholar]
 Harris, W. E. 1996, AJ, 112, 1487 [NASA ADS] [CrossRef] [Google Scholar]
 Hasani Zonoozi, A., & Haghi, H. 2010, A&A, 524, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jordi, K., Grebel, E. K., Hilker, M., et al. 2009, AJ, 137, 4586 [NASA ADS] [CrossRef] [Google Scholar]
 Klypin, A., Kravtsov, A. V., Valenzuela, O., & Prada, F. 1999, ApJ, 522, 82 [Google Scholar]
 Kroupa, P., Famaey, B., de Boer, K. S., et al. 2010, A&A, 523, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Küpper, A. H. W., & Kroupa, P. 2010, ApJ, 716, 776 [NASA ADS] [CrossRef] [Google Scholar]
 Londrillo, P., & Nipoti, C. 2009, MSAIS, 13, L89 [NASA ADS] [Google Scholar]
 Malekjani, M., Rahvar, S., & Haghi, H. 2009, ApJ, 694, 1220 [NASA ADS] [CrossRef] [Google Scholar]
 Metz, M., Kroupa, P., & Libeskind, N., 2008, ApJ, 680, 287 [NASA ADS] [CrossRef] [Google Scholar]
 Mieske, S., Hilker, M., Jordán, A., et al. 2008, A&A, 487, 921 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Milgrom, M. 1983, ApJ, 270, 365 [NASA ADS] [CrossRef] [Google Scholar]
 Milgrom, M. 1986, ApJ, 302, 617 [NASA ADS] [CrossRef] [Google Scholar]
 Milgrom, M. 1994, ApJ, 429, 540 [NASA ADS] [CrossRef] [Google Scholar]
 Milgrom, M. 1995, ApJ, 455, 439 [NASA ADS] [CrossRef] [Google Scholar]
 Milgrom, M. 2009, PRD, 80, 123536 [Google Scholar]
 Moore, B. 1996, ApJ, 461, L13 [NASA ADS] [CrossRef] [Google Scholar]
 Moore, B., Ghigna, S., Governato, F., et al. 1999, ApJ, 524, L19 [Google Scholar]
 Nipoti, C., Londrillo, P., & Ciotti, L. 2007, ApJ, 660, 256 [NASA ADS] [CrossRef] [Google Scholar]
 SanchezSalcedo, F. J., & Hernandez, X. 2007, ApJ, 667, 878 [NASA ADS] [CrossRef] [Google Scholar]
 Sanders, R. H. 2005, MNRAS, 363, 459 [NASA ADS] [Google Scholar]
 Sanders, R. H., & McGaugh, S. 2002, ARA&A, 40, 263 [NASA ADS] [CrossRef] [Google Scholar]
 Sollima, A., & Nipoti, C. 2010, MNRAS, 401, 131 [NASA ADS] [CrossRef] [Google Scholar]
 Spergel, D. N., Verde, L., Peiris, H. V., et al. 2003, ApJS, 148, 175 [NASA ADS] [CrossRef] [Google Scholar]
 Tiret, O., & Combes, F. 2007, A&A, 464, 517 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wu, X., Zhao, H. S., Famaey, B., et al. 2007, ApJ, 665, L101 [NASA ADS] [CrossRef] [Google Scholar]
 Wu, X., Famaey, B., Gentile, G., Perets, H., & Zhao, H. 2008, MNRAS, 386, 2199 [NASA ADS] [CrossRef] [Google Scholar]
 Zhao, H.S, & Famaey, B. 2006, ApJ,638, L9 [NASA ADS] [CrossRef] [Google Scholar]
 Zlosnik, T. G., Ferreira, P. G., & Starkman, G. D. 2007, PRD, 75, 044017 [Google Scholar]
All Tables
All Figures
Fig. 1
Lineofsight velocity dispersion for clusters listed in Table 1 for various masses as found by NMODY. In order to compare with the real cluster, the halfmass radii of all models are fixed. For low masses, which mean low internal accelerations, the prediction of Newtonian dynamics differs from the numerical solution. As the mass increases (at about 10^{6} M_{⊙}), the internal acceleration grows and the system enters the Newtonian regime and the numerical solutions get close to the Newtonian prediction. In order to focus on a reasonable range of mass for these clusters, we show the results in the range of 500−2 × 10^{5}. Dashed vertical lines show the expected cluster mass, calculated from the absolute luminosity by assuming a stellar masstolight ratio of M/L_{V} = 2. This value agrees with observed masstolight ratios of galactic GCs (Mieske et al. 2008). 

In the text 
Fig. 2
Top panel: the ratio of the velocity dispersion in the direction of the external field compared to the perpendicular direction for Pal 3 as an example plotted against the radius. The system is isotropic throughout the cluster. Bottom panel: the anisotropy profile, β, for Pal 3 as an example plotted against the radius. The system is isotropic throughout the cluster. 

In the text 
Fig. 3
Upper panel: radial velocity distribution of stars for a Newtonian (red triangles) and MONDian (black dots) system in which both have the same internal velocity dispersion. The smoothed line shows a Gaussian fitted to the distributions. Both distributions follow the Gaussian fit very closely. Lower panel: the relative difference between Newtonian and Mondian velocity distributions. There is a difference between the two models, but it is too small (less than 1%) to be detected observationally. 

In the text 
Fig. 4 Mean Pvalue vs. number of stars in the sample for a number of random realizations of distributions of stars following a Gaussian with a standard deviation of σ_{N} for clusters listed in Table 1. The minimum number of stars necessary to reject MOND is defined where the mean Pvalue drops below 0.05. 

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.