TRINITY: A coupled model of winds, radiation, and photoionised gas in molecular clouds. I. Methods and validation

Résumé technique : TRINITY : Un modèle couplé des vents, du rayonnement et du gaz photoionisé dans les nuages moléculaires – I. Méthodes et validation

Problème et contexte
Les nuages moléculaires géants (GMC) se dispersent dans un délai de 1 à 5 millions d'années (Myr) après l'émergence d'étoiles massives, généralement avant l'arrivée des premières supernovae à effondrement de cœur. Les relevés observationnels (par exemple, PHANGS, LEGUS) indiquent que les mécanismes de rétroaction pré-supernova — vents stellaires, pression de radiation et pression du gaz photoionisé ($P_{\text{H II}}$) — sont les principaux moteurs du nettoyage des nuages. Cependant, la dominance relative de ces canaux, leur chronologie et leur dépendance aux propriétés du nuage (masse, profil de densité, efficacité de formation stellaire) restent incertaines. Les outils semi-analytiques existants traitent souvent ces mécanismes de rétroaction de manière isolée ou s'appuient sur des profils de densité simplifiés, tandis que les simulations 3D sont coûteuses en calcul pour des relevés statistiques. Il existe un besoin d'un modèle direct efficace et auto-cohérent qui couple tous les canaux de rétroaction pré-supernova à diverses structures de nuages afin de cartographier la dominance de la rétroaction à travers l'espace des paramètres.

Méthodologie
Les auteurs présentent Trinity, un code à une dimension, sphériquement symétrique et en coquille mince, qui fait évoluer une bulle entraînée par la rétroaction et une coquille accrétée au sein d'un nuage moléculaire. Trinity succède au cadre warpfield mais introduit deux changements architecturaux critiques :

1. Pression explicite du gaz photoionisé : $P_{\text{H II}}$ est inclus explicitement dans l'équation du mouvement de la coquille, plutôt que d'être traité uniquement comme une condition aux limites ou ignoré.
2. Prescription de pilotage consciente de la phase : Le code sélectionne dynamiquement le terme de pression motrice en fonction de la phase d'évolution pour éviter le double comptage des forces :
   • Phase entraînée par l'énergie : La pression motrice est le maximum de la pression thermique de la bulle chaude ($P_b$) et de la pression du gaz photoionisé ($P_{\text{H II}}$).
   • Phase de transition : La pression motrice est le maximum de la $P_b$ restante et de la somme de $P_{\text{H II}}$ plus la pression dynamique du vent nouvellement injecté et des éjectas de supernova ($P_{\text{ram}}$).
   • Phase entraînée par la quantité de mouvement : Une fois que la bulle chaude se refroidit ($E_b \to 0$), la pression motrice est la somme $P_{\text{H II}} + P_{\text{ram}}$.

Le modèle résout un système couplé d'équations différentielles ordinaires pour le rayon de la bulle ($R_b$), la vitesse ($v_b$) et l'énergie thermique ($E_b$). Il intègre :

• Conditions initiales : Prend en charge des profils de densité uniformes, en loi de puissance par morceaux ($\rho \propto r^\alpha$) et de sphères de Bonnor–Ebert (BE).
• Modules physiques : Inclut les vents stellaires, les supernovae, la pression de radiation directe et reprocessée par la poussière, l'auto-gravité de la coquille et de l'amas, ainsi que l'évaporation conductive de la coquille vers la bulle chaude.
• Couplage rayonnement-coquille : Un calcul de structure hydrostatique de la coquille détermine la position du front d'ionisation ($R_{\text{if}}$), les fractions de rayonnement absorbées et l'épaisseur optique infrarouge de manière auto-cohérente.
• Refroidissement : Traite la transition de l'équilibre d'ionisation collisionnelle (CIE) vers des régimes non-CIE et calcule les pertes par refroidissement pour déterminer la transition énergie-quantité de mouvement.

Résultats clés
L'article valide le code contre des limites analytiques (vent de Weaver, expansion de Spitzer, chasse-neige) et des simulations warpfield, puis présente des résultats pour des nuages avec $M_{\text{cloud}} = 10^5–10^{6.5} M_\odot$, des densités de cœur $n_{\text{core}} = 10^3–10^4 \text{ cm}^{-3}$, et des efficacités de formation stellaire $\epsilon = 0.01–0.30$.

1. Impact du gaz photoionisé : L'inclusion de $P_{\text{H II}}$ modifie mesurablement la dynamique. Dans une simulation de validation de référence ($M_{\text{cloud}}=10^6 M_\odot, \epsilon=0.01$), l'inclusion de $P_{\text{H II}}$ augmente le rayon de la coquille d'environ 17 % à 10 Myr par rapport à une simulation équivalente warpfield où $P_{\text{H II}}$ est désactivé. $P_{\text{H II}}$ reste dynamiquement important même après la transition vers une évolution dominée par la quantité de mouvement.
2. Dépendance au profil de densité : La structure du nuage dicte le résultat même lorsque les populations stellaires sont fixes. Pour une $M_{\text{cloud}}$, une $n_{\text{core}}$ et un $\epsilon$ fixes :
   • Les profils homogènes ($\rho \propto r^0$) et peu pentus ($\rho \propto r^{-1}$) se réeffondrent.
   • Un profil raide ($\rho \propto r^{-2}$) permet à la coquille d'accélérer le long du gradient de densité et de continuer à se dilater.
   • Les sphères de Bonnor–Ebert se réeffondrent dans des cas où des nuages homogènes pourraient se disperser, ou se dispersent significativement plus tard.
3. Échelles de temps de dispersion : Dans une grille de modèles où la dispersion est garantie, les nuages de Bonnor–Ebert se dispersent environ 55 % plus tard que les nuages homogènes. Ce décalage systématique provient du fait que les structures BE distribuent la même masse sur un rayon physique plus grand, obligeant la coquille à parcourir une distance plus longue.
4. Dominance pré-supernova : Dans l'espace des paramètres testé, la dispersion ou le réeffondrement du nuage est déterminé par la rétroaction pré-supernova (vents, rayonnement, $P_{\text{H II}}$) avant l'arrivée de la première supernova à effondrement de cœur (généralement > 3 Myr).
5. Échappement LyC : Le code suit la fraction d'échappement du continuum de Lyman ($f_{\text{LyC}}^{\text{esc}}$) comme une « horloge » distincte du nettoyage dynamique, distinguant le moment où la rétroaction modifie la dynamique de la coquille du moment où les photons commencent à s'échapper dans le milieu interstellaire diffus.

Signification et affirmations
Les auteurs affirment que Trinity fournit un cadre hautement efficace et interprétable pour cartographier la dominance de la rétroaction à travers l'espace des paramètres des nuages. La conclusion principale est que la dominance de la rétroaction ne peut pas être déduite de la population stellaire seule. Le résultat dépend de manière critique de :

• La façon dont les canaux de rétroaction se couplent à la coquille (spécifiquement le rôle de $P_{\text{H II}}$).
• La chronologie de la transition de refroidissement de la bulle chaude.
• La structure de densité du nuage natal.

L'article affirme que l'hypothèse d'un nuage homogène canonique ou d'une unique pente de loi de puissance introduit des erreurs systématiques dans les temps d'émergence et les seuils de dispersion déduits. En offrant un outil qui couple de manière auto-cohérente les vents, le rayonnement et la photoionisation avec des profils de densité flexibles, Trinity permet l'étude des régions H II résolues et fournit des prescriptions sous-maille pour les simulations à l'échelle galactique concernant l'échappement des photons ionisants et la dispersion des nuages.

Limites
Les auteurs notent explicitement que la symétrie sphérique 1D moyenne les structures angulaires (par exemple, les écoulements champagne, les ruptures de cloques) et ne peut pas résoudre les instabilités comme Kelvin–Helmholtz ou Rayleigh–Taylor qui pourraient amplifier le refroidissement. Le traitement de $P_{\text{H II}}$ utilise un rayon unique pour le système bulle-front d'ionisation plutôt que d'évoluer $R_b$ et $R_{\text{if}}$ comme des variables dynamiques distinctes. De plus, le modèle néglige actuellement le support magnétique, les gradients de pression des rayons cosmiques et l'évolution de la fraction de couverture (coquilles fuyardes). Ceux-ci sont identifiés comme des cibles pour un développement futur.