By Apiary Staff
Introduction
For more than a century, Einstein’s General Relativity (GR) has been the cornerstone of modern cosmology. Its predictions—from the bending of starlight around the Sun to the expansion of the universe—have been confirmed countless times, and the theory underpins the Standard Model of Cosmology (ΛCDM). Yet, when we peer deeper into the cosmos, a persistent mismatch emerges: the motions of stars in galaxies, the growth of large‑scale structure, and the accelerated expansion of space all demand the existence of unseen components—dark matter and dark energy—that together account for more than 95 % of the universe’s energy budget.
Despite dozens of indirect detections, neither dark matter particles nor a fundamental dark energy field have been observed in the laboratory. This gap has motivated a vibrant community of physicists to ask a bold question: What if gravity itself behaves differently on cosmic scales? Modified or alternative gravity theories—such as Modified Newtonian Dynamics (MOND), Tensor‑Vector‑Scalar (TeVeS) gravity, emergent gravity, and massive‑gravity frameworks—attempt to rewrite the rules of gravitation, often eliminating the need for dark components or reshaping their role.
These ideas are not abstract exercises; they influence how we design astronomical surveys, interpret the data from gravitational‑wave detectors, and even inspire algorithms that govern self‑organizing AI agents. Moreover, the same principles that help a honeybee orient itself using the Earth's gravity can illuminate how alternative gravities affect the formation of structures that bees rely upon—flower fields, forests, and climate patterns. In this pillar article we dive into the most developed alternative gravity theories, examine the concrete evidence that supports or challenges them, and explore the broader implications for cosmology, technology, and conservation.
1. The Standard Model of Cosmology and Its Tensions
The ΛCDM model (Lambda‑Cold Dark Matter) rests on three pillars:
- General Relativity as the law of gravitation.
- Cold dark matter (CDM) that clusters on all scales, providing the gravitational scaffolding for galaxies.
- A cosmological constant (Λ) that drives the observed accelerated expansion.
When ΛCDM is fed with the exquisite measurements from the Planck satellite (2018 release), it yields a Hubble constant H₀ ≈ 67.4 km s⁻¹ Mpc⁻¹ and a matter density Ω_m ≈ 0.315. These numbers reproduce the angular power spectrum of the Cosmic Microwave Background (CMB) to sub‑percent precision.
However, a set of independent observations now points to tension:
| Observation | Inferred H₀ | Method | Discrepancy with Planck |
|---|---|---|---|
| Cepheid‑based distance ladder (Riess 2019) | 73.2 ± 1.3 km s⁻¹ Mpc⁻¹ | Type Ia supernovae calibrated with Cepheids | +5.8 km s⁻¹ Mpc⁻¹ (≈ 8 %) |
| Strong lensing time delays (H0LiCOW) | 73.3 ± 1.8 km s⁻¹ Mpc⁻¹ | Lensed quasars | +5.9 km s⁻¹ Mpc⁻¹ |
| Tip of the Red Giant Branch (TRGB) | 69.8 ± 1.9 km s⁻¹ Mpc⁻¹ | Stellar evolution | +2.4 km s⁻¹ Mpc⁻¹ |
These “Hubble tension” results exceed the combined statistical uncertainties, suggesting either unaccounted systematic errors or new physics beyond ΛCDM.
Beyond the Hubble constant, other discrepancies have accumulated:
- Missing satellites problem – ΛCDM predicts hundreds of dwarf galaxies around the Milky Way, yet only ~50 are observed.
- Core‑cusp problem – Simulations produce steep density cusps in dwarf galaxy centers, while observations favor shallow cores.
- Large‑scale anisotropies – The CMB exhibits a low‑ℓ quadrupole and an alignment of multipoles that are statistically unlikely under pure ΛCDM.
These tensions have spurred a renaissance of alternative gravity models, each attempting to reconcile the data without invoking exotic dark components—or at least to reduce their required abundance.
2. MOND: The Genesis of Modified Newtonian Dynamics
2.1 The Core Idea
In 1983, Mordehai Milgrom proposed Modified Newtonian Dynamics (MOND) as a phenomenological fix to the flat rotation curves of spiral galaxies. In Newtonian gravity, the circular velocity v at radius r satisfies
\[ v^2 = \frac{G M(r)}{r}, \]
where M(r) is the enclosed mass. Observationally, many galaxies show v ≈ constant (≈ 200 km s⁻¹) well beyond the luminous edge, implying \(M(r) \propto r\), which is impossible without invisible mass.
MOND replaces Newton’s second law with
\[ \mu\!\left(\frac{a}{a_0}\right) a = a_{\rm N}, \]
where a is the true acceleration, a_N the Newtonian acceleration, μ(x) an interpolating function with μ → 1 for \(x \gg 1\) (high accelerations) and μ → x for \(x \ll 1\) (low accelerations), and a₀ is a universal constant. Empirically, a₀ ≈ 1.2 × 10⁻¹⁰ m s⁻², a value intriguingly close to \(c H_0\) (the product of the speed of light and the Hubble constant).
In the deep‑MOND regime (\(a \ll a_0\)), the relation simplifies to
\[ a = \sqrt{a_0 a_{\rm N}} \quad \Rightarrow \quad v^4 = G M a_0, \]
producing the Baryonic Tully–Fisher Relation (BTFR):
\[ M_{\rm baryon} \propto v^4, \]
which is observed across five orders of magnitude in mass with a scatter of only ~0.1 dex. The BTFR is a striking success of MOND, because ΛCDM can reproduce it only after invoking a complex interplay of feedback processes.
2.2 Empirical Successes
- Galaxy rotation curves – A 1998 sample of 75 low‑surface‑brightness galaxies matched MOND predictions with χ² per degree of freedom ≈ 1.2.
- Dwarf spheroidal galaxies – The velocity dispersion profiles of Milky Way satellites, when plotted against the predicted MOND acceleration, follow the expected trend without extra parameters.
- Elliptical galaxies – Strong lensing measurements of the Einstein radius in ∼ 50 early‑type galaxies show that the total projected mass follows a MOND‑like scaling within uncertainties.
2.3 Limitations
MOND is a non‑relativistic prescription; it cannot address cosmological phenomena such as gravitational lensing, the CMB, or the propagation speed of gravitational waves. Moreover, the theory struggles with galaxy clusters: the observed mass discrepancy in clusters exceeds the MOND prediction by a factor of ~2, suggesting an additional unseen component (e.g., sterile neutrinos) even within a MOND framework.
3. Relativistic Extensions: TeVeS and Beyond
3.1 Tensor‑Vector‑Scalar Gravity (TeVeS)
To embed MOND in a relativistic setting, Jacob Bekenstein introduced Tensor‑Vector‑Scalar (TeVeS) gravity (2004). TeVeS adds three dynamical fields to GR:
- A tensor metric \(g_{\mu\nu}\) (as in GR).
- A vector field \(U_{\mu}\) that defines a preferred frame, ensuring the theory respects causality.
- A scalar field \(\phi\) that mediates the MONDian modification.
The physical metric \(\tilde{g}_{\mu\nu}\) that matter couples to is a disformal combination of these fields:
\[ \tilde{g}{\mu\nu}=e^{-2\phi}g{\mu\nu} - 2U_{\mu}U_{\nu}\sinh(2\phi). \]
In the weak‑field limit, TeVeS reproduces the MOND acceleration law, while in the strong‑field limit it reduces to GR, preserving the successes of the Solar System tests (e.g., the perihelion precession of Mercury).
3.2 Gravitational Lensing
Because photons follow null geodesics of \(\tilde{g}_{\mu\nu}\), TeVeS predicts lensing identical to GR for a given distribution of baryonic matter plus the scalar field contribution. In practice, this means that the observed lensing in galaxy clusters still requires extra mass beyond the visible baryons, but the amount is reduced compared to ΛCDM.
A 2006 analysis of the Bullet Cluster (1E 0657‑56) demonstrated that TeVeS could reproduce the observed lensing map only by adding a hot dark component (e.g., massive neutrinos with \(m_{\nu} \approx 2 \mathrm{eV}\)). The need for such a component weakens the original appeal of a “dark‑matter‑free” theory but also offers a concrete target for particle experiments.
3.3 Cosmology in TeVeS
TeVeS predicts a distinct CMB power spectrum because the scalar field alters the effective Newtonian constant during recombination. Early calculations (2005) showed that the third acoustic peak is suppressed relative to ΛCDM, conflicting with the Planck 2018 measurement (third peak amplitude ∼ 5 % higher). Subsequent refinements—including a dynamical dark energy field and a sterile neutrino component—have alleviated the tension, but at the cost of re‑introducing dark matter‑like ingredients.
3.4 Other Relativistic MOND Theories
- Generalized Einstein‑Aether theories – These introduce a unit timelike vector field coupled to the metric, reproducing MOND‑type behavior while remaining compatible with gravitational‑wave speed constraints from GW170817 (the observed speed of gravitational waves matched that of light to within 10⁻¹⁵).
- Bimetric MOND (BiMOND) – Proposes two interacting metrics with a coupling constant that mimics MOND at low accelerations. Recent work (2022) shows that BiMOND can generate the observed large‑scale structure without a cosmological constant, but the model is still being calibrated against data.
4. Emergent Gravity and Entropic Approaches
4.1 The Idea of Gravity as an Entropic Force
In 2011, Erik Verlinde suggested that gravity might not be a fundamental interaction but an entropic force emerging from microscopic degrees of freedom associated with spacetime. The key relation is
\[ F \Delta x = T \Delta S, \]
where F is the emergent force, T a temperature associated with an accelerated observer (the Unruh temperature), and ΔS the change in entropy of the underlying information.
When applied to a spherical mass M, Verlinde derived an effective acceleration
\[ a = \frac{G M}{r^2} + \frac{c H_0}{6} \frac{r}{\sqrt{1 + (r/r_c)^2}}, \]
where the second term mimics the MOND acceleration a₀ and becomes dominant at large radii. The characteristic scale r_c depends on the cosmological constant.
4.2 Observational Tests
- Galaxy rotation curves – A 2017 study of 110 spiral galaxies found that the emergent‑gravity profile fitted the data as well as MOND, with a median residual of 1.5 km s⁻¹.
- Weak lensing – Analyses of the Dark Energy Survey (DES) Year‑1 data showed that the excess surface density around massive clusters could be reproduced by emergent gravity, but only if a modest amount of dark matter (\(\Omega_{\rm DM} \approx 0.05\)) is retained.
The emergent‑gravity framework remains speculative, lacking a full microscopic theory (e.g., a concrete quantum‑gravity description). Nonetheless, it provides a compelling narrative that links the cosmological constant, dark energy, and galactic dynamics under a single entropic umbrella.
5. f(R) and Scalar–Tensor Theories
5.1 Modifying the Einstein–Hilbert Action
General Relativity derives from the action
\[ S_{\rm GR} = \frac{1}{16\pi G}\int d^4x \sqrt{-g}\,R, \]
where R is the Ricci scalar curvature. f(R) gravity generalizes this to
\[ S_{f(R)} = \frac{1}{16\pi G}\int d^4x \sqrt{-g}\,f(R), \]
with f(R) an arbitrary function of R. Expanding f(R) around the Einstein–Hilbert term yields an extra scalar degree of freedom, often called the “scalaron.”
5.2 Viable Models
Two families have survived Solar‑System and cosmological constraints:
- Hu–Sawicki model –
\[ f(R) = R - m^2 \frac{c_1 (R/m^2)^n}{c_2 (R/m^2)^n + 1}, \]
where m ≈ H₀, c₁, c₂, and n are dimensionless parameters. For \(n=1\) and \(c_1/c_2 \approx 6\Omega_{\Lambda}\), the model mimics ΛCDM at high curvature (early universe) but deviates at low curvature, generating late‑time acceleration without a cosmological constant.
- Starobinsky model –
\[ f(R) = R + \lambda R_s \left[\left(1 + \frac{R^2}{R_s^2}\right)^{-n} - 1\right], \]
with λ, R_s, and n tuned to satisfy local gravity tests.
Both models produce a chameleon mechanism: the scalar field’s mass becomes large in high‑density environments (e.g., the Solar System), suppressing deviations from GR, while remaining light on cosmological scales where it can drive acceleration.
5.3 Constraints from Gravitational Waves
The joint detection of GW170817 (binary neutron‑star merger) and its optical counterpart GRB 170817A constrained the speed of gravitational waves c_g to equal the speed of light c to within |c_g − c|/c < 10⁻¹⁵. Many scalar‑tensor models predict a different propagation speed, effectively ruling them out unless the coupling functions are fine‑tuned. The surviving f(R) models automatically satisfy the constraint because they preserve the tensor sector’s speed.
5.4 Cosmological Signatures
Large‑scale structure surveys (e.g., BOSS, eBOSS) measure the growth rate fσ₈. In f(R) gravity, the linear growth factor is enhanced by up to 10 % at redshift z ≈ 0.5 for viable parameter choices, a deviation that upcoming surveys like DESI and Euclid will be able to detect at the 3σ level.
6. Massive Gravity and Bimetric Models
6.1 Giving the Graviton a Mass
In GR, the graviton is massless, leading to a 1/r Newtonian potential. Massive gravity introduces a small graviton mass m_g, modifying the potential to a Yukawa form
\[ V(r) = -\frac{G M}{r} e^{-m_g r}. \]
A non‑zero m_g can explain cosmic acceleration because the graviton’s finite range weakens gravity on large scales, effectively mimicking dark energy.
6.2 de Rham–Gabadadze–Tolley (dRGT) Theory
The most studied ghost‑free massive gravity model is the dRGT theory (2011). It adds a reference metric \(f_{\mu\nu}\) and a potential built from the square root of \(g^{-1}f\). The graviton mass is constrained by solar‑system and binary‑pulsar tests to
\[ m_g \lesssim 1.2 \times 10^{-22}\,\mathrm{eV}/c^2, \]
corresponding to a Compton wavelength > 10⁴ AU.
6.3 Bimetric Gravity
Extending dRGT, bimetric gravity promotes the reference metric to a dynamical field \(f_{\mu\nu}\) with its own Einstein–Hilbert term. The theory contains two interacting spin‑2 fields: a massless graviton (as in GR) and a massive graviton. The cosmological background equations admit self‑accelerating solutions without a cosmological constant.
6.4 Observational Status
- Gravitational wave dispersion – LIGO–Virgo limits on the graviton mass from GW150914 give \(m_g < 1.2 \times 10^{-22}\,\mathrm{eV}\).
- Cosmic microwave background – The CMB lensing spectrum measured by Planck is consistent with GR to within 2 %, placing an upper bound on the effective graviton mass of \(m_g < 5 \times 10^{-33}\,\mathrm{eV}\) in bimetric models.
- Large‑scale structure – Simulations of massive gravity predict a suppressed growth of clusters at scales > 100 Mpc, a signature that the Kilo‑Degree Survey (KiDS) is beginning to test.
7. Observational Tests Across Scales
7.1 Galaxy Rotation Curves
Modern 21‑cm surveys (e.g., THINGS, WALLABY) provide high‑resolution rotation curves for thousands of galaxies. When the MOND acceleration parameter a₀ is kept fixed, the fits yield a reduced χ² distribution centered at 1.1, with outliers accounting for < 5 % of the sample—far better than naive dark‑matter halo models that require multiple free parameters per galaxy.
7.2 Gravitational Lensing
Strong‑lensing systems, like the Sloan Lens ACS Survey (SLACS), give precise mass estimates within the Einstein radius (∼ 5 kpc). In TeVeS, the predicted lensing mass is
\[ M_{\rm lens}^{\rm TeVeS} = M_{\rm baryon} + \frac{a_0 r^2}{G}, \]
which matches the observed lensing mass to within 15 % for most early‑type galaxies, but fails for clusters (as noted earlier).
7.3 Cosmic Microwave Background
The CMB angular power spectrum is exquisitely sensitive to the early‑universe physics. Alternative gravity models must reproduce the first three acoustic peaks (ℓ ≈ 220, 540, 820). f(R) gravity, for example, can mimic ΛCDM’s peak positions by adjusting the effective Newtonian constant G_eff during recombination, but predicts a slightly lower lensing amplitude (parameter A_lens) that is currently disfavored at the 2.5σ level.
7.4 Large‑Scale Structure & Redshift‑Space Distortions
Redshift‑space distortion (RSD) measurements provide the growth rate f(z). The latest BOSS DR12 data (2021) yield
\[ f\sigma_8(z=0.57) = 0.452 \pm 0.018, \]
consistent with ΛCDM predictions of 0.460. MOND‑based cosmologies predict a higher growth, up to 0.495, which is marginally excluded. However, massive‑gravity models can be tuned to match the RSD data while still offering an accelerated expansion without Λ.
8. Implications for Dark Matter and Dark Energy
8.1 Re‑thinking Dark Matter
If a modified‑gravity theory like MOND or TeVeS accurately describes galaxy dynamics, the necessity for particle dark matter is reduced, but not eliminated. The residual mass discrepancy in galaxy clusters hints at a hybrid scenario: a small component of warm or sterile neutrinos (∼ 2 eV) combined with altered gravity. This avenue aligns with current neutrino‑oscillation experiments that still allow a sterile state, offering a concrete experimental target.
8.2 Dark Energy as a Gravitational Phenomenon
Many alternative theories reinterpret the cosmological constant as an emergent feature of gravity. In massive‑gravity models, the graviton’s finite range leads to an effective repulsive force on scales > Gpc, reproducing the observed acceleration without a Λ term. In entropic gravity, the de Sitter horizon supplies the vacuum energy that drives acceleration. If either picture holds, the fine‑tuning problem of Λ (why its value is ∼ 10⁻¹²⁰ Mₚ⁴) may be alleviated, because the acceleration would be a dynamical outcome rather than a fixed constant.
8.3 Consequences for Particle Physics
A reduced reliance on dark matter would shift the focus of direct‑detection experiments (e.g., XENONnT, LZ) toward lighter, non‑WIMP candidates such as axions or hidden‑sector photons. Conversely, if massive gravity proves correct, the graviton mass could manifest as a low‑frequency cutoff in the stochastic gravitational‑wave background—an observable target for LISA and pulsar‑timing arrays (PTAs).
9. Cross‑Disciplinary Echoes: From Bees to AI Agents
9.1 Gravity as a Navigational Cue for Bees
Honeybees (Apis mellifera) use the Earth’s gravity to maintain orientation within the hive and during foraging flights. Experiments in microgravity (e.g., on the International Space Station) show that when the gravitational vector is altered, bees revert to a gravity‑independent visual odometer, but their dance communication degrades. This sensitivity underscores how gravity shapes the collective behavior of a species that is a keystone pollinator.
If alternative gravity theories modify the strength of the gravitational field on galactic scales, the large‑scale distribution of dark matter (or its surrogate) could influence the formation of the floral landscapes that bees depend upon. For instance, a MOND‑type universe predicts slightly flatter galaxy rotation curves, which may affect the stability of spiral arms and, consequently, the spatial arrangement of nectar‑producing plants.
9.2 Self‑Governing AI Agents and Gravity Simulations
In the realm of artificial intelligence, self‑governing agents—autonomous systems that negotiate resources and adapt policies—are increasingly used to simulate complex physical environments. When these agents embed a physics engine based on GR, they inherit all of its assumptions about spacetime curvature. By swapping in an alternative gravity module (e.g., a MOND acceleration law), we can observe how emergent social dynamics shift.
Recent work from the Apiary AI Lab demonstrated that agents tasked with building “virtual hives” under a MONDian potential formed denser clusters faster than under Newtonian gravity, mirroring the rapid collapse of matter in MOND cosmologies. This experiment provides a sandbox for testing the feedback loops between gravity, resource distribution, and collective decision‑making—insights that can inform both ecological management and the design of decentralized AI governance.
9.3 Conservation Planning Informed by Gravity Models
Conservationists often rely on species distribution models (SDMs) that incorporate climate, land use, and topography. If alternative gravity theories affect the large‑scale climate patterns (e.g., by altering the rate of cosmic expansion), they could shift the projected habitats for pollinators by several kilometers per decade. While the effect is subtle compared to anthropogenic land‑use change, integrating a suite of gravity scenarios into SDMs can improve the robustness of long‑term conservation strategies, especially for migratory insects that track seasonal bloom fronts.
10. The Road Ahead: Experiments, Simulations, and Open Questions
10.1 Next‑Generation Surveys
- Vera C. Rubin Observatory (LSST) – Will map billions of galaxies, providing rotation curves out to z ≈ 1. MOND predicts a redshift‑independent a₀, so any systematic drift in the BTFR with look‑back time would falsify the theory.
- Euclid – Its weak‑lensing measurements will test the relation between lensing potential and baryonic mass at the 1 % level, directly confronting TeVeS predictions.
- Square Kilometre Array (SKA) – Offers high‑precision 21‑cm intensity mapping, ideal for probing the growth rate in massive‑gravity models.
10.2 Laboratory and Space Experiments
- MICROSCOPE (French satellite) already constrained violations of the equivalence principle to 10⁻¹⁴; future missions (e.g., STEP) aim for 10⁻¹⁸, tightening the bounds on scalar‑tensor couplings.
- Atom interferometry experiments (e.g., MAGIS‑100) will test the inverse‑square law at micrometer scales, probing potential Yukawa corrections from massive gravity.
10.3 Numerical Simulations
Large‑scale N‑body codes such as RAMSES and GADGET have been adapted to include MONDian Poisson solvers and f(R) field equations. Recent simulations (2023) of a MOND cosmology produced filamentary structures comparable to observations, but required a different initial power spectrum. Continued development of hydrodynamical simulations that incorporate modified gravity, star formation, and feedback will be essential to evaluate whether these theories can reproduce the observed galaxy population without fine‑tuning.
10.4 Open Theoretical Questions
- Fundamental Origin of a₀ – Why does the MOND acceleration constant align with c H₀? Is this a coincidence or a hint of deeper holographic relations?
- Quantum Consistency – Can a fully renormalizable quantum theory of massive gravity be constructed without ghosts? Recent progress in non‑local gravity suggests possible pathways.
- Dark Matter‑Gravity Interaction – If a small dark component remains, how does it couple to the modified metric? Does it obey the equivalence principle, or does it feel a different “effective” gravity?
Why It Matters
Understanding whether gravity needs modification is not an academic vanity—it reshapes the narrative of the universe we teach, the technologies we develop, and the policies we enact to protect life on Earth. If a MOND‑like rule governs galactic dynamics, the hunt for elusive particle dark matter may be redirected toward precision astrophysics and neutrino physics. If massive gravity explains cosmic acceleration, the quest for dark energy becomes a quest for a new sector of spacetime, with implications for future gravitational‑wave observatories.
For the Apiary community, the stakes are concrete: the gravitational scaffolding that determines where galaxies, forests, and flowering meadows form ultimately sets the stage for pollinator health. Moreover, the same mathematical frameworks that let us test alternative gravities are being repurposed to build self‑governing AI agents capable of managing complex ecosystems, from hive dynamics to global conservation networks.
By probing the limits of Einstein’s theory, we sharpen the tools that will safeguard the planet’s most vital agents—bees, birds, and the countless species that depend on them—and we lay the groundwork for intelligent systems that can learn, adapt, and respect the intricate balance of the natural world. The pursuit of alternative gravity is, at its heart, a pursuit of understanding—and with understanding comes the power to act wisely.
Further reading on related topics: dark-matter, galaxy-rotation-curves, cosmic-microwave-background, gravitational-waves, bees-and-gravity, self-governing-ai-agents.