NISER-IUCAA New Simulations of JWST GAlaxies and Quasars(NINJA): Properties of galaxies at $5 \leq z \leq 10$

La suite de simulations hydrodynamiques cosmologiques NINJA démontre que, bien que des modèles spectraux et de poussière spécifiques puissent reproduire les fonctions de luminosité UV observées pour les galaxies à $5 \leq z \leq 10$, des dégénérescences significatives entre les propriétés de rétroaction et de poussière, ainsi que les limitations de résolution à $z > 10$, nécessitent des simulations à plus haute résolution et des observations multi-longueurs d'onde pour contraindre de manière robuste l'évolution des galaxies à haut décalage vers le rouge.

L'essentiel

Imaginez l'univers primordial comme un vaste chantier de construction sombre, quelques centaines de millions d'années seulement après le Big Bang. Le télescope spatial James Webb (JWST) est comparable à une nouvelle grue puissante et une équipe de tournage qui sont enfin arrivées pour prendre des photos haute définition des tout premiers bâtiments (galaxies) en cours de construction.

Ce document, intitulé NINJA, est un rapport d'une équipe d'astronomes qui ont construit leur propre « chantier virtuel » à l'aide de superordinateurs. Leur objectif était de vérifier si leurs modèles numériques pouvaient correspondre aux vraies photos que le JWST prend de ces galaxies anciennes.

Voici une décomposition de ce qu'ils ont fait et découvert, en utilisant des analogies simples :

1. Le chantier virtuel (Les simulations)

Les chercheurs ont créé trois univers virtuels de tailles différentes (comme construire une maquette de ville dans une boîte à chaussures, un salon et un stade). Ils ont rempli ces boîtes de matière noire et de gaz, laissant la gravité les rassembler pour former des galaxies.

• Le défi : Ils devaient s'assurer que leurs galaxies numériques ressemblaient aux vraies. Plus précisément, ils devaient correspondre à la luminosité de ces galaxies dans la lumière ultraviolette (UV), qui est la principale façon dont nous voyons les étoiles jeunes et chaudes.

2. Le problème du « filtre de poussière »

Dans le monde réel, si vous essayez de prendre une photo d'une ampoule à travers une vitre sale, la lumière paraît plus terne et plus rouge. Dans l'espace, cette « vitre sale » est la poussière cosmique.

• Le problème : L'équipe a constaté que leurs galaxies numériques étaient naturellement trop brillantes et trop bleues par rapport à ce que le JWST observe. Pour corriger cela, ils ont dû ajouter un « filtre de poussière » à leurs modèles.
• L'expérience : Ils ont essayé différents types de « recettes de poussière ». Certaines recettes supposaient que la poussière se formait selon une relation linéaire simple avec les métaux (comme mélanger de la peinture). D'autres ont essayé des recettes plus complexes où la formation de poussière changeait radicalement en fonction de la richesse en métaux de la galaxie. Ils ont également essayé différents « objectifs » (courbes d'atténuation) pour voir comment la poussière bloquait la lumière.

3. Le rapport « poussière/métaux » (L'ingrédient secret)

Pour faire correspondre leurs galaxies virtuelles aux réelles, l'équipe a dû régler un cadran appelé $\epsilon$ (epsilon). Imaginez cela comme le « bouton d'efficacité de la poussière ».

• Ce qu'ils ont découvert : Ils ont découvert que, dans l'univers primordial, les galaxies étaient beaucoup moins efficaces pour produire de la poussière que notre propre Voie lactée ne l'est aujourd'hui.
  • À un décalage vers le rouge de 5 ou 6 (très tôt), le rapport poussière/métaux n'était que d'environ 35 % de ce que nous observons dans notre voisinage local.
  • À un décalage vers le rouge de 9 ou 10 (encore plus tôt), il chutait à moins de 10 %.
• La difficulté : Le nombre exact nécessaire pour tourner le cadran dépendait fortement de quelle recette de poussière ils avaient choisie. S'ils changeaient la recette, le réglage du cadran variait d'un facteur 7 ! Cela signifie que nous ne pouvons pas être certains à 100 % de la quantité exacte de poussière existante sans davantage de données.

4. L'effet des « bébés étoiles » (Émission nébulaire)

L'équipe a réalisé qu'il manquait un ingrédient crucial : l'émission nébulaire.

• L'analogie : Imaginez un chantier de construction où les ouvriers (les étoiles) sont entourés d'un brouillard lumineux (nuages de gaz). Si vous ne comptez que la lumière des ouvriers, vous manquez la lueur du brouillard.
• Le résultat : Lorsqu'ils ont ajouté la lumière de ce « brouillard » à leurs modèles, les galaxies sont devenues plus brillantes, en particulier les plus petites et les plus faibles. Cela a permis à leurs modèles de correspondre beaucoup mieux aux observations réelles.

5. Le problème des étoiles « trop lourdes » (La fonction initiale de masse)

L'équipe a également testé ce qui se passe si l'univers primordial produisait en moyenne des étoiles « plus grosses » qu'aujourd'hui.

• L'analogie : Habituellement, une usine d'étoiles produit un mélange d'étoiles petites, moyennes et grandes (comme une boulangerie standard). Mais que se passerait-il si l'univers primordial ne cuisait que des pains géants ?
• Le résultat : S'ils supposaient que l'univers primordial produisait plus d'étoiles massives (une fonction initiale de masse « lourde en haut »), les galaxies devenaient incroyablement brillantes. Cela a aidé à mieux expliquer les galaxies les plus faibles, mais cela nécessitait encore plus de poussière pour les assombrir afin de correspondre à ce que le JWST observe.

6. Le problème « trop brillant » au bord du temps

Lorsqu'ils ont observé les toutes premières galaxies (décalage vers le rouge $z > 10$), leurs modèles ont buté contre un mur.

• Le problème : Même avec leurs meilleures recettes de poussière et hypothèses sur les étoiles, leurs galaxies virtuelles restaient trop sombres par rapport à celles que le JWST a réellement trouvées.
• La conclusion : L'article suggère que leurs modèles informatiques ne sont pas encore assez détaillés. C'est comme essayer de dessiner un portrait haute résolution avec un crayon basse résolution ; ils ont besoin d'une simulation « plus haute résolution » pour comprendre correctement ces galaxies les plus anciennes.

7. Le « rapport de Balmer » et l'« excès de couleur » (Le détective de la poussière)

L'équipe a utilisé des signatures chimiques spécifiques (comme le rapport entre deux couleurs spécifiques de lumière, H-alpha et H-beta) pour agir comme un « détective de la poussière ».

• La découverte : Ils ont constaté que la poussière autour des étoiles nouveau-nées (dans leurs « nuages de naissance ») est beaucoup plus rouge que la poussière flottant dans le reste de la galaxie.
• La divergence : Leurs modèles prédisaient que la poussière autour des étoiles et la poussière dans le reste de la galaxie devraient être quelque peu similaires. Cependant, les observations réelles suggèrent que la poussière autour des étoiles est beaucoup plus efficace pour bloquer la lumière. Cela suggère que leurs « recettes de poussière » actuelles pourraient avoir besoin d'une refonte majeure.

Résumé : Que signifie cela ?

L'équipe NINJA a réussi à construire un univers virtuel capable d'imiter la luminosité des galaxies primordiales, mais uniquement s'ils ajustent soigneusement la quantité de poussière cosmique et les types d'étoiles en formation.

• La poussière est essentielle : Même dans l'univers très primordial, la poussière se formait déjà et assombrissait la lumière, mais elle était beaucoup moins efficace qu'aujourd'hui.
• Nous avons besoin de plus de données : Parce que différentes « recettes de poussière » donnent des réponses différentes, nous avons besoin de plus d'observations (en particulier du télescope ALMA, qui observe directement la poussière) pour déterminer la recette correcte.
• Nous avons besoin d'ordinateurs plus performants : Pour comprendre les toutes premières galaxies (au-delà du décalage vers le rouge 10), leurs simulations actuelles ne sont pas assez détaillées. Ils doivent exécuter la simulation avec une résolution plus élevée pour arrêter le « pixelisation » de leurs modèles.

En bref, l'univers était un chantier de construction poussiéreux et en formation d'étoiles bien plus tôt que nous ne le pensions, mais nous cherchons encore à déterminer exactement combien de poussière se trouvait sur les vitres et à quel point les lumières étaient réellement brillantes.

Résumé technique

Résumé Technique : Simulations NINJA de Galaxies à Haut Décalage vers le Rouge (5 ≤ z ≤ 10)

Énoncé du Problème
Le télescope spatial James Webb (JWST) a fourni des observations sans précédent de galaxies à haut décalage vers le rouge ($z \gtrsim 5$), révélant des détails sur l'époque de la réionisation et la formation précoce des galaxies. Cependant, des écarts existent entre ces observations et les prédictions issues des simulations hydrodynamiques cosmologiques standard du modèle $\Lambda$CDM. Plus précisément, les observations indiquent souvent un excès de galaxies brillantes à $z > 10$ et suggèrent des systèmes plus évolués (masses stellaires et contenu en poussière plus élevés) que ce que prédisent les simulations calibrées sur des données à faible décalage vers le rouge. Un défi critique consiste à démêler les effets de l'extinction par la poussière, de la fonction initiale de masse (IMF) et des mécanismes de rétroaction, car ces paramètres sont souvent dégénérés. Les simulations actuelles peinent à reproduire simultanément la fonction de luminosité ultraviolette (UVLF), les pentes spectrales ($\beta_{UV}$), les rapports des raies de Balmer et les contraintes sur la masse de poussière à travers plusieurs décalages vers le rouge sans ajustage fin.

Méthodologie
Les auteurs présentent la suite de simulations hydrodynamiques cosmologiques NINJA (NISER-IUCAA New Simulations of JWST GAlaxies and Quasars).

• Configuration des Simulations : Les simulations utilisent le code MP-GADGET, employant une formulation de l'hydrodynamique à particules lissées (SPH) basée sur l'entropie de pression. Trois boîtes de simulation ont été exécutées avec des longueurs de côté comobiles de 50, 150 et 250 $h^{-1}$ Mpc (L50N2040, L150N2040, L250N2040), contenant chacune $2040^3$ particules de matière noire et de gaz. Cette configuration permet d'étudier une large gamme de masses d'halos tout en quantifiant les problèmes de convergence.
• Modèles Physiques : Les simulations suivent les prescriptions sous-maille de la suite ASTRID pour la formation stellaire, la rétroaction des supernovae et l'amorçage/croissance des trous noirs. La formation stellaire suit un modèle d'ISM multiphase, et les trous noirs sont amorcés dans les halos au-dessus de $5 \times 10^9 M_\odot h^{-1}$ avec une accrétion de type Bondi-Hoyle-Lyttleton subséquente et une rétroaction en double mode (thermique/cinétique).
• Post-traitement et Synthèse Spectrale :
  • Génération de SED : Les spectres intrinsèques sont générés en sommant les contributions des populations stellaires (en utilisant les modèles SSP bpass-v2.2.1) et de l'émission nébulaire (en utilisant CLOUDY).
  • Modélisation de la Poussière : L'extinction par la poussière est appliquée comme une étape de post-traitement. L'extinction visuelle ($A_V$) est dérivée de la densité de colonne d'hydrogène ($N_H$) mise à l'échelle par un paramètre de rapport poussière/métaux ($\epsilon$) et une fonction dépendante de la métallicité $f(Z_g)$.
  • Variations : L'étude explore plusieurs prescriptions :
    1. Mise à l'échelle Poussière-Métallicité : Relations linéaires vs loi de puissance brisée (double loi de puissance).
    2. Courbes d'Atténuation : Lois de Calzetti et al. (2000) vs lois d'extinction moyennes du SMC.
    3. Géométrie de la Poussière : Phase unique (ISM uniquement) vs deux phases (ISM + atténuation par « nuage de naissance » pour les étoiles < 10 Myr).
    4. IMF : IMF Chabrier standard vs une IMF Kroupa surpeuplée (KP300-TH) avec une coupure de masse supérieure à $300 M_\odot$.
• Calibration : Pour chaque décalage vers le rouge ($5 \le z \le 10$), le paramètre $\epsilon$ est ajusté par minimisation du $\chi^2$ pour correspondre à la fonction de luminosité ultraviolette (UVLF) observée.

Résultats Clés

1. Reproduction de l'UVLF : Les modèles de référence, une fois calibrés avec des valeurs appropriées de $\epsilon$, reproduisent avec succès les UVLF observées sur $5 \le z \le 10$. L'inclusion de l'émission du continuum nébulaire augmente systématiquement la luminosité UV, améliorant l'accord avec les observations, en particulier à l'extrémité faible.
2. Évolution du Rapport Poussière/Métaux : Le rapport poussière/métaux déduit diminue avec le décalage vers le rouge. Pour le modèle de référence (courbe Calzetti), le rapport à $z=5-6$ est d'environ 35 % de la valeur locale de la Voie lactée, chutant à $\le 10\%$ à $z \ge 9$. Cependant, cette normalisation varie d'un facteur d'environ 7 selon la courbe d'atténuation et la mise à l'échelle poussière-métallicité adoptées.
3. Dégénérescences et Contraintes :
   • Courbes d'Atténuation : Bien que les courbes Calzetti et SMC puissent toutes deux ajuster l'UVLF, elles prédisent des fonctions de luminosité en bande B, des luminosités H$\alpha$ et des pentes $\beta_{UV}$ différentes. La relation observée $\beta_{UV}$ vs $M_{UV}$ à $z=5-7$ favorise des courbes d'extinction plus raides (proches du SMC) par rapport à la loi Calzetti plus plate.
   • Effets de l'IMF : Une IMF surpeuplée (KP300-TH) produit plus de lumière UV par unité de taux de formation stellaire, nécessitant des rapports poussière/métaux plus élevés pour correspondre à l'UVLF. Cette IMF aide à ajuster l'UVLF à l'extrémité faible à $7 \le z \le 9$ mais implique des propriétés de poussière plus proches des valeurs locales aux époques précoces.
   • Rapports de Balmer : Les modèles incluant une atténuation par « nuage de naissance » produisent un excès de couleur nébulaire ($E(B-V)_{neb}$) plus élevé par rapport à l'excès de couleur stellaire ($E(B-V)_{star}$), cohérent avec certaines tendances à haut décalage vers le rouge, bien que le rapport $f = E(B-V)_{star}/E(B-V)_{neb}$ reste difficile à contraindre précisément avec les données actuelles.
4. Limites à Haut Décalage vers le Rouge ($z > 10$) : Les simulations sous-estiment l'UVLF à $z \ge 10$ par rapport aux observations actuelles, même en adoptant une IMF surpeuplée et en négligeant la poussière. Les auteurs attribuent cela à un manque de convergence numérique ; le nombre de galaxies dans les boîtes de simulation chute considérablement à $z \sim 15$, indiquant que des simulations à plus haute résolution sont nécessaires pour modéliser de manière robuste les galaxies à ces décalages vers le rouge.
5. Luminosité H$\alpha$ : Les simulations surestiment systématiquement la fonction de luminosité H$\alpha$ par rapport aux observations à $z \sim 5-6$. Les auteurs suggèrent que cet écart pourrait être résolu en incluant de la poussière au sein des régions HII (actuellement supposées sans poussière dans le modèle sous-maille) ou en permettant une fraction d'échappement plus élevée des photons ionisants.

Signification et Revendications
L'article revendique que la suite NINJA fournit un cadre complet pour interpréter les données du JWST en faisant varier systématiquement les paramètres de poussière, d'IMF et de rétroaction. La signification principale réside dans la démonstration que :

• La poussière est cruciale même à haut $z$ : Une extinction par la poussière non négligeable est requise pour correspondre à la forme de l'UVLF à l'extrémité brillante, impliquant une formation efficace de poussière dans les galaxies précoces.
• Des contraintes multi-observables sont nécessaires : Reproduire simultanément l'UVLF, la LF en bande B, la LF H$\alpha$ et les pentes spectrales est essentiel pour briser les dégénérescences entre les modèles de poussière et les hypothèses d'IMF.
• Limites de résolution : Les écarts actuels à $z > 10$ soulignent le besoin de simulations à plus haute résolution pour éviter les problèmes de convergence numérique.
• Perspectives Futures : Les auteurs soulignent que des contraintes indépendantes provenant d'ALMA (masse de poussière) et de la spectroscopie JWST (rapports de Balmer, ailes d'amortissement Ly$\alpha$) sont critiques pour affiner les modèles de poussière. Ils notent que leurs résultats actuels sont préliminaires concernant l'UVLF à l'extrémité faible et la luminosité H$\alpha$, qui seront abordés dans un travail futur impliquant des paramètres de rétroaction variés.

L'article conclut que si les modèles peuvent être ajustés pour correspondre à des observables spécifiques (comme l'UVLF), les paramètres physiques résultants (rapports poussière/métaux, courbes d'atténuation) varient considérablement selon les prescriptions choisies, soulignant le besoin de données multi-longueurs d'onde et multi-observables pour contraindre la physique de la formation des galaxies dans l'univers primordial.