Coupled Quintessence Dark Energy

• Coupled quintessence is an interacting dark energy model where a canonical scalar field exchanges energy with cold dark matter, altering the cosmic expansion.
• The models utilize various coupling functions—constant, field-dependent, or scale-factor parametrizations—to modify background dynamics and effective gravitational strength.
• They offer insights into the coincidence problem, phantom-divide crossing, and structure formation, with predictions tested through dynamical systems and observational constraints.

Coupled quintessence is a class of interacting dark-energy models in which the late-time acceleration of the Universe is driven by a canonical scalar field $\phi$, while cold dark matter (CDM) exchanges energy-momentum with that field through a non-gravitational interaction. In the standard Einstein-frame formulation, baryons and radiation are left uncoupled, the scalar sector has $\rho_\phi=\dot\phi^2/2+V(\phi)$ and $p_\phi=\dot\phi^2/2-V(\phi)$, and the dark sector obeys $\dot\rho_c+3H\rho_c=-Q$ and $\dot\rho_\phi+3H(\rho_\phi+p_\phi)=Q$, so total energy-momentum is conserved while the individual dark components are not (Potting et al., 2021, Wang et al., 20 Jun 2026). The model generalizes uncoupled quintessence by replacing separate conservation with a coupling term $Q$, or equivalently by allowing the dark-matter mass or effective metric to depend on $\phi$, and it has been studied as a framework for scaling behavior, coincidence-problem alleviation, effective phantom-divide crossing, and modified structure growth (Lopez-Honorez et al., 2010, Shojai et al., 2011).

1. Field-theoretic structure and background equations

A common formulation starts from Einstein gravity plus a canonical scalar with potential $V(\phi)$, with CDM coupled either through a field-dependent mass term or through a conformal factor in the CDM action. Representative actions include $S_c[\mathcal A^2(\phi)g_{\mu\nu};\psi_c]$ for conformally coupled CDM and $m_{dm}(\phi)\bar\Psi\Psi$ for field-dependent dark-matter mass, while baryons remain minimally coupled to $\rho_\phi=\dot\phi^2/2+V(\phi)$0 in order to evade local gravity constraints (Wang et al., 20 Jun 2026, Lopez-Honorez et al., 2010). In this sense, coupled quintessence is frequently described as an interacting dark-sector model rather than a modification of the visible sector.

At the homogeneous level, the coupling is encoded in a source term $\rho_\phi=\dot\phi^2/2+V(\phi)$1. For the generalized model with exponential scalar potential $\rho_\phi=\dot\phi^2/2+V(\phi)$2, the equations are

$\rho_\phi=\dot\phi^2/2+V(\phi)$3

$\rho_\phi=\dot\phi^2/2+V(\phi)$4

with interaction $\rho_\phi=\dot\phi^2/2+V(\phi)$5 in one widely used class (Potting et al., 2021). Closely related conformal models give

$\rho_\phi=\dot\phi^2/2+V(\phi)$6

with $\rho_\phi=\dot\phi^2/2+V(\phi)$7, so the usual CDM dilution law is explicitly modified (Wang et al., 20 Jun 2026). This change in the background matter evolution is one of the defining signatures of the framework.

Coupled quintessence also has a well-known scalar-tensor interpretation. In Brans–Dicke-inspired constructions, the Einstein-frame theory is equivalent to a quintessence field conformally coupled to dark matter, with interaction $\rho_\phi=\dot\phi^2/2+V(\phi)$8, unstable matter domination, a dark-energy-dominated point, and, for one branch, a mixed scaling solution relevant to the coincidence problem (Shojai et al., 2011). This establishes coupled quintessence as both a phenomenological dark-sector model and an Einstein-frame description of broader scalar-tensor theories.

2. Interaction functions and model families

The defining freedom of the subject lies in the choice of coupling function. The standard constant-coupling model corresponds to $\rho_\phi=\dot\phi^2/2+V(\phi)$9, giving $p_\phi=\dot\phi^2/2-V(\phi)$0, but many generalizations have been studied in order to alter the asymptotics, improve the matter era, or generate effective late-time dynamics beyond those of the constant-coupling case (Potting et al., 2021).

Several representative model classes appear repeatedly in the literature.

Interaction class	Representative form	Characteristic feature
Constant conformal coupling	$p_\phi=\dot\phi^2/2-V(\phi)$1	Standard coupled quintessence
Generalized field dependence	$p_\phi=\dot\phi^2/2-V(\phi)$2	Perturbed or unbounded coupling
Scale-factor parametrization	$p_\phi=\dot\phi^2/2-V(\phi)$3	Time-dependent phenomenological interaction
Dissipative coupling	$p_\phi=\dot\phi^2/2-V(\phi)$4	Warm-inflation-inspired energy transfer
Kinetic coupling	$p_\phi=\dot\phi^2/2-V(\phi)$5	Coupling through the scalar kinetic term
Curvature-modulated coupling	$p_\phi=\dot\phi^2/2-V(\phi)$6	Late-time activation by curvature suppression

In the generalized interaction model, two explicit choices are emphasized: a small perturbation around a constant,

$p_\phi=\dot\phi^2/2-V(\phi)$7

and an unbounded form,

$p_\phi=\dot\phi^2/2-V(\phi)$8

(Potting et al., 2021). In the first case, the familiar scaling accelerated solutions survive, although perturbed; in the second, true scaling solutions disappear. A distinct phenomenological route replaces the field dependence by a scale-factor function $p_\phi=\dot\phi^2/2-V(\phi)$9, with parameterizations $\dot\rho_c+3H\rho_c=-Q$0 and $\dot\rho_c+3H\rho_c=-Q$1, so that the interaction can vary explicitly over cosmic time (Costa, 2010).

Other constructions encode the interaction more directly in the CDM sector. In the “higher-order” model, the dark-matter mass is chosen as

$\dot\rho_c+3H\rho_c=-Q$2

which yields a coupling $\dot\rho_c+3H\rho_c=-Q$3 and a late-time accelerated attractor, while large $\dot\rho_c+3H\rho_c=-Q$4 suppresses the coupling and moves the theory toward uncoupled quintessence (Lopez-Honorez et al., 2010). In the linear-in-e-folds parametrization, $\dot\rho_c+3H\rho_c=-Q$5 with $\dot\rho_c+3H\rho_c=-Q$6, the coupling is a constant conformal coupling $\dot\rho_c+3H\rho_c=-Q$7, but the ansatz reconstructs an explicit potential

$\dot\rho_c+3H\rho_c=-Q$8

giving analytic control over the background evolution (Fonseca et al., 2021).

The literature has also extended coupled quintessence beyond energy exchange of the standard conformal type. Kinetic couplings link CDM to the scalar kinetic term $\dot\rho_c+3H\rho_c=-Q$9, momentum couplings modify the Euler equation without altering the CDM continuity equation, and scalar couplings to $\dot\rho_\phi+3H(\rho_\phi+p_\phi)=Q$0-form kinetic sectors generate anisotropic dark-energy dynamics rather than an isotropic perfect-fluid limit (2207.13682, Palma et al., 2023, Guarnizo et al., 2019). These variants retain the essential idea of a canonical quintessence field but enlarge the phenomenology of the dark sector.

3. Dynamical systems, scaling solutions, and cosmic sequence

A central question in coupled quintessence is whether the model yields a realistic cosmic sequence: a matter-dominated era long enough for structure formation, followed by late-time accelerated expansion. Dynamical-systems analyses have made this issue particularly precise. For the standard exponential-potential model with constant coupling, the usual dimensionless variables

$\dot\rho_\phi+3H(\rho_\phi+p_\phi)=Q$1

give the Friedmann constraint $\dot\rho_\phi+3H(\rho_\phi+p_\phi)=Q$2 and effective equation of state $\dot\rho_\phi+3H(\rho_\phi+p_\phi)=Q$3 (Potting et al., 2021). In that standard case, the system admits scaling critical points, including an accelerating one, but the key drawback recalled in the literature is that these accelerated scaling solutions are not preceded by a sufficiently long matter-dominated era.

For generalized couplings, the phase-space structure changes qualitatively. When $\dot\rho_\phi+3H(\rho_\phi+p_\phi)=Q$4, the compactified variable $\dot\rho_\phi+3H(\rho_\phi+p_\phi)=Q$5 is introduced, and the autonomous system contains a matter point $\dot\rho_\phi+3H(\rho_\phi+p_\phi)=Q$6 with $\dot\rho_\phi+3H(\rho_\phi+p_\phi)=Q$7, $\dot\rho_\phi+3H(\rho_\phi+p_\phi)=Q$8, $\dot\rho_\phi+3H(\rho_\phi+p_\phi)=Q$9, together with critical lines

$Q$0

which are scalar-field dominated (Potting et al., 2021). The striking result is that there are no true scaling solutions for this unbounded coupling. For $Q$1, the late-time attractor is a dark-energy-dominated state on the $Q$2 line, accelerating when $Q$3, with

$Q$4

For $Q$5, the attractor becomes $Q$6, a pure dark-energy state with $Q$7, but the approach is extremely slow. During that long transient, $Q$8, so the evolution mimics an accelerated scaling solution without possessing a true scaling fixed point.

A similar asymptotic structure appears in the warm-inflation-inspired model with $Q$9. There the autonomous system has the standard matter point $\phi$0, an unstable matter-like scaling point $\phi$1, and new critical lines at $\phi$2, but the attractors lie on those lines rather than on a stable mixed scaling point (Sá, 2023). For $\phi$3, the system can spend a long time in a scaling-type regime with appreciable $\phi$4, yet asymptotically it still approaches dark-energy domination. This reinforces a recurring conclusion in the subject: long-lived quasi-scaling behavior need not correspond to a true asymptotic scaling attractor.

Other parametrizations recover scaling more explicitly. In the linear-in-e-folds model, the field tracks radiation at early times, with $\phi$5, enters a matter-era scaling attractor with

$\phi$6

and then freezes at late times as the shallower exponential term dominates (Fonseca et al., 2021). This suggests that coupled quintessence can realize the radiation-to-matter-to-acceleration sequence either through true scaling phases or through slow transient regimes, depending on the interaction function.

4. Effective equation of state and phantom-divide crossing

The intrinsic scalar field in coupled quintessence is canonical, so its equation of state obeys $\phi$7. Nevertheless, several coupled models allow the effective dark-energy equation of state inferred from the background expansion to cross the phantom divide $\phi$8 without introducing a wrong-sign kinetic term (Chakraborty et al., 13 Mar 2025, Wang et al., 20 Jun 2026). This distinction between the physical scalar equation of state and the observationally reconstructed effective equation of state is now one of the most important conceptual points in the subject.

In Yukawa-type coupled models, dark matter has a field-dependent mass through

$\phi$9

so that $V(\phi)$0 rather than $V(\phi)$1 (Chakraborty et al., 13 Mar 2025). In that framework,

$V(\phi)$2

If $V(\phi)$3, then $V(\phi)$4 is shifted below $V(\phi)$5; thus $V(\phi)$6 can cross $V(\phi)$7 even though $V(\phi)$8 at all times. With suitable parameter choices, the model crosses the phantom divide around $V(\phi)$9 to $S_c[\mathcal A^2(\phi)g_{\mu\nu};\psi_c]$0, in agreement with the DESI-motivated trend discussed in that work.

The sign of the coupling can also change dynamically. In the exponential-potential, exponential-coupling CQ-EXP model, the branch $S_c[\mathcal A^2(\phi)g_{\mu\nu};\psi_c]$1 allows the interaction force to oppose the potential force at early times, so the scalar can be driven up its potential and later roll back down, making $S_c[\mathcal A^2(\phi)g_{\mu\nu};\psi_c]$2 change sign (Wang et al., 20 Jun 2026). Since $S_c[\mathcal A^2(\phi)g_{\mu\nu};\psi_c]$3, the energy transfer reverses direction during evolution, and the effective dark-energy equation of state can cross $S_c[\mathcal A^2(\phi)g_{\mu\nu};\psi_c]$4. By contrast, for $S_c[\mathcal A^2(\phi)g_{\mu\nu};\psi_c]$5, the scalar rolls monotonically, $S_c[\mathcal A^2(\phi)g_{\mu\nu};\psi_c]$6 throughout, and the effective equation of state remains above $S_c[\mathcal A^2(\phi)g_{\mu\nu};\psi_c]$7.

A different realization of sign switching appears in the $S_c[\mathcal A^2(\phi)g_{\mu\nu};\psi_c]$8CDM model with

$S_c[\mathcal A^2(\phi)g_{\mu\nu};\psi_c]$9

Because $m_{dm}(\phi)\bar\Psi\Psi$0, the coupling term changes sign when $m_{dm}(\phi)\bar\Psi\Psi$1 crosses zero, generating an early phase in which energy flows from dark matter to dark energy and a later phase in which it flows from dark energy to dark matter (Wang et al., 2 Apr 2026). The scalar itself remains canonical, with $m_{dm}(\phi)\bar\Psi\Psi$2, but the modified CDM evolution produces an effective $m_{dm}(\phi)\bar\Psi\Psi$3 that can cross the phantom divide.

SUGRA-based coupled quintessence yields an analogous mechanism. With

$m_{dm}(\phi)\bar\Psi\Psi$4

the field can begin near the potential minimum; coupling to CDM then displaces it, induces a turnaround, changes the sign of $m_{dm}(\phi)\bar\Psi\Psi$5, and produces $m_{dm}(\phi)\bar\Psi\Psi$6 crossing from below to above $m_{dm}(\phi)\bar\Psi\Psi$7 (Wang et al., 18 May 2026). This shows that phantom-divide crossing in coupled quintessence is not a unique consequence of one coupling ansatz, but a general effect of interaction-induced modification of the inferred dark-energy sector.

Not all coupled models predict eternal acceleration. In the phenomenological $m_{dm}(\phi)\bar\Psi\Psi$8 class, positive late-time coupling can drain enough energy from dark energy that the present accelerated epoch is only temporary, after which the Universe re-enters a future CDM-dominated era (Costa, 2010). This is qualitatively different from standard $m_{dm}(\phi)\bar\Psi\Psi$9CDM and from the more common coupled-quintessence asymptotics in which dark-energy domination persists.

5. Perturbations, fifth forces, and structure formation

Coupled quintessence affects not only the homogeneous background but also linear and nonlinear structure formation. In the simplest conformal models, the matter-density contrast equation acquires both modified friction and modified source terms. One representative form is

$\rho_\phi=\dot\phi^2/2+V(\phi)$00

which makes explicit that the coupling changes both the damping term and the effective gravitational strength (Patil et al., 2023). In related treatments, the fifth-force interpretation is summarized by

$\rho_\phi=\dot\phi^2/2+V(\phi)$01

so the coupling enhances or suppresses CDM growth depending on the sign and the background scalar evolution (Tarrant et al., 2011).

Kinetic couplings alter perturbations in a more distinctive way because the interaction depends on the scalar kinetic term rather than on $\rho_\phi=\dot\phi^2/2+V(\phi)$02 itself. For $\rho_\phi=\dot\phi^2/2+V(\phi)$03, the model yields early-time scaling behavior, modified CDM continuity and Euler equations, and effective friction and gravitational couplings in the sub-horizon growth equations (2207.13682). The paper reports measurable signatures in the CMB TT spectrum, the lensing potential auto-correlation, and the matter power spectrum; observationally, however, the coupling parameter is constrained to be of order $\rho_\phi=\dot\phi^2/2+V(\phi)$04, and model selection does not reveal a statistical preference over $\rho_\phi=\dot\phi^2/2+V(\phi)$05CDM.

At the nonlinear level, coupled quintessence modifies halo abundance and halo structure. Semi-analytic halo-mass-function calculations show that the effect depends strongly on the potential and on the direction of energy transfer. For the double-exponential potential, a weak positive coupling can enhance the abundance of very massive high-redshift halos, while Albrecht–Skordis and SUGRA models generally predict fewer such objects than $\rho_\phi=\dot\phi^2/2+V(\phi)$06CDM (Tarrant et al., 2011). The same study compares spherical-collapse and Newtonian prescriptions for the collapse threshold $\rho_\phi=\dot\phi^2/2+V(\phi)$07 and finds that they can differ substantially, especially at high redshift and strong coupling, although this uncertainty is secondary to parameter uncertainties and fitting-function uncertainties.

Momentum-exchange models provide yet another perturbative channel. In the Type-3 momentum coupling, the CDM continuity equation remains standard, but the Euler equation becomes

$\rho_\phi=\dot\phi^2/2+V(\phi)$08

so the coupling acts through a modified friction term and a modified effective gravitational force (Palma et al., 2023). N-body simulations of this scenario find that, for certain quintessence potentials, a positive coupling can reduce small-scale structure, enhance somewhat the large-scale power, lower halo central densities, and increase halo velocity dispersions. This marks an important distinction between energy-exchange couplings and pure momentum-exchange couplings: the former alter the background matter density itself, whereas the latter leave the continuity equation unchanged and act directly on the dark-matter trajectories.

The literature has also extended coupled-quintessence phenomenology to galactic dynamical friction. In the uncoupled limit, standard quintessence does not generate dynamical friction at galactic scales, but once the scalar couples nonminimally to matter, the scalar-mediated fifth force produces a modified wake and a nonzero dynamical-friction force (Nari et al., 15 Feb 2025). The Fornax application in that analysis shows that galactic-scale fits can imply significant departures from standard CDM, whereas cosmological limits on the coupling make the effect negligible. This illustrates a broader tension between local and cosmological constraints on dark-sector interactions.

6. Observational constraints, degeneracies, and current status

A persistent result of observational work is that coupled quintessence is difficult to infer uniquely from background data because the dark-energy equation of state and the interaction are degenerate in the expansion history. A model-independent reconstruction based on Gaussian processes shows that, for a fixed $\rho_\phi=\dot\phi^2/2+V(\phi)$09, one cannot independently determine both $\rho_\phi=\dot\phi^2/2+V(\phi)$10 and the interaction $\rho_\phi=\dot\phi^2/2+V(\phi)$11 from background data alone; a prior on one is needed to reconstruct the other (Yang, 2020). Using cosmic chronometers, BAO, Type Ia supernovae, and Ly$\rho_\phi=\dot\phi^2/2+V(\phi)$12 BAO, that analysis finds only a mild $\rho_\phi=\dot\phi^2/2+V(\phi)$13 hint of nonzero interaction at low redshift and a $\rho_\phi=\dot\phi^2/2+V(\phi)$14 hint at high redshift, with the latter driven mainly by the Ly$\rho_\phi=\dot\phi^2/2+V(\phi)$15 BAO point at $\rho_\phi=\dot\phi^2/2+V(\phi)$16. The same work notes that allowing coupling weakens the reconstructed Swampland bound $\rho_\phi=\dot\phi^2/2+V(\phi)$17 at $\rho_\phi=\dot\phi^2/2+V(\phi)$18 from $\rho_\phi=\dot\phi^2/2+V(\phi)$19 in the uncoupled case to $\rho_\phi=\dot\phi^2/2+V(\phi)$20 at $\rho_\phi=\dot\phi^2/2+V(\phi)$21 CL.

Fits to specific coupled-quintessence parametrizations generally favor small couplings, but the preferred sign can be dataset dependent. In the linear-in-e-folds model, Planck 2018 slightly prefers $\rho_\phi=\dot\phi^2/2+V(\phi)$22, corresponding to energy transfer from dark energy to dark matter, whereas KiDS-450 slightly prefers $\rho_\phi=\dot\phi^2/2+V(\phi)$23, corresponding to the opposite direction (Fonseca et al., 2021). In a broader 2023 observational analysis with CMB distance priors, BAO, Pantheon+, cosmic chronometers, megamasers, growth data, and SH0ES, the preferred interaction remains very small, $\rho_\phi=\dot\phi^2/2+V(\phi)$24, and the model behaves close to uncoupled quintessence in practice; AIC and BIC favor the coupled model over $\rho_\phi=\dot\phi^2/2+V(\phi)$25CDM, but the reduced-$\rho_\phi=\dot\phi^2/2+V(\phi)$26 criterion favors the simpler model (Patil et al., 2023). This divergence between information criteria and reduced-$\rho_\phi=\dot\phi^2/2+V(\phi)$27 is an early example of a recurring model-selection ambiguity in the literature.

DESI-era analyses have sharpened attention on effective phantom-divide crossing. A 2025 Yukawa-type study shows that a canonical coupled scalar with $\rho_\phi=\dot\phi^2/2+V(\phi)$28 can still reproduce the DESI-preferred trend in $\rho_\phi=\dot\phi^2/2+V(\phi)$29, typically crossing near $\rho_\phi=\dot\phi^2/2+V(\phi)$30 to $\rho_\phi=\dot\phi^2/2+V(\phi)$31, with several parameter choices lying inside DESI or DESI+Union3+Planck $\rho_\phi=\dot\phi^2/2+V(\phi)$32 contours (Chakraborty et al., 13 Mar 2025). A 2026 CQ-EXP likelihood analysis with Planck CMB, DESI BAO, and DES-Dovekie supernovae reports that, when $\rho_\phi=\dot\phi^2/2+V(\phi)$33 eV is fixed, both $\rho_\phi=\dot\phi^2/2+V(\phi)$34 and $\rho_\phi=\dot\phi^2/2+V(\phi)$35 branches prefer nonzero coupling, with the $\rho_\phi=\dot\phi^2/2+V(\phi)$36 branch favored more strongly and separated from zero by more than $\rho_\phi=\dot\phi^2/2+V(\phi)$37; however, allowing a free signed effective neutrino-mass parameter weakens the coupling preference and makes the $\rho_\phi=\dot\phi^2/2+V(\phi)$38 and $\rho_\phi=\dot\phi^2/2+V(\phi)$39 branches nearly identical in $\rho_\phi=\dot\phi^2/2+V(\phi)$40, leaving CQ-EXP statistically indistinguishable from $\rho_\phi=\dot\phi^2/2+V(\phi)$41CDM with the data considered there (Wang et al., 20 Jun 2026). This directly exhibits the degeneracy between dark-sector coupling and late-time neutrino-sector freedom.

The strongest DESI-era claim in the supplied literature comes from coupled SUGRA quintessence. Using DESI BAO, DES-Dovekie SNIa, and Planck CMB data, that analysis finds the coupling parameter differs from zero by more than $\rho_\phi=\dot\phi^2/2+V(\phi)$42, and the preferred branch is the one in which the interaction changes sign and the effective equation of state crosses the phantom divide (Wang et al., 18 May 2026). Yet even there the best-fit coupled SUGRA branch is described as statistically indistinguishable from the CPL parametrization, with only a very small difference in $\rho_\phi=\dot\phi^2/2+V(\phi)$43. A plausible implication is that present data may be more sensitive to the broad phenomenology of evolving dark energy than to the specific microphysical origin of that evolution.

Model-selection controversies remain active. In the 2025 Planck+DESI DR2 study of coupled quintessence with an additional warm-dark-matter sector, Deviance Information Criterion strongly prefers the quintessence model coupled with warm dark matter, while Bayesian evidence favors $\rho_\phi=\dot\phi^2/2+V(\phi)$44CDM for all dataset combinations considered (Samanta et al., 11 Sep 2025). The literature therefore supports two simultaneous conclusions: coupled quintessence is observationally viable and can reproduce several late-time signatures now under discussion, but the statistical evidence for a nonzero dark-sector coupling is highly sensitive to parametrization, auxiliary sector freedom, and the model-selection metric itself.