IllustrisTNG50 angular momentum maps: tracing the morpho-kinematic evolution of galaxies

Résumé technique : Cartes de moment angulaire d'IllustrisTNG50 : Tracer l'évolution morpho-cinématique des galaxies

Problème et motivation
Bien que les simulations cosmologiques modernes reproduisent avec succès la relation d'échelle globale entre le moment angulaire spécifique stellaire total ($j_\star$) et la masse stellaire totale ($M_\star$), connue sous le nom de relation de Fall, reproduire cette relation globale seule est insuffisant pour comprendre pleinement le couplage entre l'évolution des galaxies et l'acquisition du moment angulaire. Les observations de galaxies de type tardif ont récemment révélé une diversité morpho-cinématique inattendue au sein des disques galactiques, caractérisée par des distributions spatiales distinctes de la densité de surface du moment angulaire spécifique stellaire (sAMSD). Cette diversité suggère que le moment angulaire est redistribué dans des sous-structures spécifiques ($j_\star$-sous-structures) qui ne s'alignent pas toujours avec les classifications morphologiques traditionnelles. L'objectif principal de cette étude est de quantifier et de comprendre cette diversité morpho-cinématique au sein de simulations hydrodynamiques cosmologiques, en testant spécifiquement si les $j_\star$-sous-structures identifiées dans les observations (Pacheco-Arias et al. 2026, ci-après PA26) sont reproduites dans la suite de simulations IllustrisTNG50.

Méthodologie
Les auteurs ont analysé un échantillon d'environ 8 000 disques stellaires provenant du sous-échantillon MW/M31 d'IllustrisTNG50, couvrant des redshifts $0 \le z \le 3,5$ et des masses stellaires $10^{9,5} \le M_\star \le 10^{11,2} M_\odot$. L'échantillon a été construit en suivant les sous-halos du catalogue MW/M31 en remontant le temps, en sélectionnant uniquement ceux hébergeant des disques stables et rotationnellement supportés ($V/\sigma \ge 0,5$) avec une résolution spatiale suffisante ($R_{max} \ge 7$ kpc) et une faible contre-rotation ($f_{CR} \le 0,20$).

La méthodologie centrale impliquait :

1. Construction de la carte sAMSD : Calcul du moment angulaire spécifique stellaire ($j_{\star,k}$) pour chaque particule et projection de la densité de surface ($\iota_\star = dj_\star/dS$) sur une grille régulière (pixels de $500$ pc), lissée avec un noyau gaussien correspondant à la longueur de lissage stellaire ($\epsilon_\star \approx 300$ pc).
2. Métriques morpho-cinématiques : Quatre métriques ont été définies pour quantifier les $j_\star$-sous-structures :
   • $C_{2D}$ : Un coefficient de corrélation 2D mesurant la ressemblance de la carte sAMSD avec un modèle de disque de Freeman (indicateur de $j_\star$-anneaux).
   • $a_{max}^1$ : L'amplitude normalisée maximale du mode de Fourier $m=1$ (quantifiant l'asymétrie / les irrégularités de $j_\star$).
   • $a_{max}^2$ : L'amplitude normalisée maximale du mode de Fourier $m=2$ (quantifiant la bi-symétrie).
   • $R_{max}^2$ : La fraction d'étendue radiale où le mode $m=2$ culmine (distinguant les barres centralement concentrées des spirales étendues).
3. Classification : Un modèle de mélange gaussien (GMM) avec quatre composantes entièrement covariantes a été appliqué à l'espace de paramètres quadridimensionnel ($C_{2D}, a_{max}^1, a_{max}^2, R_{max}^2$). Cette approche probabiliste a assigné à chaque galaxie une probabilité a posteriori d'appartenir à l'un des quatre types de $j_\star$ identifiés : $j_\star$-irrégulier, $j_\star$-spirale, $j_\star$-anneau et $j_\star$-barre.

Résultats clés

• Diversité morpho-cinématique : Les simulations TNG50 présentent une diversité morpho-cinématique similaire aux observations. Les disques stellaires stockent et redistribuent le moment angulaire principalement dans les quatre sous-structures de $j_\star$ identifiées.
• Séquence évolutive : Une voie évolutive canonique a été identifiée sur la base des redshifts moyens de chaque classe et des probabilités de transition de Markov :
  • $j_\star$-irrégulier ($\bar{z} = 0,91$) $\rightarrow$ $j_\star$-spirale ($\bar{z} = 0,76$) $\rightarrow$ $j_\star$-anneau ($\bar{z} = 0,62$) $\rightarrow$ $j_\star$-barre ($\bar{z} = 0,39$).
  • Les $j_\star$-irréguliers et les $j_\star$-spirales agissent comme des états de transition, tandis que les $j_\star$-anneaux et les $j_\star$-barres agissent comme des attracteurs dynamiques ; une fois formés, les anneaux et les barres ont tendance à persister.
• Facteurs physiques :
  • Fraction de gaz ($f_{gas}$) : C'est le principal moteur séparant les classes. Une $f_{gas}$ élevée favorise les $j_\star$-irréguliers et les $j_\star$-spirales, tandis qu'une $f_{gas}$ faible favorise les $j_\star$-anneaux et les $j_\star$-barres.
  • Support de rotation ($V/\sigma$) : Une fois la $f_{gas}$ fixée, un $V/\sigma$ élevé favorise les $j_\star$-spirales dans les systèmes riches en gaz et les $j_\star$-anneaux dans les systèmes pauvres en gaz.
  • Environnement : À $z=0$, aucune corrélation forte n'a été trouvée entre les métriques morpho-cinématiques et l'environnement local de la toile cosmique (densité, distance aux filaments/nœuds) ou l'alignement du vecteur de moment angulaire avec les filaments.
• Diagramme $j_\star$-$M_\star$ : Les galaxies se regroupent dans des régions distinctes du plan $j_\star$-$M_\star$ selon leur type de $j_\star$. La relation de Fall la mieux ajustée pour l'échantillon donne une pente $\alpha = 0,80 \pm 0,01$, plus raide que certaines études antérieures de TNG50 en raison de la sélection cinématique ($V/\sigma \ge 0,5$) qui exclut les systèmes dominés par la dispersion.

Signification et conclusions
L'article établit une voie canonique pour la redistribution du moment angulaire au sein des disques évoluant de manière séculaire dans TNG50. L'étude démontre que la classification par $j_\star$-sous-structure fournit un cadre qui intègre plus efficacement l'information dynamique et morphologique que la morphologie traditionnelle seule, car les $j_\star$-sous-structures ne coïncident pas toujours avec leurs homologues morphologiques (par exemple, une barre morphologique n'implique pas nécessairement une $j_\star$-barre).

Les résultats suggèrent que l'accrétion de gaz et le feedback stellaire sont les principaux processus physiques pilotant la redistribution du moment angulaire. La prévalence observée des $j_\star$-anneaux et des $j_\star$-barres à $z=0$ est attribuée aux critères de sélection de l'échantillon (favorisant les disques stables) et à l'absence de mécanismes de régulation efficaces pour la destruction des barres dans la physique des simulations TNG50. Les auteurs notent l'absence de $j_\star$-grumeaux (clumps) dans la simulation, probablement due aux limites de résolution, ainsi qu'une surabondance de $j_\star$-barres, suggérant la nécessité de processus de feedback affinés dans les futures itérations de simulation. Enfin, l'analyse sémantique de la sAMSD lie les variations de la dynamique stellaire à la redistribution de la masse, offrant un outil pour reconstruire des histoires évolutives galactiques complètes qui vont au-delà des indicateurs cinématiques standards.