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

Resumen Técnico: Un estudio de la interacción entre Materia Oscura y Energía Oscura bajo la restricción de los datos de DESI DR2

Planteamiento del Problema
El modelo estándar $\Lambda$CDM, aunque exitoso al ajustar gran parte de los datos cosmológicos, enfrenta desafíos significativos, notablemente la tensión de Hubble ($H_0$) e indicios de una ecuación de estado (EoS) de energía oscura (DE) dinámica de los análisis recientes del Dark Energy Spectroscopic Instrument (DESI) DR2. Estas observaciones sugieren una desviación de la predicción de la constante cosmológica de $w = -1$. Los intentos previos de modelar interacciones entre Materia Oscura y Energía Oscura (DM-DE) suelen basarse en ecuaciones de continuidad fenomenológicas con términos de fuente que carecen de una derivación de teoría de campos consistente. Este trabajo aborda la necesidad de un marco riguroso de teoría de campos que imponga naturalmente la conservación de la energía mientras se investiga si las interacciones DM-DE pueden resolver las tensiones cosmológicas actuales y alinearse con las restricciones de DESI DR2.

Metodología
Los autores proponen y analizan un modelo "QCDM", una cosmología de teoría de campos que comprende $n$ campos de espín cero que interactúan. Específicamente, el estudio se centra en un sistema de dos campos: un campo de materia oscura ($\chi$) y un campo de energía oscura ($\phi$).

• Marco Teórico: La acción incluye un sector gravitacional, términos cinéticos para los campos y un potencial $V$ que consiste en potenciales individuales ($V_1(\chi)$ y $V_2(\phi)$) y un término de interacción $V_{int}(\phi, \chi) = \frac{\lambda}{2}\chi^2\phi^2$. El potencial de la materia oscura es cuadrático ($V_1 = \frac{1}{2}m_\chi^2\chi^2$), y el potencial de la energía oscura es de la forma de un bosón de pseudo-Nambu-Goldstone (PNGB) ($V_2 = \mu^4[1 + \cos(\phi/F)]$).
• Ecuaciones de Movimiento: El modelo utiliza las ecuaciones de Klein-Gordon para los campos de fondo y perturbaciones lineales en un universo FLRW plano. El término de interacción modifica las ecuaciones de continuidad, introduciendo un término de fuente $Q_i$ que acopla la evolución de los campos.
• Implementación Numérica: El modelo se implementa en el solver de Boltzmann CLASS. Los autores evolucionan las ecuaciones de fondo y de perturbación desde $a_{ini} = 10^{-14}$ hasta $a_0 = 1$. Para manejar las oscilaciones rápidas del campo de materia oscura, se introduce una redefinición de la densidad de energía de la materia oscura y nuevas variables adimensionales.
• Análisis de Datos: El estudio confronta el modelo con datos cosmológicos actuales, incluyendo mediciones de CMB de Planck 2018, BAO de DESI-DR2 y conjuntos de datos de supernovas (PantheonPlus+SH0ES, Union3 y DESY5).
• Inferencia Estadística: Se realiza un muestreo de Monte Carlo por Cadenas de Markov (MCMC) utilizando Cobaya para restringir los parámetros ($\log \mu^4, F, \phi_{ini}, \log \lambda$) junto con los parámetros estándar de $\Lambda$CDM. Los autores emplean XGBoost y SHAP (explicaciones aditivas de Shapley) para identificar qué parámetros de entrada impulsan más fuertemente las predicciones del modelo hacia la región de $2\sigma$ de las restricciones de DESI sobre $(w_0, w_a)$. La selección del modelo se evalúa utilizando el Criterio de Información de Akaike (AIC).

Contribuciones Clave y Resultados

1. Clasificación de Regímenes: El estudio identifica dos regímenes distintos basados en la fuerza de interacción $\lambda$:

   • Acoplamiento Fuerte ($\lambda \gtrsim 10^{-2}$): En este régimen, la interacción DM-DE transforma el comportamiento de la DE de una fase temprana de "descongelamiento" (thawing) a una fase tardía de "escalamiento-congelamiento" (scaling-freezing). Los autores demuestran que este comportamiento de escalamiento-congelamiento es incompatible con los datos de DESI DR2, que favorecen un descongelamiento o una evolución suave. Consecuentemente, el régimen de acoplamiento fuerte queda efectivamente descartado.
   • Acoplamiento Débil ($\lambda \lesssim 10^{-2}$): Este régimen produce una evolución de la EoS consistente con los datos de DESI. La interacción induce una transición de un comportamiento de descongelamiento a uno de escalamiento-congelamiento, pero para un acoplamiento débil, el comportamiento tardío permanece cerca de $w \approx -1$ con una evolución suave, evitando el cruce de la división fantasma (phantom divide).

2. Restricciones sobre la Fuerza de Interacción: Utilizando el análisis SHAP y ajustes MCMC, los autores derivan límites superiores en la fuerza de interacción $\lambda$. El conjunto de datos PantheonPlus impone la restricción más fuerte: $\log \lambda < -2.78$. Otras combinaciones (por ejemplo, con Union3 o DESY5) arrojan límites aún más estrictos ($\log \lambda < -3.26$).

3. Restricciones de Parámetros Cosmológicos:

   • Ecuación de Estado: Los puntos de mejor ajuste para el régimen de acoplamiento débil caen en el cuarto cuadrante del plano $(w_0, w_a)$ ($w_0 < 0, w_a < 0$), lo que indica una preferencia por una quintaesencia evolutiva sin cruce fantasma. Los valores derivados permanecen cerca de $\Lambda$CDM ($w_0 \simeq -1, w_a \simeq 0$) pero permiten una evolución suave.
   • Tensión de Hubble: El modelo proporciona solo una mitigación modesta de la tensión de $H_0$ al nivel de confianza del 95%. Los valores de $H_0$ inferidos (por ejemplo, $67.6 - 69.4$ km/s/Mpc dependiendo del conjunto de datos) no resuelven completamente la discrepancia con las mediciones locales ($H_0 \approx 73$ km/s/Mpc).
   • Crecimiento de Estructuras: El parámetro $S_8$ derivado es consistente con los resultados recientes de KiDS-Legacy, sin mostrar una tensión significativa.

4. Comparación de Modelos: El análisis del Criterio de Información de Akaike (AIC) arroja un $\Delta \text{AIC} > 10$ para todas las combinaciones de datos, lo que indica que los datos actuales favorecen decisivamente al modelo estándar $\Lambda$CDM sobre el modelo QCDM interactuante, a pesar de la capacidad de este último para acomodar una energía oscura evolutiva.

Significado
El artículo demuestra que una descripción de teoría de campos completa de la interacción DM-DE, que impone automáticamente la conservación de la energía, puede reproducir una ecuación de estado de energía oscura dinámica consistente con las observaciones de DESI DR2, siempre que la interacción sea débil. Aunque el modelo no resuelve decisivamente la tensión de Hubble ni supera al modelo $\Lambda$CDM en términos de ajuste estadístico (AIC), logra descartar escenarios de acoplamiento fuerte que conducirían a un comportamiento de escalamiento-congelamiento. El trabajo establece que las interacciones débiles DM-DE son una extensión viable, aunque restringida, al modelo estándar que permite una quintaesencia evolutiva sin invocar términos de fuente fenomenológicos.