As I understand it, Milgrom's MOND model can be interpreted in one way by saying that Newton's Law of Gravitation on its own is insufficient to predict the gravitational acceleration when acceleration is very small.
A second interpretation is that Newton's Second Law (force = time rate of change in momentum, in the absence of other unbalanced forces) breaks down at very small accelerations. I have not found an understandable account of this interpretation, so I will focus on the first interpretation.
Milgrom, 1983 suggested a new fundamental constant with dimensions of acceleration and whose suggested magnitude (from fitting to observed astronomical data) is $|a_0| \approx 2*10e-8 cm/s^2$. For example this is the magnitude of the Newtonian gravitational acceleration $a=GM/r^2$ caused by a 3Kg cannonball, acting on a pea at a distance of 1m.
For a given source mass and reacting massive test particle, three regimes, based on the appropriate formula for acceleration $a_r$ at distance $r$, can be defined:-
(1) High acceleration Newtonian regime where $|a_r|>>|a_0$| and acceleration is well-described by Newton's Law: $ a_r = F_N/m = -GM/r^2$. The $-$ sign indicates that the acceleration is directed in the opposite direction to that in which distance $r$ from source to receiver increases.
(2) Transitional regime where some kind of transitional formula applies. Within this regime lies the transition distance $r_{t}=\sqrt{GM/a_0}$ where, for given $M$, the formulae for regimes (1) and (3) give the same acceleration : $ -GM/r^2 = a_r = -\sqrt{a_0 GM}/r $
(3) Low acceleration Deep MOND regime where where $|a_r|<<|a_0|$ and acceleration is given by:- $F_N/m = -GM/r^2 = a_r|a_r/a_0|$ and thus $a_r = -\sqrt{a_0 GM}/r.$
This seems fairly straightforward to me in the case of an isolated test particle some distance from a static concentrated source mass. For given $GM$ the appropriate acceleration regime is a variable function of source:target distance $r$. As the test particle moves away from the source, acceleration falls of initially as $\propto 1/r^2$ and eventually as $\propto 1/r$, with a gradual transition between the two regimes. And using the circular orbit centripetal force relation ($-V_R^2/r = a$) it is not hard to see how a fixed central mass $M$ can generate the sort of flat $V_R$ vs $r$ rotation curve interpreted (from observed doppler shifts) in the outer parts of galactic disks:-
$$\frac{V_R^2}{r} =-a_r = \frac{\sqrt{a_0 GM}}{r} \implies V_R^2 = \sqrt{a_0 GM} \implies V_R = constant $$
Also, over a range of galactic masses $M_i$, we derive $ {V_{iR}}^4 \propto M_i $ and, if luminosity $L_i \propto M_i$ then the observed Tully-Fisher relation $ {V_{iR}}^4 \propto L_i $ is obtained (as MOND was designed to do).
But I have a problem with this. While the model seems clear in the case of an isolated test particle far from the gravitational source I do not understand how it applies to a test particle in the interior of a large massive object such as a star, planet, moon, or asteroid. At any point inside a spherically symmetrical example of such an object (apart from a very small region at the centre) the net gravitational acceleration from all the object mass closer to the object centre, will be much larger than the given value of $a_0$. (The effect of object mass further out from the centre will net to zero in accord with Newton's Shell Theorem).
One way that a test particle might respond differently to gravitons from the object and to gravitons from a far distance source might be if some characteristic of the gravitons of either or both sources acted as an indicator of distance (or age, or divergence which both depend on distance). But Milgrom and others, e.g. Sanders & McGaugh, 2002, have emphasized that a distance dependency is not supported by the data. Moreover, how would the reaction mechanism "know" to respond differently to gravitons coming from one distant source e.g. a star of 2 solar masses compared to to gravitons coming from two distant (but close to each other) sources of same total mass (e.g. a binary pair of stars with each star having 1 solar mass).
I know that Milgram and other MOND supporters do not strongly defend their model in terms of any particular physical mechanism but rather point to its claimed success in predicting galactic rotation curves and the potential insights that might be gleaned from that. Maybe modern (astro)physics simply does not reject a model simply because of lack of even an intuitive mechanism. Maybe there is a plausible explanation of how massive particles in a star can react differently to gravitons from remote objects. Or maybe I am missing something obvious in the MOND model.
What do you think?