Modified Newtonian Dynamics (MOND)

• Modified Newtonian Dynamics (MOND) is a phenomenological framework that modifies the relation between observed baryonic acceleration and actual dynamical acceleration below a characteristic acceleration, providing an alternative or complement to the nonbaryonic dark-matter interpretation of galactic phenomena.
• MOND explains several galaxy-scale regularities, including flat rotation curves, the baryonic Tully–Fisher relation, and the Radial Acceleration Relation, through a modification of Newtonian dynamics, which become crucial at low accelerations.
• MOND’s predictive power is showcased in its ability to delineate galactic dynamics based on baryonic mass, aggregates substantial observational support, yet faces challenges such as with non-Newtonian dynamics, infinite escape energies, and relativistic generalization.

Modified Newtonian Dynamics (MOND), also called Milgromian Dynamics, is a phenomenological framework proposed as an alternative or complement to the nonbaryonic dark-matter interpretation of galactic mass discrepancies. It posits that the relation between the acceleration generated by observed baryons, $g_N$, and the actual dynamical acceleration, $g$, changes below a characteristic acceleration $a_0\simeq 1.2\times10^{-10}\ {\rm m\,s^{-2}}$. In its standard algebraic form,

$\mu\!\left(\frac{g}{a_0}\right)g=g_N,$

where $\mu(x)\to1$ for $x\gg1$ and $\mu(x)\to x$ for $x\ll1$. The deep-MOND limit is therefore $g\simeq\sqrt{a_0g_N}$. MOND provides a compact account of several galaxy-scale regularities, particularly flat rotation curves, the baryonic Tully–Fisher relation, and the Radial Acceleration Relation, but it remains incomplete as a relativistic and cosmological theory.

1. Physical motivation and foundational postulates

The principal empirical motivation for MOND is the mass-discrepancy problem. Once most of a galaxy’s visible mass has been enclosed, Newtonian gravity predicts a Keplerian decline,

$V(r)\simeq \sqrt{\frac{GM}{r}},$

whereas spiral-galaxy rotation curves generally remain approximately flat. Under the conventional interpretation, this requires extended dark-matter halos. MOND instead attributes the discrepancy to a modification of dynamics or gravity at sufficiently low acceleration.

The controlling variable in MOND is acceleration rather than total mass, size, luminosity, or morphological class. A compact, massive galaxy can remain approximately Newtonian if its internal acceleration is high, while a low-surface-density system can exhibit a large apparent mass discrepancy even when its total mass is modest. The empirical transition occurs near $g$0, whereas the deep-MOND regime is characterized by $g$1.

The basic limiting laws are

$g$2

and

$g$3

The deep-MOND equations are associated with space-time scale invariance under

$g$4

under which velocities remain invariant while accelerations scale as $g$5. In this limit, $g$6 and $g$7 appear through the combination $g$8.

A related phenomenological regularity is the close numerical association between $g$9 and cosmological scales:

$a_0\simeq 1.2\times10^{-10}\ {\rm m\,s^{-2}}$0

up to convention-dependent numerical factors. This coincidence motivates several proposed cosmological, Machian, quantum-gravitational, and infrared interpretations, but no unique derivation of $a_0\simeq 1.2\times10^{-10}\ {\rm m\,s^{-2}}$1 has been established.

The algebraic MOND prescription is exact only in certain symmetric configurations and, in modified-inertia formulations, for particular orbital situations. In general modified-gravity theories, the relation between $a_0\simeq 1.2\times10^{-10}\ {\rm m\,s^{-2}}$2 and $a_0\simeq 1.2\times10^{-10}\ {\rm m\,s^{-2}}$3 can include a divergence-free curl field. Consequently, the algebraic law is a phenomenological summary rather than a complete many-body theory.

2. Interpolating functions and characteristic predictions

The function $a_0\simeq 1.2\times10^{-10}\ {\rm m\,s^{-2}}$4 interpolates between the Newtonian and deep-MOND regimes. Frequently used choices include the simple function,

$a_0\simeq 1.2\times10^{-10}\ {\rm m\,s^{-2}}$5

and the standard function,

$a_0\simeq 1.2\times10^{-10}\ {\rm m\,s^{-2}}$6

The corresponding inverse formulation is

$a_0\simeq 1.2\times10^{-10}\ {\rm m\,s^{-2}}$7

with $a_0\simeq 1.2\times10^{-10}\ {\rm m\,s^{-2}}$8 for $a_0\simeq 1.2\times10^{-10}\ {\rm m\,s^{-2}}$9 and $\mu\!\left(\frac{g}{a_0}\right)g=g_N,$0 for $\mu\!\left(\frac{g}{a_0}\right)g=g_N,$1. The interpolation function is not uniquely specified by the asymptotic limits. Galaxy rotation-curve studies often favor a transition resembling the simple function, whereas Solar-System tests require a more rapid return to Newtonian behavior at high acceleration. The simple function, used without modification, approaches unity too slowly and can conflict with planetary constraints (Gentile et al., 2010).

For an isolated spherical mass $\mu\!\left(\frac{g}{a_0}\right)g=g_N,$2,

$\mu\!\left(\frac{g}{a_0}\right)g=g_N,$3

The deep-MOND relation then becomes

$\mu\!\left(\frac{g}{a_0}\right)g=g_N,$4

For circular motion, $\mu\!\left(\frac{g}{a_0}\right)g=g_N,$5, so

$\mu\!\left(\frac{g}{a_0}\right)g=g_N,$6

This is the baryonic Tully–Fisher relation (BTFR):

$\mu\!\left(\frac{g}{a_0}\right)g=g_N,$7

where $\mu\!\left(\frac{g}{a_0}\right)g=g_N,$8 is the baryonic mass and $\mu\!\left(\frac{g}{a_0}\right)g=g_N,$9 is the asymptotic flat rotation speed. The asymptotic velocity depends on total baryonic mass but not on galaxy size or baryonic surface density at fixed mass.

MOND also predicts that the onset of the mass discrepancy occurs when the relevant acceleration crosses $\mu(x)\to1$0:

$\mu(x)\to1$1

The detailed transition-region rotation curve remains dependent on the interpolation function, the mass distribution, and the specific MOND formulation.

The deep-MOND scaling also implies characteristic relations for pressure-supported systems. For suitable isolated, virialized systems,

$\mu(x)\to1$2

although exact coefficients depend on morphology, anisotropy, constituent masses, and whether the theory modifies gravity or inertia. In a smooth deep-MOND modified-gravity system, a virial relation gives

$\mu(x)\to1$3

under the assumptions specified by Milgrom (Milgrom, 2012).

MOND defines a characteristic surface-density scale,

$\mu(x)\to1$4

Quasi-isothermal systems with mean surface density below or near this scale are expected to exhibit deep-MOND behavior. The apparent halo surface density obtained by reinterpreting MOND in Newtonian terms is likewise expected to be of order $\mu(x)\to1$5 or below.

3. Galactic phenomenology and observational tests

The strongest empirical domain of MOND is the radial dynamics of galaxies. Given a baryonic distribution, MOND predicts the Newtonian acceleration $\mu(x)\to1$6 and maps it to the observed acceleration through $\mu(x)\to1$7 or $\mu(x)\to1$8. The model generally adjusts the stellar mass-to-light ratio, while gas distributions are often directly constrained by observations.

The Radial Acceleration Relation (RAR) connects the observed acceleration,

$\mu(x)\to1$9

with the Newtonian acceleration generated by baryons, $x\gg1$0. Its limiting behavior is

$x\gg1$1

at high acceleration and

$x\gg1$2

at low acceleration. The RAR is regarded as MOND’s strongest observational success because it organizes thousands of measurements across many galaxies with a small intrinsic scatter. For the SPARC sample, one analysis reported an intrinsic scatter of $x\gg1$3 dex under a basic error model, falling to $x\gg1$4 dex when correlated galaxy-wide nuisance parameters are fitted simultaneously (Desmond, 27 May 2025).

MOND predicts a systematic dependence of inner rotation-curve shape on baryonic surface density. High-surface-brightness galaxies reach high internal accelerations and tend to have rapidly rising inner curves. Low-surface-brightness galaxies can remain in the deep-MOND regime over much of their extent and therefore have more slowly rising curves. At fixed baryonic mass, both classes approach the same asymptotic velocity through the BTFR.

This behavior is particularly significant in gas-dominated low-surface-brightness galaxies, where uncertainties in stellar mass-to-light ratios are reduced. MOND analyses also emphasize the correspondence between baryonic features and rotation-curve features, commonly summarized as Renzo’s rule. The inferred discrepancy is correlated with local baryonic acceleration rather than merely with radius or total luminosity.

High-resolution rotation curves from the THINGS survey provided fits for 12 nearby galaxies spanning a wide luminosity range. The simple interpolating function performed better overall than the standard function, and the fitted acceleration scale was

$x\gg1$5

close to earlier determinations. Allowing distances to vary within their observational uncertainties produced excellent-quality fits for approximately 75% of the sample. The main exception was NGC 3198, for which non-circular motions, a possible bar, and distance uncertainties were discussed as possible sources of tension (Gentile et al., 2010).

The Milky Way provides a more restrictive test because its baryonic morphology and kinematic data can be modeled in detail. An analysis using seven bulge morphologies, five disk morphologies, and two gas models found that the standard interpolating function generally required

$x\gg1$6

inconsistent with the value inferred from external galaxies. The simple function performed better and was compatible with the external-galaxy range for a subset of baryonic models, although not for most of the 70 combinations considered (Iocco et al., 2015).

Early-type galaxies provide an additional test. Modeling approximately 9046 ellipticals and lenticulars with spherical Sérsic profiles, constant anisotropy, and MOND gravity showed that about 92% could be fitted under the adopted assumptions. The inferred stellar mass-to-light correction was close to unity, with modest radial anisotropy. MOND also produced a fundamental-plane-like relation, although its predicted coefficient of the velocity-dispersion term remained generally larger than observed (Cardone et al., 2010).

Groups extend MOND tests to accelerations far below those commonly probed by galaxy rotation curves. For a selected sample of 56 intermediate-richness groups, the deep-MOND virial relation

$x\gg1$7

gave a median $x\gg1$8-band mass-to-light ratio of approximately $x\gg1$9 in solar units, compared with a Newtonian dynamical median near $\mu(x)\to x$0. The groups followed the deep-MOND acceleration branch down to accelerations of order $\mu(x)\to x$1, subject to uncertainties in virialization, membership, anisotropy, distances, and the external-field effect (Milgrom, 2018).

4. Nonrelativistic formulations and the external-field effect

The algebraic MOND relation does not by itself guarantee momentum conservation in general many-body systems. Two principal modified-gravity formulations provide action-based alternatives.

AQUAL, or the A-quadratic Lagrangian theory, modifies the Poisson equation:

$\mu(x)\to x$2

In spherical symmetry this reproduces the algebraic MOND relation. In the deep-MOND regime, the action contains a term proportional to $\mu(x)\to x$3, producing the nonlinear operator

$\mu(x)\to x$4

QUMOND, or quasi-linear MOND, first solves the Newtonian Poisson equation,

$\mu(x)\to x$5

and then solves

$\mu(x)\to x$6

Its computational structure is closer to two ordinary Poisson solves. A Newtonian-equivalent phantom density can be defined by

$\mu(x)\to x$7

This density is a bookkeeping device rather than a required material component, and it can be non-spherical, displaced from the baryons, or locally negative.

The external-field effect (EFE) follows from the nonlinearity of MOND. A system’s internal dynamics can depend on an approximately uniform external acceleration, violating the strong equivalence principle. In the deep-MOND limit of AQUAL, when the external field dominates the internal field,

$\mu(x)\to x$8

the internal dynamics becomes approximately Newtonian with

$\mu(x)\to x$9

together with anisotropic corrections. Thus a system with very low internal acceleration can behave quasi-Newtonian if embedded in a sufficiently strong external field.

The EFE is relevant to dwarf satellites, globular clusters, wide binaries, the outer Milky Way, galaxy groups, and tidal systems. It can also regulate escape speeds because an isolated deep-MOND potential is logarithmic and formally has infinite escape energy, whereas an external field eventually restores quasi-Newtonian behavior.

Modified-inertia MOND constitutes a distinct class. In this approach, the gravitational field may remain Newtonian while the kinetic response depends on acceleration or, more generally, on the full particle trajectory. Such theories can reproduce the circular-orbit MOND relation but may be nonlocal in time and have difficulty with relativistic generalization. They can also differ from modified-gravity theories in their treatment of noncircular motions and the EFE.

5. Relativistic completions and proposed theoretical origins

The algebraic MOND law is nonrelativistic. A relativistic completion must reproduce MOND in the weak-field, quasistatic, low-acceleration regime; recover general relativity at high acceleration; predict gravitational lensing; satisfy Solar-System and gravitational-wave constraints; and provide a viable cosmology.

TeVeS, the Tensor–Vector–Scalar theory, supplements an Einstein metric with scalar and vector fields and uses a disformal physical metric. This allows the potential affecting massive particles to contribute appropriately to light deflection. Generalized Einstein–Aether theories use a dynamical unit timelike vector field, while BIMOND introduces two metrics coupled through the difference of their connections. Other proposals include scalar–tensor models, dipolar dark matter, nonlocal metric theories, and theories with additional massive fields.

A separate class of proposals attempts to derive the deep-MOND asymptote from modified metric actions. One construction uses an $x\ll1$0-type action with

$x\ll1$1

and a source-dependent length $x\ll1$2. In the weak-field, far-from-source regime, this yields

$x\ll1$3

but it does not derive a complete interpolation function or an all-scale theory (Mendoza et al., 2012).

Another recent proposal introduces an infrared deformation of the Einstein–Hilbert action,

$x\ll1$4

with an acceleration invariant constructed from the lapse. In the quasistatic limit, the deformation is designed to reduce to AQUAL, while its high-acceleration limit becomes negligible. The deep-MOND action is proportional to $x\ll1$5 and is associated with a ten-parameter conformal symmetry isomorphic to $x\ll1$6. The proposal relates $x\ll1$7 to an infrared de Sitter length $x\ll1$8 through

$x\ll1$9

but a fully covariant completion, perturbation theory, stability analysis, and dynamical derivation of the infrared vacuum remain unresolved (Singh, 7 Jan 2026).

Other speculative approaches interpret MOND as modified inertia arising from cosmological or quantum-gravitational effects. One proposal associates the low-acceleration regime with de Sitter and Unruh temperatures, suggesting an environmentally dependent ratio of inertial to gravitational mass. Another uses Sciama-type gravitoelectric induction and cosmic frame-dragging, with the scale

$g\simeq\sqrt{a_0g_N}$0

arising from the relation between the cosmic horizon and an observer’s Rindler horizon. These constructions derive the simple interpolation function only heuristically and do not yet provide complete relativistic field theories (Smolin, 2017, Uruena, 7 Oct 2026).

A phenomenological inverse-Yukawa model instead supplements Newtonian gravity with a scalar potential

$g\simeq\sqrt{a_0g_N}$1

designed to be negligible in the Solar System, attractive on galactic and cluster scales, and repulsive on cosmological scales. Such proposals seek to connect MOND-like galaxy dynamics with Friedmann–Robertson–Walker cosmology, but they remain heuristic and do not supply a covariant action, stress-energy tensor, or consistent treatment of a spatially varying cosmological term (Falcon, 2010).

6. Limitations, controversies, and scientific status

Clusters remain the most persistent astrophysical difficulty. MOND reduces the cluster mass discrepancy but generally leaves a residual factor of roughly $g\simeq\sqrt{a_0g_N}$2–$g\simeq\sqrt{a_0g_N}$3, especially in central regions. Cluster dynamics, X-ray gas, and lensing require more gravitational influence than the observed baryons provide. Proposed additions include undetected baryons, ordinary or sterile neutrinos, and additional scalar or vector fields, but no universally accepted solution has emerged.

The Bullet Cluster is particularly challenging because the strongest lensing features are associated with the collisionless galaxies rather than the dominant X-ray gas. MOND phantom density can be displaced from the baryons, so the system is not an automatic logical refutation of every MOND theory, but reproducing its lensing map generally requires additional mass or a more elaborate relativistic completion.

Solar-System constraints generate a tension with galaxy phenomenology. In AQUAL and QUMOND, the Milky Way’s external field induces a Solar-System quadrupole. Cassini tracking constrains its amplitude to

$g\simeq\sqrt{a_0g_N}$4

whereas some galaxy-favored interpolation functions predict substantially larger values. The interpolation-function parameters preferred by Cassini and by the SPARC RAR have been reported to differ at approximately the $g\simeq\sqrt{a_0g_N}$5 level (Desmond, 27 May 2025). This may indicate a need for screening, frequency-dependent dynamics, modified inertia, or a different relativistic formulation.

Non-radial dynamics is another unresolved issue. A MOND model that fits the radial rotation curve need not reproduce the vertical force of a galactic disk, the kinematics of radial stellar orbits, or the dynamics of bars and warps. In the Milky Way, comparisons have found that vertical velocity dispersions can favor a dark-matter model over MOND, although the result depends on baryonic structure and the specific MOND implementation.

Wide binaries provide a potentially sensitive test because their internal accelerations can enter the MOND regime while remaining in the Solar neighborhood. Some analyses report approximately Newtonian behavior, whereas others report a MOND-like enhancement. The disagreement has been attributed to unresolved companions, contamination, eccentricity distributions, sample selection, and error modeling. Consequently, the wide-binary test remains unsettled rather than constituting a decisive confirmation or refutation.

MOND’s cosmological status is weaker than its galaxy-scale status. The original framework contains no Friedmann equation, no complete perturbation theory, and no unique prediction for the CMB, matter power spectrum, or structure growth. Relativistic extensions can reproduce portions of the cosmological phenomenology, but often introduce additional fields or unseen matter. A pure baryon-only MOND cosmology historically encounters difficulty with the third acoustic peak, cluster formation, and the detailed growth of structure.

The theoretical status is also contested. MOND can be formally reduced to Newtonian gravity in the high-acceleration limit, but the interpolation function is not derived from a universally accepted fundamental theory. A philosophical analysis characterizes this as a failure of “reduction-wise justification”: formal limiting reduction is present, but downward theoretical grounding of the interpolation function and upward embedding in a unified all-scale mathematical structure are absent (Antoniou et al., 2 Dec 2025).

MOND therefore occupies an intermediate scientific status. Its galaxy-scale phenomenology is highly constrained and empirically successful: the BTFR, RAR, acceleration-linked onset of discrepancies, low-surface-brightness behavior, and many rotation-curve shapes follow from a common acceleration scale. Its difficulties involve clusters, lensing, Solar-System constraints, non-radial dynamics, wide binaries, and relativistic cosmology. The principal unresolved question is whether the MOND regularities represent a fundamental modification of gravity or inertia, an emergent property of a more complicated dark sector, or a regularity produced by baryon–halo coupling within $g\simeq\sqrt{a_0g_N}$6CDM.