Exponential $f(R)$ cosmology with massive neutrinos as a dynamical dark energy framework

Résumé Technique : Cosmologie $f(R)$ exponentielle avec neutrinos massifs

Énoncé du Problème
Le modèle cosmologique standard $\Lambda$CDM, bien que fructueux, fait face à des défis théoriques et observationnels persistants. Les problèmes théoriques incluent les problèmes d'ajustement fin et de coïncidence associés à la constante cosmologique $\Lambda$, ainsi que l'absence d'une particule de matière noire confirmée. Sur le plan observationnel, des tensions significatives existent, notamment la « tension de Hubble », un écart de $\sim 5\sigma$ entre les mesures locales de la distance (échelle de distance) ($H_0 \approx 73,04$ km s$^{-1}$ Mpc$^{-1}$) et les inférences de l'Univers primordial issues du Fond Diffus Cosmologique (CMB) au sein du cadre $\Lambda$CDM ($H_0 \approx 67,36$ km s$^{-1}$ Mpc$^{-1}$). De plus, des données récentes sur les Oscillations Acoustiques des Baryons (BAO) provenant de la collaboration DESI suggèrent des déviations par rapport à $\Lambda$CDM, laissant présager une énergie noire dynamique.

Une tension secondaire concerne la masse totale des neutrinos ($\sum m_\nu$). Si les expériences d'oscillation des neutrinos établissent une limite inférieure de $\sum m_\nu \gtrsim 0,06$ eV, les données cosmologiques au sein du modèle $\Lambda$CDM produisent souvent des limites supérieures inférieures à ce seuil (par exemple, $\sum m_\nu < 0,064$ eV selon DESI). Les modèles phénoménologiques tels que $w_0w_a$CDM peuvent relaxer ces limites mais manquent de fondement théorique solide. Ce travail examine si une modification de la gravité, motivée théoriquement, spécifiquement la gravité $f(R)$ exponentielle, combinée aux neutrinos massifs, peut simultanément répondre à la tension de Hubble et aux contraintes de masse des neutrinos tout en maintenant la cohérence avec le modèle standard aux hautes redshifts.

Méthodologie
Les auteurs utilisent un modèle de gravité $f(R)$ exponentielle généralisé défini par le Lagrangien :
$$f(R) = R - 2\Lambda \left\{ 1 - \exp\left[ -\beta \left( \frac{R}{2\Lambda} \right)^\alpha \right] \right\}$$
où $\alpha$ et $\beta$ sont des paramètres sans dimension contrôlant la transition de la Relativité Générale (RG) à haute courbure vers la gravité modifiée à faible courbure. Le modèle inclut un secteur de matière composé de matière noire froide, de baryons, de rayonnement et de neutrinos massifs. La densité d'énergie des neutrinos est calculée via une intégrale d'espace des phases sur la distribution de Fermi–Dirac, tenant compte de la transition des régimes relativistes aux régimes non relativistes.

La dynamique cosmologique est dérivée des équations de Friedmann modifiées dans une métrique FLRW spatialement plate. Pour faciliter l'analyse numérique, les auteurs projettent les paramètres $f(R)$ vers des paramètres $\Lambda$CDM équivalents ($H_0^*, \Omega_i^*$) en utilisant un paramètre de Hubble sans dimension $E(z) = H(z)/H_0^*$. Le système d'équations différentielles couplées pour $E(z)$ et le scalaire de Ricci sans dimension $\mathcal{R}(z)$ est intégré numériquement à partir d'une condition initiale à haute redshift (où le modèle retrouve $\Lambda$CDM) jusqu'à $z=0$.

Une analyse de Monte Carlo par chaînes de Markov (MCMC) bayésienne est réalisée en utilisant le package emcee pour contraindre l'espace des paramètres $(H_0^*, \alpha, \beta, \omega_b^*, \omega_{bc}^*, \omega_\nu^*)$. L'analyse utilise quatre ensembles de données observationnelles :

1. Chronomètres Cosmiques (CC) : 31 points de données mesurant $H(z)$.
2. DESI DR2 BAO : Mesures de $D_M(z)/r_d$ et $D_H(z)/r_d$ à travers six compartiments de redshift.
3. CMB : Vraisemblances compressées de Planck 2018 pour l'échelle acoustique $\theta_*$ et les densités physiques $\omega_b, \omega_{bc}$.
4. Pantheon+ (Pan+) : 1550 modules de distance de Supernovae de type Ia.

Les résultats sont comparés au modèle standard $\Lambda$CDM et au modèle phénoménologique $w_0w_a$CDM en utilisant le critère d'information d'Akaike (AIC) et le critère d'information bayésien (BIC) pour évaluer la préférence statistique.

Contributions Clés et Résultats
L'étude fournit la première contrainte conjointe sur les paramètres du modèle $f(R)$ exponentiel généralisé ($\alpha, \beta$) et la masse totale des neutrinos ($\sum m_\nu$) en utilisant les dernières données DESI DR2.

• Viabilité du Modèle : Le modèle $f(R)$ exponentiel reste cohérent avec les observations actuelles. Les conditions de stabilité $f_R > 0$ et $f_{RR} > 0$ sont satisfaites sur l'ensemble de l'espace des paramètres a posteriori et de la plage de redshift, confirmant l'absence d'instabilités de type fantôme ou tachyonique.
• Tension de Hubble : Le modèle atténue partiellement la tension de Hubble. En incluant les supernovae Pantheon+, la valeur de $H_0$ inférée pour le modèle $f(R)$ est de $70,89 \pm 0,19$ km s$^{-1}$ Mpc$^{-1}$. Cette valeur réduit la tension avec les mesures locales par rapport à $\Lambda$CDM ($71,43 \pm 0,13$) et $w_0w_a$CDM ($71,13 \pm 1,20$), bien qu'un écart de $2,3\sigma$ subsiste entre les valeurs $f(R)$ et $\Lambda$CDM dans la combinaison complète des données.
• Contraintes sur la Masse des Neutrinos : L'inclusion de neutrinos massifs dans le cadre $f(R)$ conduit à une légère relaxation de la limite supérieure sur $\sum m_\nu$ par rapport à $\Lambda$CDM. Avec les données Pantheon+, la limite supérieure à $1\sigma$ est $\sum m_\nu < 0,0290$ eV pour $f(R)$, contre $< 0,0126$ eV pour $\Lambda$CDM et $< 0,0613$ eV pour $w_0w_a$CDM. Bien qu'il s'agisse d'une relaxation, la limite reste inférieure à la limite inférieure d'oscillation ($\sim 0,06$ eV), indiquant que la tension n'est pas totalement résolue.
• Préférence Statistique :
  • Sans les données Pantheon+, $\Lambda$CDM est statistiquement préféré.
  • Avec l'inclusion des données Pantheon+, $w_0w_a$CDM ($\Delta$AIC = 43,75) et $f(R)$ exponentiel ($\Delta$AIC = 20,46) montrent respectivement une preuve « très forte » et « forte » sur $\Lambda$CDM selon les métriques AIC et BIC.
• Dégénérescences de Paramètres : L'analyse révèle de fortes dégénérescences entre les paramètres de modification ($\alpha, \beta$) et les paramètres cosmologiques ($H_0, \Omega_m$), indiquant que la gravité modifiée peut partiellement compenser les changements de densité de matière et d'histoire de l'expansion. La distribution a posteriori pour $\beta$ est non gaussienne, reflétant le rôle non linéaire de la correction exponentielle.

Signification et Revendications
L'article affirme que l'interaction entre la gravité modifiée et la physique des neutrinos offre un cadre viable et théoriquement robuste pour investiguer les tensions cosmologiques. Contrairement aux paramétrisations phénoménologiques de l'énergie noire, le modèle $f(R)$ exponentiel fournit une base théorique solide ancrée dans les modifications du Lagrangien gravitationnel.

Les auteurs soulignent que, bien que le modèle ne résolve pas entièrement la tension de Hubble ou la discrépance de la masse des neutrinos, il démontre qu'une approche combinée — modifiant le secteur gravitationnel tout en incluant des neutrinos massifs — peut améliorer l'ajustement statistique aux données observationnelles et relaxer les contraintes sur la masse des neutrinos plus efficacement que $\Lambda$CDM seul. Ce travail souligne l'importance des jeux de données à bas redshift (comme les SNe Ia) pour discriminer entre les modèles alternatifs et la cosmologie standard. Les auteurs concluent que les investigations futures doivent inclure les effets perturbatifs (par exemple, le lentillage gravitationnel faible, les distorsions dans l'espace des redshifts) pour évaluer pleinement la viabilité du modèle, car l'analyse actuelle est restreinte à la cosmologie de fond.