Empirical estimates of how massive galaxies can be in {\Lambda}CDM

Résumé technique : Estimations empiriques de la masse maximale des galaxies dans le modèle ΛCDM

Énoncé du problème
Les observations récentes du télescope spatial James Webb (JWST) ont révélé une population de galaxies inattendument massives à des décalages vers le rouge élevés ($6 \lesssim z \lesssim 10$), avec des masses stellaires déduites dépassant $M_* \sim 10^{10} M_\odot$. Ces découvertes ont suscité un débat concernant leur compatibilité avec le modèle cosmologique standard $\Lambda$CDM. Des études antérieures suggéraient que de telles galaxies pourraient nécessiter des efficacités de formation d'étoiles approchant ou dépassant le budget baryonique cosmique disponible, remettant potentiellement en cause l'assemblage hiérarchique des structures prédit par $\Lambda$CDM. Cependant, ces comparaisons reposent souvent sur les fonctions de masse stellaire observées (GSMF) sans tenir pleinement compte des incertitudes systématiques, telles que le biais d'Eddington résultant d'erreurs aléatoires dans les estimations de masse stellaire et la dispersion intrinsèque de la relation masse étoile-masse halo (SHMR). Cet article examine si la tension apparente entre les observations du JWST et $\Lambda$CDM persiste lorsque ces effets statistiques et observationnels sont rigoureusement intégrés.

Méthodologie
Les auteurs utilisent un cadre empirique basé sur les statistiques des valeurs extrêmes (EVS) pour estimer les limites probabilistes des galaxies et des halos de matière noire les plus massifs en fonction du décalage vers le rouge et de la surface d'enquête.

1. Cadre EVS : Au lieu de s'appuyer sur des modèles spécifiques de formation galactique, l'étude utilise la GSMF observée et la fonction de masse des halos (HMF) pour dériver la distribution de probabilité des objets les plus extrêmes. La fonction de répartition cumulative (CDF) de la masse stellaire maximale, $M_{*,\text{max}}$, est dérivée de la CDF de la GSMF, $F_*(M_*)$, élevée à la puissance du nombre de galaxies indépendantes, $N$, dans un volume d'enquête donné.
2. Correction des incertitudes : L'analyse distingue trois niveaux de GSMF :
   • Observée ($\phi_{*,\text{obs}}$) : Les données brutes.
   • Vraie ($\phi_{*,\text{true}}$) : Déconvoluée pour éliminer les erreurs aléatoires (biais d'Eddington) dans les estimations de masse stellaire.
   • Intrinsèque ($\phi_{*,\text{int}}$) : Déconvoluée davantage pour tenir compte de la dispersion intrinsèque de la SHMR.
     Les auteurs modélisent ces erreurs à l'aide de distributions gaussiennes, avec des plages d'incertitude augmentant avec le décalage vers le rouge (jusqu'à $\sim 0,6$ dex à $z \sim 7-8$).
3. Contraintes sur les halos : L'étude calcule la masse baryonique théorique maximale disponible dans les halos les plus massifs ($M_{\text{bar,max}} = f_{\text{bar}} M_{\text{vir,max}}$) en utilisant sept HMF différentes (par exemple, Despali et al. 2016, Sheth & Tormen 1999) et les paramètres cosmologiques de Planck.
4. Lien Galaxie-Halo : Pour comparer les galaxies et les halos, les auteurs explorent trois scénarios de corrélation :
   • Stochastique : En supposant que la galaxie la plus massive et le halo le plus massif sont des variables aléatoires indépendantes.
   • Dispersion corrélée positivement (PCS) : En supposant une correspondance d'ordre de rang où le halo le plus massif héberge la galaxie la plus massive.
   • Rang inversé : Le halo le plus massif héberge la galaxie la moins massive (utilisé comme borne inférieure).
5. Réservoirs baryoniques : L'étude estime la masse totale de gaz froid ($M_{\text{cold}}$) disponible pour la formation d'étoiles en modélisant les flux froids et les écoulements de refroidissement radiatif, en distinguant les modes d'accrétion chaude et froide.

Contributions et résultats clés

• Limites empiriques sur la masse des galaxies : Les auteurs fournissent des prédictions EVS pour les galaxies les plus massives à travers diverses surfaces d'enquête (de $10 \text{ arcmin}^2$ à l'ensemble du ciel). Pour une enquête couvrant tout le ciel à $z \sim 0$, la galaxie la plus massive attendue à $5\sigma$ a une masse $M_* \sim 7 \times 10^{12} M_\odot$, diminuant jusqu'à $\sim 10^{10} M_\odot$ à $z \sim 16$. Ces estimations sont globalement cohérentes avec les observations actuelles du JWST lorsque la surface d'enquête est prise en compte.
• Impact du biais d'Eddington : L'étude démontre que les erreurs aléatoires dans les estimations de masse stellaire gonflent considérablement l'abondance observée des galaxies massives. Lors de la correction de ces erreurs (passage de $\phi_{*,\text{obs}}$ à $\phi_{*,\text{true}}$), la masse stellaire maximale déduite diminue d'environ $\sim 0,4-0,5$ dex à $z \sim 5$ et presque $\sim 1$ dex à $z \sim 16$. L'inclusion de la dispersion intrinsèque de la SHMR réduit encore la masse maximale d'environ $\sim 0,1-0,2$ dex à $z \sim 5$ et $\sim 0,5-0,6$ dex à $z \sim 16$.
• Résolution de la tension de la « galaxie impossible » : En comparant la masse stellaire maximale intrinsèque ($M_{*,\text{max|int}}$) à la masse baryonique maximale disponible ($M_{\text{bar,max}}$), la tension apparente disparaît largement. Sous l'hypothèse d'une correspondance parfaite d'ordre de rang (PCS), aucune galaxie n'est trouvée dépassant la limite baryonique théorique ($M_* < 0,16 \times M_{\text{vir}}$) à aucun décalage vers le rouge une fois les incertitudes de mesure incluses.
• Efficacité de formation d'étoiles (SFE) :
  • Aux décalages vers le rouge intermédiaires ($2 \lesssim z \lesssim 6$), les galaxies les plus massives approchent le réservoir total de gaz froid ($M_{*,\text{max|int}} \sim M_{\text{cold}}$), suggérant des efficacités de formation d'étoiles proches du maximum ($\epsilon \sim 1$) durant cette époque, cohérent avec le pic des galaxies formant des étoiles submillimétriques/poussiéreuses.
  • Aux décalages vers le rouge plus élevés ($z \gtrsim 8$), l'efficacité déduite chute considérablement ($\epsilon \ll 1$) lors de la prise en compte des incertitudes, indiquant que les halos ne convertissent pas les baryons avec l'efficacité extrême précédemment suggérée.
• Luminosité UV et poussière : L'étude modélise les luminosités UV en tenant compte de l'atténuation par la poussière et de la fraction de gaz déjà convertie en étoiles. Bien que le modèle reproduise la luminosité UV observée pour de petites surfaces d'enquête à haut décalage vers le rouge, il prédit que les galaxies les plus brillantes dans de plus grandes surfaces d'enquête devraient être plus faibles que certaines détections actuelles du JWST, suggérant une surestimation potentielle des luminosités ou des effets de sélection dans les plus grands champs.
• Paramètres cosmologiques : Les variations des paramètres cosmologiques (dans les incertitudes de Planck) ont un impact négligeable sur les prédictions EVS pour les halos les plus massifs, indiquant que la tension n'est pas entraînée par des incertitudes sur les paramètres cosmologiques.

Importance et affirmations
L'article affirme que le conflit apparent entre les observations du JWST de galaxies massives à haut décalage vers le rouge et le modèle $\Lambda$CDM est largement un artefact d'effets statistiques et observationnels, spécifiquement le biais d'Eddington et la dispersion intrinsèque de la SHMR. En appliquant un cadre EVS entièrement empirique qui tient rigoureusement compte de ces incertitudes et des volumes d'enquête, les auteurs concluent que les observations actuelles restent globalement cohérentes avec $\Lambda$CDM.

L'étude souligne que la surface d'enquête est un paramètre fondamental dans la définition des extrêmes attendus des populations galactiques ; les enquêtes de petite surface sont intrinsèquement biaisées vers des extrêmes de masse inférieurs et ne peuvent sonder la pleine queue de la distribution. Les auteurs affirment que le problème de la « galaxie impossible » découle probablement d'incertitudes systématiques dans les estimations de masse stellaire plutôt que d'un échec du modèle cosmologique standard ou de la nécessité d'une nouvelle physique (comme l'énergie noire précoce). Ils suggèrent en outre que les halos les plus massifs à $2 \lesssim z \lesssim 6$ peuvent en effet héberger des galaxies opérant à des efficacités de formation d'étoiles proches du maximum, tandis qu'à des décalages vers le rouge plus élevés, l'efficacité est inférieure à ce qui avait été précédemment revendiqué.