arxiv:astro-ph/ v1 27 Nov 2005

Transcription

1 Astronomy & Astrophysics manuscript no. tian April 25, 28 (DOI: will be inserted by hand later) Roche Lobe Shapes in MOND-like Modified Gravity HongSheng Zhao,2, LanLan Tian 3,4 arxiv:astro-ph/5754v 27 Nov 25 National Astronomical Observatories, Chinese Academy of Sciences, Beijing 2, PRC and 2 SUPA, School of Physics and Astronomy, University of St Andrews, KY6 9SS, Fife, UK, 3 National Astronomical Observatories, Chinese Academy of Sciences, Nanjing 28, PRC and 4 University of Victoria, Department of Physics and Astronomy, Victoria, British Columbia, V8P 5C2, Canada Abstract. We consider how to break the near degeneracy between dark matter and baryonic MOdified Newtonian Dynamics (MOND). We show that the Roche Lobes of a two-body baryonic system (e.g., a globular cluster orbiting a host galaxy) are sensitive to modifications of the law of gravity. We generalise the analytical results obtained in the deep-mond limit by Zhao (25, astro-ph/573), consider in the framework of a general non-newtonian law g R ζ. We give analytical expressions for the inner Lagrange point and Robe lobe axis ratios. The Roche lobe volume is proven to scale linearly with the baryonic mass ratio which applies to any modification function µ(g) or power-law function ζ. The lobes are squashed, with the aspect ratio varying with µ(g) due to the anisotropic dilation effect in MOND-like gravity. We generalise these results for extended mass distribution. Precise measurement of Roche lobe shapes for different satellite-host separations could constrain the law of gravity. Key words. dark matter galaxy kinematics and dynamics gravitation galaxies: dwarf globular clusters. Introduction The law of gravity is still uncertain experimentally on very small scales and very large scales. The running controversy over dark matter or MOND is an example. The observed acceleration g, traced by gas clouds and stars in outer parts of steady state galaxies, is weak, often weaker than a characteristic scale, a 8 cmsec 2. Yet the Newtonian gravity g N from these baryonic matter is often a factor of ten lower than the observed g. This is normally the justification to invoke massive halos of dark matter particles on galaxy scales. However, the alternative gravity theory of Modified Newtonian Dynamics (MOND) (Milgrom 983) has been doing very well in fitting kinematic data on galaxy scales, often better than the standard cold dark matter theory. In the MOND theory the gravity g(r) at distance R from a baryonic point mass M is boosted from the Newtonian value g N (R) = GM/R 2 by a factor /µ > such that far away from the mass point we have µ g ( ) ( ) N/a g R, () g/a a R GM This reproduces the flat rotation curves in bright disk galaxies at radii beyond R = a since V 2 gr R cst. The predictive power of this 2-year-old classical theory with virtually no free parameters (Bekenstein & Milgrom 984) is recently highlighted by the astonishingly good fits to contemporary kinematic data of a wide variety of high and low surface brightness spiral and elliptical galaxies; even the fine details of the ups and downs of velocity curves are elegantly reproduced without fine tuning of the baryonic model (Sanders & McGaugh 22, Milgrom & Sanders 23). For a long time a problem has been that the gravitational lensing in this non-relativistic theory is illly-posed, and it misses a factor of two in the bending angle (e.g., Qin, Wu & Zou 995). However, this empirical MOND has now a respectable relativistic field theory formulation (called TeVeS by Bekenstein 24), which passes standard tests to check General Relativity, and allows for rigourous modeling of Hubble expansion and gravitational lensing. This has generated wide interests, and many are examining the consequences of modifications to gravity (Skordis, Mota, Send offprint requests to: H. Zhao, hz4@st-and.ac.uk Outstanding Young Overseas Scholar and PPARC Advanced Fellow

2 2 HongSheng Zhao,2, LanLan Tian 3,4 : Roche Lobe Shapes in MOND-like Modified Gravity Ferreira et al. 25, Hao & Akhoury 25, Chiu, Ko & Tian 25, Zhao, Bacon, Taylor et al. 25, Pointecouteau & Silk 25, Ciotti & Binney 24, Baumgardt, Grebel, Kroupa 25, Read & Moore 25, Famaey & Binney 25, Zhao & Famaey 25). Zhao (25) has shown that the shape of the Roche lobes in the deep-mond limit is more squashed than in the Newtonian case. Here we study the shape of the Roche lobe of a two-body system in a MOND-like gravity. We show that Roche lobes shape varies with the assumed law of gravity, and could be used to differentiate among laws of gravity. The results are first shown for point masses and are generalised to extended mass distribution. 2. MOND-like gravity In a modified gravity, the gravity g and gravitational potential Φ follow a Poisson-like equation [ µg ] = 4πG N M i δ (R R i ), g = Φ, (2) i= for an isolated system of N baryonic point masses M,..., M N at distances R i. Here µ is the modification factor, and is a function of the gravity g; the Poisson s equation for the Newtonian gravity is recovered if µ =. This modification satisfies the usual conservations of total energy and angular momentum for any geometry of baryonic distribution (see Bekenstein & Milgrom 984). A very broad class of modified gravity models could be parametrized by a modification function given as follows, [ µ(g) = g ( ) ] kn n a = +, g N g N (3) which makes µ(g) an implicit function of the gravitational strength g = g. The Newtonian gravity corresponds to models with k = n =. Exotic gravity can be achieved by letting, e.g., (k, n) = (, ) or (k, n) = ( 3 2, )2. Conventional MOND gravity corresponds to models with k = /2. Such a model goes from g g N = GMR 2 in strong gravity to g g N a = GMa R in weak gravity around a point mass M. The sharpness of the transition GM is controlled the parameter n. The location of the transition zone for a point mass M is given by R a. For more compact notations later on we shall introduce a few auxilary functions, related to µ in a spherical potential. We define (g) d lnµ d lng, We also define + = d lng N d lng. (4) ( ) d lnr(g) ζ(g), ζ ζ + = d lng d ln Ω2 d lnr, (5) where Ω(R) = g R is the angular frequency for a circular orbit of radius R. In writing down R(g) we treat the spherical radius R as a function of the gravitational field strength g through the relation MR 2 g N = µg. Note in the case of the field around a point mass, and assume non-negative n and k, (g) is related to the parameter ζ(g) by ζ 2 = d lng d lnr 2 = { near a point mass = k far from a point mass. (6) Conventional MOND models have the property that and ζ 2 when the gravity is strong and and ζ when the gravity is weak compared to the acceleration scale a. The result in this paper applies to any function (g) or ζ(g). In general, we shall call the above gravity g the MOND-like Gravity or Modified Gravity. Fig. shows how ζ and change as the MONDian gravity decreases away from a point mass. These MOND models have a narrow range of ζ 2 or. so that the gravity g cst at large distances as in, e.g., (Mannheim 997). 2 so that g R 5 at small distances as in, e.g., the theory of Qin et al. 25 with three extra nanometer-thick spatial dimensions.

3 HongSheng Zhao,2, LanLan Tian 3,4 : Roche Lobe Shapes in MOND-like Modified Gravity 3 2.e Log R/R, 2 Log R/R.e2.e Log R/R,. 2 Log R/R Fig.. Upper panels: shows ζ (red circles) and (black diamonds) as functions of rescaled orbital distance R/R for a range of modification functions (cf. eq. 3). Left panels: for point mass with (k, n) = (/2, ) (smallest symbols), (k, n) = (/2, 8) (bigger symbols); typical MOND models have ζ and in the range between these two curves. Right panels: for an extended Hernquist profile mass with a scale b = /2 in a stronger-than-mond gravity with (k, n) = (9/, ). Lower panels: same parameters as in upper panels, except for showing the circular velocity curves V cir (R) = gr in modified gravity (red circle) and V cir (R) = g N R in Newtonian gravity (blue cross); both in units of (GM a ) /4. 3. Roche Lobes of a two-body system in a MOND-like gravity A binary system is specified by two baryonic masses M and m separated by a distance D o along, say, the z-axis. The two bodies could be a galaxy and a satellite. The MOND-like gravity distribution is determined by eq. 2 where N = 2 for a binary. Define the inner Lagrange point of the binary at a distance r and D o r from the masses m and M. There are only 3 independent dimensionless quantities in this system, r /(D o r ), m/m and a typical value for the parameter ζ. This suggests that the Roche radius might be expressable in general as ( ) B(ζ) r m =, (7) D o A(ζ)M where A and B are dimensionless functions with dimensionless arguments; the exact form requires detailed modeling of the gravitational potential. We set up the coordinate system such that the origin is on the smaller mass m, at a distance D o from the mass M on the z-axis. We shall treat the gravity due to the smaller mass m as a perturbation to that of the bigger mass M. First consider the case without a satellite m =. At a distance R far way from the isolated baryonic mass M, the gravitational potential can in general be approximated as spherical with Φ = g(r)dr, R = (z D o ) 2 + y 2 + x 2 (8) where the gravity g(r) at radius R from the spherical galaxy centre is related to the Newtonian gravity GMR/R 3 by µg = g N = GM R 3 [x, y, z D o], g = Φ (R). (9) Taylor-expand the three components of the gravity in the vicinity of the origin (x, y, z) = (,, ), we have Φ [x, y, D o ]Ω 2 + [,, z] dg(d o) dd o, Ω 2 g(d o) D o. ()

4 4 HongSheng Zhao,2, LanLan Tian 3,4 : Roche Lobe Shapes in MOND-like Modified Gravity.8 z/s y/s, x/s Fig.2. shows one octant of the 3-dimensional Roche Lobe of a satellite of mass m rotating in the xz plane around ( /3 m a mass M at a distance D o on the z-axis. The sizes are scaled by a length s M(D o)) Do. The curves are for ζ = d log g d log D o = 2 (Newtonian, thick black) and ζ = (Deep-MOND, red). A slimer Roche lobe for ζ = /5 (thin green) and a fatter Roche lobe for ζ = 5 (thin blue) are also made for comparison. The gravity field is such that if a test particle is put on circular orbit with an orbital radius R = D o around M, then the rotation angular frequency must be Ω in order to balance gravity with centrifugal force. Integrate the above potential gradient, we find Φ (R) Φ (D o ) + [ x 2 + y 2 z 2 ζ ] Ω zd oω 2 () in the vicinity of the origin (,, ), where according to the definitions of ζ and we have ζ = d log g(d o) d log D o. Let a low-mass satellite with m/m rotate around the mass M (fixed) with an angular velocity Ωŷ, then particles in the corotating frame conserve the Jacobi energy with an effective (triaxial) potential Φ e (x, y, z) Φ + Φ x2 + (z D o ) 2 Ω 2, 2 where Φ (x, y, z) is the perturbation due to the satellite. Now we consider the effect of the self-gravity of a satellite mass distribution inside a weak external galaxy field. Sufficiently far away from the satellite, the gravity is dominated by the spatially slow-varying gravity field of the galaxy g(r), and the galaxy s gravity determines the slow-varying MOND dielectric index µ. The perturbation due to the self-gravity of the satellite of total mass m is almost Newtonian except in an anisotropic medium (see eq. of Milgrom 986). To the first order in m the perturbation in potential is given by Φ (x, y, z) = Gm r, where m is the modified inertia of the satellite given by Gm Gm = g(d o) g N (D o ) = Ω2 D o GMDo 2, and r is the effective distance from the centre of the satellite and is given by r = z 2 + (y 2 + x 2 ), where = + is the MOND dilation factor. The combined potential Φ(x, y, z) = Φ (R) + Φ (x, y, z) clearly has axial symmetry around the z-axis, and has two centres separated by distance D o along the z-axis. The above (2) (3) (4) (5) (6)

5 HongSheng Zhao,2, LanLan Tian 3,4 : Roche Lobe Shapes in MOND-like Modified Gravity 5 formulation approximates the potential of, e.g., the Milky Way galaxy of baryonic mass M plus a satellite of baryonic mass m. We also note that the density ρ = along the z-axis for the gravity field g = Φ of the above combined potential; this can be verified by evaluating the modified Poisson s equation ρ = 4πG (gµ). The effective potential Φ e (x, y, z), if Taylor expanded near the secondary (,, ), has an triaxial shape with [ ] y 2 ζ z 2 Φ e cst + Ω 2 D 2 o 2D 2 o D o m/m z2 + (y 2 + x 2 ), (7) where the terms to the first order of z and second order to x are canceled due to force balance. The inner or outer Lagrangian points is then calculated from the saddle point of the effective potential where = Φ [ e(,, z) z = Ω 2 D o r ζ ± m D 2 ] o z=±r D o M r 2. (8) This definition of the Lagrange radius r reduces to ( ) ( ) m M r 3 = ζ Do 3. (9) So inside the Lagrange radius r the average density of the satellite equals k-times the average density inside the orbital radius D o of the satellite. The equation can also be written as ( ) r m 3 =, ζ = + ζ. (2) D o ζ M Note that the masses here m and M are true baryonic masses of the binary, not the modified inertia masses. Eq. (9) can be further massaged into Ω Gm = ζ Ω, (2) r 3 where Ω takes the meaning of the angular frequency of circulation around the secondary at its Robe radius. Hence the tidal radius of the secondary scales linearly with the terminal velocity of the secondary, and the internal period is within a factor of ζ of the period of the secondary s orbit. Finally note that near the tidal radius, the gravity due to the secondary and primary has the ratio Gm ( m ) r 2g(D o) = ζ 3 M, which confirms that near the tidal radius r it is valid to approximate the self-gravity of a low-mass secondary (m/m ) as a perturbation to the slow-varying gravity g(d o ) of the primary. The shape of the Roche Lobe is specified by contour of equal effective potential with the contour passing the Lagrange points (,, ±r ). The effective potential can be redefined as a dimensionless function Ψ e (x, y, z) after massaging zero point and a prefactor. Substitute (m/m) in favor of r /D o using eq. (2), we found that the points (x, y, z) along the Roche lobe satisfy Ψ e = 3ζ r 2 + y 2 ζ z 2 2r 2 ζ r =. (23) z2 + y 2 + x2 Clearly the z-axis intersects the Roche Lobe at (,, ±r ), where ( ) r m 3 =, D o ζ M which defines the first axis of the Roche Lobe. Draw a line parallel to the x-axis direction, it intersects the Roche lobe at (±r 2,, ). Solving equation 23) we find 3r 2 =, (25) 2r which defines the second axis of the Roche Lobe. Likewise the third axis of the Roche Lobe is defined by the intersection points (, ±r 3, ) with a line parallel to the rotational y-axis. We find 2r [ 3 = (u 2 ) 3 (u + ) 3 (u ) 3 3r 2 (22) (24) ], u + ζ (26) by simply solving equation (23), which reduces to an essentially cubic equation ( ) 2 r 2 4 r3 =. r 3 27 ζ r 2 (27)

6 6 HongSheng Zhao,2, LanLan Tian 3,4 : Roche Lobe Shapes in MOND-like Modified Gravity.8 r r2 r dlog g/dlog R Fig.3. shows the Roche lobe semi-axis sizes r, r 2, r 3 of a satellite m at orbital distance R = D o as functions of ( /3 the gravity power-law index ζ = d log g m d log D o. Here the sizes are rescaled by s M(D o)) Do. Also shown are the arithematic mean (r +r 2 +r 3 )/3 (thicker dashed curve), and the geometric mean size r r 2 r 3 /3 (thin dashed curve). 4. Generalisation and Discussion 4.. Extended spherical primary mass Although we refer to our primary mass as a mass point, our result applies exactly for a spatially extended spherical mass distribution of the primary M(R), e.g., a galaxy, if we simply make the substitution M M(R). As a result the ζ relation changes to 2 d lnm d lnr = ζ. (29) As a specific example, a host galaxy can be approximated by a spherical Hernquist mass profile with ( ) 2 R GM M(R) = M, R, (3) R + y h R a with the scalelength parameter y h = /2, where we adopted notations of Zhao, Bacon, Taylor et al. (25). The Newtonian gravity of the model is g N = GM(R) R 2 = a (y h + y) 2, y R R. (3) Note g N = a /( + b) 2 = 4a /9 < a at R = R and g N = a /b 2 = 4a is finite at R =. We give the Hernquist model a stronger-than-mond gravity with (k, n) = (9/, ) so that the model has a rising terminal velocity curve. As shown by Fig.b this model has ζ = near the primary, and ζ = /5 far away Prolate or oblate? Interestingly, the second and the third (i.e., rotation) axes of the Roche lobe are always very similar in length. This is 4 4 because equation (27) contains a factor 27 ζ = 27( +3) for any < <. In fact that r 2 = [,.2,.3,.45,.55,.65], (32) r 3 (28)

7 HongSheng Zhao,2, LanLan Tian 3,4 : Roche Lobe Shapes in MOND-like Modified Gravity 7.8 r2/r r3/r dlog g/dlog R Fig.4. similar to the previous one, but shows the aspect ratio of the Roche lobe semi-axis sizes r 2 /r, r 3 /r as functions of ζ = d log g d log R. for ζ = [,.5,, 2, 4, ], where ζ =, 2 are the deep-mond, and deep-newtonian cases. The first axis is often substantially longer than the other two, but not always. We find for point masses r = [, 3, 2.2,.5,.95, ], for ζ = [,.5,, 2, 5, ]. (33) r 2 The generally triaxial Roche lobe for point masses is nearly prolate for < ζ < 5, and it becomes nearly spherical for ζ = 4 6, and becomes nearly oblate only in exotic situations where the power-law index ζ = d lng/d lnr is harder than 5. For extended mass, the ratio can be calculated from r r 2 = 3 2 = 3 2 d log g N d log g, together with a relation between g and g N (cf. eq. 3) and g N = GM(R) R. 2 Overall the Roche lobe axis ratio changes by a significant amount for different assumptions of the laws of gravity and the mass profile (cf. Fig. 5), e.g., the MONDian Roche lobe (k = /2) at large distance is significantly more squashed than the Newtonian case, and a modified gravity with k = 9/ at large distance is even more squashed. We speculate that it might be feasible to differentiate these laws with a precise measurement of a Roche lobe. One advantage of using axis ratio is that it does not require precise knowledge of the satellite baryonic mass and its actual distance and absolute size; none of these factors enter equations (25) and (26). In particular we note that in Newtonian gravity =, hence r r 2 = 3 2 holds for any mass distribution at any distance; it would be interesting for testing non-newtonian gravity to quantify the deviation from this ratio in observed Roche lobe. (34) 4.3. Size-mass and size-distance relations In general, the Roche radii of two-body systems (whether stellar binaries or satellite-host galactic binaries) provide a measure of the relative weights of the binary in baryons in any MOND-like gravity. It is interesting that the Roche Lobe size in any MOND-like gravity has the same simple scaling relation m r 3 / M(D o) D 3 o = ζ = 2 d log V 2 cir d log R. The above relation holds for an extended primary as well. As in the Newtonian case the relevant mass is the unmodified true baryonic mass of the binary, the difference is only in the prefactor, which depends on the M and D o. For point (35)

8 8 HongSheng Zhao,2, LanLan Tian 3,4 : Roche Lobe Shapes in MOND-like Modified Gravity.7.5 r2/r r3/r.7.5 r2/r r3/r Log R/R, Log R/R Fig.5. showing the Roche lobe axis ratios r 2 /r, r 3 /r as functions of R/R with R = D o being the orbital distance of the secondary. Left panel: Models are as in Fig.a and Fig.c, with (k, n) = (/2, ) (solid) and (k, n) = (/2, 8) (dashed). Right panel: as in Fig.b and Fig.d, with (k, n) = (9/, ) (solid) and (k, n) = (/4, ) (dashed) around a Hernquist profile mass. r.8.8 r r23 r Log R/R, 2 Log R/R Fig.6. showing the Roche lobe semi-axis sizes r, r 2 r 3 as functions of R/R. Here the sizes are rescaled by s ( /3 m M(R)) R with R = Do being the orbital distance of the secondary. Left panel: models as in Fig.a and Fig.c, with (k, n) = (/2, ) (solid) and (k, n) = (/2, 8) (dashed). Right panel: models as in Fig.b and Fig.d, with (k, n) = (9/, ) (solid) and (k, n) = (/4, ) (dashed) around a Hernquist profile mass. masses the prefactor varies from ζ = 3 in the region with Keplerian rotation V cir R /2 (Binney & Tremaine 987) to ζ = 2 for the flat rotation curve region. A similar relation holds for r 2, Θ IR r ( 2 m ) /3 2 =, D o M 3 2 ζ 3 (36) where ζ = 2 d log V 2 cir d log R, = d log g N d log g. Here Θ IR takes the meaning of a rescaled Roche lobe intermediate semi-axis. The size of the lobe varies with the distance from the primary. The Roche Lobe sizes in terms of r and r 2 r 3 r 2 can be read off from Fig. 6 for different assumptions of the µ-function and mass profile; in practise we can treat the Roche lobe as prolate. (37) 5. An application: Limiting radii of outer satellites of MW and M3 As an illustration of the usuage of the Roche lobe, we compare the predicted Roche sizes with the observed limiting radii of outer satellites of galaxies. Slightly different from Zhao (25), the sample consists of all globular clusters and dwarf galaxies beyond 2 kpc (outside the disks) of both the Milky Way and Andromeda. Only Andromeda satellites with reliable 3-dimensional distances to the Andromeda centre are used (McConnachie et al. 24); tidal radius and distance of the newly discovered AndIX are from Harbeck et al. (25). The limiting angular radii θ lim,min are taken primarily from Harris (996) and Mateo (998). Most systems are consistent with no flattening in projection. For a few significant flattened systems we take the limiting radii on the minor axis. This radius should approximate the

9 HongSheng Zhao,2, LanLan Tian 3,4 : Roche Lobe Shapes in MOND-like Modified Gravity 9 intermediate semi-axis r 2 of a nearly prolate Roche lobe if Roche lobes are filled; a prolate object appears the same width on the minor axis at any inclination angle. Hence we have r 2 = F /3 D D θ lim,min, F, D o D D o where the factor F takes into the facts that (i) the Roche volume is smaller at the pericenter of an eccentric orbit of a satellite, and that (ii) the volume filling factor of the pericentric Roche lobe can still be less than unity if star counts drop below detection limit before reaching the edge of the lobe (e.g., Grebel et al. 2). Several distances are involved here, D o is the present orbital distance, D is the present distance of the Sun from the satellite, and D is the Sun s distance to the centre of the host galaxy (here the Milky Way or the Andromeda). The rhs is insensitive to errors on distance D; for the outer Milky Way satellites we have D D o, and for Andromeda satellites we have D D. Assuming mass-tracing-light, we expect that the satellite-to-host mass ratio is the ratio q of observable stars within a factor of two, allowing for some differences in age and metallicity. So we can define F 2 m Mq.5 2. Multiplying a factor (m/m) /3 to both sides of the above equation 38, we have (38) (39) r ( 2 m ) /3 /3 F Θ min Θ min, D o M (4) where 3 Θ min D D D D o θ lim,min q /3, is an observable, taking the meaning of the rescaled limiting minor axis size of a satellite. If MOND-like mass-tracinglight models are correct, substitute in eqs. 24,and 25 we expect that Θ min Θ IR = 2 2/6 =.374, 3 2 ζ 3 3 where we used the fact that all the baryons are contained well inside 2kpc, and ζ = 2 d log V 2 cir d log R 2, = d log g N d log g (4) (42) 2, (43) Note that we do not have freedom in chosing the law of the gravity, or ζ or : they are all determined by the nearly flat circular velocity curves in the Milky Way and Andromeda, which are inferred from the fairly constant velocity dispersions of the outer satellites and other tracers of the outer halo such as PNe. Satellites with known proper motions often suggest only mildly eccentric orbits: pericentric distances are rarely smaller than /2 of the present distances, and almost never smaller than /5 of the present distance (Dinescu et al. 999, 24 and references therein). Outer globulars on long period orbits also have time to fill the instantaneous Roche lobe by two-body relaxation (Bellazzini 24). Globular clusters as a class generally exhibit extra-tidal stars under deep observations; a sharp break of the star count profile and/or sometimes two-dimensional confirmation of unvirialised structures outside the King radius have been shown for globular clusters in the Milky Way, Andromeda, and other galaxies (Leon, Meylan & Combes 2, Lehmann & Scholz 997, Sohn et al. 23, Harris et al. 22). To be generous we assume the instantaneous Roche lobe has a filling factor F., or F /3 F /3 5. So equation 42 predicts that the majority of satellites should lie in a band with dashed line boundaries as marked in Fig. 7. In reality we observe a number of outliers to the inequality 42 and the expectation band. Overall we observe a very large scatter of Θ min for satellites in similar gravitational field with similar expected values of Θ IR (cf. Fig. 7). In summary, the shapes of Roche lobes of a binary in a class of gravity theory including MOND are given analytically. The volume of the Roche lobe is proven to scale with the true mass ratio of the binary. The axis ratios of the Roche lobe are independent of the mass ratio, but are sensitive to the modifications to the law of gravity. Observed limiting radii of outer satellites of the Milky Way have a large scatter, with some outliers by a factor of 4 above or factor below predictions of our basic models for the Roche lobes in MOND-like mass-tracing-light gravity theories. Acknowledgements. We acknowledge travel and publication support from Chinese NSF grant to HSZ. 3 We drop the minor factor F /

10 HongSheng Zhao,2, LanLan Tian 3,4 : Roche Lobe Shapes in MOND-like Modified Gravity 2.5 Observation Expectation Fig.7. shows expectations vs. observations of the rescaled limiting minor axis radii Θ min of outer satellites (beyond 2kpc) of the Milky Way and Andromeda. The horizontal axis is the value of Θ min predicted from mass-tracinglight models assuming a satellite on a circular orbit filling its Roche lobe. The expectation band is shown by two slanted dashed lines; the expectation-observation equality line (thick diagonal dashed line) should be an upper limit for satellites are on eccentric orbits and not filling their Roche lobes. There are several outliers. To be conservative, L we adopt a relatively low luminosity for the Milky Way and Andromeda so that q = sat 4 L. To be accurate we ( ] 2 /4 r adopt a slightly falling circular velocity curve V cir = V [ + kpc) (cf. Wilkinson & Evans 999) for the Milky Way and Andromeda. References Baumgardt H., Grebel E.K., Kroupa P., 25, MNRAS, 359, L Bekenstein J. 24, Phys. Rev. D. 7, 359 Bekenstein J. & Milgrom M. 984, ApJ, 286, 7, Bellazzini, M. 24, MNRAS, 347, 9 Binney J. & Tremaine S., 987, Galaxy dynamics, PUP Brada R., Milgrom M. 2, ApJ, 54, 556 Chiu M., Ko C.-M., & Tian Y. 25, ApJ in press, astro-ph/57332 Ciotti L., Binney, J. 24, MNRAS, 35, 285 Dinescu, D. I., Girard, T. M., and van Altena, W. F. 999, AJ, 7, 792 Dinescu D., Keeney B.A., Majewski S.R., Girard T.M., 24, AJ, 28, 687, Famaey B., & Binney J. 25, MNRAS, in press, astro-ph/56723 Grebel E., Odenkirchen M., Harbeck, D., 2, AAS, 2, 467 Hao J. & Akhoury R. 25, astro-ph/543 Harbeck D., Gallagher J.S., Grebel E.K., Koch A., Zucker D.B. 25, ApJ, 623, 59 Harris, W. E., 996, AJ, 2, 487 Harris W.E., Harris G.L.H., Holland S.T., McLaughlin D.E. 22, AJ, 24, 435 Leon S., Meylan G., Combes F., 2, A&A 359, 97 Lehmann I. & Scholz R.D. 997, A&A, 32, 776 Mannheim P., 997, ApJ, 479, 659 Mateo M., 998, ARA&A, 36, 435 McConnachie A.W., Irwin, M., Ibata R., Ferguson A., Lewis G., Tavir N. 24, MNRAS, submitted, Milgrom M. 983, ApJ, 27, 365 Milgrom M., 986, ApJ, 32, 67 Milgrom M., Sanders R. H. 23, ApJ, 599, L25, Pointecouteau E., & Silk J., 25, MNRAS, in press, astro-ph/557 Qin B., Peng U., Silk J. 25, Phy. Rev. Lett. submitted (astro-ph/58572)

11 HongSheng Zhao,2, LanLan Tian 3,4 : Roche Lobe Shapes in MOND-like Modified Gravity Qin B., Wu X.P., Zou, L.Z., et al. 995, A&A, 296, 264 Read J.I., Moore B. 25, MNRAS, 36, 97 Sanders, R., McGaugh S. 22, ARA&A, 4, 263, Skordis C., Mota D.F., Ferreira P.G., & Boehm C., astro-ph/5559 Sohn Y. et al. 23 AJ, 26, 83 Wilkinson M., & Evans N.W. 999, MNRAS, 3, 645 Zhao H.S, 25, A&A Letters, in press (astro-ph/573) Zhao H.S., Bacon D., Taylor A.N., Horne K.D., 25, MNRAS submitted (astro-ph/5959) Zhao H.S, & Famaey B. 25, ApJ Letters, submitted