CMBComp: A Simple and Accurate Compressed CMB Likelihood for Dark Energy, Curvature, and Massive Neutrinos

Résumé technique de « CMBComp : Un likelihood du CMB compressé, simple et précis pour l'énergie noire, la courbure et les neutrinos massifs »

Énoncé du problème
Le likelihood complet du Fond Diffus Cosmologique (CMB) est un outil extrêmement coûteux en termes de calcul et structurellement complexe, ce qui rend son déploiement difficile dans les analyses exploratoires, particulièrement pour les cosmologies phénoménologiques tardives. Cette limitation entrave les balayages rapides à travers de vastes espaces de modèles et les analyses conjointes avec des sondes géométriques à bas redshift, telles que les Oscillations Acoustiques des Baryons (BAO). Bien que l'information dominante du CMB pour l'inférence tardive soit conceptuellement comprise comme étant encodée dans un petit nombre de quantités géométriques (la taille angulaire de l'horizon sonore, la distance comobile jusqu'à la surface de dernière diffusion, et les paramètres de densité physique), un likelihood compressé, pratique, adapté aux modèles et précis, était nécessaire pour remplacer le likelihood complet sans sacrifier la précision.

Méthodologie
Les auteurs présentent CMBComp, un ensemble de likelihoods du CMB compressés, gaussiens et compacts, dérivés du jeu de données combiné SPA (SPT-3G D1 + ACT DR6 + anisotropies primaires Planck PR3 + lentillage CMB Planck PR4/NPIPE). La stratégie de compression repose sur le fait que pour les extensions du modèle standard lisses, pilotées par l'énergie noire et la géométrie, le CMB agit principalement comme une ancre géométrique à haut redshift.

Le cadre construit cinq likelihoods compressés distincts, chacun étant adapté à un espace de modèle cosmologique spécifique. Les observables géométriques centrales sont le paramètre de décalage $R$ et l'échelle angulaire acoustique $\ell_a$. Selon le modèle, ces paramètres sont augmentés d'un ou deux paramètres de densité physique pour former le vecteur de données $\mathbf{v}$. Le likelihood est défini par :
$$-2 \ln \mathcal{L}_{\text{comp}} = (\mathbf{v}_{\text{th}} - \bar{\mathbf{v}})^\top \mathbf{C}^{-1} (\mathbf{v}_{\text{th}} - \bar{\mathbf{v}})$$
où $\bar{\mathbf{v}}$ et $\mathbf{C}$ sont le vecteur moyen et la matrice de covariance extraits des chaînes SPA pour le modèle spécifique, et $\mathbf{v}_{\text{th}}$ est la prédiction théorique calculée via CAMB.

Les cinq compressions spécifiques sont :

1. CMB-3 : Une compression à trois paramètres $(R, \ell_a, \Omega_b h^2)$ pour le modèle $\Lambda$CDM spatialement plat.
2. CMB-3w : Le même vecteur à trois paramètres appliqué au modèle d'énergie noire dynamique $w_0 w_a$CDM, avec une moyenne et une covariance ré-extraites des chaînes $w_0 w_a$CDM pour tenir compte de l'expansion de l'espace des paramètres.
3. CMB-4$\nu$ : Une compression à quatre paramètres $(R, \ell_a, \Omega_b h^2, \Omega_c h^2)$ pour le modèle à neutrinos massifs $\nu\Lambda$CDM. L'inclusion de $\Omega_c h^2$ est nécessaire pour suivre correctement le bilan de la matière lorsque la somme des masses des neutrinos ($\sum m_\nu$) varie, évitant ainsi le traitement incorrect des neutrinos massifs comme de la matière froide lors de la recombinaison.
4. CMB-4k : Une compression à quatre paramètres $(R, \ell_a, \Omega_b h^2, \Omega_\Lambda)$ pour le modèle non plat $o\Lambda$CDM. Crucialement, pour ce modèle, le redshift de la dernière diffusion $z_*$ est évalué de manière auto-cohérente à partir de CAMB à chaque étape MCMC plutôt que d'être fixé, afin de rendre compte des décalages induits par la courbure dans la cartographie de l'horizon sonore.
5. CMB-4$\nu w$ : Une compression à quatre paramètres $(R, \ell_a, \Omega_b h^2, \Omega_c h^2)$ pour le modèle conjoint de neutrinos massifs et d'énergie noire dynamique $\nu w_0 w_a$CDM, avec des statistiques ré-extraites des chaînes conjointes correspondantes.

Principales contributions

• Likelihoods adaptés aux modèles : L'article fournit cinq own likelihoods compressés spécifiques, prêts à l'emploi (vecteurs moyens et matrices de covariance) dérivés de la génération la plus récente de données CMB (SPA).
• Indépendance d'implémentation : Les auteurs fournissent une implémentation Python autonome et légère (CMBComp) qui ne nécessite pas de solveur de Boltzmann complet (comme CAMB) ni de cadre d'inférence complexe (comme Cobaya ou CosmoMC) pour fonctionner. Elle inclut un calculateur d'évolution de fond personnalisé et un échantillonneur Metropolis-Hastings adaptatif.
• Validation : Les likelihoods compressés sont rigoureusement validés en les combinant avec le likelihood BAO de DESI DR2 et en comparant les posteriors marginalisés résultants aux chaînes de référence générées à l'aide du likelihood CMB SPA complet.

Résultats
La validation démontre que les likelihoods compressés reproduisent les contraintes marginalisées et les formes de contours 2D de l'analyse CMB complète avec une grande précision à travers les cinq classes de modèles :

• $\Lambda$CDM : Le posterior compressé montre un léger décalage résiduel ($\sim 0,6\sigma$) dans $H_0$ et $\Omega_m$ par rapport au CMB complet. Les auteurs attribuent cela à l'incapacité de la compression à 3 paramètres à modéliser pleinement les effets de lentillage du CMB, qui déplacent préférentiellement la constante de Hubble vers le bas dans le likelihood complet.
• $w_0 w_a$CDM : L'accord est excellent. Les paramètres supplémentaires d'énergie noire absorbent les biais observés dans le $\Lambda$CDM, entraînant des décalages négligeables ($< 0,3\sigma$) pour tous les paramètres.
• $\nu\Lambda$CDM : Les paramètres géométriques ($\Omega_m, H_0$) sont reproduits avec des erreurs $< 0,12\sigma$. La limite supérieure de la masse des neutrinos est légèrement plus faible ($\sim 13\%$) que la référence CMB complète, reflétant une réduction du levier sur $\sum m_\nu$ via la relation indirecte du bilan de la matière, mais la structure qualitative est préservée.
• $o\Lambda$CDM : L'accord est essentiellement exact, avec tous les décalages de moyenne inférieurs à $0,13\sigma$. Le calcul auto-cohérent de $z_*$ et l'inclusion de $\Omega_\Lambda$ éliminent la fuite géométrique observée dans le cas $\Lambda$CDM.
• $\nu w_0 w_a$CDM : Les paramètres d'énergie noire et géométriques concordent à $0,09\sigma$. La contrainte sur la masse des neutrinos est quelque peu plus large ($\sim 14\%$) que la référence en raison des directions de dégénérescence supplémentaires ouvertes par les paramètres d'énergie noire, mais la structure du posterior est fidèlement reproduite.

Signification et affirmations
L'article affirme que CMBComp est particulièrement bien adapté à l'inférence cosmologique dans les modèles qui modifient l'histoire de l'expansion tardive. Il permet d'incorporer des contraintes CMB précises dans de nouvelles analyses avec un surcoût de calcul minimal et sans dépendre d'un pipeline basé sur un solveur de Boltzmann complet.

Les auteurs soulignent que ces compressions sont valides pour les modèles où les effets sont confinés à l'histoire de l'expansion tardive et à l'évolution du fond, préservant ainsi la physique standard de l'Univers primordial. Ils avertissent explicitement que les modèles introduisant une nouvelle physique précoce, modifiant la recombinaison ou altérant de manière significative le secteur des perturbations primordiales, devraient toujours utiliser le likelihood CMB complet. Ce travail fournit un outil pratique et accessible aux chercheurs, étudiants et éducateurs pour réaliser une inférence cosmologique de haute précision à l'aide d'un cadre simplifié, transparent et informatiquement efficace.