Mathematical Models Of Interacting Dark Energy And Dark Matter

Two independent lines of evidence forced twentieth-century cosmology to accept components beyond ordinary matter: galactic rotation curves and cluster dynamics that only make sense with a dominant unseen mass (dark matter), and, at the end of the 1990s, supernova distance measurements showing that cosmic expansion is speeding up rather than slowing down (Riess et al., 1998; Perlmutter et al., 1999), which requires a component with negative pressure (dark energy). The ΛCDM model absorbs both into general relativity with minimal assumptions: dark matter as cold, pressureless dust, and dark energy as a constant Λ, with no interaction of any kind between them beyond shared gravity.

This economy comes at a conceptual cost. ΛCDM gives no reason why ρ_dm and ρ_de should be within an order of magnitude of one another at the present epoch — the coincidence problem — given that their ratio would otherwise sweep through many orders of magnitude over cosmic time. More recently, a second anomaly has sharpened this discomfort: the Hubble constant inferred from the CMB, Hâ‚€ = 67.4 ± 0.5 km/s/Mpc (Planck Collaboration, 2020), sits in persistent tension with the value measured from the local distance ladder, Hâ‚€ = 73.0 ± 1.0 km/s/Mpc (Riess et al., 2022). Interacting dark energy (IDE) models respond to both puzzles by proposing that dark matter and dark energy are not mutually blind but exchange energy directly, governed by a coupling function Q built into the cosmological continuity equations.

What follows sets up the mathematical machinery behind this proposal (Section 2), reviews the principal ways Q has been constructed in the literature (Section 3), tabulates what current data say about the dark sector (Section 4), and discusses what a working coupling would buy us in terms of the coincidence problem and the Hubble tension (Section 5), before closing in Section 6.

2. Setting Up the Dark-Sector Conservation Laws

2.1 Background geometry

Cosmology adopts, as a working idealization, a spatially flat, homogeneous, isotropic background — the FLRW metric, Eq. (1):

ds² = −dt² + a(t)² [ dr² + r²(dθ² + sin²θ dφ²) ]    (1)

Inserted into the Einstein field equations, this symmetry collapses the ten independent components of the field equations down to the two ordinary differential equations for the scale factor a(t) given in Eq. (2):

H² = (8πG/3) Σáµ¢ ρáµ¢ ,      ḁ/a = −(4πG/3) Σáµ¢ (ρáµ¢ + 3páµ¢)         (2)

with H = ā/a the Hubble parameter and the sum running over every energy component present — radiation, baryons, dark matter, and dark energy.

2.2 Where the coupling enters

General covariance demands only that the total stress-energy tensor be conserved; it says nothing about whether individual matter species must be separately conserved. IDE models exploit exactly this freedom, letting dark matter and dark energy trade energy while their sum still obeys the combined conservation law of Eq. (3):

ρ̇_dm + 3Hρ_dm = Q ,      ρ̇_de + 3H(1+w)ρ_de = −Q   (3)

Here w = p_de/ρ_de is the dark energy equation of state, and Q is the interaction kernel whose sign fixes the direction of flow: positive Q drains dark energy into dark matter, negative Q does the reverse. The entire model-building problem in this literature reduces to choosing a defensible functional form for Q, since the underlying microphysics of neither sector is known.

3. How the Coupling Has Been Built: Four Approaches

3.1 A scalar field that talks to matter directly

Wetterich (1988) first proposed that a light quintessence scalar φ could couple to dark matter if the dark matter particle mass itself depends on φ, m(φ); Amendola (2000) developed this into the now-standard “coupled quintessence” framework, with Q = βκφ̇ρ_dm for dimensionless coupling constant β and κ = √(8πG). Because the coupling descends from an explicit Lagrangian, this construction remains the reference point against which purely phenomenological alternatives are measured.

3.2 Guessing Q directly from dimensional analysis

A parallel and more pragmatic tradition skips the Lagrangian altogether and simply writes down plausible forms for Q built from H and the densities themselves — Q = 3Hβρ_dm, Q = 3Hβρ_de, or their sum. Wei and Zhang (2007) study the sign-changeable versions of these couplings, and Zimdahl (2012), in a systematic review, organizes the resulting late-time attractor solutions by which density the coupling is tied to, showing that the choice materially affects whether and how the coincidence problem gets resolved.

Tamanini (2015) raises an important caution here: introducing Q only at the background level, while leaving the perturbed Einstein equations untouched, is not generally self-consistent. A coupling of this kind is expected to seed an effective dissipative pressure once perturbations are included, with observable consequences for structure growth — meaning many of the phenomenological Q's in circulation are, strictly, incomplete models until their perturbative sector is worked out.

3.3 Deriving the coupling instead of assuming it

Ludwick (2018) takes a different route: rather than postulating Q, a non-minimal coupling ξφ²R between the dark energy scalar and the Ricci scalar R is introduced — a term that is in fact required for the field's renormalizability in curved spacetime. An effective energy exchange between the dark sectors then falls out automatically from this curvature term, without any additional assumption, and astrophysical data can be used to bound the resulting coupling strength and an associated dark matter mass scale.

3.4 Two fluids, or one field changing character?

Santos (2015) reframes the problem entirely with the Continuous Tower of Scalar Fields model: instead of two separate fluids exchanging energy, a continuum of scalar fields is postulated that behave like dark energy early on and gradually take on dark-matter-like behaviour later. Under this picture, Q is not an energy flux between fixed species but a bookkeeping device for the rate at which fields cross over from one dynamical regime to the other — a genuinely different resolution of the coincidence problem.

At smaller, galactic scales, Nari and Roshan (2025) show that Amendola's coupled quintessence produces a measurable dynamical friction force on massive bodies moving through a dark matter halo — an effect absent for ordinary, uncoupled quintessence — opening up an astrophysical test of the coupling that is independent of background cosmology altogether.

4. What the Data Currently Say

Table 1 collects representative constraints from CMB, supernova, BAO, and local-distance-ladder measurements bearing on dark-sector interaction models.

Figure 1. Numerically integrated dark-sector density parameters, Ω_dm(z) and Ω_de(z), comparing an uncoupled trajectory (c = 0) against a coupled trajectory (c = 0.15), with w = −0.95 held fixed; illustrative only.

CONCLUSION

Four distinct strategies — Lagrangian-derived coupled quintessence, phenomenological density-dependent couplings, curvature-induced coupling, and multi-field tower constructions — all converge on the same underlying mathematical device: a source term Q linking the dark matter and dark energy conservation equations. None has yet unseated ΛCDM as the preferred fit to precision data, but each offers a concrete, falsifiable mechanism for addressing the coincidence problem, and several offer at least a partial handle on the Hubble tension. Sharper constraints from upcoming CMB, weak-lensing, and galaxy-survey data, combined with the kind of perturbative consistency checks Tamanini's work calls for, should determine within the coming decade whether any dark-sector coupling is actually realized in nature.

REFERENCES

1. Amendola, L. (2000). Coupled quintessence. Physical Review D, 62(4), 043511.
2. Caldwell, R. R., Davé, R., & Steinhardt, P. J. (1998). Cosmological imprint of an energy component with general equation of state. Physical Review Letters, 80(8), 1582–1585.
3. Copeland, E. J., Sahni, V., & Tsujikawa, S. (2006). Dynamics of dark energy. International Journal of Modern Physics D, 15(11), 1753–1935.
4. Ludwick, K. J. (2018). Possible couplings of dark matter. In Redefining Standard Model Cosmology. IntechOpen. arXiv:1809.09971.
5. Nari, N., & Roshan, M. (2025). Dynamical friction by coupled dark energy. arXiv:2502.10676.
6. Perlmutter, S., et al. (1999). Measurements of Ω and Λ from 42 high-redshift supernovae. The Astrophysical Journal, 517(2), 565–586.
7. Planck Collaboration. (2020). Planck 2018 results. VI. Cosmological parameters. Astronomy & Astrophysics, 641, A6.
8. Riess, A. G., et al. (1998). Observational evidence from supernovae for an accelerating universe and a cosmological constant. The Astronomical Journal, 116(3), 1009–1038.
9. Riess, A. G., et al. (2022). A comprehensive measurement of the local value of the Hubble constant with 1 km/s/Mpc uncertainty from the Hubble Space Telescope and the SH0ES team. The Astrophysical Journal Letters, 934(1), L7.
10. Santos, P. (2015). The continuous tower of scalar fields as a system of interacting dark matter–dark energy. arXiv:1508.04809.
11. Tamanini, N. (2015). On phenomenological models of dark energy interacting with dark matter. Physical Review D, 92(4), 043524. arXiv:1504.07397.
12. Wang, B., Abdalla, E., Atrio-Barandela, F., & Pavón, D. (2016). Dark matter and dark energy interactions: Theoretical challenges, cosmological implications, and observational signatures. Reports on Progress in Physics, 79(9), 096901.
13. Wei, H., & Zhang, S. N. (2007). Cosmological constraints on the sign-changeable interactions. Physical Review D, 76(6), 063003.
14. Wetterich, C. (1988). Cosmology and the fate of dilatation symmetry. Nuclear Physics B, 302(4), 668–696.
15. Zhang, X. (2005). Coupled quintessence in a power-law case and the cosmic coincidence problem. Modern Physics Letters A, 20(31), 2575–2582. arXiv:hep-ph/0410292.
16. Zimdahl, W. (2012). Models of interacting dark energy. arXiv:1204.5892.