Introduction

Galaxy rotation curves present one of the most striking puzzles in modern physics. When astronomers measure how fast stars orbit around galactic centers, they consistently find velocities that remain nearly constant far beyond the visible edge of the galaxy, where Newtonian gravity predicts they should fall off as the square root of distance. This discrepancy appears across galaxies of all sizes and types, from dwarf spheroidals to giant ellipticals, always becoming significant when the gravitational acceleration drops below a characteristic scale of approximately \(1.2 \times 10^{-10}~\text{m/s}^2\)—about ten billion times weaker than Earth’s surface gravity.

The standard explanation invokes dark matter: an invisible substance that surrounds galaxies in extended halos, providing the additional gravitational pull needed to keep stars moving fast. Dark matter models successfully reproduce many observations, but they require fine-tuning to explain certain regularities in the data. Most notably, the baryonic Tully-Fisher relation shows that the asymptotic rotation velocity \(v\) of a galaxy depends only on its visible (baryonic) mass \(M\) according to \(v^4 \propto M\), with little scatter. If dark matter halos had independent properties from the visible galaxy, we would expect this relation to show much more variation.

Milgrom  proposed an alternative interpretation: perhaps gravity itself behaves differently at very low accelerations. His Modified Newtonian Dynamics (MOND) hypothesis posits that when the Newtonian gravitational acceleration \(g_N = GM/r^2\) falls below a critical value \(\azero\), the actual acceleration experienced by a test particle transitions to \(g \approx \sqrt{g_N \cdot \azero}\). This simple modification, with \(\azero \approx 1.2 \times 10^{-10}~\text{m/s}^2\), automatically explains flat rotation curves, the Tully-Fisher relation, and numerous other galactic scaling relations . However, MOND has traditionally lacked a theoretical foundation—the acceleration scale \(\azero\) must be fitted to data rather than predicted from first principles.

The purpose of this paper is to provide a theoretical foundation for the MOND acceleration scale. We show that MOND phenomenology is consistent with the thermodynamics of de Sitter space, with the acceleration scale \(\azero = c\Hz/(2\pi)\) following from physically motivated assumptions about entropy dynamics. This connection between cosmological horizons and galactic dynamics, if correct, would reflect a deep relationship between large-scale and small-scale gravitational phenomena.

Our starting point is Jacobson’s result  that Einstein’s equations can be derived from thermodynamic principles. When one requires that the generalized entropy \(S_{\text{gen}} = A/(4\ell_P^2) + S_{\text{bulk}}\) remains stationary under small variations of the spacetime geometry, the Einstein equations follow as a consistency condition. This suggests that gravity is not a fundamental force but an emergent phenomenon arising from the equilibrium of quantum information distributed across spacetime.

Jacobson’s original derivation assumed flat spacetime at large distances. But our universe is not flat—it is described by de Sitter space, with an accelerating expansion driven by dark energy. The de Sitter geometry introduces a cosmological horizon at distance \(c/\Hz\) from any observer, beyond which events are causally disconnected. Like a black hole horizon, this cosmological horizon has thermodynamic properties: it radiates at the Gibbons-Hawking temperature \(\TdS = \hbar\Hz/(2\pi k_B)\) and carries entropy proportional to its area.

The central point of this paper is that the de Sitter horizon changes the entropy accounting. In flat spacetime, entropy obeys an area law: the maximum entropy in a region scales with its boundary area. But in de Sitter space, the thermal bath at temperature \(\TdS\) contributes volume-law entropy \(S_{\text{vol}} = s_\Lambda V\) throughout the interior. When a mass \(M\) is placed in this thermal bath, it “displaces” entropy \(\Delta S = Mc^2/\TdS\) from the surrounding volume. We postulate that if the available volume entropy cannot accommodate this displacement, local thermodynamic equilibrium fails, and the gravitational dynamics must change. This failure mechanism is physically motivated but not rigorously derived from first principles.

This failure of equilibrium occurs precisely when the Newtonian acceleration falls below \(\azero = c\Hz/(2\pi)\). At higher accelerations, the local geometry can absorb the entropy displacement through small curvature adjustments, and standard general relativity applies. At lower accelerations, the entropy deficit spreads non-locally through the de Sitter volume, producing the characteristic MOND behavior \(g \propto \sqrt{g_N}\). The transition is not put in by hand but emerges from the competition between area-law and volume-law contributions to the generalized entropy.

We develop this argument in stages. 2 establishes the thermodynamic properties of de Sitter space, including the Gibbons-Hawking temperature and the volume entropy density. 3 demonstrates that the total volume entropy saturates the Bousso covariant entropy bound, resolving any apparent conflict with holography. 4 derives the MOND acceleration scale from entropy equilibrium and shows how the Tully-Fisher relation emerges automatically. 5 provides the \(\epsilon\)-expansion that quantifies when general relativity versus MOND applies, showing that solar system tests are safely in the GR regime.

Having established the non-relativistic framework, we then present its covariant completion. 6 develops Entanglement-Elastic Gravity (EEG), a fully relativistic field theory in which an elastic displacement field \(\psi\) encodes the entropy strain caused by matter. The EEG action yields modified Einstein equations with an elastic stress tensor that reduces to the MOND Poisson equation in the Newtonian limit. The theory is ghost-free, preserves gravitational wave propagation at speed \(c\) (consistent with GW170817), and provides a complete framework for analyzing strong-field phenomena.

7 presents observational signatures. Already-confirmed results include flat rotation curves and the Tully-Fisher relation. Expected signatures—weak lensing slip \(|\Phi - \Psi|/|\Phi| \sim 15\%\) beyond 50 kpc, growth rate suppression \(\Delta f\sigma_8 \sim -0.03\) at \(z = 0.8\), and anomalous acceleration in wide binary stars detectable with Gaia DR4/DR5—are qualitative expectations from the modified gravity structure; quantitative derivations of the latter two are left for future work. 8 summarizes our findings and discusses open problems.

Series context. This paper belongs to the Quantum-Geometric Correspondence series. The canonical core paper  presents the complete axiomatic framework, develops gravitational decoherence from the same entropy principles, and establishes the information-theoretic bounds; the holographic dark-energy paper derives holographic dark energy; the entanglement-decay paper presents the unique entanglement-decoherence prediction. Each paper is self-contained but cross-references the others for extended discussion.

De Sitter Thermodynamics

The cosmological constant \(\Lambda\) is not merely a parameter in Einstein’s equations—it fundamentally alters the causal structure of spacetime, creating a horizon that radiates like a black hole turned inside out. This section establishes the thermodynamic properties of de Sitter space that will underpin our derivation of MOND phenomenology.

The Static Patch and Cosmological Horizon

De Sitter space is the maximally symmetric solution to Einstein’s equations with positive cosmological constant, \(G_{\mu\nu} + \Lambda g_{\mu\nu} = 0\). Unlike Minkowski space, de Sitter has a characteristic length scale set by \(\Lambda\) that determines how far any observer can see. In static coordinates centered on a freely falling observer, the metric takes the form \[\begin{equation} ds^2 = -\left(1 - \frac{r^2}{\RH^2}\right)c^2 dt^2 + \frac{dr^2}{1 - r^2/\RH^2} + r^2 d\Omega^2, \label{eq:desitter-metric} \end{equation}\] where the de Sitter radius is \[\begin{equation} \RH = \sqrt{\frac{3}{\Lambda}} = \frac{c}{\Hz}. \label{eq:horizon-radius} \end{equation}\]

The significance of \(\RH\) becomes apparent from the metric: the coefficient \(g_{tt}\) vanishes at \(r = \RH\), and the coefficient \(g_{rr}\) diverges. This marks a coordinate horizon—not a singularity, but a boundary beyond which events are causally disconnected from the central observer. Light signals sent from \(r > \RH\) can never reach \(r = 0\), just as signals sent into a black hole can never escape. For our universe with \(\Hz \approx 70~\text{km/s/Mpc}\), the horizon lies at \(\RH \approx 14\) billion light-years, defining the edge of the observable universe.

The region \(r < \RH\) is called the static patch because it admits a timelike Killing vector \(\partial/\partial t\). This symmetry allows us to define a conserved energy and construct a meaningful thermodynamics. The static patch contains a finite proper volume \(V_H = (4\pi/3)\RH^3\), bounded by the cosmological horizon. Unlike the exterior of a black hole, where freely falling observers cross the horizon inward, the de Sitter horizon is a future boundary—all observers eventually find their signals redshifted into oblivion as they approach it.

The Gibbons-Hawking Temperature

A central result of semiclassical gravity is that horizons radiate. Hawking showed that black holes emit thermal radiation at a temperature inversely proportional to their mass; Gibbons and Hawking  extended this result to cosmological horizons. The de Sitter horizon bathes its interior in thermal radiation at temperature \[\begin{equation} \TdS = \frac{\hbar \Hz}{2\pi k_B} \approx 2.76 \times 10^{-30}~\text{K}. \label{eq:GH-temperature} \end{equation}\] This extraordinarily low temperature—ten billion times colder than the cosmic microwave background—might seem irrelevant to practical physics. But as we shall see, it sets the acceleration scale below which gravitational dynamics fundamentally change.

The Gibbons-Hawking temperature can be derived by three independent methods, all yielding the same result. The first method uses surface gravity: the quantity \(\kappa_{\rm dS}\) that measures the “strength” of the horizon is defined geometrically by the Killing equation \(\nabla_\mu(k_\nu k^\nu) = -2\kappa_{\rm dS} k_\mu\). For the de Sitter metric, direct calculation gives \(\kappa_{\rm dS} = c\Hz\), and the temperature follows from the universal relation \(T = \hbar\kappa_{\rm dS}/(2\pi ck_B)\). (Note: the symbol \(\kappa\) without subscript is reserved for the elastic modulus \(\kappa \equiv c^4/(8\pi G\azero)\) introduced in Section 6.)

The second method uses Euclidean continuation. Under the analytic continuation \(t \to -i\tau\), the Lorentzian metric becomes Riemannian, but a conical singularity appears at \(r = \RH\) unless \(\tau\) is periodic with a specific period. Requiring regularity uniquely fixes this period to be \(\beta = 2\pi/\Hz\), which corresponds to temperature \(T = \hbar/(2\pi k_B\beta) = \hbar\Hz/(2\pi k_B)\).

The third method appeals to the Unruh effect. A static observer at any radius \(r < \RH\) must accelerate to resist the cosmic expansion. As \(r \to \RH\), this proper acceleration diverges, and the Unruh radiation seen by the accelerating observer approaches the Gibbons-Hawking temperature after accounting for gravitational redshift. That three completely different approaches—geometric, topological, and kinematic—yield the identical answer strongly suggests that \(\TdS\) is a fundamental property of de Sitter space.

The First Law and Horizon Entropy

The cosmological horizon is not merely a kinematic boundary but a genuine thermodynamic system. Like a black hole, it carries entropy proportional to its area—an entropy that in our universe reaches the staggering value \(S_{dS} = A_H/(4\lP^2) \approx 10^{122}\) in natural units. This is the largest entropy of any system in the observable universe, dwarfing even the combined entropy of all black holes.

The thermodynamic nature of the horizon is expressed by the first law \[\begin{equation} dE_{dS} = \TdS\, dS_{dS}, \label{eq:first-law} \end{equation}\] where \(E_{dS}\) is the vacuum energy contained within the static patch. This equation states that adding energy to the interior increases the horizon area, just as feeding a black hole makes it grow. The cosmological horizon acts as a heat reservoir at temperature \(\TdS\), in thermal equilibrium with the interior.

The Characteristic Energy and Acceleration Scales

What physical scale does the Gibbons-Hawking temperature set? Converting temperature to energy gives \[\begin{equation} E_{dS} = k_B\TdS = \frac{\hbar\Hz}{2\pi} \approx 3.6 \times 10^{-52}~\text{J}. \end{equation}\] This is an absurdly small energy—roughly the kinetic energy of a single atom moving at \(10^{-12}\) meters per second. No laboratory experiment could ever detect individual quanta at this scale.

Yet when we ask what acceleration this energy corresponds to, a familiar scale appears. Using natural gravitational units defined by the Planck mass \(\mP = \sqrt{\hbar c/G}\) and Planck length \(\lP = \sqrt{\hbar G/c^3}\), we find \[\begin{equation} \azero = \frac{E_{dS}}{\mP\lP} = \frac{c\Hz}{2\pi} \approx 1.08 \times 10^{-10}~\text{m/s}^2. \end{equation}\] This is the MOND acceleration scale, the value below which galaxy rotation curves deviate from Newtonian predictions. The coincidence \(\azero \sim c\Hz\) has puzzled physicists since Milgrom first noted it in 1983; here we see it is not a coincidence at all, but a direct consequence of de Sitter thermodynamics. The Gibbons-Hawking temperature of our universe determines the acceleration scale at which gravity changes character.

Volume-Law Entropy in de Sitter Space

Holography has taught us that gravitational systems carry entropy proportional to their boundary area, not their volume. A black hole’s entropy scales as \(r^2\), not \(r^3\), regardless of what fell in to form it. This area law is often taken as a fundamental principle, suggesting that the degrees of freedom in any gravitational system are somehow “painted” on surfaces rather than distributed through space.

De Sitter space challenges this picture. The thermal radiation from the cosmological horizon fills the interior with entropy that scales with volume, not area. This section shows that this volume-law contribution is not only consistent with holography but required by it—the volume entropy exactly saturates the Bousso covariant entropy bound at the cosmological scale.

The de Sitter Vacuum as a Thermal State

To understand the origin of volume entropy, we must examine how quantum field theory behaves when the spacetime itself has a horizon. Consider the vacuum state of a quantum field in de Sitter space. The full Hilbert space naturally factorizes as \(\mathcal{H}_{\text{full}} = \mathcal{H}_{\text{in}} \otimes \mathcal{H}_{\text{out}}\), where \(\mathcal{H}_{\text{in}}\) contains modes inside the cosmological horizon and \(\mathcal{H}_{\text{out}}\) contains modes beyond it.

What we call “the vacuum” from a global perspective is not a vacuum at all from the perspective of an observer confined to the static patch. The global vacuum state is entangled across the horizon: \[\begin{equation} |0\rangle = \sum_n e^{-\pi E_n/\hbar\Hz} |n\rangle_{\text{in}} \otimes |n\rangle_{\text{out}}. \end{equation}\] When the observer traces over the inaccessible super-horizon modes, the remaining state is thermal: \[\begin{equation} \rho_{\text{in}} = \text{Tr}_{\text{out}}|0\rangle\langle 0| = \frac{1}{Z}\sum_n e^{-2\pi E_n/\hbar\Hz}|n\rangle_{\text{in}}\langle n|. \end{equation}\] This is a Boltzmann distribution at the Gibbons-Hawking temperature \(\TdS = \hbar\Hz/(2\pi k_B)\). The de Sitter horizon does not merely radiate; it holds its interior in thermal equilibrium with itself.

The Volume Entropy Density

A thermal state at temperature \(T\) carrying energy density \(\rho\) has entropy density \(s = \rho/T\). For de Sitter space, the vacuum energy density is \[\begin{equation} \rho_\Lambda = \frac{\Lambda c^4}{8\pi G} = \frac{3c^2\Hz^2}{8\pi G}, \end{equation}\] where the second equality uses \(\Lambda = 3\Hz^2/c^2\) for a flat universe approaching de Sitter at late times. Dividing by the Gibbons-Hawking temperature gives the entropy density \[\begin{equation} s_\Lambda= \frac{\rho_\Lambda}{\TdS} = \frac{3c^2\Hz^2}{8\pi G} \cdot \frac{2\pi k_B}{\hbar\Hz} = \frac{3c^2\Hz k_B}{4G\hbar}. \label{eq:sLambda} \end{equation}\]

In numerical terms, \(s_\Lambda\approx 2.2 \times 10^{43}~k_B\,\text{m}^{-3}\)—an enormous entropy density by ordinary standards. A cubic meter of de Sitter vacuum carries more entropy than most astrophysical systems. Yet because this entropy is spread uniformly through space, it is invisible to local experiments; it manifests only when we consider gravitational dynamics on scales where significant volume is enclosed.

Total Entropy: Area Plus Volume

A spherical region of radius \(r\) inside the static patch therefore carries entropy from two sources: the standard area-law contribution from entanglement across its boundary, and the volume-law contribution from the thermal bath: \[\begin{equation} S_{dS}(r) = \frac{A(r)c^3}{4G\hbar} + s_\Lambda V(r) = \frac{\pi r^2 c^3}{G\hbar} + \frac{\pi c^2\Hz k_Br^3}{G\hbar}. \label{eq:total-entropy} \end{equation}\] The area term dominates at small scales; the volume term dominates at large scales. The crossover occurs when \(r \sim c/\Hz = \RH\)—precisely at the cosmological horizon.

Consistency with Holography

The volume-law entropy raises an apparent paradox: for sufficiently large regions, \(S_{\text{vol}}> \Sarea\), seemingly violating the holographic principle that entropy cannot exceed \(A/(4\lP^2)\). This contradiction is resolved by recognizing that the naive statement of holography is incorrect.

The correct formulation is the Bousso covariant entropy bound . For any surface \(B\), construct a light sheet—a null hypersurface generated by non-expanding null geodesics orthogonal to \(B\). The bound states that the entropy on this light sheet satisfies \(S \leq A(B)/(4\lP^2)\). Crucially, the bound applies to light sheets, not arbitrary spatial volumes. A spatial volume can exceed the area bound; a light sheet cannot.

For de Sitter space, an exact saturation occurs when we apply the Bousso bound at the cosmological horizon itself. The total volume entropy within the static patch equals \[\begin{equation} S_{\text{vol}}^{\text{total}} = s_\Lambda\cdot V_H = \frac{3c^2\Hz k_B}{4G\hbar} \cdot \frac{4\pi c^3}{3\Hz^3} = \frac{\pi c^5k_B}{G\hbar\Hz^2}, \end{equation}\] while the horizon entropy is \[\begin{equation} S_H = \frac{A_H}{4\lP^2} = \frac{4\pi\RH^2 \cdot c^3}{4G\hbar} = \frac{\pi c^5}{G\hbar\Hz^2}. \end{equation}\] These are identical (in natural units where \(k_B= 1\)). The volume entropy exactly saturates the Bousso bound—it does not violate holography because it is the holographic entropy of the horizon, expressed in bulk variables. The thermal entropy filling the interior is the same entropy that lives on the boundary, just described in a different language.

This saturation has physical content: it means that de Sitter space is a maximally entropic configuration, unable to hold more information without collapsing into a state with a different cosmological constant. Any matter placed inside must “displace” entropy from this thermal bath, and when the displacement becomes significant, the gravitational dynamics change. This is the origin of MOND phenomenology, as we develop in the next section.

Emergence of the MOND Acceleration Scale

Having established that de Sitter space carries extensive volume entropy at the Gibbons-Hawking temperature, we now present the central result of this paper: the MOND acceleration scale emerges as a consequence of entropy equilibrium arguments, with the characteristic scale determined entirely by cosmological parameters.

Entropy Displacement by Matter

The first step is to understand how matter interacts with the de Sitter thermal bath. Consider placing a mass \(M\) at rest in a region of de Sitter space. The mass represents bound energy \(E = Mc^2\) that has been “removed” from the thermal background. By the first law of thermodynamics, removing energy \(dE\) from a system at temperature \(T\) reduces its entropy by \(dS = dE/T\). The mass therefore creates an entropy displacement \[\begin{equation} \Delta S(M) = \frac{Mc^2}{\TdS} = \frac{2\pi k_BMc^2}{\hbar\Hz}. \label{eq:entropy-displacement} \end{equation}\] This is not a small effect. The de Sitter temperature \(\TdS \approx 2.76 \times 10^{-30}\) K is extraordinarily low, so even modest masses create enormous entropy deficits relative to the thermal background.

To develop physical intuition, consider a Milky Way-sized galaxy with baryonic mass \(M \approx 10^{11} M_\odot\). Substituting into Eq. \(\eqref{eq:entropy-displacement}\) yields \(\Delta S \approx 10^{99}k_B\)—an incomprehensibly large number. This entropy debt cannot simply vanish; it must be accommodated by rearrangements in the surrounding de Sitter geometry. The question is whether the available entropy budget is sufficient.

The Entropy Budget

The de Sitter volume entropy, derived in Section 3, fills space with an entropy density \(s_\Lambda= 3c^2\Hz k_B/(4G\hbar)\). Within a spherical region of radius \(r\), the available entropy is \[\begin{equation} S_{\text{available}}(r) = s_\Lambda\cdot \frac{4\pi r^3}{3} = \frac{\pi c^2 \Hz k_Br^3}{G\hbar}. \end{equation}\] This grows as \(r^3\), so sufficiently large regions can accommodate any finite entropy displacement. The crucial question is: how large must the region be?

Local thermodynamic equilibrium requires that the entropy displaced by the mass can be absorbed within the region where the mass exerts significant gravitational influence. If we define the “gravitational influence region” as the volume where the Newtonian potential exceeds some threshold, equilibrium demands \[\begin{equation} S_{\text{available}}(r) \geq S_{\text{displaced}}(M). \end{equation}\] When this inequality is satisfied, small curvature adjustments can restore equilibrium, and standard general relativity applies. When it fails—when the entropy deficit exceeds what the local volume can provide—the excess must spread non-locally, fundamentally altering the gravitational dynamics.

The Critical Radius

The boundary between these regimes occurs at the critical radius \(r_c\) where the available and displaced entropies are exactly equal: \[\begin{equation} S_{\text{displaced}}(M) = S_{\text{available}}(r_c). \end{equation}\] Substituting the expressions derived above: \[\begin{equation} \frac{2\pi k_BMc^2}{\hbar\Hz} = \frac{\pi c^2 \Hz k_Br_c^3}{G\hbar}. \end{equation}\] The factors of \(\pi\), \(k_B\), \(c^2\), and \(\hbar\) cancel, leaving a simple result: \[\begin{equation} r_c= \left(\frac{2GM}{\Hz^2}\right)^{1/3}. \label{eq:critical-radius} \end{equation}\]

This critical radius has a natural physical interpretation: it marks the distance at which the de Sitter expansion velocity \(v = \Hz r\) equals the escape velocity \(v_{\text{esc}} = \sqrt{2GM/r}\) from the central mass. Beyond \(r_c\), the cosmic expansion dominates over the local gravitational binding, and the mass can no longer be treated as an isolated system in flat spacetime.

The numerical values illustrate the relevant scales. For the Sun (\(M = M_\odot\)), the critical radius is \(r_c\approx 394\) light-years—far beyond the solar system but well within the Milky Way, confirming that solar-system dynamics lie safely in the Newtonian regime. For a typical galaxy (\(M = 10^{11} M_\odot\)), the formula gives \(r_c\approx 560\) kpc, which exceeds the optical radius where flat rotation curves are observed (typically 10–50 kpc). This means the outer rotation-curve region lies well within \(r_c\): the system is in the entropy-equilibrium regime but at low enough acceleration (\(g_N\lesssim \azero\)) to show MOND dynamics, since \(r_c\) marks the boundary of Newtonian dominance, not the onset of MOND. For a galaxy cluster (\(M = 10^{14} M_\odot\)), the critical radius extends to \(r_c\approx 5.6\) Mpc. These scales are broadly consistent with where non-Newtonian behavior is observed, though the quantitative correspondence requires the full \(\epsilon\)-expansion of Section 5.

The Universal Acceleration Scale

The critical radius depends on the mass \(M\), so different objects have different values of \(r_c\). However, the acceleration at the critical radius reveals a universal scale that depends only on cosmological parameters.

The Newtonian gravitational acceleration at radius \(r_c\) is \[\begin{equation} a_c = \frac{GM}{r_c^2} = \frac{GM}{\left(2GM/\Hz^2\right)^{2/3}} = \frac{(GM)^{1/3}\Hz^{4/3}}{2^{2/3}}. \end{equation}\] This still depends on \(M\), so \(a_c\) itself is not universal. But the characteristic scale that governs the transition emerges from the intrinsic properties of the de Sitter thermal bath, independent of any particular mass.

The de Sitter horizon radiates at temperature \(\TdS = \hbar\Hz/(2\pi k_B)\), which sets a natural energy scale \(E_{\text{dS}} = k_B\TdS = \hbar\Hz/(2\pi)\). To convert this energy to an acceleration, we use the Planck mass \(\mP = \sqrt{\hbar c/G}\) and Planck length \(\lP = \sqrt{\hbar G/c^3}\) as the natural units of mass and length: \[\begin{equation} \azero = \frac{E_{\text{dS}}}{\mP \lP} = \frac{\hbar\Hz/(2\pi)}{\sqrt{\hbar c/G} \cdot \sqrt{\hbar G/c^3}} = \frac{c\Hz}{2\pi}. \label{eq:a0-derivation} \end{equation}\]

This is our central result: the MOND acceleration scale is determined entirely by fundamental constants and the present-day Hubble parameter. Two logically distinct steps enter here: the entropy-budget argument establishes that \(cH\) is the relevant acceleration scale, while the specific coefficient \(1/(2\pi)\) is the Gibbons-Hawking temperature assignment \(\TdS = \hbar\Hz/(2\pi k_B)\)—a dimensional identification rather than an independent derivation, in the same spirit in which Bekenstein separated the black-hole area law from the value of its coefficient.

Numerically, using \(c = 2.998 \times 10^8\) m/s and \(\Hz = 2.268 \times 10^{-18}\) s\(^{-1}\) (corresponding to \(H_0 = 70\) km/s/Mpc): \[\begin{equation} \azero = \frac{(2.998 \times 10^8)(2.268 \times 10^{-18})}{2\pi} \approx 1.08 \times 10^{-10}~\text{m/s}^2. \end{equation}\] The observed MOND scale from galaxy rotation curve fits is \(\azero^{\text{obs}} \approx (1.2 \pm 0.2) \times 10^{-10}~\text{m/s}^2\) . The agreement to within 10% holds for a parameter-free prediction connecting cosmology to galactic dynamics.

Modified Gravitational Dynamics

Having identified the critical acceleration scale, we now derive how gravitational dynamics change when \(g_N< \azero\). The entropy equilibrium condition \(\delta\Sgen = 0\) must include the volume entropy contribution when local equilibrium fails.

The physical picture is as follows. At high accelerations (\(g_N\gg \azero\)), the entropy displaced by a mass is easily accommodated within a small region, and the standard Einstein equations apply. At low accelerations (\(g_N\ll \azero\)), the entropy deficit cannot be localized; it spreads through the de Sitter volume, creating a non-local “strain” in the entropy distribution. This strain generates an additional contribution to the gravitational field that supplements the Newtonian term.

The mathematical expression of this physics is the modified Poisson equation: \[\begin{equation} \nabla \cdot \left[\mu\left(\frac{|\nabla\Phi|}{\azero}\right)\nabla\Phi\right] = 4\pi G\rho, \label{eq:modified-poisson} \end{equation}\] where \(\Phi\) is the gravitational potential, \(\rho\) is the matter density, and \(\mu(x)\) is an interpolation function satisfying

• \(\mu(x) \to 1\) for \(x \gg 1\) (Newtonian regime), and

• \(\mu(x) \to x\) for \(x \ll 1\) (deep MOND regime).

The form of the interpolation function encodes how the transition between regimes occurs. The entropy arguments specify the limiting behaviors but not the detailed shape of \(\mu(x)\). Phenomenologically, the standard ansatz \(\mu(x) = x/\sqrt{1 + x^2}\) provides excellent fits to rotation curve data . Deriving \(\mu(x)\) rigorously from the microscopic entropy dynamics remains an open problem—this is a significant gap in the theoretical framework that we acknowledge explicitly.

The effective gravitational field \(g_{\text{eff}}= -\nabla\Phi\) takes different forms in the two regimes. In the Newtonian regime where \(g_N\gg \azero\): \[\begin{equation} g_{\text{eff}}= g_N= \frac{GM}{r^2}, \end{equation}\] recovering Newton’s inverse-square law exactly. In the deep MOND regime where \(g_N\ll \azero\): \[\begin{equation} g_{\text{eff}}= \sqrt{g_N\cdot \azero} = \sqrt{\frac{GM\azero}{r^2}} = \frac{\sqrt{GM\azero}}{r}. \end{equation}\] The acceleration now falls as \(1/r\) rather than \(1/r^2\), which directly produces flat rotation curves: if \(v^2/r = g_{\text{eff}} \propto 1/r\), then \(v = \text{constant}\).

The Baryonic Tully-Fisher Relation

The most celebrated success of MOND phenomenology is the baryonic Tully-Fisher relation, which relates the asymptotic rotation velocity of a disk galaxy to its total baryonic mass. Our framework predicts this relation with no free parameters.

For a test particle in circular orbit at radius \(r\) in the deep MOND regime, the centripetal acceleration equals the gravitational field: \[\begin{equation} \frac{v^2}{r} = g_{\text{eff}}= \sqrt{\frac{GM\azero}{r^2}}. \end{equation}\] Squaring both sides and rearranging: \[\begin{equation} \frac{v^4}{r^2} = \frac{GM\azero}{r^2}. \end{equation}\] The factors of \(r^2\) cancel, yielding the mass-independent relation \[\begin{equation} v^4 = GM\azero = GM \cdot \frac{c\Hz}{2\pi}. \label{eq:tully-fisher} \end{equation}\]

This is the baryonic Tully-Fisher relation with normalization completely fixed by cosmology. The relation \(v^4 \propto M\) holds exactly in the deep MOND limit, with slope 4 on a log-log plot of velocity versus mass. The observed slope is \(3.98 \pm 0.12\) , in excellent agreement with the prediction.

To verify the normalization, consider a galaxy with baryonic mass \(M = 10^{11} M_\odot = 2 \times 10^{41}\) kg. Substituting into Eq. \(\eqref{eq:tully-fisher}\): \[\begin{align} v^4 &= (6.67 \times 10^{-11}~\text{m}^3\text{kg}^{-1}\text{s}^{-2})(2 \times 10^{41}~\text{kg})(1.08 \times 10^{-10}~\text{m/s}^2) \nonumber \\ &\approx 1.5 \times 10^{21}~\text{m}^4/\text{s}^4. \end{align}\] Taking the fourth root: \(v \approx 200\) km/s, matching the observed rotation velocities of Milky Way-sized galaxies.

The tight scatter observed in the Tully-Fisher relation now has a natural explanation. In dark matter models, the relation arises from a correlation between visible and dark matter that must be imposed by hand. In our framework, the relation is fundamental: it follows directly from the thermodynamic properties of de Sitter space, with no room for scatter from varying dark matter halo properties.

Regime Separation: GR versus MOND

A theory that modifies gravity at low accelerations must explain why no deviations from general relativity are observed in the solar system, where precision tests constrain any departure to parts per billion. This section develops a systematic framework for understanding when GR applies exactly and when MOND corrections become significant.

The Control Parameter

The natural dimensionless parameter that governs the transition between regimes is the ratio of the local Newtonian acceleration to the MOND scale: \[\begin{equation} \epsilon(r) \equiv \frac{g_N(r)}{\azero} = \frac{GM}{r^2\azero} = \frac{2\pi GM}{r^2 c\Hz}. \label{eq:epsilon} \end{equation}\] When \(\epsilon \gg 1\), Newtonian gravity dominates and cosmological effects are negligible; when \(\epsilon \lesssim 1\), the volume entropy becomes dynamically important and MOND corrections emerge.

To see this quantitatively, consider how the different contributions to the generalized entropy scale. The entropy equilibrium condition \(\delta\Sgen = 0\) involves three terms: the area entropy \(\delta\Sarea\), the bulk matter entropy \(\delta\Sbulk\), and the de Sitter volume entropy \(\delta S_{\text{vol}}\). For a spherically symmetric perturbation at radius \(r\), dimensional analysis gives \[\begin{align} \delta\Sarea &\sim \frac{r}{G\hbar}\,\delta r, \\ \delta\Sbulk &\sim \frac{c^2 r^2 \epsilon}{G\hbar}\,\delta r, \\ \delta S_{\text{vol}}&\sim \frac{c^2\Hz r^2}{G\hbar}\,\delta r. \end{align}\] The ratio of volume to bulk contributions scales as \(\delta S_{\text{vol}}/\delta\Sbulk \sim 1/\epsilon\). This motivates organizing the equilibrium condition as an expansion in inverse powers of \(\epsilon\): \[\begin{equation} \delta\Sgen = \underbrace{\delta\Sarea + \delta\Sbulk}_{O(1)} + \underbrace{\epsilon^{-1}\delta S_{\text{vol}}^{(1)}}_{O(\epsilon^{-1})} + O(\epsilon^{-2}). \label{eq:epsilon-expansion} \end{equation}\]

The Three Dynamical Regimes

The expansion reveals three distinct physical regimes, each with characteristic gravitational dynamics.

In the high-acceleration regime where \(\epsilon \gg 1\), the \(1/\epsilon\) terms are negligible. The equilibrium condition reduces to \(\delta\Sarea + \delta\Sbulk = 0\), which is precisely Jacobson’s derivation of Einstein’s equations. General relativity emerges exactly, with corrections suppressed by \(\azero/g_N\). This is why solar system tests see pure GR: at Earth’s surface where \(g_N\approx 10~\text{m/s}^2\), we have \(\epsilon \approx 10^{11}\), and any MOND correction is at the level of \(10^{-11}\)—far below measurement precision.

In the transition regime where \(\epsilon \sim 1\), all three contributions compete comparably. The equilibrium condition \(\delta\Sarea + \delta\Sbulk + \delta S_{\text{vol}}= 0\) yields interpolated dynamics described by \[\begin{equation} g_{\text{eff}}\cdot \mu\left(\frac{g_{\text{eff}}}{\azero}\right) = g_N, \end{equation}\] where \(\mu(x)\) is an interpolation function that smoothly connects the two limiting behaviors.

In the low-acceleration regime where \(\epsilon \ll 1\), the volume term dominates and the equilibrium condition gives the deep MOND result \(g_{\text{eff}}= \sqrt{g_N\azero}\). The effective gravitational field now falls as \(1/r\) rather than \(1/r^2\), producing flat rotation curves.

Precision Tests in the Solar System

The framework must explain why no MOND effects appear in precision solar system tests. The answer lies in the enormous values of \(\epsilon\) throughout the inner solar system. At Earth’s surface, where \(g_N\approx 10~\text{m/s}^2\), the ratio is \(\epsilon \approx 10^{11}\), suppressing any MOND correction to the level of \(10^{-11}\). At Mercury’s perihelion, \(\epsilon \approx 4 \times 10^8\), predicting deviations at the \(10^{-9}\) level—safely below the precision of perihelion precession measurements. Even at Voyager 2’s current distance of 120 AU, where \(g_N\approx 4 \times 10^{-6}~\text{m/s}^2\), we have \(\epsilon \approx 10^4\), keeping corrections at \(10^{-4}\)—below the precision of deep-space tracking.

The critical radius \(r_c\) at which \(\epsilon = 1\) marks where the inverse-square law begins to fail. For the Sun, this radius is \[\begin{equation} r_c= \left(\frac{2G\Msun}{\Hz^2}\right)^{1/3} \approx 394~\text{light years}. \end{equation}\] The entire solar system lies well within this radius, explaining why Newtonian gravity appears exact. For a typical galaxy with \(M = 10^{11}\Msun\), the formula gives \(r_c\approx 560\) kpc, which is larger than the region probed by rotation curves (typically \(<50\) kpc). The rotation-curve regime lies inside \(r_c\) but at low enough local acceleration (\(g_N\lesssim \azero\)) to exhibit MOND dynamics; the \(\epsilon\)-expansion of Section 5 shows this is precisely where volume-entropy corrections dominate. The critical radius is thus a necessary but not sufficient condition: MOND effects appear when \(\epsilon \lesssim 1\), which occurs within \(r_c\) for sufficiently diffuse mass distributions.

The Spectrum of Gravitational Environments

Astrophysical systems span an enormous range in \(\epsilon\), from \(10^{13}\) in binary pulsars to \(0.01\) in cosmic voids. The framework makes specific predictions for each environment: GR should be exact to parts in \(10^{7}\)–\(10^{13}\) in the solar system and pulsars; near-Newtonian behavior with small corrections should appear in galactic bulges (\(\epsilon \sim 10\)); full MOND dynamics should operate in galaxy outskirts and halos (\(\epsilon \lesssim 1\)); and deep MOND with possible higher-order corrections should characterize galaxy clusters (\(\epsilon \sim 0.3\)) and void regions (\(\epsilon \sim 0.01\)). No tuning is required to achieve this hierarchy—it emerges automatically from the single scale \(\azero = c\Hz/(2\pi)\).

Cosmological Evolution of the MOND Scale

The robust statement of this framework is that \(\azero\) is locked to a cosmological horizon rate, \(\azero = cH/(2\pi)\) (the horizon-clock lock; see the Falsification Criteria in Section 7). Whether \(\azero\) then evolves with cosmic time depends on which horizon rate enters: in the instantaneous-horizon reading the scale was larger in the past, whereas the alternative event-horizon reading gives a constant \(\azero\). Current high-redshift rotation curves favour an approximately constant \(\azero\) over \(z = 0\)–\(2\), so we present the evolving case below as the signature that would distinguish the instantaneous reading rather than as an established prediction. In that reading, since \(\azero = cH(z)/(2\pi)\) and the instantaneous Hubble parameter changes as the universe expands, \[\begin{equation} \azero(z) = \frac{cH(z)}{2\pi} = \azero(0) \cdot \sqrt{\Omega_m(1+z)^3 + \Omega_\Lambda}. \label{eq:a0-z} \end{equation}\] At redshift \(z = 1\), the MOND scale was about 1.78 times larger than today; at \(z = 3\), about 4.54 times larger; at recombination (\(z \sim 1100\)), roughly \(2.4 \times 10^{4}\) times larger (including the radiation contribution \(\Omega_r \approx 9.4 \times 10^{-5}\)).

This evolution has direct observational consequences. At higher redshift, with a larger \(\azero\), more systems fall into the Newtonian regime. Rotation curve anomalies should appear at smaller radii for high-redshift galaxies, and extremely high-redshift systems may appear entirely Newtonian. This is a testable prediction: surveys with JWST and Euclid can measure rotation curves at \(z \sim 1\)–\(3\) and check whether the transition radius scales as predicted.

Consistency with the Cosmic Microwave Background

The framework’s most stringent test comes from the cosmic microwave background. At recombination, \(\azero(z=1100) \approx 2.5 \times 10^{-6}~\text{m/s}^2\)—much larger than the present-day value. The critical radius for a proto-cluster shrinks to \(\sim\!1\) kpc, while the acoustic oscillations occur on scales of 150 Mpc. The entire CMB physics operates deep in the Newtonian regime, where our framework reduces exactly to standard gravity.

This resolves a long-standing objection to MOND: how can a theory that modifies gravity explain the CMB acoustic peaks without dark matter? The answer is that the CMB probes accelerations far above \(\azero(z_{\text{rec}})\), where the framework predicts no modification. The acoustic peaks are consistent with standard \(\Lambda\)CDM precisely because that is what our framework predicts at early times.

The Galaxy Cluster Tension

Galaxy clusters present a known challenge for MOND: even after accounting for MOND dynamics, clusters show a mass discrepancy of factor 2–5 between the dynamically inferred mass and the observed baryonic mass. This residual discrepancy has been used to argue that MOND cannot fully replace dark matter.

Within our framework, clusters occupy the transition regime with \(\epsilon \sim 0.1\)–1, where the \(\epsilon\)-expansion requires higher-order terms. The leading MOND correction scales as \(\epsilon^{-1}\), but subleading corrections proportional to \(\epsilon^{-2}\) and higher become non-negligible when \(\epsilon\) approaches unity. These higher-order terms are sensitive to the density gradient \(\nabla^2\rho/\rho\), not just the local density, making them significant in the extended, non-uniform environment of galaxy clusters.

Several additional effects may contribute to the cluster tension. Galaxy clusters are triaxial, not spherical, and the departure from spherical symmetry affects the entropy equilibrium in ways not captured by our spherically symmetric analysis. The intracluster medium at \(T \sim 10^8\) K contributes its own thermodynamics that may modify the effective temperature governing entropy balance. Neutrinos with \(\sum m_\nu \sim 0.1\) eV contribute 1–5% of the cluster mass, partially closing the gap. A complete resolution requires detailed calculations of all these effects, which remains an important direction for future work.

Covariant Completion: Entanglement-Elastic Gravity

The preceding sections developed a non-relativistic framework in which MOND phenomenology emerges from the thermodynamics of de Sitter space. While this successfully explains galactic dynamics and predicts the acceleration scale \(\azero = c\Hz/(2\pi)\), a complete theory requires fully covariant field equations that reduce to the modified Poisson equation in the appropriate limit. This section sketches Entanglement-Elastic Gravity (EEG), a candidate relativistic extension. We present it as a schematic proposal, not a finished theory: as we make explicit in Section 6.5, the quadratic elastic action used here is a linearized caricature, and a consistent covariant completion that reproduces the deep-MOND limit requires a nonlinear (Bekenstein–Milgrom/AQUAL-type) kinetic functional . What we set out below is the structure of such a theory, the elastic-modulus matching it must satisfy, and the open problems a full completion must solve.

Conceptual Foundation

The key insight underlying EEG is that spacetime exhibits two distinct but interacting forms of quantum entanglement, each contributing to the gravitational dynamics in different regimes.

The first form is local entanglement, which obeys an area law. Across any infinitesimal surface in spacetime, quantum degrees of freedom form short-range Bell pairs, contributing entropy proportional to the surface area: \(S_{\text{area}} = A/(4\ell_P^2)\). This is the entanglement that Jacobson  showed gives rise to Einstein’s equations when one demands that the generalized entropy remain stationary. The area-law network encodes the local curvature of spacetime and is described mathematically by the metric tensor \(g_{\mu\nu}\).

The second form is global entanglement, which obeys a volume law. The de Sitter horizon at distance \(c/\Hz\) creates a thermal bath that fills the interior with long-range entanglement at density \(s_\Lambda= 3c^2\Hz k_B/(4G\hbar)\). This volume-law network has no analog in flat spacetime—it is a genuinely cosmological contribution that arises from the finite size of the observable universe.

Matter interacts with both networks simultaneously. When a mass \(M\) is placed in spacetime, it requires area-law equilibrium on surrounding surfaces (producing local curvature) while also displacing volume-law entropy \(\Delta S \sim Mc^2/\TdS\) from the thermal background (creating global strain). At high accelerations where \(g \gg \azero\), the area-law contribution dominates, and standard general relativity applies. At low accelerations where \(g \lesssim \azero\), the volume-law strain becomes significant, producing the MOND-like modifications we observe in galactic dynamics. The critical scale \(\azero = c\Hz/(2\pi)\) marks the crossover where both contributions are comparable.

To describe this physics mathematically, we introduce a scalar field \(\psi(x)\) that quantifies the local displacement of volume-law entanglement. Where matter is present, \(\psi\) is large; in empty regions far from matter, \(\psi\) approaches zero. The gradients of \(\psi\) represent the “strain” in the entanglement network, and this strain contributes to the gravitational field equations through an elastic stress tensor.

Microscopic Model

Before developing the field equations, it is helpful to have a concrete microscopic picture of how entanglement elasticity arises. We present a toy model that captures the essential physics while remaining agnostic about the ultraviolet completion.

Consider a closed spacelike surface \(\mathcal{H}\) of area \(A\). We discretize this surface into \(N = A/\ell_P^2\) Planck-sized plaquettes, each carrying a pair of “horizon qubits” \((q_i^{(L)}, q_i^{(R)})\). In the vacuum state, these qubits are maximally entangled across the surface: \[\begin{equation} \ket{\Psi_{\text{area}}} = \bigotimes_{i=1}^{N} \ket{\Phi^+}_i, \quad \text{where} \quad \ket{\Phi^+}_i = \frac{\ket{00}_i + \ket{11}_i}{\sqrt{2}}. \end{equation}\] This state has von Neumann entropy \(S = N \log 2 = A/(4\ell_P^2)\), reproducing the Bekenstein-Hawking formula.

Injecting a mass \(M\) at the center of this surface disrupts the entanglement pattern. The mass effectively “removes” a number of Bell pairs proportional to its rest energy: \[\begin{equation} n = \frac{2\pi Mc}{\hbar \Hz}. \end{equation}\] The removed pairs create an entanglement deficit that spreads through the surrounding space. If we coarse-grain over many plaquettes, this deficit can be described by a continuous displacement field \(\psi\) whose gradient measures the local strain. The elastic energy associated with this strain takes the standard quadratic form: \[\begin{equation} E_{\text{el}} = \frac{\kappa}{2}\int (\nabla\psi)^2 \, d^3x, \end{equation}\] where dimensional analysis combined with the requirement that the elastic force reproduce MOND phenomenology fixes the modulus to \[\begin{equation} \kappa = \frac{c^4}{8\pi G \azero}. \end{equation}\]

This toy model has obvious limitations—it ignores gauge constraints, assumes instantaneous equilibration, and leaves the scrambling dynamics unspecified. Nevertheless, it achieves three important goals: it shows how area-law and volume-law entanglement can coexist in a single quantum state, it reproduces the entropy displacement formula underlying MOND, and it predicts the correct elastic modulus without fine-tuning. These features make it a useful guide for constructing the covariant theory.

Effective Action

The covariant formulation of EEG follows from an action principle that enforces entanglement equilibrium. The total action consists of three parts: the Einstein-Hilbert term for gravity, the matter Lagrangian, and an elastic sector for the displacement field \(\psi\): \[\begin{align} S_{\text{eff}} &= \frac{1}{16\pi G}\int R\sqrt{-g}\,d^4x + \int \mathcal{L}_{\text{matter}}\sqrt{-g}\,d^4x \nonumber\\ &\quad + \int\left[-\frac{c^4}{8\pi G \azero}(\partial\psi)^2 - \lambda\left(\psi - \frac{8\pi G}{\azero c^2}\rho_{\text{rest}}\right)\right]\sqrt{-g}\,d^4x. \label{eq:Seff} \end{align}\] Here \(\lambda(x)\) is a Lagrange multiplier that enforces the constraint relating \(\psi\) to the matter density. The coefficient \(c^4/(8\pi G\azero)\) in the kinetic term is fixed by requiring consistency with the non-relativistic limit.

Several features of this action merit comment. First, the kinetic term for \(\psi\) has the standard sign, ensuring that the theory is free of ghosts. The Hamiltonian density \[\begin{equation} \mathcal{H}_\psi = \frac{c^4}{8\pi G \azero}\left[(\partial_t\psi)^2 + (\nabla\psi)^2\right] \end{equation}\] is manifestly non-negative, so the elastic sector is stable both classically and quantum mechanically. Second, the constraint term couples \(\psi\) algebraically to the matter density, not dynamically. This means that \(\psi\) does not propagate independently—it is determined instantaneously by the matter distribution, similar to the Newtonian potential in non-relativistic gravity. Third, the action reduces to standard general relativity when \(\azero \to 0\) (equivalently, when \(\Hz \to 0\)), confirming that EEG is a genuine extension rather than a replacement.

Field Equations

Varying the action with respect to \(g_{\mu\nu}\) yields the modified Einstein equations: \[\begin{equation} G_{\mu\nu} + \Lambda g_{\mu\nu} = 8\pi G\left(T_{\mu\nu}^{\text{matter}} + T_{\mu\nu}^{\text{elastic}}\right), \label{eq:Einstein-EEG} \end{equation}\] where the elastic stress tensor takes the form \[\begin{equation} T_{\mu\nu}^{\text{elastic}} = \frac{c^4}{4\pi G \azero}\left(\partial_\mu\psi \partial_\nu\psi - \frac{1}{2}g_{\mu\nu}(\partial\psi)^2\right). \label{eq:elastic-stress} \end{equation}\] This has the same structure as the stress tensor for a massless scalar field, but with \(\psi\) constrained to track the matter distribution rather than evolving freely.

Varying with respect to \(\psi\) and \(\lambda\) yields the constraint equations: \[\begin{align} \nabla^2\psi &= \frac{4\pi G \azero}{c^4}\lambda(x), \\ \psi &= \frac{8\pi G}{\azero c^2}\rho_{\text{rest}}(x). \end{align}\] Combining these determines \(\lambda\) in terms of known quantities, completing the system.

Toward the MOND Limit

The Newtonian limit is where the schematic action of this section is deliberately incomplete, and it is worth being explicit about why. In the weak-field, slow-motion limit (\(g_{00} \approx -(1 + 2\Phi/c^2)\), static fields), the time–time component of Eq. \(\eqref{eq:Einstein-EEG}\) reads \[\begin{equation} \nabla^2\Phi = 4\pi G\rho + \frac{c^4}{\azero}(\nabla\psi)^2. \end{equation}\] With \(\psi\) algebraically slaved to the rest-mass density by the constraint \(\psi = 8\pi G\rho/(\azero c^2)\), the elastic source \((\nabla\psi)^2 \propto (\nabla\rho)^2\) depends on gradients of the matter distribution, not on the gravitational field \(\nabla\Phi\). It therefore does not reproduce the MOND structure \[\begin{equation} \nabla\cdot\!\left[\mu\!\left(\frac{|\nabla\Phi|}{\azero}\right)\nabla\Phi\right] = 4\pi G\rho, \end{equation}\] in which the modification is controlled by the field strength \(|\nabla\Phi|/\azero\). Recovering this structure requires promoting \(\psi\) to a genuinely dynamical field with a nonlinear kinetic functional \(\mathcal{F}\!\left((\partial\psi)^2/\azero^2\right)\) whose deep-MOND limit \(\mathcal{F}(z)\propto z^{3/2}\) yields the \(\sqrt{\azero\,g_N}\) behaviour—exactly the route taken in the Bekenstein–Milgrom AQUAL theory  and its relativistic descendants (TeVeS , the Skordis–Złośnik model ). The quadratic action written above is the linearization of such a functional about the high-acceleration regime, and is presented only to fix notation and the elastic modulus.

What EEG fixes without the nonlinear completion is the acceleration scale and the elastic modulus: \(\azero = c\Hz/(2\pi)\) from de Sitter thermodynamics (Section 4) and \(\kappa = c^4/(8\pi G\azero)\) from the strain energy of Section 6. Constructing the nonlinear functional explicitly, and then verifying ghost-freedom and the gravitational-wave speed in the resulting theory, is left to future work.

Gravitational Waves

A crucial test of any modified gravity theory is consistency with gravitational wave observations. The detection of GW170817 and its electromagnetic counterpart GRB 170817A  established that gravitational waves travel at the speed of light to within one part in \(10^{15}\). Many modified gravity theories, including some MOND extensions, predict deviations from \(c\) and are therefore ruled out.

A scalar elastic sector of the kind proposed here is, by construction, compatible with this bound: because \(\psi\) enters only through an elastic stress tensor that sources curvature but does not mix with the transverse-traceless metric kinetic term, tensor perturbations propagate on null geodesics of the background metric at speed \(c\). We note this as a design requirement that any nonlinear completion (Section 6.5) must preserve, rather than as an established property of a finished theory—unlike vector–tensor MOND extensions, where the GW speed must be checked explicitly .

Comparison with Alternative Approaches

EEG offers several advantages over both dark matter models and phenomenological MOND. Compared to dark matter, EEG requires no new particles, no fine-tuning of halo profiles, and no “cosmic conspiracy” to explain the tight Tully-Fisher relation—the scaling emerges automatically from the thermodynamic origin of \(\azero\). Compared to phenomenological MOND, EEG provides a complete relativistic framework with a microscopic foundation in quantum entanglement, natural incorporation of the cosmological constant, and compatibility with gravitational wave constraints.

The theory does have limitations. The interpolation function \(\mu(x)\) is not derived from first principles but must be chosen phenomenologically. The strong-field regime (black holes, neutron stars) has not been fully analyzed. Cosmological perturbation theory, including predictions for the CMB and large-scale structure, requires numerical implementation in a Boltzmann code. These are topics for future work.

Summary

Entanglement-Elastic Gravity sketches a candidate covariant completion of the thermodynamic MOND framework developed in earlier sections. It introduces a scalar displacement field \(\psi\) that encodes the strain in the cosmological entanglement network caused by the presence of matter, contributing an elastic stress tensor to Einstein’s equations. The proposal is designed to be ghost-free and to leave the gravitational-wave speed equal to \(c\); reproducing the modified Poisson equation in the Newtonian limit, however, requires replacing the linearized quadratic action used here with a nonlinear (AQUAL-type) kinetic functional, which we leave to future work (Section 6.5). What is fixed already, independent of that completion, is the acceleration scale \(\azero = c\Hz/(2\pi)\) and the elastic modulus \(\kappa = c^4/(8\pi G\azero)\), from the same de Sitter thermodynamics that motivated the non-relativistic framework, requiring no new free parameters.

Testable Predictions

A theoretical framework earns credibility through predictions that can be tested and potentially falsified. This section presents the observational consequences of entanglement-elastic gravity, ranging from already-verified galactic phenomenology to novel signatures that distinguish this framework from alternatives.

Parameter-Free Predictions Already Confirmed

The framework makes two central predictions with no adjustable parameters, both of which match observations.

The first is the MOND acceleration scale itself. From de Sitter thermodynamics, we predict \[\begin{equation} \azero = \frac{c\Hz}{2\pi} \approx 1.08 \times 10^{-10}~\text{m/s}^2. \end{equation}\] The observed value from galaxy rotation curve fits is \(\azero^{\text{obs}} = (1.2 \pm 0.2) \times 10^{-10}~\text{m/s}^2\) . The prediction agrees to within 10% for a quantity spanning 11 orders of magnitude below terrestrial gravity, with no free parameters.

The second is the baryonic Tully-Fisher relation. In the deep MOND regime, \(v^4 = GM\azero\), predicting a slope of exactly 4 on a log-log plot of velocity versus mass, with normalization fixed by the same \(\azero\). The observed slope is \(3.98 \pm 0.12\) , and the scatter in the relation is consistent with measurement error alone—no intrinsic scatter from varying “dark matter halo” properties. This tight correlation, which requires fine-tuning in dark matter models, emerges automatically from thermodynamic principles.

Predictions for Precision Tests

The framework predicts that deviations from GR scale as \(\azero/g_N\) in any high-acceleration environment. This means lunar laser ranging, with \(\epsilon \sim 10^{11}\), should show GR to parts in \(10^{11}\); Mercury’s perihelion precession, with \(\epsilon \sim 10^8\), should match GR to parts in \(10^8\); and binary pulsar timing, with \(\epsilon \sim 10^{13}\), should be the most stringent test of pure GR. All current precision tests are consistent with these predictions. Importantly, the framework predicts that no anomalous effects should appear in these high-\(\epsilon\) systems—any claimed detection of MOND-like effects in the inner solar system would falsify the framework.

Cosmological Evolution of the Acceleration Scale

The robust prediction of this framework is that \(\azero\) is locked to a cosmological horizon frequency: \(\azero = cH/(2\pi)\), with the same horizon rate \(H\) fixing the cosmic gravitational decoherence rate through \(\azero/(c\,\Gamma_{\rm cosmic}) = 2/g(1)\) (the horizon-clock lock; see the Falsification Criteria below). A separate, less certain question is which horizon rate enters \(H\), and hence whether \(\azero\) evolves with redshift. In the instantaneous-horizon reading, \(\azero(z) = c H(z)/(2\pi)\) grows in the past, so galaxies at high redshift would show MOND effects at smaller radii than their present-day counterparts; in the asymptotic event-horizon reading, \(\azero\) is constant. We retain the evolving case below because it would distinguish the instantaneous reading, but note that current high-redshift rotation curves favour an approximately constant \(\azero\).

At redshift \(z = 1\), the MOND scale was approximately 1.78 times larger. A galaxy with fixed baryonic mass should therefore have had an asymptotic rotation velocity about 19% higher than the same galaxy would have today. At \(z = 2\), with \(\azero\) about 3.0 times larger, many galaxies should appear nearly Newtonian. These predictions can be tested with JWST and Euclid observations of high-redshift rotation curves; a measured \(\azero(z)\propto H(z)\) would single out the instantaneous-horizon reading over both phenomenological MOND and the asymptotic-horizon reading (both of which give constant \(\azero\)).

Priority targets. The transition regime at \(z \sim 2\)–2.5 provides optimal leverage. Galaxies zC406690 (\(z = 2.196\)) and BX482 (\(z = 2.26\)) lie precisely where \(\azero(z)\) is large enough to suppress MOND effects substantially but not so large that signals vanish. Dwarf galaxies at \(z > 2\) should appear completely Newtonian, with no evidence of the flat rotation curves seen in local dwarfs. JWST observations during the 2027–2030 window can achieve 3–5\(\sigma\) sensitivity to the predicted evolution; ELT commissioning in 2030+ will provide definitive \(>10\sigma\) measurements.

CMB Consistency

At recombination (\(z \approx 1100\)), the MOND scale was \(\azero(z_{\rm rec}) \approx 2.5 \times 10^{-6}~\text{m/s}^2\) (computed as \(\azero \times \sqrt{\Omega_m(1+z)^3 + \Omega_\Lambda + \Omega_r(1+z)^4}\) with \(\Omega_r \approx 9.4 \times 10^{-5}\)), vastly larger than any acceleration relevant to CMB physics. All cosmological perturbations lie deep in the Newtonian regime, where the elastic corrections are negligible. The framework therefore predicts that CMB acoustic peaks, lensing, and polarization should match \(\Lambda\)CDM predictions identically. No modification to the standard CMB analysis is required.

One subtle prediction concerns the integrated Sachs-Wolfe (ISW) effect. Modified gravity in the late universe can alter the growth of structure, leading to ISW suppression at the 10–20% level compared to \(\Lambda\)CDM. This signature lies within cosmic variance for current datasets and is not a strong discriminator, but future cross-correlations with large-scale structure may provide a consistency check. In short, CMB observations constrain the framework at recombination—where it predicts standard physics—rather than testing the MOND regime.

Falsification Criteria

The framework makes strong predictions that can be falsified. If precision measurements establish that \(\azero \neq c\Hz/(2\pi)\) at greater than 20% discrepancy, the thermodynamic derivation would be wrong. If MOND-like effects are detected in high-\(\epsilon\) systems at levels exceeding \(\azero/g_N\)—for example, anomalous perihelion precession at parts in \(10^6\) rather than \(10^9\)—the regime separation would be violated. Most distinctively, the framework ties the MOND scale to a cosmological horizon, \(\azero = cH/(2\pi)\). Whether the relevant \(H\) is the instantaneous apparent-horizon rate \(H(z)\) (giving the evolving \(\azero(z)\propto H(z)\) used above) or the asymptotic event-horizon rate \(H_\Lambda\) (giving a constant \(\azero\)) is an open question we do not settle here; current high-redshift rotation curves in fact favour an approximately constant \(\azero\) over \(z=0\)–\(2\). A measured \(\azero(z)\) tracking \(H(z)\) would confirm the instantaneous reading; a constant \(\azero\) would select the asymptotic (event-horizon) reading. The latter does not falsify the \(\azero=cH/(2\pi)\) derivation—which holds with \(H=H_\Lambda\)—nor the rest of the framework (dark energy, decoherence, the horizon-clock lock). The horizon-clock paper shows that \(\azero\) and the cosmic gravitational decoherence rate are the same de Sitter horizon frequency, locked at \(\azero/(c\,\Gamma_{\rm cosmic})=2/g(1)\) independent of which horizon is chosen; that locked ratio, not a particular \(\azero(z)\) law, is the robust prediction.

The galaxy cluster tension represents an incomplete rather than falsified prediction. If detailed calculations of higher-order corrections, non-spherical geometry, hot gas effects, and neutrino contributions fail to close the factor 2–5 gap, this would indicate that additional physics is needed—though not necessarily that the framework is wrong at the galactic level.

Non-Equilibrium Cluster Dynamics

One mechanism addresses the cluster tension partially: thermodynamic equilibrium timescales in galaxy clusters are much longer than dynamical times. The elastic stress tensor derives from vacuum entanglement entropy, which requires thermal equilibrium to saturate at the predicted level. For a cluster with mass \(M_{\text{cl}} \sim 10^{15}~M_\odot\) and characteristic scale \(R \sim 1~\text{Mpc}\), the equilibration timescale is estimated to be \[\begin{equation} \tau_{\text{eq}} \gg \tau_{\text{dyn}} \sim \frac{R}{v_{\text{circ}}} \sim 1~\text{Gyr}, \end{equation}\] where \(v_{\text{circ}} \sim 10^3\) km/s is the cluster velocity dispersion. A rough estimate gives \(\tau_{\text{eq}} \sim 100\)–300 Gyr, though a complete derivation of the equilibration timescale from the EEG equations of motion has not been carried out and should be treated as an order-of-magnitude guide. What matters is that \(\tau_{\text{eq}} \gg \tau_{\text{dyn}}\), so clusters are far from equilibrium.

Clusters are therefore out of equilibrium: the elastic stress has not had time to fully develop since formation at \(z \sim 1\)–2. This explains why clusters show only a factor \(\sim 2\) MOND boost rather than the factor \(\sim 3.5\) expected in equilibrium. The framework predicts a testable correlation: relaxed, cool-core clusters that have had longer to equilibrate should show larger MOND boosts than recently merged or dynamically disturbed clusters. This can be tested with X-ray morphology classifications combined with weak lensing mass estimates, providing a novel signature distinct from both \(\Lambda\)CDM and phenomenological MOND.

The Bullet Cluster Challenge

The Bullet Cluster (1E 0657-558) presents the most cited evidence against MOND-like theories. In this merging system at \(z = 0.296\), weak lensing reveals that the gravitational mass is spatially offset from the dominant baryonic component (hot X-ray gas), apparently requiring invisible matter that passed through the collision while gas was shocked and left behind.

This framework addresses the Bullet Cluster through several mechanisms. First, the collision occurred approximately 150 Myr ago, vastly shorter than the entanglement equilibration timescale \(\tau_{\text{eq}} \sim 100\)–300 Gyr. The elastic stress tensor, which provides the MOND enhancement, requires equilibrium to develop fully. In a violent merger, pre-existing equilibrium is disrupted, and the system has not had time to re-equilibrate. Second, at \(z = 0.296\), the MOND scale was \(\azero(0.296) \approx 1.3 \times \azero(0)\), slightly reducing the expected MOND enhancement compared to local clusters.

Quantitatively, the Bullet Cluster requires approximately 5–6 times the baryonic mass to explain lensing observations. This framework, combined with standard contributions (massive neutrinos at 5–15%, hydrostatic bias at 10–30%, triaxial geometry at 10–20%), can account for approximately 3–4 times the baryonic mass. A factor of \(\sim 2\) residual discrepancy remains.

We assess this honestly: the Bullet Cluster challenges but does not definitively falsify the framework. The discrepancy is smaller than often claimed (factor \(\sim 2\), not \(\sim 6\)), and specific predictions distinguish this framework from both \(\Lambda\)CDM and standard MOND. In particular, we predict that the mass discrepancy should decrease with cluster relaxation state—the Bullet Cluster, as a recent violent merger, should show larger discrepancy than relaxed systems. This correlation is testable and would distinguish the non-equilibrium entanglement mechanism from particle dark matter.

Comparison with Alternatives

Several approaches have been proposed to explain MOND phenomenology. Standard MOND treats \(\azero\) as an empirical constant fitted to data, with no explanation for its value or cosmic evolution. Tensor-vector-scalar (TeVeS) theories  provide a relativistic completion but require multiple fields and free parameters, and typically assume constant \(\azero\); more recent relativistic MOND theories  reproduce the CMB and matter power spectra at the cost of additional fields. Verlinde’s emergent gravity  derives an acceleration scale from holographic principles but does not specify its redshift evolution.

This framework is unique in predicting \(\azero = cH/(2\pi)\) from first principles, with zero free parameters and an explicit tie to the cosmological horizon rate. The robust prediction is the horizon-clock lock \(\azero/(c\,\Gamma_{\rm cosmic}) = 2/g(1)\), which holds in either horizon reading. A redshift evolution \(\azero(z) \propto H(z)\) would additionally select the instantaneous-horizon reading over phenomenological MOND, though current high-redshift data favour an approximately constant \(\azero\); next-generation surveys can settle which reading applies.

Implications for Dark Matter Searches

If this framework correctly explains galactic dynamics, the “missing mass” in galaxies is a gravitational effect arising from cosmological entropy, not a particle physics phenomenon. Direct detection experiments would find no WIMP dark matter, and collider searches would produce no dark matter candidates, because none exist at the relevant scales. This does not preclude the existence of dark matter particles—sterile neutrinos, axions, or other relics may well exist—but suggests they are not responsible for galaxy rotation curves. The cluster tension may require residual dark matter at the 1–10% level, possibly from massive neutrinos with \(\sum m_\nu \sim 0.1\)–0.5 eV.

Conclusions

This paper has developed a theoretical framework for MOND phenomenology based on the thermodynamics of de Sitter space. The central physical picture is that the cosmological horizon acts as a thermal reservoir, filling the observable universe with volume-law entropy that competes with the standard area-law contribution to gravitational dynamics. This framework rests on physically motivated assumptions—particularly that local entropy equilibrium fails at low accelerations—that are not rigorously derived from first principles.

When matter is placed in this thermal bath, it displaces entropy from the background, creating a thermodynamic strain that contributes to the gravitational field. At high accelerations, where the local curvature can easily accommodate this displacement, standard general relativity applies. At low accelerations, where the entropy deficit spreads non-locally through the de Sitter volume, modified dynamics emerge. The transition occurs at the acceleration scale \(\azero = c\Hz/(2\pi) \approx 1.08 \times 10^{-10}~\text{m/s}^2\), matching the observed MOND scale with no adjustable parameters.

The framework makes concrete, falsifiable predictions. It explains why no MOND effects appear in the solar system (where \(\epsilon \gg 1\)) while producing flat rotation curves and the Tully-Fisher relation in galaxies (where \(\epsilon \lesssim 1\)). It predicts that the MOND scale evolves with cosmic time as \(\azero(z) \propto H(z)\), so high-redshift galaxies should show MOND effects at smaller radii. It preserves the success of standard cosmology at the CMB epoch, where all relevant scales lie deep in the Newtonian regime.

Relation to the Broader Framework

This paper is part of a series developing quantum-geometric correspondence—the thesis that quantum mechanics and general relativity are complementary descriptions of the same underlying reality. The canonical core paper  presents the complete axiomatic framework and develops gravitational decoherence from entanglement equilibrium arguments; the holographic dark-energy paper shows how holographic dark energy is consistent with thermodynamic principles. The present work extends these ideas to de Sitter backgrounds, where the cosmological horizon introduces new thermodynamic effects absent in asymptotically flat spacetime.

The volume entropy is a background-dependent phenomenon that appears when \(\Lambda > 0\), not a modification of the universal axioms. In the limit \(\Lambda \to 0\), the de Sitter temperature vanishes, the volume entropy disappears, and standard general relativity is recovered exactly. The MOND-like modifications are thus a cosmological effect, arising from the finite size of the observable universe rather than from new fundamental physics at short distances.

Open Problems

Several important questions remain for future work. The interpolation function \(\mu(x)\) that governs the transition between Newtonian and deep MOND regimes is introduced phenomenologically; deriving it rigorously from the underlying entropy dynamics would strengthen the framework considerably. The galaxy cluster tension—a factor 2–5 mass discrepancy even after MOND corrections—likely involves higher-order terms in the \(\epsilon\)-expansion, non-spherical geometry, and hot gas thermodynamics, but detailed calculations are needed.

The framework as developed here assumes quasi-static configurations. Extending it to dynamical situations—galaxy mergers, structure formation, time-dependent gravitational fields—requires a time-dependent entropy analysis that has not been attempted. The volume entropy is fundamentally a semiclassical concept; a full quantum gravity treatment might confirm the volume-law from explicit entanglement calculations, reveal the microscopic origin of \(s_\Lambda\), and suggest quantum corrections that could resolve residual tensions.

Experimental Outlook

The next decade offers unprecedented opportunities to test this framework. JWST and Euclid can measure rotation curves at \(z \sim 1\)–3, directly testing whether \(\azero(z) \propto H(z)\) as predicted. Gaia DR4 and DR5 will provide precision dynamics for wide binary stars and the outer Milky Way, probing the transition regime. DESI and future spectroscopic surveys can measure the growth rate of structure, testing whether the predicted suppression \(\Delta f\sigma_8 \approx -0.03\) at \(z = 0.8\) is observed. Gravitational wave observations with next-generation detectors can probe strong-field predictions of the covariant extension.

Final Remarks

The connection between MOND phenomenology and de Sitter thermodynamics provides a compelling framework for addressing a 40-year-old puzzle: why does the MOND scale coincide with \(c\Hz\)? Within this framework, both arise from the same physics—the thermodynamic properties of the cosmological horizon. The acceleration \(\azero = c\Hz/(2\pi)\) emerges from the framework with no free parameters, connecting cosmology to galaxy dynamics through the universal language of entropy. Whether this picture is correct remains to be determined by observation.

This framework does not eliminate dark matter as a possibility, but suggests that whatever “dark” component exists is not responsible for the regularities observed in galactic dynamics. The tight baryonic Tully-Fisher relation, the universal acceleration scale, the flat rotation curves extending far beyond visible matter—all are consistent with this thermodynamic framework, requiring no new particles and no free parameters beyond those already present in cosmology. Whether this picture survives confrontation with data remains to be seen, but the framework is precise enough to be falsified and predictive enough to be tested.

Appendices

Derivation Details

Surface Gravity of de Sitter Horizon

From the metric: \[\begin{equation} ds^2 = -f(r)c^2 dt^2 + f(r)^{-1}dr^2 + r^2 d\Omega^2 \end{equation}\] with \(f(r) = 1 - r^2/\RH^2\).

The Killing vector is \(k^\mu = (1, 0, 0, 0)\). The acceleration of a static observer at radius \(r\): \[\begin{equation} a^\mu = k^\nu\nabla_\nu k^\mu = \frac{c^2 r/\RH^2}{\sqrt{f(r)}}\hat{r} \end{equation}\]

The surface gravity (denoted \(\kappa_{\rm dS}\) to distinguish it from the elastic modulus \(\kappa = c^4/(8\pi G\azero)\) of Section 6): \[\begin{equation} \kappa_{\rm dS} = \lim_{r \to \RH} \sqrt{f(r)} \cdot a = \frac{c^2}{\RH} = c\Hz \end{equation}\]

Volume Entropy from Partition Function

The partition function for a scalar field in de Sitter: \[\begin{equation} \ln Z = -\sum_n \ln(1 - e^{-\beta\omega_n}) \end{equation}\]

In the continuum limit with density of states \(D(\omega)\): \[\begin{equation} \ln Z = -\int_0^\infty d\omega\, D(\omega)\ln(1 - e^{-\beta\omega}) \end{equation}\]

For de Sitter, \(D(\omega) \approx \omega^2 V\) for low frequencies, giving: \[\begin{equation} \ln Z \approx \frac{\pi^2}{90}\frac{V}{(\beta\hbar c)^3} \end{equation}\]

The entropy: \[\begin{equation} S = \beta^2\frac{\partial}{\partial\beta}(\beta^{-1}\ln Z) = \frac{2\pi^2}{45}g_*k_B\left(\frac{k_BT}{\hbar c}\right)^3 V \end{equation}\]

For \(g_* \sim O(1)\) and \(T = \TdS\), this gives the volume entropy \(S_{\text{vol}}= s_\Lambda V\).

Modified Poisson Equation from \(\delta\Sgen = 0\)

Starting from \(\delta\Sgen = 0\) with volume term: \[\begin{equation} \frac{\delta A}{4G\hbar} + \delta\Sbulk + s_\Lambda\delta V = 0 \end{equation}\]

Using:

• \(\delta A/(4G\hbar) \approx -\ell^{d+1}G_{00}/(G\hbar)\) (Raychaudhuri)

• \(\delta\Sbulk \approx \ell^{d+1}T_{00}/T\) (first law)

• \(\delta V = \ell^3\) for perturbation scale \(\ell\)

For a localized mass, the equilibrium condition in the low-acceleration limit becomes: \[\begin{equation} \nabla^2\Phi = 4\pi G\rho + \lambda_{\text{MOND}}\nabla \cdot \left(\frac{\nabla\Phi}{|\nabla\Phi|}\sqrt{|\nabla\Phi|}\right) \end{equation}\] where \(\lambda_{\text{MOND}} = \sqrt{\azero}\).

This can be rewritten in the standard MOND form: \[\begin{equation} \nabla \cdot \left[\mu\left(\frac{|\nabla\Phi|}{\azero}\right)\nabla\Phi\right] = 4\pi G\rho \end{equation}\]

Response Kernel from Modular Hamiltonian

For the vacuum state restricted to a causal diamond \(D\) in de Sitter, the modular Hamiltonian is: \[\begin{equation} K = 2\pi\int_D \frac{(\ell^2 - |\mathbf{x}|^2)}{2\ell}T_{00}(\mathbf{x})\,d^3x + K_{\text{horizon}} \end{equation}\]

The response kernel arises from the fluctuation-dissipation relation: \[\begin{equation} \chi(x, x') = \frac{1}{\TdS}\langle\delta K(x)\delta\rho_M(x')\rangle \end{equation}\]

From the de Sitter Green’s function \(G(r) \propto e^{-r/\RH}/r\): \[\begin{equation} \chi(r) = \chi_0 \cdot \frac{\RH}{r} \cdot e^{-r/\RH} \end{equation}\]

The normalization \(\chi_0\) is fixed by: \[\begin{equation} \int d^3x\,\chi(x) = \frac{2\pi k_Bc^2}{\hbar\Hz} \end{equation}\]

Numerical Constants and Values

Fundamental Constants

Cosmological Parameters

Derived Quantities

System-Specific Values

Comparison with Observations

Agreement is within 10% with no adjustable parameters.

← Back to Research