The Entangling of Supernova Feedback Impacts with Coarsening Simulation Resolution

Résumé technique : L'enchevêtrement des impacts du feedback des supernovas avec la résolution de simulation de coarsening (grossissement)

Énoncé du problème
Le feedback des supernovas (SN) est largement reconnu comme le mécanisme principal régulant la formation stellaire (SF) et pilotant les écoulements galactiques. Cependant, les canaux spécifiques par lesquels les SN parviennent à ces effets macroscopiques restent obscurs. Alors que les simulations à haute résolution peuvent résoudre explicitement l'expansion thermique et le refroidissement radiatif des rémanents de SN, la plupart des simulations de formation de galaxies opèrent à des résolutions plus faibles. Ces modèles à basse résolution reposent sur des approximations « subgrid » (sous-maillage) qui supposent que les conditions du milieu interstellaire (ISM) au sein d'une seule cellule de gaz sont homogènes et dérivables des données disponibles de la simulation. Le problème central abordé est de savoir si les simulations à résolution grossière peuvent distinguer de manière auto-cohérente les impacts physiques distincts des SN se produisant dans différents environnements (par exemple, régions diffuses vs régions denses) ou si la moyenne inhérente aux grilles grossières occulte ces processus physiques critiques, conduisant à des estimations erronées des rayons de refroidissement et des taux de formation stellaire.

Méthodologie
Les auteurs analysent deux simulations indépendantes à ultra-haute résolution de galaxies naines isolées : LYRA (résolution de masse de gaz $m_{gas} \sim 4 M_\odot$) et RIGEL ($m_{gas} \sim 1 M_\odot$). Les deux simulations utilisent le code d'hydrodynamique à maillage mobile AREPO et partagent des conditions initiales similaires (faible fraction baryonique, $M_{vir} = 2 \times 10^{10} M_\odot$). Crucialement, les deux simulations excluent le feedback radiatif des jeunes étoiles pour isoler les effets de l'injection d'énergie par SN et du refroidissement radiatif.

L'étude se concentre sur l'environnement immédiat d'événements SN individuels, défini par la masse de gaz $M_{env}$ la plus proche de l'étoile progénitrice. Les auteurs :

1. Ont caractérisé les environnements locaux : Mesuré la densité locale du gaz ($\rho$) et le changement d'énergie spécifique ($\Delta e$) pour les événements SN dans les deux simulations.
2. Ont analysé les effets de résolution : Ont systématiquement « grossi » les données à haute résolution en faisant varier l'échelle de masse de l'environnement local ($M_{env}$) de $20 M_\odot$ (fin) à $200 M_\odot$ et $2000 M_\odot$ (grossier). Cela simule le processus de moyennage des simulations à plus basse résolution.
3. Ont corrélé avec les effets globaux : Comparé les propriétés des SN aux taux d'écoulement galactique (flux d'énergie à des hauteurs spécifiques au-dessus du disque) et aux taux de formation stellaire (SFR) locaux avant et après les événements SN.
4. Ont validé avec des modèles analytiques : Comparé les données de simulation à des solutions analytiques pour l'évolution de l'onde de choc de la SN (phase de Sedov-Taylor, estimations du rayon de refroidissement selon Kim & Ostriker 2015).

Contributions clés et résultats

• Distribution bimodale des environnements de SN : L'analyse révèle une distribution bimodale distincte des densités locales où les SN se produisent.

  • Mode de faible densité : Les SN se produisant dans des régions diffuses ($n \sim 10^{-4} - 10^{-1} \text{ cm}^{-3}$). Ces événements présentent de grands changements de l'énergie spécifique ($\Delta e$), dépassant souvent la vitesse de libération du disque. Ils sont corrélés temporellement avec les écoulements galactiques et sont capables de s'étendre à des échelles de kiloparsecs avant que les pertes radiatives ne dominent.
  • Mode de haute densité : Les SN se produisant dans des régions denses et de formation stellaire ($n \sim 10^4 \text{ cm}^{-3}$). Ces événements montrent de plus petits changements de $\Delta e$ en raison d'un refroidissement radiatif efficace. Ils perturbent efficacement les amas de formation stellaire locaux à des échelles de parsecs mais ne pilotent pas d'écoulements significatifs.

• Canaux de feedback distincts : Le papier établit que ces deux modes pilotent des comportements émergents séparés :

  • Les SN de faible densité sont les principaux moteurs des écoulements galactiques (injection d'énergie dans le milieu circumgalactique).
  • Les SN de haute densité sont les principaux moteurs de la suppression directe de la formation stellaire (disruption des nuages de gaz denses).
  • Crucialement, les SN de faible densité ne suppriment pas de manière significative la formation stellaire à leur proximité immédiate (car le gaz est trop diffus pour former des étoiles), et les SN de haute densité ne pilotent pas efficacement les écoulements (car l'énergie est rayonnée avant que la coque ne puisse s'échapper du disque).

• Le seuil de résolution : La distinction entre ces deux modes n'est préservée qu'à des résolutions très élevées ($M_{env} \lesssim 100 M_\odot$).

  • Lorsque $M_{env}$ augmente jusqu'à $200 M_\odot$, la bimodalité de la densité commence à s'estomper.
  • À $M_{env} = 2000 M_\odot$ (typique de nombreuses simulations cosmologiques), les deux modes deviennent indiscernables. Le processus de moyennage efface le contraste de densité, rendant impossible de déterminer à partir des données grossières seules si une SN va piloter un écoulement ou supprimer la formation stellaire.

• Mauvaise estimation des rayons de refroidissement : Les estimations de densité à résolution grossière ($\rho_{2000}$) échouent à prédire le véritable rayon de refroidissement ($r_{cool}$) d'une SN.

  • Pour une densité grossière donnée, le véritable rayon de refroidissement (dérivé de la densité à fine échelle) peut varier de plusieurs ordres de grandeur.
  • Les estimations grossières ont tendance à sous-estimer $r_{cool}$ dans les environnements de faible densité et ne parviennent pas à capturer la nature bimodale de l'effet dans les environnements de haute densité. Cela compromet les modèles subgrid qui reposent sur la densité locale pour décider s'ils doivent injecter de l'énergie thermique ou du moment.

• Seuils de formation stellaire : Lorsque l'on effectue la moyenne sur des échelles de masse grossières ($2000 M_\odot$), des régions avec des densités moyennes aussi basses que $10^{-2} \text{ cm}^{-3}$ peuvent présenter des taux de formation stellaire non nuls. Cela suggère que les seuils de densité élevés ($\rho_{th} \sim 10^3 - 10^4 \text{ cm}^{-3}$) utilisés dans les simulations à haute résolution pour mimer les nuages moléculaires peuvent être artificiellement restrictifs lorsqu'ils sont appliqués à des cellules de résolution grossière contenant des phases mixtes.

Signification et affirmations
Les auteurs concluent que les effets macroscopiques du feedback des SN ne sont pas uniformes mais dépendent fortement de la structure à petite échelle de l'ISM. Le papier affirme que :

1. L'auto-cohérence nécessite une haute résolution : Les simulations visant à modéliser de manière auto-cohérente les doubles rôles du feedback des SN (écoulements et suppression de la SF) doivent résoudre des échelles de masse de $\lesssim 100 M_\odot$ par cellule. En dessous de ce seuil, les canaux physiques distincts du feedback sont enchevêtrés et perdus.
2. Limites des modèles subgrid : Les paramètres couramment utilisés dans les modèles subgrid à résolution grossière, tels que le rayon de refroidissement des SN et les seuils de densité de SF, pourraient nécessiter une réévaluation. Ces modèles supposent souvent des conditions de l'ISM homogènes au sein d'une cellule, une hypothèse que ce papier démontre mener à des erreurs significatives dans la prédiction des résultats des SN.
3. Stratégies de modélisation alternatives : Pour les simulations qui ne peuvent atteindre la résolution requise, les auteurs suggèrent que les modèles efficaces doivent mimer les effets émergents des résultats à haute résolution plutôt que de compter uniquement sur les variables hydrodynamiques locales. Ils mettent en évidence les modèles de « vent découplé » (qui restreignent les écoulements au gaz de faible densité) et les approches multi-fluides comme des voies potentielles pour capturer ces effets sans résoudre explicitement la longueur de Jeans des amas denses.

Le papier ne prétend pas que les simulations à basse résolution sont inutiles, mais qu'elles ne peuvent pas déterminer de manière auto-cohérente le canal de feedback spécifique (écoulement vs suppression) d'un événement SN donné en se basant uniquement sur des informations de résolution grossière.