Exponential quintessence with momentum coupling to dark matter

We present updated constraints on an interacting dark energy - dark matter model with pure momentum transfer, where dark energy is in the form of a quintessence scalar field with an exponential potential. We run a suite of MCMC analyses using the DESI DR2 BAO measurements, in combination with CMB data from Planck and supernovae data from DESY5. In contrast to the standard case of uncoupled quintessence, we find that values for the potential’s slope parameter $λ\geq \sqrt{2}$, which are conjectured by string theory scenarios, are not excluded. If $λ$ is fixed to such a value, we find that the data favour the negative coupling branch of the model, which is the branch exhibiting late-time growth suppression. We also derive 95% upper limits on the sum of the neutrino masses, finding $\sum m_ν< 0.06$ eV ($\sum m_ν< 0.16$ eV) when $λ$ is fixed (varied). Our results motivate further studies on dynamical dark energy models that obey string theory bounds and can be constrained with cosmological observations.

💡 Research Summary

The authors investigate a dark‑energy–dark‑matter interaction model in which only momentum, not energy, is transferred between the two sectors. Dark energy is modeled as a quintessence scalar field ϕ with an exponential potential V(ϕ)=V₀ exp(−λϕ/Mₚₗ), a form that frequently appears in supergravity and string‑theoretic constructions. The interaction is introduced through a term β Z² in the Lagrangian, where Z≡u^μ∇_μϕ couples the field gradient to the dark‑matter four‑velocity. The dimensionless coupling constant β is allowed to be negative or positive, but theoretical considerations (avoidance of ghost instabilities) restrict β to the range −2 ≤ β < 0.5. The potential slope λ is of particular interest because string‑theory arguments often require λ ≥ √2.

The background cosmology retains the standard continuity equation for cold dark matter, while the Euler equation for the dark‑matter velocity acquires β‑dependent terms. Consequently, β > 0 enhances the growth of matter perturbations, whereas β < 0 suppresses it, offering a possible route to alleviate the current S₈ tension. The scalar‑field perturbation equation is also modified, but the Einstein equations remain unchanged.

To confront the model with data, the authors modify the CLASS Boltzmann solver and use the Cobaya framework for Markov‑Chain Monte‑Carlo sampling. The data sets comprise Planck 2018 temperature and polarization spectra (low‑ℓ, high‑ℓ TTTEEE, and lensing), the DESI Data Release 2 BAO measurements, and the DES‑Year 5 Type Ia supernova sample. Standard cosmological parameters (ω_b, ω_cdm, θ_s, A_s, n_s, τ_reio) are varied together with β, λ, and the sum of neutrino masses Σm_ν. Priors are flat: −2 ≤ β < 0.5, 0 ≤ λ ≤ 2.1, and Σm_ν ≥ 0 (with a fixed value of 0.06 eV in some runs).

The analysis proceeds in three stages. First, with Σm_ν fixed at 0.06 eV, the authors find that the data do not exclude λ ≥ √2 when β is free, in contrast to uncoupled quintessence (β = 0) where such large λ values are ruled out. The posterior shows a strong λ–β degeneracy; β remains essentially unconstrained, allowing the string‑theory motivated region of parameter space.

Second, fixing λ to 1.5 (≈√2) and allowing β to vary, the posterior shifts decisively toward negative β values (β ≈ −0.8 ± 0.9). This negative‑coupling branch predicts a suppression of late‑time growth, yielding a lower σ₈ (≈0.78) and a modestly higher H₀ (≈66.9 km s⁻¹ Mpc⁻¹) compared with the uncoupled case. The result suggests that, if the exponential potential slope is set by theoretical considerations, the data favor a momentum‑exchange that damps structure formation.

Third, the authors let Σm_ν vary. With λ fixed at 1.5, they obtain a 95 % upper limit Σm_ν < 0.06 eV, essentially at the lower bound from neutrino‑oscillation experiments. When λ is also allowed to vary, the neutrino mass bound relaxes to Σm_ν < 0.16 eV (95 % C.L.). The relaxation is driven primarily by the λ–β degeneracy: a larger λ can mimic the effect of massive neutrinos on the matter power spectrum, weakening the neutrino constraint. Nonetheless, both coupled and uncoupled quintessence models produce neutrino limits comparable to those from the w₀w_a parametrisation, though they do not reproduce the positive‑mass peak seen in the more flexible CPL analyses.

A χ² comparison shows that all dynamical dark‑energy models improve over ΛCDM (Δχ² ≈ −12 to −20), with the CPL (w₀w_a) model giving the largest improvement. The coupled quintessence with β < 0 also yields a modest reduction in σ₈ while slightly raising H₀, hinting at a simultaneous alleviation of the H₀ and S₈ tensions.

In summary, the paper demonstrates that an interacting dark‑energy model with pure momentum transfer can accommodate exponential potentials satisfying string‑theory bounds (λ ≥ √2) and remain fully compatible with current cosmological observations. The coupling parameter β is weakly constrained but shows a preference for negative values when λ is fixed to a theoretically motivated value, leading to late‑time growth suppression. Neutrino mass limits are sensitive to the λ–β degeneracy, tightening when λ is fixed and loosening when it is free. The authors conclude that future high‑precision large‑scale‑structure surveys (e.g., Euclid, LSST) and next‑generation CMB experiments (CMB‑S₄) will be able to break these degeneracies, providing sharper tests of momentum‑exchange interactions and their implications for fundamental physics.