---
title: Dark Matter–Dark Energy Interactions
url: https://www.emergentmind.com/topics/dark-matter-dark-energy-interactions
type: topic
---

# Dark Matter–Dark Energy Interactions

Dark matter–dark energy (DM–DE) interactions are hypothesized non-gravitational couplings between the two dominant components of the cosmic energy budget. Such models have been developed to address theoretical puzzles like the coincidence problem and to alleviate empirical tensions in cosmology, such as those related to the Hubble constant ($H_0$) and the amplitude of structure growth ($S_8$ and $\sigma_8$). The literature has established a landscape of phenomenological, field-theoretical, and simulation-calibrated frameworks to study these interactions, with increasing observational constraints.

## 1. Theoretical Frameworks and Interaction Forms

Models of DM–DE interactions are constructed by modifying the standard cosmological conservation equations to include a direct energy–momentum exchange term $Q$:
\[
\dot\rho_c + 3H\rho_c = +Q, \qquad \dot\rho_{\text{de}} + 3H(1+w)\rho_{\text{de}} = -Q,
\]
where $\rho_c$ and $\rho_{\text{de}}$ are the dark matter and dark energy densities, $w$ is the DE equation of state, and $Q$ encodes the interaction rate [1204.5892][2402.00819].

**Canonical functional forms** for $Q$ found in the literature include:
- $Q = 3H\xi\rho_c$ (proportional to DM density)
- $Q = 3H\xi\rho_{\text{de}}$ (proportional to DE density)
- $Q = 3H\xi(\rho_c + \rho_{\text{de}})$ (linear in both densities)
- $Q = 3H\xi\,\rho_c^\alpha\,\rho_{\text{de}}^\beta$ (bilinear, generalized)
Field-theoretic models often generate $Q$ via a coupling between a quintessence field and DM, leading to terms like $Q \propto \dot\phi\rho_c$ or more complex structures involving both conformal and disformal couplings [1709.04882][2402.00819].

Stability considerations dictate the sign and allowed magnitude of $\xi$ for a given $w$ (typically $w \approx -1$ for accelerated expansion). For example, stability of early-time perturbations for $Q \propto \rho_{\text{de}}$ with $w \approx -1$ requires $\xi < 0$ [2002.06127].

## 2. Dynamical Implications and Motivations

Non-gravitational coupling between DM and DE profoundly modifies cosmic expansion dynamics. A key motivation is the **“coincidence problem”**—the unexplained near-equality of DM and DE densities today despite their different redshift scalings in $\Lambda$CDM. Many interaction models produce a scaling solution $r \equiv \rho_c/\rho_{\text{de}} = r_0 a^{-\xi}$, with $\xi < 3$ slowing the decay of $r$ and potentially yielding $r(a)\approx\text{constant}$ over cosmological timescales [1204.5892][1603.08299].

In models like $Q \propto H\rho_{\text{de}}$, the energy flow at late times can raise the Hubble rate inferred from CMB, bringing it closer to local $H_0$ measurements and thus partially alleviating the $H_0$ tension [2002.06127][2402.00819]. Conversely, $Q \propto H\rho_c$ tends to suppress the amplitude of structure growth at late times, addressing the $S_8$ or $\sigma_8$ tension between CMB and weak-lensing/galaxy-clustering data [2105.09249][2402.00819].

For specific interaction forms:
- Bilinear models $Q = \gamma\rho_c^\alpha\rho_{\text{de}}^\beta$ can be motivated by analogy with decay and annihilation rates in particle physics, and their parameter space can interpolate between various phenomenological scenarios [1911.04520].
- Field-theoretic constructions, including Chameleon-type and axion-monodromy models, can naturally produce late-time DM–DE coupling through the dependence of DM mass or couplings on scalar fields [1709.04882][1605.00996][2403.14247].

## 3. Linear Perturbations, Stability, and Structure Growth

DM–DE coupling modifies not only background dynamics but also the evolution of cosmic perturbations. In synchronous or Newtonian gauge, additional source and friction terms proportional to $Q$ appear in the equations for DM and DE density contrasts ($\delta_c$, $\delta_{\text{de}}$) and velocity divergences ($\theta_c$, $\theta_{\text{de}}$) [2002.06127][1911.04520][2105.09249]:
\[
\dot\delta_c = -\theta_c - \frac12\dot h + Q~\text{terms},\qquad
\dot\delta_{\text{de}} = -(1+w)\left(\theta_{\text{de}} + \frac12\dot h\right) + Q~\text{terms}.
\]
Correctly specifying the perturbative $Q^\mu$ is crucial for avoiding non-adiabatic and large-scale instabilities, especially at early times [1504.07397][1911.04520].

Phenomenologically, $Q>0$ (DE $\rightarrow$ DM) can suppress the late-time growth rate $f\sigma_8(z)$, aligning theory with the observed low clustering amplitude; $Q<0$ (DM $\rightarrow$ DE) may instead enhance growth unless tightly constrained. Best-fit values of parameters like $f\sigma_8(z)$ in interacting models typically show improved compatibility with redshift-space distortion and weak-lensing datasets [1911.04520][2105.09249][2010.10823].

## 4. Observational Constraints and Model Selection

Cosmological constraints on DM–DE interactions are derived from multi-probe analyses employing:
- CMB anisotropies (Planck 2015/2018, WMAP9)
- BAO (6dFGS, SDSS/BOSS/CMASS/DR14, DESI DR2)
- SNe Ia (Pantheon, Pantheon+, JLA, Union)
- Growth/RSD ($f\sigma_8$)
- Local $H_0$ direct measurements (SH0ES/R19)
- Cosmic chronometers, lookback time, lensing, and cluster counts

Representative model constraints include:
- For $Q = \xi H \rho_{\text{de}}$, joint CMB+R19+BAO+Pantheon fits favor $\xi = -0.179_{-0.074}^{+0.090}$ with $H_0 = 69.82_{-0.76}^{+0.63}$ km/s/Mpc, reducing the Hubble tension to $2.5\sigma$ [2002.06127].
- For $Q = \eta H \rho_{\text{dm}}$, $\eta = 0.00378 \pm 0.0012$ is obtained from CMB+BAO+SNIa, with the 1D PDF of cosmic void ellipticity shifting at the $\sim3\sigma$ level [2010.10823].
- Power-law interaction forms $Q \propto \rho_c^\alpha\rho_{\text{de}}^\beta$ yield best-fit $\alpha = 0.37_{-0.11}^{+0.25}$, $\beta = 1.19_{-1.19}^{+0.11}$ for $w < -1$ [1911.04520].
- Direct field-theoretic analyses (Chameleon, conformal, and disformal couplings) report best-fit $\beta \approx 0.26-0.34$ at $>2\sigma$ when DESI DR2 and SH0ES are included, with preference driven primarily by the late ISW effect [2403.14247].
- Current constraints generally require $|\xi| \lesssim 0.03$–$0.16$ at 95% C.L. for most phenomenological forms [2402.00819][1404.1620][1709.04882].

Statistical evidence for or against interaction is mixed. While some datasets or combinations can provide $\sim2$–$4\sigma$ preference for nonzero coupling when including local $H_0$ or specific weak-lensing datasets, combined global analyses using the full array of data typically yield inconclusive or weak (statistically insignificant) evidence, with Bayesian model comparison failing to favor interacting models over $\Lambda$CDM [2002.06127][2105.09249][2402.00819].

| Model/Class            | Key Coupling(s)         | Best-fit/Constraint        | Tension Addressed   |
|------------------------|------------------------|----------------------------|---------------------|
| $Q=\xi H\rho_{\text{de}}$   | $\xi$                 | $-0.179^{+0.090}_{-0.074}$ (all probes) | $H_0$              |
| $Q=\gamma\rho_c^\alpha\rho_{\text{de}}^\beta$ | $\alpha,\beta,\gamma$      | $\alpha\sim0.3$–$0.6$, $\beta\sim0$–$1$ (phantom) | $f\sigma_8$, $H_0$ |
| Chameleon/Conformal    | $\beta,\alpha$         | $\beta = 0.26^{+0.09}_{-0.06}$       | ISW anomaly, $w(z)$  |
| Power-law ($a^\beta$)  | $\xi_c,\beta$          | $\xi_c=0.34^{+0.14}_{-0.19}$, $\beta=-0.51^{+0.48}_{-0.41}$ | Early/late $\rho_c$  |

## 5. Nonlinear Regime, Simulations, and Physical Novelty

Large-scale N-body simulations for interacting models require careful modifications:
- Time-dependent DM mass evolution
- Drag forces from energy–momentum exchange
- Modified Poisson equation (fifth-force terms)
- Clustering signatures: changes in halo mass function, power spectrum, and void ellipticity

These simulations have shown:
- Modest changes in the high-mass end of the halo mass function for $|\xi| \sim 0.1$ [2402.00819]
- Matter power spectrum modifications with non-trivial scale dependence, requiring customized Halofit-like prescriptions or full N-body to match nonlinear observables

More speculative work has proposed that mutual DM–DE interactions may exhibit *chaotic attractor behavior* under general thermodynamic open-system principles, associated with the emergence of novel macroscopic order parameters and possibly observable fractal-like properties in cosmological structure [2304.06614].

## 6. Advanced Models and Fundamental Constraints

Field-theory-based models allow for a theoretically consistent embedding of DM–DE couplings:
- Chameleon models use conformal coupling of DM to a scalar field $\phi$, with action $S_\text{DM}[f(\phi)g_{\mu\nu},\psi]$, motivated by string compactifications. Current $\beta$ constraints are compatible with a fifth force exclusively in the dark sector [2403.14247].
- Monodromy, axion, and multi-scalar constructions provide radiatively stable mechanisms for mixing DE and ultra-light axion DM, with associated signatures in $w(z)$ oscillations and structure growth [1605.00996].
- “Holographic” and Ricci-scalar-cutoff DE models exploit connections between IR physics and the vacuum energy, introducing further covariant interaction possibilities [1206.1073][1009.6198][0901.1215].

First-principles kinetic theory, e.g., via Boltzmann equations for a Yukawa coupling between DM and DE scalars, predicts physically allowed microphysical $Q$ to be vastly smaller than phenomenological fits ($Q\sim10^{-100}$–$10^{-80}$ eV$^5$) [1909.10890], suggesting that empirically meaningful interactions are likely emergent or effective rather than fundamental.

## 7. Current Status, Limitations, and Future Prospects

Current cosmological data provide model-dependent upper bounds on DM–DE couplings at the $|\xi| \sim 10^{-2}$–$10^{-1}$ level, with no robust global detection across all datasets; some specific extensions or data combinations suggest mild preference but remain statistically inconclusive [2002.06127][1709.04882][2402.00819][2403.14247]. Constraints degrade rapidly if only background expansion is used—growth and lensing data are crucial.

Major open difficulties include:
- Maintaining stability and avoiding non-adiabatic instabilities at the perturbative level [1504.07397]
- Ambiguity in defining the energy–momentum transfer at both background and perturbation order, unless a Lagrangian approach is specified
- Non-trivial degeneracies with other cosmological parameters and the broad phenomenological flexibility permitted by current precision

Forthcoming probes will decisively test much of the allowed parameter space:
- ESA/Euclid, LSST, DESI for weak lensing, large-scale structure, and BAO
- CMB Stage-4 (lensing, ISW)
- 21cm intensity mapping, gravitational wave standard sirens, and fast radio bursts for geometry–growth combination

The minimum coupling viable with $\Lambda$CDM is expected to be probed at the $|\xi|\sim10^{-3}$ level, with nonlinear clustering statistics and cosmic void shapes providing independent constraints [2010.10823]. Laboratory and astrophysical searches for fifth forces in the dark sector may become relevant if evidence continues to mount for nonzero DM–DE interaction [2403.14247]. The ongoing search for theoretical consistency, data-driven phenomenology, and robust observable signatures remains central in the field.

---

**Key references:**  
[2002.06127], [2402.00819], [1911.04520], [1404.1620], [2010.10823], [1204.5892], [1506.02518], [2105.09249], [1709.04882], [1009.6198], [0901.1215], [2403.14247], [2304.06614], [1603.08299], [1504.07397], [1206.1073], [1605.00996], [1909.10890], [2410.02261].

Source: https://www.emergentmind.com/topics/dark-matter-dark-energy-interactions