Bound or blown: the fate of hot gas in galaxy groups

Résumé technique : Lié ou expulsé : le sort du gaz chaud dans les groupes de galaxies

Problème et motivation
L'état thermodynamique du gaz chaud dans les groupes de galaxies demeure une incertitude critique pour la compréhension de la formation des galaxies et de la structure à grande échelle (LSS) de l'Univers. Bien que les simulations hydrodynamiques cosmologiques aient progressé de manière significative dans la modélisation de la formation stellaire et du feedback des noyaux actifs de galaxies (AGN), l'impact précis du feedback sur le milieu intra-groupe (IGrM) fait toujours l'objet de débats. Les observations suggèrent que le gaz dans les groupes est moins dense et plus étendu que ce que prédit un simple effondrement gravitationnel, impliquant une expulsion importante de baryons. Cependant, des contraintes récentes issues de mesures de l'effet Sunyaev-Zel'dovich cinétique (kSZ) favorisent des modèles de feedback extrêmement éjectifs, tandis que d'autres analyses, incluant les observations d'XRISM, suggèrent des dispersions de vitesse plus faibles qui pourraient indiquer des modèles trop éjectifs. Il est nécessaire de comparer directement les prédictions des simulations avec les données observationnelles tout en tenant compte rigoureusement des effets de sélection et des biais de mesure afin de déterminer la force correcte du feedback des AGN dans les halos de faible masse ($M_{500c} \sim 10^{13}–10^{14} M_\odot$).

Méthodologie
Cette étude emploie une approche de modélisation directe (forward-modelling) pour comparer l'échantillon de l'XMM-Newton Group AGN Project (X-GAP) à la suite de simulations hydrodynamiques cosmologiques FLAMINGO. La méthodologie se déroule selon plusieurs étapes clés :

1. Sélection de l'échantillon et analogues : Les auteurs utilisent la suite de simulations FLAMINGO, qui inclut diverses prescriptions de feedback des AGN allant d'un faible à un extrême éjection (paramétrées par des décalages dans la calibration de la fraction de gaz, $f_{gas} \pm N\sigma$, et des modèles de jets). Pour créer des analogues réalistes de l'échantillon X-GAP, les auteurs appliquent la fonction de sélection de X-GAP (basée sur le flux X, le décalage vers le rouge et la dispersion de vitesse) aux données de simulation. Cette sélection tient compte de la complétude et de la pureté spécifiques de l'échantillon X-GAP dérivé du catalogue AXES.
2. Modélisation directe des observations : Une innovation cruciale est la génération d'observations fictives (mock observations) de bout en bout pour XMM-Newton. Partant de la distribution 3D des particules de gaz chaud dans les simulations, les auteurs :
   • Génèrent des listes de photons idéalisées à l'aide du logiciel pyXSIM, en incorporant les températures, les métallicités et les abondances individuelles des particules.
   • Simulent les effets instrumentaux réalistes (Fonction d'étalement du point [PSF], variations de l'aire effective et fond diffus) en utilisant le package xmm_simulator pour produire des listes d'événements fictives.
   • Analysent ces données fictives en utilisant exactement les mêmes pipelines d'imagerie et d'analyse spectrale (incluant le masquage des sources ponctuelles, la modélisation du fond et l'ajustement spectral) appliqués aux données réelles de X-GAP.
3. Extraction des observables : À partir des données fictives analysées, les auteurs dérivent les luminosités X (noyau exclu) ($L_X$), les masses de gaz ($M_{gas}$), les températures ($T$) et les dispersions de vitesse ($\sigma_v$). Ils tiennent également compte de la variance cosmique en générant plusieurs cônes de lumière indépendants et les incertitudes de la fonction de sélection.
4. Comparaison statistique : L'étude compare les distributions de ces observables entre les données X-GAP et les différents modèles FLAMINGO. Crucialement, les auteurs utilisent une approche de copule gaussienne pour combiner plusieurs observables (normalisations des relations d'échelle, température moyenne, comptage des groupes et dispersion de vitesse) tout en tenant compte de leur covariance, plutôt que de les traiter comme des contraintes indépendantes.

Principales contributions

• Modélisation directe de bout en bout : L'article établit un cadre rigoureux pour comparer les simulations aux observations en générant des données fictives XMM-Newton qui incluent les effets instrumentaux et en les analysant de manière identique aux données réelles. Cela élimine les biais découlant de la comparaison directe de quantités de simulation « vraies » à des quantités « observées ».
• Contraintes multi-observables : En combinant la normalisation des relations d'échelle $L_X-T$ et $M_{gas}-T$ avec le comptage des groupes, la température moyenne et la dispersion de vitesse, l'étude fournit un discriminateur plus robuste pour les modèles de feedback qu'une seule observable isolée.
• Quantification des biais systématiques : Le travail quantifie explicitement l'impact de la variance cosmique et des incertitudes de la fonction de sélection, démontant que le comptage des groupes seul est un faible discriminateur en raison de fluctuations importantes (>20 %), tandis que les propriétés thermodynamiques offrent des contraintes plus serrées.

Résultats
L'analyse produit les conclusions suivantes concernant la force du feedback des AGN dans les groupes de galaxies :

• Modèle le mieux ajusté : Le modèle FLAMINGO avec une force de feedback intermédiaire ($f_{gas} - 2\sigma$) offre le meilleur accord avec les données X-GAP, présentant une tension globale de seulement $0,8\sigma$. Ce modèle produit une fraction de gaz typique d'environ $6,1\%$ pour les groupes sélectionnés.
• Modèles exclus : Le scénario de feedback éjectif le plus extrême ($f_{gas} - 8\sigma$) est exclu à plus de $4\sigma$ (spécifiquement $4,4\sigma$). Ce modèle prédit des fractions de gaz aussi basses que $\sim 3,8\%$ et est incompatible avec les relations d'échelle et les distributions de température observées.
• Fiduciaire vs Données : La calibration fiduciaire de FLAMINGO ($L1\_m8$) présente une tension de $3,1\sigma$, indiquant que la force de feedback standard est probablement trop faible pour reproduire le contenu de gaz chaud observé dans les groupes.
• Relations d'échelle : La normalisation des relations $L_X-T$ et $M_{gas}-T$ sert de principal discriminateur. Un feedback plus fort déplace les systèmes vers des luminosités plus faibles et des températures plus élevées à masse fixe, les éloignant de la distribution X-GAP.
• Variance cosmique : Le nombre de groupes détectés est fortement influencé par la variance cosmique, ce qui en fait un critère peu performant pour distinguer les modèles de feedback dans l'échantillon X-GAP à faible redshift et faible volume.

Signification et affirmations
L'article affirme que les propriétés thermodynamiques des groupes de galaxies favorisent une force de feedback des AGN supérieure à la calibration fiduciaire de FLAMINGO, mais nettement plus faible que les modèles les plus éjectifs actuellement proposés. Ce résultat souligne la nécessité de combiner la modélisation directe avec des contraintes multi-observables pour sonder avec précision le sort des baryons chauds dans les halos de faible masse.

Les auteurs notent que leurs conclusions sont en tension avec les mesures récentes du kSZ qui favorisent les modèles les plus éjectifs ($f_{gas} - 8\sigma$). Ils suggèrent que cette divergence peut provenir de différences dans les plages de masse et de redshift sondées, des effets de sélection, ou des échelles spécifiques auxquelles le feedback opère. L'étude conclut que, bien que les simulations actuelles puissent reproduire la fonction de masse stellaire des galaxies, leurs prédictions pour le contenu de gaz chaud dans les groupes de galaxies restent sensibles aux détails de l'implémentation du feedback, et que le « sort » du gaz chaud dans ces systèmes est probablement intermédiaire entre les extrêmes testés.