Testing the inference of kinematics from mock JWST NIRSpec/MSA observations of TNG50 galaxies at $z\sim2-6$

Résumé technique : Test de l'inférence de la cinématique à partir d'observations simulées de JWST NIRSpec/MSA de galaxies TNG50 à $z \sim 2-6$

Énoncé du problème
Comprendre l'évolution dynamique des galaxies nécessite de mesurer les vitesses de rotation ($v$) et les dispersions de vitesse ($\sigma$) du gaz ionisé à travers les époques cosmiques. Bien que les simulations hydrodynamiques comme IllustrisTNG fournissent des données spatiales et cinématiques en 6 dimensions, elles manquent de solutions explicites de transfert de rayonnement pour des configurations observationnelles spécifiques. Inversement, les études observationnelles utilisant des instruments tels que le NIRSpec MSA de JWST reposent sur la modélisation directe pour inférer la cinématique à partir de spectres 2D. Ces modèles supposent généralement que les galaxies sont des disques tournants fins et axisymétriques avec des profils d'intensité de type Sérsic et des courbes de vitesse en arc tangente. Cependant, les galaxies à haut redshift ($z > 2$) présentent souvent des morphologies complexes (sphéroïdales, parsemées de amas/clumps) et des mouvements non circulaires (flux entrants/sortants) qui peuvent violer ces hypothèses. Il existe un manque critique de validation statistique concernant la capacité des modèles paramétriques simples à récupérer précisément la cinématique intrinsèque pour les galaxies typiques à haut redshift, particulièrement compte tenu de la résolution spatiale limitée et de la nature par fentes des observations MSA par rapport aux unités de champ intégral (IFU).

Méthodologie
Les auteurs utilisent la simulation cosmologique haute résolution TNG50 pour générer une chaîne de traitement complète d'observations simulées pour plus de 4 000 galaxies à formation stellaire ($M_* = 10^8 - 10^{11,5} M_\odot$) aux redshifts $z = 2, 3, 4, 6$.

1. Mesures intrinsèques : Les auteurs définissent les quantités cinématiques « intrinsèques » ($v_{max}$ et $\sigma_e$) au sein du rayon effectif ($r_e$) du gaz émettant l'H$\alpha$. Ils testent explicitement la sensibilité de ces valeurs intrinsèques aux choix de mesure, incluant la taille de l'ouverture, le binnage spatial (0,5 vs 5 ckpc) et les schémas de pondération (non pondéré vs pondéré par le taux de formation stellaire [SFR]).
2. Observations simulées :
   • Imagerie : Des images simulées NIRCam sont générées en projetant les cellules de gaz en formation stellaire, en convertissant les taux de formation stellaire (SFR) en luminosité H$\alpha$ via la relation de Kennicutt (1998), et en convoluant avec les PSF de NIRCam. Le bruit est ajouté pour imiter des expositions profondes de type JADES.
   • Spectroscopie : Des spectres 2D simulés NIRSpec/MSA sont créés en projetant les cellules de gaz selon un axe aléatoire, en assignant des vitesses de ligne de visée, et en convoluant avec la fonction de diffusion de la ligne (LSF) et la PSF de NIRSpec. Un bruit gaussien est appliqué pour atteindre un rapport signal sur bruit (S/N) comparable aux données JWST les plus profondes (S/N $\sim$ 50).
3. Pipeline d'inférence :
   • Morphologie : Les images simulées NIRCam sont ajustées avec pysersic pour dériver les priors morphologiques (rayon effectif $r_e$, indice de Sérsic $n_s$, rapport d'aspect $q$, angle de position PA).
   • Cinématique : Les spectres NIRSpec simulés sont ajustés en utilisant msafit, supposant un modèle de disque fin en rotation avec un profil de vitesse en arc tangente et une dispersion de vitesse constante. Les paramètres morphologiques dérivés de l'imagerie sont utilisés comme priors forts. L'ajustement emploie une minimisation $\chi^2$ pondérée, incluant une mise à l'échelle du flux non paramétrique ligne par ligne pour tenir compte des distributions de lumière irrégulières (clumpy).

Contributions clés et résultats
L'étude compare les paramètres cinématiques optimaux dérivés des spectres simulés par rapport aux valeurs intrinsèques de la simulation afin de quantifier les biais observationnels.

• Précision de la récupération : Pour les sources alignées et résolues à $z=3$, les vitesses de rotation ($v$) et les dispersions de vitesse ($\sigma$) inférées concordent généralement avec les mesures intrinsèques à un facteur de $\sim 2$ et $\sim 1,5$, respectivement. La dispersion dans $v$ est plus grande que celle de $\sigma$.
• Biais systématiques :
  • Dépendance morphologique : La récupération est robuste pour les systèmes lisses et disqués. Cependant, pour les systèmes intrinsèquement allongés ou sphéroïdaux (grands rapports d'axes intrinsèques $C/A$), $v$ tend à être surestimée et $\sigma$ sous-estimée d'un facteur de plus de 2. Cela se produit car l'hypothèse du disque fin force le modèle à interpréter l'élargissement non rotationnel comme un effet d'inclinaison.
  • Évolution avec le redshift : La dispersion de la précision de récupération augmente à des redshifts plus élevés ($z > 3$) à mesure que les galaxies s'écartent davantage de l'hypothèse du disque fin.
  • Individuel vs Population : Bien que la cinématique de galaxies individuelles ($v/\sigma$) puisse être mal classée par un facteur de 2 à 3 (par exemple, un système supporté par la dispersion apparaissant comme supporté par la rotation), les tendances de population sont récupérées de manière robuste. L'évolution faible de $\sigma$ et de $v/\sigma$ avec le redshift, ainsi que la dépendance de $\sigma$ vis-à-vis du SFR global, sont préservées dans les données inférées par simulation.
• Relations d'échelle : La relation inférée $\sigma$-SFR s'aligne bien avec les modèles théoriques de « transport de masse + feedback » (Krumholz et al. 2018), suggérant qu'en dépit des incertitudes individuelles, la modélisation cinématique peut contraindre de manière fiable les moteurs physiques de la turbulence (feedback vs gravité) dans les populations de galaxies.
• Modes d'échec : Environ 5 à 10 % des sources échouent à l'ajustement au modèle de disque fin, souvent en raison de morphologies irrégulières, de systèmes à masse stellaire élevée présentant des gradients de vitesse centraux abrupts non capturés par les profils en arc tangente, ou de sources où la $\sigma$ intrinsèque est inférieure à la limite inférieure du modèle (10 km s$^{-1}$).

Signification et affirmations
L'article affirme fournir le premier cadre complet permettant de comparer les observations NIRSpec MSA aux simulations cosmologiques pour quantifier les biais de mesure de la cinématique du gaz ionisé. Les auteurs concluent que si les modèles de disques fins simples sont insuffisants pour déterminer la « discuité » ($v/\sigma$) précise de galaxies individuelles à haut redshift (en raison des grandes incertitudes), ils sont efficaces pour récupérer les tendances statistiques de population.

Le travail souligne que les inférences observationnelles de $v/\sigma$ pour des systèmes individuels à $z \ge 2$ doivent être traitées avec prudence, car elles peuvent être biaisées par l'hypothèse de géométrie de disque fin et par la présence de mouvements non circulaires. Les auteurs soutiennent que pour améliorer les contraintes sur l'état dynamique des galaxies primitives, les travaux futurs doivent dépasser la spectroscopie par fente 2D pour s'orienter vers des observations IFU 3D ou des stratégies de déplacement de fente (slit-stepping), et développer des modèles cinématiques dédiés capables de prendre en compte les morphologies 3D complexes et l'aspect granulaire (clumpiness) inhérents aux galaxies à haut redshift. L'étude souligne que les limitations actuelles ne sont pas simplement instrumentales, mais découlent du décalage entre les modèles paramétriques simples et la complexité intrinsèque de la physique de la formation des galaxies.