Cosmology with one galaxy: An analytic formula relating $Ω_{\rm m}$ with galaxy properties

Résumé technique : Cosmologie avec une seule galaxie

Énoncé du problème
Les analyses cosmologiques standard découpent généralement la formation des galaxies de l'inférence des paramètres cosmologiques, en s'appuyant sur des statistiques au niveau de la population telles que le regroupement, le lentillage gravitationnel ou les abondances d'halos pour contraindre le paramètre de densité de matière, $\Omega_m$. Bien que des études récentes utilisant l'apprentissage automatique sur la suite de simulations CAMELS aient démontré que les propriétés individuelles des galaxies peuvent encoder des informations cosmologiques, l'origine physique de ce signal est restée obscure, et les approches précédentes reposaient sur des réseaux de neurones « boîte noire » manquant d'interprétabilité physique. Ce travail aborde le défi de dériver une relation analytique physiquement transparente reliant directement $\Omega_m$ à des observables intrinsèques de galaxies uniques, évitant ainsi le besoin de statistiques de regroupement à grande échelle.

Méthodologie
Les auteurs utilisent le projet CAMELS (Cosmology and Astrophysics with MachinE Learning Simulations), spécifiquement les ensembles de données Latin Hypercube (LH) et l'ensemble de données étendu TNG-SB28. Ces suites contiennent plus de 16 000 simulations hydrodynamiques faisant varier les paramètres cosmologiques ($\Omega_m$, $\sigma_8$, $\Omega_b$, $H_0$, etc.) et les paramètres astrophysiques sous-maille (efficacités de rétroaction des supernovae et des noyaux actifs de galaxie).

La méthodologie centrale emploie la Régression Symbolique (en utilisant l'algorithme PySR) pour découvrir des formes fonctionnelles analytiques compactes reliant les propriétés des galaxies à $\Omega_m$.

1. Préparation des données : L'ensemble d'entraînement est composé de galaxies échantillonnées à partir des simulations IllustrisTNG à $z=0$. Des variables sans dimension clés sont construites, notamment la métallicité stellaire ($Z_\star$), la fraction de masse stellaire ($f_\star$), la vitesse dynamique effective ($V_{\text{eff}}$) et un paramètre de compacité ($R_{\text{compact}}$).
2. Recherche symbolique : L'algorithme recherche une expression minimisant l'erreur quadratique moyenne entre le $\Omega_m$ prédit et le vrai, sous réserve de pénalités de parcimonie et de contraintes sur la complexité des opérateurs (par exemple, en limitant la profondeur d'imbrication et les plages d'exposants).
3. Validation : La formule analytique résultante, entraînée exclusivement sur IllustrisTNG, est testée contre des suites de simulations indépendantes (ASTRID, SIMBA, Swift-EAGLE) et à travers différents décalages vers le rouge ($z=0$ à $z=3$).
4. Espace de paramètres étendu : Le cadre est davantage testé sur l'ensemble de données TNG-SB28, qui fait varier 28 paramètres simultanément, y compris $\Omega_b$ et $H_0$, pour évaluer la robustesse face à des variations cosmologiques et de rétroaction plus larges.

Contributions clés

1. Premier traceur analytique de galaxie unique : L'article présente la première formule analytique physiquement interprétable reliant $\Omega_m$ directement aux propriétés individuelles des galaxies. L'expression dérivée prend la forme :
   $$\Omega_m \approx \ln \left[ \sigma \left( \frac{Z_\star}{M_\star} \frac{\Omega_b V_{\text{eff}}^3}{G H_0} \frac{\Omega_b}{k} \right) + c_0 \right] - \frac{a_0}{R_{\text{compact}}}$$
   où $\sigma(x)$ est une fonction sigmoïde.
2. Robustesse inter-suites : La forme fonctionnelle dérivée d'IllustrisTNG récupère avec succès $\Omega_m$ dans les simulations ASTRID, SIMBA et Swift-EAGLE avec seulement un recalibrage mineur des coefficients numériques ($k, c_0, a_0$). Cela suggère que le couplage physique sous-jacent est universel à travers différentes implémentations numériques de la physique de formation des galaxies.
3. Interprétation physique : Les auteurs fournissent une interprétation physique transparente de la variable latente. Le terme à l'intérieur de la sigmoïde encode un couplage baryon-cosmologie médié par la relation Masse-Métallicité (MZR). Il représente l'interplay entre la profondeur du potentiel gravitationnel (déterminée par $V_{\text{eff}}$ et la cosmologie) et l'efficacité de rétention des baryons et d'enrichissement chimique régulée par la rétroaction.
4. Évolution en décalage vers le rouge : Les auteurs dérivent une extension de la formule dépendante du décalage vers le rouge en utilisant des lois d'échelle en puissance de $(1+z)$, qui maintient une précision prédictive jusqu'à $z=3$ pour IllustrisTNG et ASTRID.
5. Généralisation à un espace de paramètres de haute dimension : Dans la suite TNG-SB28 (où $\Omega_b$ et $H_0$ varient), la forme analytique s'adapte pour inclure une dépendance explicite au rapport du paramètre de Hubble ($H_0/H_{0,\text{fid}}$) et un terme d'efficacité de rétention affiné $E(Z_\star, f_\star)$, cohérent avec un cadre de régulateur de gaz.

Résultats

• Performance prédictive : À $z=0$, le modèle atteint des coefficients de détermination ($R^2$) de 0,77 pour IllustrisTNG et 0,73 pour ASTRID. SIMBA et Swift-EAGLE montrent une performance légèrement dégradée ($R^2 \approx 0,6-0,7$) mais présentent toujours une forte corrélation monotone, indiquant que le signal n'est pas un artefact de surajustement à une simulation spécifique.
• Comportement de saturation : Dans l'ensemble de données multidimensionnel TNG-SB28, le prédicteur présente un régime de saturation à haut $\Omega_m$ ($\Omega_m \gtrsim 0,43$). Les auteurs attribuent cela à de fortes dégénérescences entre la physique de la rétroaction et la cosmologie dans l'espace de paramètres de haute dimension, où les changements de cosmologie sont de plus en plus compensés par les processus baryoniques, créant un plafond d'information.
• Mécanisme physique : L'analyse confirme que le « genou » de la relation $\Omega_m$–métallicité correspond à la transition entre l'enrichissement dominé par la cosmologie (faible métallicité) et l'enrichissement saturé par la rétroaction (forte métallicité). La pente de cette transition est régie par la prescription de rétroaction, expliquant pourquoi différentes simulations nécessitent des calibrations de coefficients différentes.

Signification et affirmations
L'article affirme établir un pont direct et interprétable entre la physique des galaxies à petite échelle et la cosmologie à grande échelle. En démontrant que $\Omega_m$ est localement encodé dans le cycle baryonique des galaxies individuelles, ce travail ouvre une voie complémentaire pour l'inférence cosmologique qui ne repose pas sur les statistiques traditionnelles de la structure à grande échelle.

Les auteurs soulignent que cette approche ne vise pas à remplacer les sondes standard mais à fournir un cadre novateur pour comprendre comment l'information cosmologique survit au traitement baryonique. La découverte d'une relation analytique compacte suggère que l'empreinte de la densité de matière sur la formation des galaxies est robuste et régie par des principes fondamentaux de rétention des baryons et de profondeur du potentiel gravitationnel. Le travail met en évidence que, si l'apprentissage automatique peut identifier des corrélations, la régression symbolique est essentielle pour révéler les lois physiques sous-jacentes régissant la connexion galaxie-cosmologie.