A study of dark matter-dark energy interaction under the DESI DR2 data constraint

Technical Summary: A Study of Dark Matter-Dark Energy Interaction under the DESI DR2 Data Constraint

Problem Statement
The standard $\Lambda$CDM model, while successful in fitting much cosmological data, faces significant challenges, most notably the Hubble tension ($H_0$) and indications of a dynamical dark energy (DE) equation of state (EoS) from recent Dark Energy Spectroscopic Instrument (DESI) DR2 analyses. These observations suggest a deviation from the cosmological constant prediction of $w = -1$. Previous attempts to model Dark Matter-Dark Energy (DM-DE) interactions often rely on phenomenological continuity equations with source terms that lack a consistent field-theoretic derivation. This work addresses the need for a rigorous field-theoretic framework that naturally enforces energy conservation while investigating whether DM-DE interactions can resolve current cosmological tensions and align with DESI DR2 constraints.

Methodology
The authors propose and analyze a "QCDM" model, a field-theoretic cosmology comprising $n$ interacting spin-zero fields. Specifically, the study focuses on a two-field system: a DM field ($\chi$) and a DE field ($\phi$).

• Theoretical Framework: The action includes a gravitational sector, kinetic terms for the fields, and a potential $V$ consisting of individual potentials ($V_1(\chi)$ and $V_2(\phi)$) and an interaction term $V_{int}(\phi, \chi) = \frac{\lambda}{2}\chi^2\phi^2$. The DM potential is quadratic ($V_1 = \frac{1}{2}m_\chi^2\chi^2$), and the DE potential is a pseudo-Nambu-Goldstone boson (PNGB) form ($V_2 = \mu^4[1 + \cos(\phi/F)]$).
• Equations of Motion: The model utilizes the Klein-Gordon equations for the background fields and linear perturbations in a flat FLRW universe. The interaction term modifies the continuity equations, introducing a source term $Q_i$ that couples the evolution of the fields.
• Numerical Implementation: The model is implemented in the Boltzmann solver CLASS. The authors evolve background and perturbation equations from $a_{ini} = 10^{-14}$ to $a_0 = 1$. To handle rapid oscillations of the DM field, a redefinition of the DM energy density and new dimensionless variables are introduced.
• Data Analysis: The study confronts the model with current cosmological data, including Planck 2018 CMB measurements, DESI-DR2 BAO, and supernova datasets (PantheonPlus+SH0ES, Union3, and DESY5).
• Statistical Inference: Markov Chain Monte Carlo (MCMC) sampling is performed using Cobaya to constrain parameters ($\log \mu^4, F, \phi_{ini}, \log \lambda$) alongside standard $\Lambda$CDM parameters. The authors employ XGBoost and SHAP (SHapley Additive exPlanations) to identify which input parameters most strongly drive the model's predictions toward the $2\sigma$ region of DESI constraints on $(w_0, w_a)$. Model selection is evaluated using the Akaike Information Criterion (AIC).

Key Contributions and Results

1. Regime Classification: The study identifies two distinct regimes based on the interaction strength $\lambda$:

   • Strong Coupling ($\lambda \gtrsim 10^{-2}$): In this regime, the DM-DE interaction transforms the DE behavior from an early-time "thawing" phase to a late-time "scaling-freezing" phase. The authors demonstrate that this scaling-freezing behavior is incompatible with DESI DR2 data, which favors thawing or mild evolution. Consequently, the strong coupling regime is effectively ruled out.
   • Weak Coupling ($\lambda \lesssim 10^{-2}$): This regime yields an EoS evolution consistent with DESI data. The interaction induces a transition from thawing to scaling-freezing behavior, but for weak coupling, the late-time behavior remains close to $w \approx -1$ with mild evolution, avoiding the phantom divide crossing.

2. Constraints on Interaction Strength: Using SHAP analysis and MCMC fits, the authors derive upper limits on the interaction strength $\lambda$. The PantheonPlus dataset imposes the strongest constraint: $\log \lambda < -2.78$. Other combinations (e.g., with Union3 or DESY5) yield even tighter bounds ($\log \lambda < -3.26$).

3. Cosmological Parameter Constraints:

   • Equation of State: The best-fit points for the weak coupling regime fall in the fourth quadrant of the $(w_0, w_a)$ plane ($w_0 < 0, w_a < 0$), indicating a preference for evolving quintessence without phantom crossing. The derived values remain close to $\Lambda$CDM ($w_0 \simeq -1, w_a \simeq 0$) but allow for mild evolution.
   • Hubble Tension: The model provides only a modest alleviation of the $H_0$ tension at the 95% confidence level. The inferred $H_0$ values (e.g., $67.6 - 69.4$ km/s/Mpc depending on the dataset) do not fully resolve the discrepancy with local measurements ($H_0 \approx 73$ km/s/Mpc).
   • Structure Growth: The derived $S_8$ parameter is consistent with recent KiDS-Legacy results, showing no significant tension.

4. Model Comparison: The Akaike Information Criterion (AIC) analysis yields $\Delta \text{AIC} > 10$ for all data combinations, indicating that the current data decisively favors the standard $\Lambda$CDM model over the interacting QCDM model, despite the latter's ability to accommodate evolving dark energy.

Significance
The paper demonstrates that a fully field-theoretic description of DM-DE interaction, which automatically enforces energy conservation, can reproduce a dynamical dark energy equation of state consistent with DESI DR2 observations, provided the interaction is weak. While the model does not decisively resolve the Hubble tension or outperform $\Lambda$CDM in terms of statistical fit (AIC), it successfully rules out strong coupling scenarios that would lead to scaling-freezing behavior. The work establishes that weak DM-DE interactions are a viable, albeit constrained, extension to the standard model that permits evolving quintessence without invoking phenomenological source terms.