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

Technische Samenvatting: TRINITY: Een gekoppeld model van winden, straling en foto-ioniseerd gas in moleculaire wolken – I. Methoden en validatie

Probleem en Context
Gigantische moleculaire wolken (GMW's) verspreiden zich binnen 1–5 Myr na het ontstaan van massieve sterren, doorgaans voordat de eerste kerninstortings-supernova's arriveren. Observatiesurveys (bijv. PHANGS, LEGUS) geven aan dat feedbackmechanismen voorafgaand aan supernova's – sterrenwinden, stralingsdruk en de druk van foto-ioniseerd gas ($P_{\text{H II}}$) – de primaire drijvers zijn van het zuiveren van wolken. De relatieve dominantie van deze kanalen, hun timing en hun afhankelijkheid van wolkkenmerken (massa, dichtheidsprofiel, sterrenvormingsefficiëntie) blijven echter onzeker. Bestaande semi-analytische hulpmiddelen behandelen deze feedbackmechanismen vaak geïsoleerd of vertrouwen op vereenvoudigde dichtheidsprofielen, terwijl 3D-simulaties voor statistische surveys computationeel duur zijn. Er is behoefte aan een efficiënt, zelfconsistent voorwaarts model dat alle feedbackkanalen voorafgaand aan supernova's koppelt aan diverse wolkstructuren om de feedback-dominantie over de parameter ruimte in kaart te brengen.

Methodologie
De auteurs introduceren Trinity, een eendimensionale, sferisch symmetrische dunne-schaalcode die een feedback-gedreven bel en opgeveegde schaal evolueert binnen een moleculaire wolk. Trinity volgt het warpfield-kader, maar introduceert twee kritieke architecturale wijzigingen:

1. Expliciete Druk van Foto-ioniseerd Gas: $P_{\text{H II}}$ wordt expliciet opgenomen in de bewegingsvergelijking van de schaal, in plaats van uitsluitend als randvoorwaarde te worden behandeld of genegeerd.
2. Fase-bewuste Drijvende Voorschrift: De code selecteert dynamisch de drijvende drukterm op basis van de evolutionaire fase om dubbeltelling van krachten te voorkomen:
   • Energie-gedreven fase: De drijvende druk is het maximum van de thermische druk van de hete bel ($P_b$) en de druk van foto-ioniseerd gas ($P_{\text{H II}}$).
   • Overgangsfase: De drijvende druk is het maximum van de resterende $P_b$ en de som van $P_{\text{H II}}$ plus de ram-druk van nieuw geïnjecteerde wind en supernova-uitwerpselen ($P_{\text{ram}}$).
   • Momentum-gedreven fase: Zodra de hete bel afkoelt ($E_b \to 0$), is de drijvende druk de som $P_{\text{H II}} + P_{\text{ram}}$.

Het model lost een gekoppeld systeem van gewone differentiaalvergelijkingen op voor de belstraal ($R_b$), snelheid ($v_b$) en thermische energie ($E_b$). Het omvat:

• Initiële Voorwaarden: Ondersteunt uniforme, stuksgewijze machts-wet ($\rho \propto r^\alpha$) en Bonnor–Ebert (BE) bol-dichtheidsprofielen.
• Fysica-modules: Omvat sterrenwinden, supernova's, directe en door stof herverwerkte stralingsdruk, zelfzwaartekracht van de schaal en de cluster, en geleidend verdamping van de schaal naar de hete bel.
• Koppeling Straling-Schaal: Een hydrostatische schaal-structuurberekening bepaalt de positie van de ionisatiefront ($R_{\text{if}}$), geabsorbeerde stralingsfracties en de infrarode optische diepte op een zelfconsistente manier.
• Koeling: Behandelt de overgang van collisionele ionisatie-evenwicht (CIE) naar niet-CIE-regimes en berekent koelingsverliezen om de overgang van energie naar momentum te bepalen.

Belangrijkste Resultaten
Het artikel valideert de code tegen analytische limieten (Weaver-wind, Spitzer-uitbreiding, sneeuwploeg) en warpfield-runs, en presenteert vervolgens resultaten voor wolken met $M_{\text{cloud}} = 10^5–10^{6.5} M_\odot$, kern-dichtheden $n_{\text{core}} = 10^3–10^4 \text{ cm}^{-3}$ en sterrenvormingsefficiënties $\epsilon = 0.01–0.30$.

1. Invloed van Foto-ioniseerd Gas: Het opnemen van $P_{\text{H II}}$ verandert de dynamica meetbaar. In een referentie-validatierun ($M_{\text{cloud}}=10^6 M_\odot, \epsilon=0.01$) verhoogt het opnemen van $P_{\text{H II}}$ de straal van de schaal met $\sim 17\%$ op 10 Myr vergeleken met een warpfield-equivalente run waarbij $P_{\text{H II}}$ is uitgeschakeld. $P_{\text{H II}}$ blijft dynamisch belangrijk zelfs na de overgang naar momentum-gedomineerde evolutie.
2. Afhankelijkheid van Dichtheidsprofiel: De wolkstructuur bepaalt het resultaat, zelfs wanneer sterpopulaties vaststaan. Voor een vaste $M_{\text{cloud}}$, $n_{\text{core}}$ en $\epsilon$:
   • Homogene ($\rho \propto r^0$) en ondiepe ($\rho \propto r^{-1}$) profielen her-collaberen.
   • Een steil profiel ($\rho \propto r^{-2}$) laat de schaal versnellen langs het dichtheidsgradiënt en doorgaan met uitbreiden.
   • Bonnor–Ebert bollen her-collaberen in gevallen waar homogene wolken zouden verspreiden, of verspreiden aanzienlijk later.
3. Verspreidingstijdschalen: In een rooster van modellen waarbij verspreiding gegarandeerd is, verspreiden Bonnor–Ebert-wolken zich ongeveer 55% later dan homogene wolken. Deze systematische verschuiving ontstaat omdat BE-structuren dezelfde massa over een grotere fysieke straal verdelen, waardoor de schaal verder moet reizen.
4. Dominantie voorafgaand aan Supernova: In de geteste parameter ruimte wordt de verspreiding of her-collaps van de wolk bepaald door feedback voorafgaand aan supernova's (winden, straling, $P_{\text{H II}}$) voordat de eerste kerninstortings-supernova arriveert (doorgaans $>3$ Myr).
5. LyC-ontsnapping: De code volgt het ontsnappingspercentage van het Lyman-continuüm (LyC) ($f_{\text{LyC}}^{\text{esc}}$) als een aparte "klok" dan dynamische zuivering, waardoor het onderscheid wordt gemaakt tussen het tijdstip waarop feedback de schaaldynamica verandert en het tijdstip waarop fotonen beginnen te lekken in het diffuse ISM.

Betekenis en Beweringen
De auteurs beweren dat Trinity een zeer efficiënt, interpreteerbaar kader biedt voor het in kaart brengen van feedback-dominantie over de parameter ruimte van wolken. De belangrijkste conclusie is dat feedback-dominantie niet kan worden afgeleid uit de sterpopulatie alleen. Het resultaat is kritiek afhankelijk van:

• Hoe feedbackkanalen koppelen aan de schaal (specifiek de rol van $P_{\text{H II}}$).
• De timing van de overgang van afkoeling van de hete bel.
• De dichtheidsstructuur van de geboortewolk.

Het artikel stelt dat het aannemen van een canonieke homogene wolk of een enkele machts-wet-helling systematische fouten introduceert in afgeleide ontstaanstijden en verspreidingsdrempels. Door een hulpmiddel aan te bieden dat winden, straling en foto-ionisatie op een zelfconsistente manier koppelt met flexibele dichtheidsprofielen, stelt Trinity het bestuderen van opgeloste H II-regio's mogelijk en biedt het sub-grid voorschriften voor simulaties op galaxie-schaal met betrekking tot lekken van ioniserende fotonen en wolkverspreiding.

Beperkingen
De auteurs merken expliciet op dat de 1D sferische symmetrie middelt over hoekstructuren (bijv. champagne-stromen, blister-uitbraken) en instabiliteiten zoals Kelvin–Helmholtz of Rayleigh–Taylor die de koeling kunnen versterken, niet kan oplossen. De behandeling van $P_{\text{H II}}$ maakt gebruik van één enkele straal voor het systeem van bel-ionisatiefront in plaats van $R_b$ en $R_{\text{if}}$ te evolueren als aparte dynamische variabelen. Bovendien negeert het model momenteel magnetische ondersteuning, stralingsdrukgradiënten van kosmische straling en de evolutie van het dekkingspercentage (lekende schalen). Deze worden geïdentificeerd als doelen voor toekomstige ontwikkeling.