We show that Modified Newtonian Dynamics (MOND) phenomenology emerges naturally from the thermodynamics of de Sitter space within the framework of quantum-geometric duality. The de Sitter cosmological horizon establishes a thermal bath at the Gibbons-Hawking temperature \(T_{\text{dS}}= \hbar H_0/(2\pi k_B)\), which generates extensive (volume-law) entropy in addition to the standard area-law contribution. When local entanglement equilibrium—the principle that generalized entropy is stationary under allowed variations—fails at low accelerations, the volume entropy becomes dynamically relevant, producing MOND-like corrections to Newtonian gravity.

We demonstrate that the characteristic acceleration scale emerges as \(a_0= cH_0/(2\pi) \approx 1.1 \times 10^{-10}~\text{m/s}^2\), matching observed MOND values within 10% with no adjustable parameters. We establish clear criteria for regime separation: GR remains exact when \(\epsilon = g_N/a_0\gg 1\) (solar system, pulsars), while MOND corrections become significant when \(\epsilon \lesssim 1\) (galaxy outskirts). The volume entropy saturates the Bousso covariant entropy bound, resolving any apparent conflict with holography.

We present Entanglement-Elastic Gravity (EEG), a covariant field theory that extends this framework to relativistic regimes. The EEG action introduces an elastic displacement field \(\psi\) encoding entanglement strain, yielding modified Einstein equations with an elastic stress tensor. The theory is ghost-free, preserves gravitational wave speed \(c\), and reduces to the modified Poisson equation in the Newtonian limit. Predictions: Flat rotation curves, Tully-Fisher relation \(v^4 = GMa_0\), weak lensing slip \(|\Phi - \Psi|/|\Phi| \approx 15\%\) beyond 50 kpc, and growth rate suppression \(\Delta f\sigma_8 \approx -0.03\) testable with Euclid and DESI.

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 remarkably 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 \(a_0\), the actual acceleration experienced by a test particle transitions to \(g \approx \sqrt{g_N \cdot a_0}\). This simple modification, with \(a_0\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 \(a_0\) must be fitted to data rather than predicted from first principles.

The purpose of this paper is to provide that theoretical foundation. We show that the MOND phenomenology emerges naturally from the thermodynamics of de Sitter space, with the acceleration scale \(a_0= cH_0/(2\pi)\) predicted directly from the Hubble constant. This is not merely a numerical coincidence: it reflects a deep connection between cosmological horizons and the information-theoretic structure of gravity.

Our starting point is Jacobson’s remarkable discovery  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/H_0\) 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 \(T_{\text{dS}}= \hbar H_0/(2\pi k_B)\) and carries entropy proportional to its area.

The key insight of this paper is that the de Sitter horizon fundamentally 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 \(T_{\text{dS}}\) 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/T_{\text{dS}}\) from the surrounding volume. If the available volume entropy cannot accommodate this displacement, local thermodynamic equilibrium fails, and the gravitational dynamics must change.

This failure of equilibrium occurs precisely when the Newtonian acceleration falls below \(a_0= cH_0/(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 predictions, including the weak lensing slip parameter \(|\Phi - \Psi|/|\Phi| \approx 15\%\) beyond 50 kpc, growth rate suppression \(\Delta f\sigma_8 \approx -0.03\) at \(z = 0.8\), and anomalous acceleration in wide binary stars detectable with Gaia DR4/DR5. 8 summarizes our findings and discusses open problems.

Series context. This is Paper IV of the Quantum-Geometric Duality series. Paper A develops gravitational decoherence from the same entropy principles; Paper B derives holographic dark energy; Paper C presents the complete axiomatic framework; Paper H establishes information-theoretic bounds. 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}{R_H^2}\right)c^2 dt^2 + \frac{dr^2}{1 - r^2/R_H^2} + r^2 d\Omega^2, \label{eq:desitter-metric} \end{equation}\] where the de Sitter radius is \[\begin{equation} R_H= \sqrt{\frac{3}{\Lambda}} = \frac{c}{H_0}. \label{eq:horizon-radius} \end{equation}\]

The significance of \(R_H\) becomes apparent from the metric: the coefficient \(g_{tt}\) vanishes at \(r = R_H\), 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 > R_H\) can never reach \(r = 0\), just as signals sent into a black hole can never escape. For our universe with \(H_0\approx 70~\text{km/s/Mpc}\), the horizon lies at \(R_H\approx 14\) billion light-years, defining the edge of the observable universe.

The region \(r < R_H\) 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)R_H^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

Perhaps the most profound discovery in 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} T_{\text{dS}}= \frac{\hbar H_0}{2\pi k_B} \approx 2.8 \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\) that measures the “strength” of the horizon is defined geometrically by the Killing equation \(\nabla_\mu(k_\nu k^\nu) = -2\kappa k_\mu\). For the de Sitter metric, direct calculation gives \(\kappa = cH_0\), and the temperature follows from the universal relation \(T = \hbar\kappa/(2\pi ck_B)\).

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 = R_H\) unless \(\tau\) is periodic with a specific period. Requiring regularity uniquely fixes this period to be \(\beta = 2\pi/H_0\), which corresponds to temperature \(T = \hbar/(2\pi k_B\beta) = \hbar H_0/(2\pi k_B)\).

The third method appeals to the Unruh effect. A static observer at any radius \(r < R_H\) must accelerate to resist the cosmic expansion. As \(r \to R_H\), 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 \(T_{\text{dS}}\) 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\ell_P^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} = T_{\text{dS}}\, 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 \(T_{\text{dS}}\), 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_BT_{\text{dS}}= \frac{\hbar H_0}{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, something remarkable happens. Using natural gravitational units defined by the Planck mass \(m_P= \sqrt{\hbar c/G}\) and Planck length \(\ell_P= \sqrt{\hbar G/c^3}\), we find \[\begin{equation} a_0= \frac{E_{dS}}{m_P\ell_P} = \frac{cH_0}{2\pi} \approx 1.1 \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 \(a_0\sim cH_0\) 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.

The crucial insight is that 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 H_0} |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 H_0}|n\rangle_{\text{in}}\langle n|. \end{equation}\] This is a Boltzmann distribution at precisely the Gibbons-Hawking temperature \(T_{\text{dS}}= \hbar H_0/(2\pi k_B)\). The de Sitter horizon is not just radiating—it is maintaining 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^2H_0^2}{8\pi G}, \end{equation}\] where the second equality uses \(\Lambda = 3H_0^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}{T_{\text{dS}}} = \frac{3c^2H_0^2}{8\pi G} \cdot \frac{2\pi k_B}{\hbar H_0} = \frac{3c^2H_0k_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)}{4G\hbar} + s_\Lambda V(r) = \frac{\pi r^2 c^3}{G\hbar} + \frac{\pi c^2H_0k_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/H_0= R_H\)—precisely at the cosmological horizon.

Consistency with Holography

The volume-law entropy raises an apparent paradox: for sufficiently large regions, \(S_{\text{vol}}> S_{\text{area}}\), seemingly violating the holographic principle that entropy cannot exceed \(A/(4\ell_P^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\ell_P^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, something remarkable happens 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^2H_0k_B}{4G\hbar} \cdot \frac{4\pi c^3}{3H_0^3} = \frac{\pi c^5k_B}{G\hbar H_0^2}, \end{equation}\] while the horizon entropy is \[\begin{equation} S_H = \frac{A_H}{4\ell_P^2} = \frac{4\pi R_H^2 \cdot c^3}{4G\hbar} = \frac{\pi c^5}{G\hbar H_0^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 derive the central result of this paper: the MOND acceleration scale emerges as a consequence of entropy equilibrium, with no free 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}{T_{\text{dS}}} = \frac{2\pi k_BMc^2}{\hbar H_0}. \label{eq:entropy-displacement} \end{equation}\] This is not a small effect. The de Sitter temperature \(T_{\text{dS}}\approx 2.7 \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^2H_0k_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 H_0k_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 H_0} = \frac{\pi c^2 H_0k_Br_c^3}{G\hbar}. \end{equation}\] The factors of \(\pi\), \(k_B\), \(c^2\), and \(\hbar\) cancel, leaving a remarkably simple result: \[\begin{equation} r_c= \left(\frac{2GM}{H_0^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 = H_0r\) 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 are striking. For the Sun (\(M = M_\odot\)), the critical radius is \(r_c\approx 700\) light-years—far beyond the solar system but well within the Milky Way. For a typical galaxy (\(M = 10^{11} M_\odot\)), the critical radius is approximately 100 kpc, comparable to the observed extent of flat rotation curves. For a galaxy cluster (\(M = 10^{14} M_\odot\)), the critical radius extends to several Mpc. These scales are precisely where the transition from Newtonian to MOND-like behavior is observed.

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/H_0^2\right)^{2/3}} = \frac{(GM)^{1/3}H_0^{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 \(T_{\text{dS}}= \hbar H_0/(2\pi k_B)\), which sets a natural energy scale \(E_{\text{dS}} = k_BT_{\text{dS}}= \hbar H_0/(2\pi)\). To convert this energy to an acceleration, we use the Planck mass \(m_P= \sqrt{\hbar c/G}\) and Planck length \(\ell_P= \sqrt{\hbar G/c^3}\) as the natural units of mass and length: \[\begin{equation} a_0= \frac{E_{\text{dS}}}{m_P\ell_P} = \frac{\hbar H_0/(2\pi)}{\sqrt{\hbar c/G} \cdot \sqrt{\hbar G/c^3}} = \frac{cH_0}{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.

Numerically, using \(c = 3 \times 10^8\) m/s and \(H_0= 2.3 \times 10^{-18}\) s\(^{-1}\) (corresponding to \(H_0 = 70\) km/s/Mpc): \[\begin{equation} a_0= \frac{(3 \times 10^8)(2.3 \times 10^{-18})}{2\pi} \approx 1.1 \times 10^{-10}~\text{m/s}^2. \end{equation}\] The observed MOND scale from galaxy rotation curve fits is \(a_0^{\text{obs}} \approx (1.2 \pm 0.2) \times 10^{-10}~\text{m/s}^2\) . The agreement to within 10% is remarkable 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< a_0\). The entropy equilibrium condition \(\delta S_{\text{gen}}= 0\) must include the volume entropy contribution when local equilibrium fails.

The physical picture is as follows. At high accelerations (\(g_N\gg a_0\)), 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 a_0\)), 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|}{a_0}\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.

The effective gravitational field \(g_{\text{eff}}= -\nabla\Phi\) takes different forms in the two regimes. In the Newtonian regime where \(g_N\gg a_0\): \[\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 a_0\): \[\begin{equation} g_{\text{eff}}= \sqrt{g_N\cdot a_0} = \sqrt{\frac{GMa_0}{r^2}} = \frac{\sqrt{GMa_0}}{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{GMa_0}{r^2}}. \end{equation}\] Squaring both sides and rearranging: \[\begin{equation} \frac{v^4}{r^2} = \frac{GMa_0}{r^2}. \end{equation}\] The factors of \(r^2\) cancel, yielding the mass-independent relation \[\begin{equation} v^4 = GMa_0= GM \cdot \frac{cH_0}{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.1 \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)}{a_0} = \frac{GM}{r^2a_0} = \frac{2\pi GM}{r^2 cH_0}. \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 S_{\text{gen}}= 0\) involves three terms: the area entropy \(\delta S_{\text{area}}\), the bulk matter entropy \(\delta S_{\text{bulk}}\), and the de Sitter volume entropy \(\delta S_{\text{vol}}\). For a spherically symmetric perturbation at radius \(r\), dimensional analysis gives \[\begin{align} \delta S_{\text{area}}&\sim \frac{r}{G\hbar}\,\delta r, \\ \delta S_{\text{bulk}}&\sim \frac{c^2 r^2 \epsilon}{G\hbar}\,\delta r, \\ \delta S_{\text{vol}}&\sim \frac{c^2H_0r^2}{G\hbar}\,\delta r. \end{align}\] The ratio of volume to bulk contributions scales as \(\delta S_{\text{vol}}/\delta S_{\text{bulk}}\sim 1/\epsilon\). This motivates organizing the equilibrium condition as an expansion in inverse powers of \(\epsilon\): \[\begin{equation} \delta S_{\text{gen}}= \underbrace{\delta S_{\text{area}}+ \delta S_{\text{bulk}}}_{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 S_{\text{area}}+ \delta S_{\text{bulk}}= 0\), which is precisely Jacobson’s derivation of Einstein’s equations. General relativity emerges exactly, with corrections suppressed by \(a_0/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 S_{\text{area}}+ \delta S_{\text{bulk}}+ \delta S_{\text{vol}}= 0\) yields interpolated dynamics described by \[\begin{equation} g_{\text{eff}}\cdot \mu\left(\frac{g_{\text{eff}}}{a_0}\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_Na_0}\). 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{2GM_\odot}{H_0^2}\right)^{1/3} \approx 700~\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}M_\odot\), the critical radius is \(r_c\approx 100\) kpc, beyond the optical radius but within the region probed by rotation curves and satellite galaxy dynamics. This is precisely where MOND phenomenology is observed.

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 \(a_0= cH_0/(2\pi)\).

Cosmological Evolution of the MOND Scale

A crucial prediction distinguishes this framework from phenomenological MOND: the acceleration scale evolves with cosmic time. Since \(a_0= cH/(2\pi)\) and the Hubble parameter changes as the universe expands, the MOND scale was larger in the past. Specifically, \[\begin{equation} a_0(z) = \frac{cH(z)}{2\pi} = a_0(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.7 times larger than today; at \(z = 3\), about 3.7 times larger; at recombination (\(z \sim 1100\)), roughly 5700 times larger.

This evolution has profound consequences. At higher redshift, with a larger \(a_0\), 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, \(a_0(z=1100) \approx 6 \times 10^{-7}~\text{m/s}^2\)—much larger than the present-day value. The critical radius for a proto-cluster shrinks to roughly 3 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 \(a_0(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 \(a_0= cH_0/(2\pi)\), a complete theory requires fully covariant field equations that reduce to the modified Poisson equation in the appropriate limit. This section develops Entanglement-Elastic Gravity (EEG), a relativistic extension that achieves this goal while maintaining consistency with gravitational wave observations and avoiding the pathologies that plague many modified gravity theories.

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/H_0\) creates a thermal bath that fills the interior with long-range entanglement at density \(s_\Lambda= 3c^2H_0k_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/T_{\text{dS}}\) from the thermal background (creating global strain). At high accelerations where \(g \gg a_0\), the area-law contribution dominates, and standard general relativity applies. At low accelerations where \(g \lesssim a_0\), the volume-law strain becomes significant, producing the MOND-like modifications we observe in galactic dynamics. The critical scale \(a_0= cH_0/(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 H_0}. \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 a_0}. \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 a_0}(\partial\psi)^2 - \lambda\left(\psi - \frac{8\pi G}{a_0c^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 Ga_0)\) 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 a_0}\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 \(a_0\to 0\) (equivalently, when \(H_0\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 a_0}\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 a_0}{c^4}\lambda(x), \\ \psi &= \frac{8\pi G}{a_0c^2}\rho_{\text{rest}}(x). \end{align}\] Combining these determines \(\lambda\) in terms of known quantities, completing the system.

Recovery of MOND in the Newtonian Limit

To verify that EEG reduces to the modified Poisson equation derived earlier, consider the weak-field, slow-motion limit where \(g_{00} \approx -(1 + 2\Phi/c^2)\) and time derivatives can be neglected. The time-time component of Eq. \(\eqref{eq:Einstein-EEG}\) becomes \[\begin{equation} \nabla^2\Phi = 4\pi G\rho + \frac{c^4}{a_0}(\nabla\psi)^2. \end{equation}\] For a spherically symmetric, static matter distribution, the constraint \(\psi = 8\pi G\rho/(a_0c^2)\) can be solved to give \(\nabla\psi \propto \nabla\Phi\). Substituting this relation and simplifying yields \[\begin{equation} \nabla^2\Phi + \frac{1}{a_0}\nabla \cdot (\nabla\Phi \cdot \nabla\Phi) = 4\pi G\rho, \end{equation}\] which is precisely the modified Poisson equation of Section 4. The covariant theory therefore reproduces the non-relativistic phenomenology by construction.

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.

EEG passes this test. The elastic sector modifies only the scalar part of the gravitational field—it contributes to the energy density and pressure that source curvature but does not alter the propagation of tensor perturbations. Gravitational waves in EEG travel on null geodesics of the background metric at exactly speed \(c\), in agreement with observations.

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 \(a_0\). 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 provides a covariant completion of the thermodynamic MOND framework developed in earlier sections. The theory introduces a scalar displacement field \(\psi\) that encodes the strain in the cosmological entanglement network caused by the presence of matter. This field contributes an elastic stress tensor to Einstein’s equations, producing MOND-like modifications at low accelerations while preserving general relativity at high accelerations. The theory is ghost-free, predicts gravitational wave speed equal to \(c\), and reduces to the modified Poisson equation in the Newtonian limit. The acceleration scale \(a_0= cH_0/(2\pi)\) emerges from the same thermodynamic arguments 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} a_0= \frac{cH_0}{2\pi} \approx 1.08 \times 10^{-10}~\text{m/s}^2. \end{equation}\] The observed value from galaxy rotation curve fits is \(a_0^{\text{obs}} = (1.2 \pm 0.1) \times 10^{-10}~\text{m/s}^2\). Agreement to within 10% for a quantity spanning 11 orders of magnitude below terrestrial gravity, with no free parameters, is remarkable.

The second is the baryonic Tully-Fisher relation. In the deep MOND regime, \(v^4 = GMa_0\), predicting a slope of exactly 4 on a log-log plot of velocity versus mass, with normalization fixed by the same \(a_0\). 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 \(a_0/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.

Novel Predictions: Cosmological Evolution

The most distinctive prediction is the redshift evolution of the MOND scale. Since \(a_0(z) = c H(z)/(2\pi)\) and the Hubble parameter was larger in the past, galaxies at high redshift should show MOND effects at smaller radii than their present-day counterparts.

At redshift \(z = 1\), the MOND scale was approximately 1.7 times larger. A galaxy with fixed baryonic mass should therefore have had an asymptotic rotation velocity about 14% higher than the same galaxy would have today. At \(z = 2\), with \(a_0\) about 2.8 times larger, many galaxies should appear nearly Newtonian. These predictions can be tested with JWST and Euclid observations of high-redshift rotation curves, providing a clear discriminator between this framework and phenomenological MOND (which assumes constant \(a_0\)).

The CMB provides another consistency check. At recombination, \(a_0\approx 6 \times 10^{-7}~\text{m/s}^2\), so large that all CMB-relevant scales lie deep in the Newtonian regime. The acoustic peaks should match \(\Lambda\)CDM predictions exactly, which they do. No modification to the standard CMB analysis is required.

Falsification Criteria

The framework makes strong predictions that can be falsified. If precision measurements establish that \(a_0\neq cH_0/(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 \(a_0/g_N\)—for example, anomalous perihelion precession at parts in \(10^6\) rather than \(10^9\)—the regime separation would be violated. Most distinctively, if high-redshift galaxies show the same MOND scale as local galaxies rather than the predicted \(a_0(z) \propto H(z)\), the entire framework would be falsified.

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.

Comparison with Alternatives

Several approaches have been proposed to explain MOND phenomenology. Standard MOND treats \(a_0\) 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 \(a_0\). Verlinde’s emergent gravity derives an acceleration scale from holographic principles but does not specify its redshift evolution.

This framework is unique in predicting \(a_0= cH/(2\pi)\) from first principles, with zero free parameters and explicit dependence on \(H(z)\). The redshift evolution is the key discriminator: phenomenological MOND predicts constant \(a_0\), while we predict \(a_0(z) \propto H(z)\). Next-generation surveys can distinguish these scenarios.

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 foundation for MOND phenomenology based on the thermodynamics of de Sitter space. The central insight 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.

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 \(a_0= cH_0/(2\pi) \approx 1.1 \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 \(a_0(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 duality—the thesis that quantum mechanics and general relativity are dual descriptions of the same underlying reality. Paper A derives gravitational decoherence from entanglement equilibrium; Paper B shows how holographic dark energy emerges; Paper C presents the complete axiomatic framework. 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 \(a_0(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 emergence of MOND phenomenology from de Sitter thermodynamics provides a compelling answer to a 40-year-old puzzle: why does the MOND scale coincide with \(cH_0\)? The answer is that both arise from the same physics—the thermodynamic properties of the cosmological horizon. The acceleration \(a_0= cH_0/(2\pi)\) is not a coincidence but a prediction, connecting cosmology to galaxy dynamics through the universal language of entropy.

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 emerge naturally from thermodynamic principles, with no new particles and no free parameters. 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/R_H^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/R_H^2}{\sqrt{f(r)}}\hat{r} \end{equation}\]

The surface gravity: \[\begin{equation} \kappa = \lim_{r \to R_H} \sqrt{f(r)} \cdot a = \frac{c^2}{R_H} = cH_0 \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 = T_{\text{dS}}\), this gives the volume entropy \(S_{\text{vol}}= s_\Lambda V\).

Modified Poisson Equation from \(\delta S_{\text{gen}}= 0\)

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

Using:

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

• \(\delta S_{\text{bulk}}\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{a_0}\).

This can be rewritten in the standard MOND form: \[\begin{equation} \nabla \cdot \left[\mu\left(\frac{|\nabla\Phi|}{a_0}\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}{T_{\text{dS}}}\langle\delta K(x)\delta\rho_M(x')\rangle \end{equation}\]

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

The normalization \(\chi_0\) is fixed by: \[\begin{equation} \int d^3x\,\chi(x) = \frac{2\pi k_Bc^2}{\hbar H_0} \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