Hint of dark matter-dark energy interaction in DESI DR2 and current cosmological dataset?

1. Problem Statement

The standard $\Lambda$CDM cosmological model faces persistent tensions, most notably the Hubble tension (discrepancy in $H_0$ measurements) and the $S_8$ tension (discrepancy in matter clustering). Furthermore, recent data from the Dark Energy Spectroscopic Instrument (DESI) suggests a possible departure from a pure cosmological constant, hinting at dynamical dark energy and a potential crossing of the "phantom divide" ($w < -1$) at redshifts $z \sim 0.5$.

Theoretical frameworks often propose interactions between Dark Matter (DM) and Dark Energy (DE) to resolve these issues. However, previous constraints on such interactions have typically yielded only upper bounds or required complex mechanisms (like sign-switching couplings) to fit data. This paper investigates a specific, theoretically motivated interaction derived from string theory compactifications (the Chameleon model) to see if current datasets, particularly the new DESI DR2 release, provide evidence for a non-zero coupling strength without ad-hoc modifications.

2. Methodology

Theoretical Model: Chameleon Dark Energy

The authors utilize a Chameleon model where a quintessential scalar field ($\phi$) couples to dark matter fermions via a Yukawa-type interaction:
$$f(\phi) = \exp\left(\frac{\beta \phi}{\sqrt{8\pi M_{Pl}}}\right)$$

• Mechanism: The scalar field evolves under an effective potential $V_{eff}(\phi)$, which is the sum of a self-interaction potential (runaway type, $V(\phi) \propto \phi^{-\alpha}$) and an interaction term sourced by the dark matter density ($\rho_{dm}$).
• Dynamics: In high-density environments, the field acquires a large mass (screening mechanism), suppressing fifth forces. In the low-density cosmic background, the field evolves.
• Key Assumption: The authors impose a boundary condition that the scalar field must settle at the minimum of the effective potential at the present epoch ($z=0$) with vanishing kinetic energy. This ensures the current dark energy density is purely potential-driven, leading to an effective equation of state $w_{eff} \to -1$ today.
• Phantom Behavior: While the true scalar field satisfies the Null Energy Condition ($w_\phi \geq -1$), an observer unaware of the interaction would infer an effective equation of state ($w_{eff}$) that can cross into the phantom regime ($w_{eff} < -1$) at earlier times due to the energy exchange between DM and DE.

Numerical Implementation

• Code Modification: The authors modified the Boltzmann code CLASS to include the modified background Friedmann equations and perturbation equations (in synchronous gauge) for the Chameleon model.
• Shooting Algorithm: A critical innovation is the use of a shooting algorithm to solve the Klein-Gordon equation. Instead of using approximate analytical solutions, the code iteratively adjusts the initial scalar field value ($\phi_i$) at high redshift to ensure the field evolves precisely to the analytically derived minimum ($\phi_0$) at $z=0$ with zero velocity. This ensures high numerical precision in the present-day energy density calculation.
• Perturbations: The code includes modified perturbation equations for the scalar field and dark matter, accounting for the coupling term $\beta$ in the evolution of density contrasts ($\delta_{dm}$) and metric perturbations.

Statistical Analysis

• Datasets: The study performs a Markov Chain Monte Carlo (MCMC) analysis using MontePython-v3 combined with the modified CLASS. The datasets include:
  • Planck 2018 (CMB TT, TE, EE, low-$\ell$).
  • SDSS BAO (BOSS, 6dFGS, DR7).
  • DESI DR2 BAO (New release with high-precision anisotropic constraints).
  • Pantheon+ (Type Ia Supernovae).
  • SH0ES (Local distance ladder/Cepheid variables).
• Parameters: The analysis constrains the standard 6 $\Lambda$CDM parameters plus the Chameleon parameters: the potential steepness $\alpha$ and the coupling strength $\beta$.
• Model Selection: The Akaike Information Criterion (AIC) is used to quantify the preference for the Chameleon model over $\Lambda$CDM, penalizing for the extra free parameter.

3. Key Contributions

1. Precise Numerical Solution: Unlike previous works relying on approximate analytical solutions, this paper employs a shooting algorithm to numerically solve the Klein-Gordon equation, ensuring the field strictly satisfies the present-day boundary conditions (minimum potential, zero velocity).
2. DESI DR2 Integration: This is one of the first comprehensive studies to incorporate the DESI DR2 BAO data into constraints on interacting dark sector models.
3. Sign-Stable Interaction: The model demonstrates that a non-zero coupling can be constrained without invoking "sign-switching" mechanisms (where energy flow reverses direction), relying solely on the natural dynamics of string-motivated compactification.
4. Phantom Crossing without Instability: The model naturally produces an effective phantom crossing ($w_{eff} < -1$) at $z \sim 0.5$ while maintaining theoretical stability (no ghost fields) and satisfying the Null Energy Condition for the fundamental field.

4. Results

Constraints on Coupling Strength ($\beta$)

The analysis reveals a mild preference for a non-zero coupling strength ($\beta > 0$), particularly when specific datasets are combined:

• Planck + SDSS BAO + Pantheon+: $\beta \approx 0.197$ (weak preference).
• Planck + DESI DR2 BAO + Pantheon+: $\beta \approx 0.256^{+0.092}_{-0.063}$. The inclusion of DESI DR2 increases the inferred interaction strength.
• Planck + SDSS BAO + Pantheon+ + SH0ES: $\beta \approx 0.335^{+0.081}_{-0.054}$. This combination yields the highest preference for interaction.

Statistical Evidence (AIC and $\Delta\chi^2$)

• $\Delta\chi^2_{min}$: The Chameleon model reduces the minimum chi-squared compared to $\Lambda$CDM.
  • With DESI DR2: $\Delta\chi^2_{min} = -4.75$.
  • With SH0ES: $\Delta\chi^2_{min} = -6.41$.
• AIC Analysis:
  • DESI DR2 combination: $\Delta AIC = -0.75$ (Weak evidence in favor).
  • SH0ES combination: $\Delta AIC = -2.41$ (Positive evidence in favor).
  • Note: Models with $\Delta AIC < 0$ are preferred; values between -2 and -6 indicate positive evidence.

Cosmological Implications

• Phantom Crossing: The effective equation of state $w_{eff}$ crosses the phantom divide ($w = -1$) at $z \sim 0.5$, aligning with the trend observed in DESI data.
• Hubble Constant: The model shifts the inferred $H_0$ slightly higher (e.g., $H_0 \sim 68.7$ with DESI, $69.4$ with SH0ES), helping to alleviate the Hubble tension, though not fully resolving it to the SH0ES value.
• $S_8$ Tension: The model predicts a slightly lower $S_8$ value compared to standard $\Lambda$CDM fits, which is consistent with the direction needed to resolve the $S_8$ tension, though the authors note that non-linear structure formation (requiring N-body simulations) is needed for a definitive conclusion.
• Limitation: Despite the phantom crossing at $z \sim 0.5$, the model predicts $w_{eff} \to -1$ at $z=0$ due to the field settling at the potential minimum. This conflicts with DESI's indication that $w_{de} > -1$ today, limiting the statistical significance of the model's fit to the DESI data compared to more flexible $w_0w_a$CDM models.

5. Significance and Future Outlook

• Theoretical Viability: The study validates that string-theory motivated interactions can naturally explain observational hints of dynamical dark energy and phantom crossing without fine-tuned sign-switching.
• DESI Impact: The results suggest that DESI DR2 data provides a "hint" of dark sector interaction, pushing the preferred coupling strength higher than previous constraints.
• Future Work: The authors emphasize that the current analysis is limited to linear scales. The non-zero coupling implies a "fifth force" that affects dark matter and baryons differently. Future work will utilize N-body simulations to study the impact on non-linear structure formation, halo profiles, and the differential infall of baryons vs. dark matter, which could provide stronger statistical evidence for or against this model.

In summary, the paper provides a robust numerical framework for testing interacting dark energy, finding that current data (especially DESI DR2 and SH0ES) mildly favors a non-zero interaction strength ($\beta \sim 0.3$) over the standard $\Lambda$CDM model, offering a potential pathway to resolve cosmological tensions.