Impact of initial mass function on the chemical evolution of high-redshift galaxies

Titre : Impact de la fonction de masse initiale sur l'évolution chimique des galaxies à haut redshift

1. Problématique et Contexte

Les observations récentes du télescope spatial James Webb (JWST) ont mis en évidence une relation invariante entre la masse stellaire, la métallicité et le taux de formation d'étoiles (SFR) dans les galaxies jusqu'à un redshift $z \sim 8$, suivie d'une rupture à des redshifts plus élevés. Comprendre la physique sous-jacente à ces corrélations est crucial pour modéliser l'évolution des premières galaxies.

Le problème central abordé par cet article est l'incertitude concernant la Fonction de Masse Initiale (IMF) dans l'univers primordial. Bien que l'IMF soit souvent supposée universelle (de type Kroupa), des preuves observationnelles et théoriques suggèrent qu'elle pourrait être "plus lourde" (top-heavy) à haut redshift, favorisant la formation d'étoiles très massives ($\gtrsim 100 M_\odot$). Ces étoiles massives sont les principaux moteurs de l'enrichissement chimique via les vents stellaires et les explosions de supernovae (SN), notamment les supernovae à instabilité de paires (PISNe). L'objectif est de déterminer comment la limite supérieure de l'IMF ($m_{max}$) et les paramètres des vents galactiques influencent l'évolution chimique et l'histoire de la formation d'étoiles pour reproduire les observations du JWST.

2. Méthodologie

Les auteurs utilisent le code semi-analytique de évolution galactique a-sloth, appliqué aux arbres de fusion issus d'une simulation cosmologique à haute résolution (Ishiyama et al. 2016) dans un univers $\Lambda$CDM standard.

Mises à jour et améliorations du modèle :

• Formation stellaire : La prescription a été mise à jour pour mieux capturer la formation d'étoiles "explosive" (bursty) et la taille des galaxies. Le rayon caractéristique $R_\star$ est désormais défini comme $2.5 R_{50}$ (rayon contenant 50% de la masse stellaire) plutôt que le rayon de l'halo, ce qui réduit la taille effective des régions de formation stellaire et améliore l'accord avec les observations.
• Évolution stellaire et rendements : Le modèle intègre des tracks d'évolution stellaire complets (code PARSEC) couvrant une large gamme de métallicités ($Z \sim 10^{-11} - 0.03$) et de masses ($2 - 600 M_\odot$). Les rendements en métaux incluent désormais non seulement les supernovae (effondrement de cœur, PISNe, SN de type Ia), mais aussi les vents stellaires, qui contribuent significativement à l'enrichissement chimique.
• IMF variable : L'étude explore une IMF de type Kroupa avec une limite supérieure de masse ($m_{max}$) variable entre $100$ et $600 M_\odot$, tout en maintenant la pente de la partie haute masse fixe ($\alpha = 2.3$).
• Rétroaction (Feedback) : Les vents galactiques sont modélisés avec une efficacité dépendante de la masse de l'halo, contrôlée par deux paramètres libres : la pente $\alpha_{out}$ et l'échelle de masse $M_{out0}$.

Paramétrisation et Contraintes :

• Espace des paramètres : 165 simulations ont été exécutées en variant $\alpha_{out}$, $M_{out0}$ et $m_{max}$.
• Contraintes observationnelles : Les modèles sont comparés à trois ensembles de données via une fonction de vraisemblance ($L_{all}$) :
  1. L'histoire de la formation d'étoiles cosmique (CSFH) et la densité de masse stellaire (CSMD).
  2. La relation Masse-Métallicité-SFR (MZSFR), incluant la relation fondamentale métallicité (FMR) et la relation masse-métallicité (MZR) observées par le JWST ($z \sim 4-10$).
  3. La relation Masse Étoile-Masse Halo (SHMR).

3. Contributions Clés

• Intégration des vents stellaires et des PISNe : L'étude démontre que l'inclusion des rendements en métaux provenant des vents stellaires et des PISNe est essentielle pour reproduire les abondances chimiques observées.
• Contrainte rigoureuse sur l'IMF : En combinant l'histoire de la formation d'étoiles et la relation MZSFR, l'article impose des contraintes strictes sur la limite supérieure de l'IMF, montrant que les modèles standards ($m_{max} \le 100 M_\odot$) échouent à reproduire simultanément les deux observables.
• Modélisation des galaxies faibles : Le modèle prédit une population importante de galaxies faibles, peu massives et pauvres en métaux qui dominent le budget de formation d'étoiles à très haut redshift, mais qui sont invisibles aux sondages actuels.

4. Résultats Principaux

A. Nécessité d'étoiles très massives ($m_{max} \gtrsim 200 M_\odot$)

• Les modèles avec $m_{max} = 100 M_\odot$ ne peuvent pas reproduire la relation MZSFR observée à $z \gtrsim 4$ ; ils produisent trop peu de métaux.
• Seul un modèle avec une limite supérieure de masse $m_{max} \gtrsim 200 M_\odot$ (idéalement $250 M_\odot$) permet de générer suffisamment de métaux via les PISNe pour correspondre aux observations, tout en respectant l'histoire de la formation d'étoiles cosmique.
• Le modèle optimal trouvé est : $\alpha_{out} = 0$ (efficacité de vent constante), $M_{out0} = 3 \times 10^9 M_\odot$, et $m_{max} = 250 M_\odot$. Ce modèle atteint un accord à moins de $1\sigma$ avec les données.

B. Rôle des Supernovae à Instabilité de Paires (PISNe)

• L'enrichissement chimique requis pour expliquer les observations du JWST implique une contribution significative des PISNe provenant d'étoiles de Population II (métalliques) à $z \lesssim 15$.
• Le taux de PISNe prédit pour le modèle optimal atteint un pic de $\sim 10$ événements par an par degré carré à $z \sim 5$, dominé par les étoiles de Population II. Sans étoiles très massives, ce taux chuterait drastiquement.

C. Population de galaxies faibles et réionisation

• Les galaxies détectables par le JWST (luminieuses en UV) ne représentent qu'une fraction mineure de la formation d'étoiles totale.
• À $z \gtrsim 12$, les galaxies faibles ($SFR < 0.3 M_\odot yr^{-1}$) contribuent à plus de 90% de la densité de taux de formation d'étoiles (SFRD).
• Ces galaxies, bien que peu lumineuses, sont majoritairement pauvres en métaux ($Z \lesssim 0.1 Z_\odot$) et constituent le réservoir principal pour la formation de binaires d'étoiles massives.

D. Implications pour les ondes gravitationnelles

• La forte contribution des galaxies faibles et pauvres en métaux suggère que le taux de fusion de trous noirs binaires à haut redshift est sous-estimé si l'on ne considère que les galaxies lumineuses.
• Les auteurs prévoient que ces galaxies pauvres en métaux dominent le SFRD à $z \gtrsim 8$, offrant un canal de formation privilégié pour les fusions de trous noirs détectables par les futurs détecteurs d'ondes gravitationnelles (3ème génération).

5. Signification et Conclusion

Cette étude fournit une preuve théorique robuste reliant les observations du JWST à la physique des étoiles massives. Elle conclut que :

1. L'IMF dans l'univers primordial ($z \gtrsim 4$) doit être top-heavy, permettant la formation d'étoiles jusqu'à $\sim 250 M_\odot$.
2. Les PISNe jouent un rôle central, non seulement dans l'enrichissement chimique rapide de l'univers, mais aussi dans la génération d'événements transitoires électromagnétiques et de sources d'ondes gravitationnelles.
3. L'histoire de la formation d'étoiles cosmique est dominée par une population de galaxies faibles et pauvres en métaux, invisibles aux sondages actuels mais cruciales pour la réionisation et l'évolution chimique.

Ces résultats ouvrent la voie à de nouvelles observations spectroscopiques pour détecter les signatures des PISNe et des étoiles très massives, et suggèrent que les futurs détecteurs d'ondes gravitationnelles pourraient contraindre l'évolution des galaxies à des époques inaccessibles par la lumière directe.