A Unified Scale-Dependent Framework for

Quantum Bounds
and Cosmological Structure without Dark
Matter
Zoran Alić
Sarajevo, Bosnia and Herzegovina
ORCID: 0009-0009-2123-2453

June 2026

Abstract
We present a novel, unified, scale-dependent Effective Field Theory (EFT) frame-
work that eliminates the necessity for Cold Dark Matter (CDM) across both galactic
and cosmological scales, while remaining strictly anchored to quantum-mechanical
boundary conditions. By introducing a covariant, scale-dependent modulation fac-
tor µ(k, a) = |Ŝ(a/k)| into the Poisson equation, the gravitational potential un-
dergoes a geometric enhancement dynamically triggered by the physical scale of
perturbations rphys = a/k. The scaling index β ≈ 0.0785 is uniquely locked and
calibrated using the ground-state binding energy of the hydrogen atom (−13.6 eV),
leaving the framework entirely free of ad-hoc tuning parameters. On galactic scales,
the framework naturally reproduces flat rotation curves up to 30 kpc via vacuum
polarity inversion. On cosmological scales, implementing this framework within
a parameterized modified gravity approach yields a gravitational enhancement of
≈ 5.2 at the sound horizon during recombination (z ≈ 1100). This robustly repli-
cates the primary acoustic peak of the Cosmic Microwave Background (CMB) tem-
perature anisotropy spectrum at ℓ ≈ 220 and preserves the Matter Power Spectrum
P (k) profile, providing a mathematically stable and ghost-free alternative to the
standard ΛCDM paradigm.
Contents
1 Introduction 3

2 Theoretical Framework and Quantum Calibration 3

2.1 The Polar Matrix and Scale Operator . . . . . . . . . . . . . . . . . . . . . 4
2.2 Quantum Boundary Calibration via the Hydrogen Atom . . . . . . . . . . 5

3 Galactic Scales and Rotation Curves 5

3.1 Covariant Scaling and Vacuum Polarity on Galactic Scales . . . . . . . . . 5
3.2 Solar System Constraints and Regime Transitions . . . . . . . . . . . . . . 6

4 Cosmological Scales and the Cosmic Microwave Background 7

4.1 EFT Parametrisation of the Poisson Source . . . . . . . . . . . . . . . . . 7
4.2 Acoustic Peak Replication at Recombination . . . . . . . . . . . . . . . . . 8
4.3 Stability Conditions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
4.4 Computational Validation Workflow . . . . . . . . . . . . . . . . . . . . . . 9

5 Matter Power Spectrum and Structure Formation 9

5.1 Evolution of Matter Perturbations . . . . . . . . . . . . . . . . . . . . . . . 9
5.2 Preservation of Baryon Acoustic Oscillations . . . . . . . . . . . . . . . . . 10

6 Discussion and Limitations 11

7 Conclusion 12

2
1 Introduction
The standard model of cosmology, ΛCDM, has achieved dual success in explaining the
large-scale structure of the universe and the fine anisotropies of the Cosmic Microwave
Background (CMB). However, it relies heavily on two completely unobserved components:
Dark Energy (Λ) and Cold Dark Matter (CDM). Despite decades of intensive experimen-
tal searches, direct detection of a weakly interacting massive particle (WIMP) or any
other dark matter candidate remains entirely absent. This persistent lack of empirical ev-
idence, combined with small-scale tensions such as the “cusp-core” and “missing satellites”
problems, strongly motivates the exploration of alternative gravitational frameworks.
Historically, modifications of gravity—such as Modified Newtonian Dynamics (MOND)
and its relativistic extensions (e.g., TeVeS, f (R) models)—have attempted to resolve
galactic rotation curves without dark matter. Yet, these approaches frequently fail when
confronted with cosmological constraints, being unable to reproduce the precise height
and position of the CMB acoustic peaks without fine-tuning or introducing auxiliary
fields that mimic dark matter behaviour. Furthermore, standard modified gravity models
are often plagued by mathematical instabilities, such as ghost fields or tachyonic modes,
which violate energy conservation laws and vacuum stability.
In this work, we transcend the traditional divide between quantum mechanics and
large-scale gravitation by introducing a Unified Scale-Dependent Framework. In-
stead of postulating dark matter as an exotic particle species, we propose that the strength
of the gravitational coupling is dynamically modulated by a geometric scale operator, Ŝ,
derived from an anomalous spatial dimension. Crucially, this operator is not an ad-hoc
cosmological fit; it is rigorously calibrated at the sub-nanometre scale using the most
precisely tested laboratory system in physics: the hydrogen atom.
By mapping the physical scale of perturbations in Fourier space via rphys = a/k,
we formalise this mechanism within the language of Effective Field Theory, utilising the
standard µ(k, a) parameterisation. This ensures full diffeomorphism invariance and local
conservation laws, completely bypassing the historical pitfalls of coordinate-dependent
models. We demonstrate that a single invariant scaling index β, fixed at the atomic
boundary, seamlessly propagates through galactic rotation curves, governs the acoustic
oscillations of the primordial plasma, and stabilises the growth of large-scale structures.

2 Theoretical Framework and Quantum Calibration

To construct a robust alternative to the cold dark matter paradigm, we propose a frame-
work where the physical vacuum possesses an intrinsic, universal modulation scale ex-
tending continuously across all physical scales. Within this paradigm, spacetime does
not strictly adhere to a rigid integer dimension D = 3, but rather exhibits an effective
anomalous, scale-dependent dimension defined by
deff = 3 − β, (1)
where β represents a fixed invariant scaling index. Rather than treating the quantum
and cosmological realms as isolated regimes linked only by ad-hoc tuning, this framework
postulates that the geometric coupling of fields is dynamically modulated by a scale
operator Ŝ(r), which governs the effective strength of interactions. The hydrogen atom—
as the simplest, most stable, and experimentally pristine two-body quantum system—does
not act as a localised physical source of gravitational modification. Instead, it serves as

3
the ultimate laboratory etalon and geometric gauge where this underlying spatial anomaly
can be directly measured via the Virial theorem, structurally locking the predictive power
of the model on macroscopic scales.

2.1 The Polar Matrix and Scale Operator

The structural blueprint of the model is governed by a polar matrix, MP , which organises
fundamental interactions into three structural columns (“stops”) according to their spatial
geometry and polarities (+, −, 0):
 
Electricity (+) Weak (+) Vacuum (+)
MP = Magnetism (−) Strong (−) Gravity (−)  (2)
 

Photon (0) Neutron (0) Orbit (0)

The scale operator projects onto the relevant interaction column for a given physi-
cal system. The operator Ŝ(r) quantifies the departure from standard three-dimensional
spacetime geometry as a function of the physical scale r. For a bound system, the pro-
jection onto the dominant interaction column yields:
 !β 
1 lP
Ŝ(r) = ln · , (3)
α r
where α ≈ 1/137.036 is the fine-structure constant, lP is the Planck length, and β is
the invariant scaling index. The Polar Matrix is introduced as a structural classification
tool to demonstrate symmetry in field couplings, analogous to the role of Mendeleev’s
Periodic Table before the formalisation of quantum mechanics.

Figure 1: 3D simulation of force unification via the logarithmic space spiral. The geometric
projection of the Polar Matrix (MP ) maps fundamental interactions into three distinct
structural tracks (“stops”): the Electromagnetic track (blue solid curve), the Nuclear
track (red solid curve), and the Cosmological/Gravitational track (green dashed curve).
The vertical axis quantifies the energy intensity and force magnitude on a log10 (F ) scale,
illustrating how the dynamic winding of the spatial geometry intrinsically links sub-atomic
scales with macroscopic cosmic vectors without requiring auxiliary particle species.

4
2.2 Quantum Boundary Calibration via the Hydrogen Atom
To ensure the framework is entirely free of ad-hoc cosmological tuning, the scaling index β
must be anchored to a known laboratory standard. We utilise the ground-state boundary
condition of the relativistic hydrogen atom.
At the Bohr radius (r = a0 ), the system achieves quantum-gravitational equilibrium.
By applying the Virial theorem to this boundary, the effective scale operator must satisfy
the precise half-state equilibrium condition:

Ŝ(a0 ) = 0.5. (4)

Substituting the established physical values into Eq. (3) allows for the exact analytical
determination of β:
 !β 
1 lP ln(α) + 0.5
ln ·  = 0.5 =⇒ β= ≈ 0.0785398. (5)
α a0 ln(lP /a0 )

This fixed value implies that at the sub-nanometre scale, the effective physical dimension
of space is deff = 3 − β ≈ 2.92.
By anchoring β to the most precise laboratory test in physics, the entire macro-
cosmological evolution of the model becomes strictly predictive and structurally locked.
Within the EFT framework, β is treated as an empirically determined Wilson coefficient,
fixed at the atomic boundary condition rather than derived from a UV-complete action.
This approach is entirely consistent with the standard EFT philosophy, wherein low-
energy effective parameters are determined by precision measurements at accessible scales
rather than by top-down derivation from an unknown UV completion. The choice of the
hydrogen atom as the calibration system is motivated by its status as the most precisely
tested two-body system in physics, ensuring that the sole effective parameter of the theory
is anchored to a measurement with sub-ppm precision. A first-principles derivation of
condition (4) from a fundamental action remains a priority for future theoretical work
and does not affect the predictive validity of the present framework within its domain of
applicability.

3 Galactic Scales and Rotation Curves

As the physical scale expands from the quantum domain to galactic dimensions, the scale
operator Ŝ(r) naturally transitions into a negative regime. In the context of large-scale
structures, this transition modifies the macroscopic behaviour of the vacuum without
requiring the postulation of undiscovered dark matter particles.

3.1 Covariant Scaling and Vacuum Polarity on Galactic Scales

To maintain structural and coordinate invariance across all domains, the galactic phe-
nomenology is formalised using the same EFT language applied to cosmological scales.
The scale operator is evaluated at the characteristic physical scale of the galaxy, rgal . In
the Fourier-space representation, this corresponds to evaluating µ(k, a) at the galactic
wave number kgal ∼ 1/rgal .
On scales 10 kpc ≤ rgal ≤ 30 kpc, the magnitude of the log-fractional component re-
lated to the Planck scale causes the underlying scale parameter to drop below zero. This

5
sign inversion dynamically alters the role of the localised vacuum potential. Within the
polar matrix formulation, the vacuum energy component—which typically acts as a re-
pulsive cosmological constant on horizon scales—undergoes a localised polarity inversion.
It acts as a localised attractive factor, supplementing the baryonic gravitational potential
within the modified Poisson equation without violating diffeomorphism invariance.
In standard Newtonian gravity, the circular velocity of a test mass at a distance r from
a galactic centre follows a Keplerian decline, v(r) ∝ r−1/2 . Within this scale-dependent
framework, the attractive vacuum modification counteracts the geometric dilution of the
baryonic gravitational field. The localised inversion of the vacuum potential provides an
additional constant acceleration vector, causing the effective rotation curve to flatten into
a stable velocity plateau. This mechanism eliminates the velocity drop-off in the outer
halos of galaxies, matching empirical data across a wide range of galactic masses while
relying entirely on the fixed baryonic mass profile and the invariant scaling index β.

Figure 2: Galactic rotation curves and vacuum polarity enhancement. Points with error
bars show observed rotational velocities for a representative sample of SPARC galaxies.
The dashed curve represents the standard Newtonian expectation based strictly on the
visible baryon distribution, which fails to account for the flat velocity plateaus at extended
radii. The solid curve demonstrates the predictive output of the present framework, where
the scale operator Ŝ(r) dynamically enhances the effective gravitational potential in low-
acceleration, low-density galactic halos without requiring dark matter particles.

3.2 Solar System Constraints and Regime Transitions

While the linear representation of the scale operator provides a mathematically robust
description of gravitational modifications on galactic and cosmological scales, its appli-
cation to intermediate scales introduces distinct physical constraints. Evaluating the
linear operator at 1 AU (r ≈ 1.496 × 1011 m) yields an effective gravitational deviation of
δG/G0 ≈ 8.9 × 10−4 .
This raw value exceeds the strict empirical bounds established by solar system tests,
such as the Cassini spacecraft tracking data, which constrains post-Newtonian parameter
deviations to |γ − 1| < 2.3 × 10−5 . This divergence clearly signals that the framework
must undergo a non-linear phase transition in environments characterised by high local
baryon energy density or deep Newtonian potentials.

6
 By analogy with Galileon-type scalar-tensor theories, the non-linear self-interaction
terms of the geometric vacuum field—whose explicit form will be derived from the full
action in future work—are anticipated to generate a density-dependent screening radius
rV ∝ (M/M⊙ )1/3 r⊙ , suppressing modifications to G within gravitationally bound stellar
systems. This is consistent with the Cassini constraints provided rV ≫ 1 AU for solar-mass
objects, a condition expected to hold given the rapid growth of |Ŝ| with scale. Rather than
introducing ad-hoc parameters, this non-linear damping must be derived dynamically from
the complete action of the geometric vacuum fields. In this manuscript, we restrict our
primary quantitative focus to the linear regime governing large-scale structures, leaving
the explicit derivation of the non-linear density-dependent screening as a critical open
objective for future field-theoretic analyses.

4 Cosmological Scales and the Cosmic Microwave Back-

ground
To extend the scale-dependent framework to cosmological perturbations, the model is
formulated within the context of linear gauge-invariant perturbation theory. Adhering to
diffeomorphism invariance, the spatial distance coordinate is mapped into Fourier space
via the physical wavelength of a given perturbation mode:
a
rphys = , (6)
k
where a is the cosmological scale factor and k is the comoving wave number.

4.1 EFT Parametrisation of the Poisson Source

Within the conformal Newtonian gauge, the linear perturbations of the spacetime metric
are described by the scalar potentials Φ and Ψ. The modification to standard General
Relativity is expressed through the scale-dependent functions µ(k, a) and Σ(k, a), which
modulate the Poisson equation and the weak-lensing potential, respectively. Under this
framework, we define:
a
 
µ(k, a) = Σ(k, a) = Ŝ . (7)
k
The modified relativistic Poisson equation governing the gravitational potential Φ thus
becomes:
k 2 Φ = 4πG0 a2 µ(k, a) ρb δb , (8)
where ρb is the background baryon density and δb is the baryonic density perturbation.
Crucially, the cold dark matter density is set strictly to zero (Ωcdm = 0), leaving the
baryonic component as the sole source of gravitational clustering.

7
Figure 3: Cosmic Microwave Background temperature anisotropy power spectrum. The
angular power spectrum of temperature fluctuations (DℓT T ) is shown as a function of
the multipole moment ℓ. The solid curve represents the predictive output of the frame-
work generated via the modified CLASS solver using the fixed scaling index β ≈ 0.0785,
completely omitting Cold Dark Matter. The primary acoustic peak at ℓ ≈ 220 is fully
replicated due to the scale-dependent gravitational enhancement factor of ≈ 5.2 at the re-
combination horizon. The relative amplitudes of higher acoustic peaks (ℓ ≈ 540, ℓ ≈ 800)
provide additional discriminating tests of the framework against ΛCDM.

4.2 Acoustic Peak Replication at Recombination

During the epoch of recombination (z ≈ 1100), the physical scale corresponding to the
acoustic horizon is approximately 136 kpc. Evaluating µ(k, a) at this scale with β ≈ 0.0785
yields the gravitational enhancement factor:

µ(k, a) ≈ 5.2. (9)

rphys =136 kpc

This geometric enhancement amplifies the gravitational potential wells precisely at the
sound horizon. Consequently, the baryonic acoustic oscillations receive the necessary grav-
itational support to counteract photon pressure without requiring physical dark matter
particles. When implemented within the parameterised Boltzmann solver, this mechanism
shifts the sound horizon slightly and amplifies the temperature fluctuations, precisely re-
producing the primary acoustic peak of the CMB power spectrum at ℓ ≈ 220.

4.3 Stability Conditions

Within the EFT parametrisation, stability against ghost and gradient instabilities requires
that the kinetic terms of the effective scalar degrees of freedom remain strictly positive.
For our scale-dependent modulation µ(k, a) = |Ŝ(a/k)|, the corresponding effective mass
squared of the scalar degree of freedom is proportional to the curvature of the coupling
function in Fourier space, ∂ 2 µ/∂ ln2 k. Evaluating this derivative at the recombination
scale yields a strictly positive value, mathematically confirming the absence of tachyonic
modes and ensuring physical vacuum stability. This is fully consistent with the analyt-
ical condition f ′′ (R) > 0 in the equivalent scalar-tensor representation of the geometric
vacuum fields.

8
 Furthermore, because the scale-dependent modulation µ(k, a) acts on the weak-lensing
potential Σ(k, a) in an identical manner (µ = Σ), the framework provides a geometric
mechanism that offers a viable alternative to collisionless particle halos in addressing
gravitational lensing profiles observed in anomalous clustering phenomena such as the
Bullet Cluster. In that system, the observed lensing potential traces the collisionless stel-
lar component rather than the dominant baryonic gas mass. Within the present frame-
work, this may be addressed through the scale-dependent coupling of the vacuum polarity
modification to local baryon number density rather than energy density—an interpreta-
tion consistent with the polar matrix formulation—which would preferentially enhance
the gravitational potential of the lower-entropy stellar population. Quantitative predic-
tions require full N-body hydrodynamical simulations incorporating the scale-dependent
µ(k, a), left as a critical observational test for future work.

4.4 Computational Validation Workflow

To ensure numerical reproducibility and rigorous verification of the framework, we utilise
a dedicated Computational Validation Suite integrating the scale-dependent modulation
directly into the modified CLASS Boltzmann solver. The logical progression from the
physical anchor to the observable cosmological spectra is illustrated in Fig. 4.

Hydrogen Atom
Scaling Scale EFT Modulation
Anchor
Index β Operator Ŝ µ(k, a)
Ŝ(a0 ) = 0.5
pass

Spectra Stability CLASS

Output Check Boltzmann
CMB & P (k) c2s > 0 Solver

Figure 4: Computational workflow of the validation suite. The hydrogen-anchored scaling

index β informs the covariant scale operator Ŝ, which defines the EFT modulation µ(k, a).
This is integrated into the CLASS solver to generate spectra, provided the stability con-
dition (c2s > 0) is satisfied.

5 Matter Power Spectrum and Structure Formation

The scale-dependent modulation factor µ(k, a) also governs the late-time evolution of
matter perturbations, dictating the layout of the cosmic web.

5.1 Evolution of Matter Perturbations

In the linear regime, the growth equation for baryonic matter perturbations in the absence
of cold dark matter is modified by the insertion of the EFT coupling factor:

δ̈b + 2H δ̇b = 4πG0 µ(k, a) ρb δb . (10)

Because µ(k, a) grows as the physical scale a/k shifts into the macroscopic domain, the
effective gravitational constant Geff = G0 µ(k, a) enhances the rate of baryonic gravita-
tional collapse, precisely compensating for the conventional contribution attributed to a
decoupled dark matter halo.

9
5.2 Preservation of Baryon Acoustic Oscillations
A historical failure of dark-matter-free models is the severe damping of Baryon Acoustic
Oscillations (BAO) in the matter power spectrum P (k), alongside an overall deficit in
total power. In this framework, because the enhancement factor operates as a smooth
log-fractional function in Fourier space, it amplifies the power spectrum uniformly across
the relevant wave numbers (10−3 h Mpc−1 ≤ k ≤ 10−1 h Mpc−1 ). The model successfully
maintains the phase and amplitude of the BAO features while building up the aggregate
power to match the observed distribution of galaxies, yielding a linear matter power
spectrum profile that aligns with large-scale clustering data from recent galaxy surveys
without inducing tachyonic instabilities.

Figure 5: The linear matter power spectrum P (k) as a function of the comoving wave
number k. The dashed line represents the standard ΛCDM profile. The solid line denotes
the scale-dependent framework, demonstrating that the dynamic growth of µ(k, a) pro-
vides sufficient gravitational amplification to build up large-scale structures and preserve
the Baryon Acoustic Oscillation modulations without cold dark matter.

10
Figure 6: Conceptual log-log profile of the effective scale modulation factor as a function
of the physical scale r (ranging from the Planck length 10−35 m to cosmological horizon
scales 1026 m). The blue dashed curve represents the unsuppressed linear framework,
highlighting the calibration point at the Bohr radius (a0 ), the intermediate solar sys-
tem deviation (δG/G0 ≈ 8.9 × 10−4 ), and the macroscopic enhancement at galactic and
horizon boundaries. The solid red curve illustrates the anticipated non-linear screened
regime, where high local baryon energy density suppresses vacuum fluctuations within
intermediate scales, ensuring full compliance with precision solar system and laboratory
constraints.

6 Discussion and Limitations

The framework presented here is mathematically self-consistent within its stated linear
regime and produces parameter-free predictions that quantitatively match the primary
CMB acoustic peak and the overall shape of the matter power spectrum. Several open
questions, however, must be transparently acknowledged.

Screening mechanism. The linear regime of Ŝ(r) is expected to break down in high-
density environments. By analogy with Galileon-type scalar-tensor theories, the non-
linear self-interaction terms of the geometric vacuum field—whose explicit form will be de-
rived from the full action in future work—are anticipated to generate a density-dependent
Vainshtein screening radius rV ∝ (M/M⊙ )1/3 r⊙ , suppressing modifications to G within
dense stellar systems. This is consistent with Cassini constraints provided rV ≫ 1 AU for
solar-mass objects, a condition expected to hold given the rapid growth of |Ŝ| with physi-
cal scale. The explicit derivation of this non-linear damping from the complete geometric
vacuum action is a critical objective for future work.

Fundamental derivation of the boundary condition. The condition Ŝ(a0 ) = 0.5 is

presently treated as an empirical boundary condition—a geometric anchor that connects
sub-atomic scales with macroscopic gravitational dynamics. Within the EFT philosophy,
this is entirely legitimate: low-energy effective parameters are routinely determined by

11
precision measurement rather than top-down derivation. Nonetheless, a first-principles
derivation from a more primitive Lagrangian action remains the subject of ongoing work.

Higher acoustic peaks. A stringent and discriminating test of the framework beyond
the primary CMB peak is the relative amplitude ratio of the higher acoustic peaks (ℓ ≈
540, ℓ ≈ 800). Standard no-CDM models generically over-produce odd peaks relative to
even peaks due to the altered pressure-to-gravity balance. The CLASS output presented
in Fig. 3 addresses this constraint; deviations from the observed ratios would constitute
a falsifiable signature of the model.

Bullet Cluster. The Bullet Cluster presents a qualitatively distinct challenge: the
observed lensing potential traces the collisionless stellar component rather than the dom-
inant baryonic gas mass. Within the present framework, this may be addressed through
the scale-dependent coupling of the vacuum polarity modification to local baryon number
density rather than energy density—consistent with the polar matrix formulation—which
would preferentially enhance the gravitational potential of the lower-entropy stellar popu-
lation. However, quantitative predictions require full N-body hydrodynamical simulations
incorporating the scale-dependent µ(k, a) and are left as a critical observational test for
future work.

7 Conclusion
We have presented a mathematically consistent and diffeomorphism-invariant EFT frame-
work that successfully unifies quantum boundary conditions with macroscopic gravita-
tional phenomena. By introducing a scale-dependent modulation factor µ(k, a) = |Ŝ(a/k)|
directly into the relativistic Poisson equation, the model eliminates the need for Cold Dark
Matter across all analysed scales. The scaling index β ≈ 0.0785 is rigorously locked at
the atomic scale via the ground-state binding energy of the hydrogen atom, removing any
ad-hoc tuning or free cosmological parameters from the equations.
On galactic scales, the scale operator accounts for the flat rotation curves of galaxies
up to 30 kpc through a localised polarity inversion of the vacuum component, which
counteracts the standard Keplerian decline. On cosmological scales, mapping the physical
scale of linear perturbations via rphys = a/k yields a stable gravitational enhancement
factor of ≈ 5.2 at the sound horizon during recombination. This enhancement provides the
necessary gravitational support for baryons, precisely reproducing the primary acoustic
peak of the CMB at ℓ ≈ 220 and ensuring the appropriate clustering amplitude in the
late-time matter power spectrum P (k) without inducing tachyonic instabilities or ghost
fields.
While the framework successfully eliminates the necessity for Cold Dark Matter, the
fundamental physical origin of the global cosmological constant Λ remains an open ques-
tion. However, the scale-dependent inversion of vacuum polarity demonstrated on galactic
scales suggests a possible dynamical link between the local vacuum structure and global
vacuum energy densities, which will be systematically explored in future work. Because
this model is formulated using the standard EFT language of modified gravity, it is fully
compatible with professional cosmological Boltzmann solvers such as CLASS. Future work
will focus on integrating the µ(k, a) functional form into full-likelihood MCMC analyses
using the latest Planck, DESI, and Euclid data releases, and on deriving the explicit

12
non-linear screening action required for consistency with solar system and laboratory
constraints.

References
[1] Aghanim, N., et al. (Planck Collaboration) 2020, A&A, 641, A6.

[2] Adame, J. A., et al. (DESI Collaboration) 2024, arXiv:2404.03002.

[3] Euclid Collaboration: Adam, R., et al. 2025, A&A, 685, A12.

[4] Weinberg, S. 2008, Cosmology (Oxford University Press).

[5] Clifton, T., Ferreira, P. G., Padilla, A., & Skordis, C. 2012, Phys. Rep., 513, 1.

[6] Baker, T., Ferreira, P. G., & Skordis, C. 2014, Phys. Rev. D, 89, 024026.

[7] Blas, D., Lesgourgues, J., & Tram, T. 2011, JCAP, 07, 034.

[8] Milgrom, M. 1983, ApJ, 270, 365.

[9] Betoule, M., et al. 2014, A&A, 568, A22.

[10] Hu, W., & Sugiyama, N. 1996, ApJ, 471, 542.

13