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

Sintesi Tecnica: TRINITY: Un modello accoppiato di venti, radiazione e gas fotoionizzato nelle nubi molecolari – I. Metodi e validazione

Problema e Contesto
Le nubi molecolari giganti (GMC) si disperdono entro 1–5 Myr dalla comparsa di stelle massicce, tipicamente prima dell'arrivo delle prime supernove a collasso del nucleo. I sondaggi osservativi (ad es. PHANGS, LEGUS) indicano che i meccanismi di feedback pre-supernova – venti stellari, pressione di radiazione e pressione del gas fotoionizzato ($P_{\text{H II}}$) – sono i principali responsabili della dispersione delle nubi. Tuttavia, la dominanza relativa di questi canali, la loro tempistica e la loro dipendenza dalle proprietà della nube (massa, profilo di densità, efficienza di formazione stellare) rimangono incerte. Gli strumenti semi-analitici esistenti spesso trattano questi meccanismi di feedback in isolamento o si basano su profili di densità semplificati, mentre le simulazioni 3D sono computazionalmente costose per sondaggi statistici. È necessario un modello forward efficiente e coerente che accoppi tutti i canali di feedback pre-supernova a strutture di nube diverse per mappare la dominanza del feedback attraverso lo spazio dei parametri.

Metodologia
Gli autori presentano Trinity, un codice a guscio sottile unidimensionale e sfericamente simmetrico che evolve una bolla guidata dal feedback e un guscio spazzato via all'interno di una nube molecolare. Trinity succede al framework warpfield ma introduce due cambiamenti architetturali critici:

1. Pressione Esplicita del Gas Fotoionizzato: $P_{\text{H II}}$ è inclusa esplicitamente nell'equazione del moto del guscio, anziché essere trattata esclusivamente come condizione al contorno o ignorata.
2. Prescrizione di Guida Consapevole della Fase: Il codice seleziona dinamicamente il termine di pressione di guida in base alla fase evolutiva per evitare il doppio conteggio delle forze:
   • Fase guidata dall'energia: La pressione di guida è il massimo della pressione termica della bolla calda ($P_b$) e della pressione del gas fotoionizzato ($P_{\text{H II}}$).
   • Fase di transizione: La pressione di guida è il massimo della $P_b$ residua e della somma di $P_{\text{H II}}$ più la pressione dinamica del vento appena iniettato e degli ejecta delle supernove ($P_{\text{ram}}$).
   • Fase guidata dalla quantità di moto: Una volta che la bolla calda si raffredda ($E_b \to 0$), la pressione di guida è la somma $P_{\text{H II}} + P_{\text{ram}}$.

Il modello risolve un sistema accoppiato di equazioni differenziali ordinarie per il raggio della bolla ($R_b$), la velocità ($v_b$) e l'energia termica ($E_b$). Incorpora:

• Condizioni Iniziali: Supporta profili di densità uniformi, a legge di potenza a tratti ($\rho \propto r^\alpha$) e sfere di Bonnor–Ebert (BE).
• Moduli Fisici: Include venti stellari, supernove, pressione di radiazione diretta e riprocessata dalla polvere, autogravità del guscio e dell'ammasso, ed evaporazione conduttiva dal guscio verso la bolla calda.
• Accoppiamento Radiazione-Guscio: Un calcolo della struttura idrostatica del guscio determina la posizione del fronte di ionizzazione ($R_{\text{if}}$), le frazioni di radiazione assorbita e la profondità ottica infrarossa in modo coerente.
• Raffreddamento: Tratta la transizione dall'equilibrio di ionizzazione collisionale (CIE) a regimi non-CIE e calcola le perdite per raffreddamento per determinare la transizione da energia a quantità di moto.

Risultati Chiave
Il documento valida il codice contro limiti analitici (vento di Weaver, espansione di Spitzer, snowplow) ed esecuzioni warpfield, quindi presenta risultati per nubi con $M_{\text{cloud}} = 10^5–10^{6.5} M_\odot$, densità del nucleo $n_{\text{core}} = 10^3–10^4 \text{ cm}^{-3}$ ed efficienze di formazione stellare $\epsilon = 0.01–0.30$.

1. Impatto del Gas Fotoionizzato: L'inclusione di $P_{\text{H II}}$ altera in modo misurabile la dinamica. In una esecuzione di validazione fiduciale ($M_{\text{cloud}}=10^6 M_\odot, \epsilon=0.01$), l'inclusione di $P_{\text{H II}}$ aumenta il raggio del guscio di $\sim 17\%$ a 10 Myr rispetto a un'esecuzione equivalente warpfield in cui $P_{\text{H II}}$ è disabilitata. $P_{\text{H II}}$ rimane dinamicamente importante anche dopo la transizione all'evoluzione dominata dalla quantità di moto.
2. Dipendenza dal Profilo di Densità: La struttura della nube determina l'esito anche quando le popolazioni stellari sono fisse. Per una $M_{\text{cloud}}$, $n_{\text{core}}$ e $\epsilon$ fisse:
   • I profili omogenei ($\rho \propto r^0$) e poco ripidi ($\rho \propto r^{-1}$) collassano nuovamente.
   • Un profilo ripido ($\rho \propto r^{-2}$) permette al guscio di accelerare lungo il gradiente di densità e continuare ad espandersi.
   • Le sfere di Bonnor–Ebert collassano nuovamente in casi in cui nubi omogenee potrebbero disperdersi, o si disperdono significativamente più tardi.
3. Scale Temporali di Dispersione: In una griglia di modelli in cui la dispersione è garantita, le nubi di Bonnor–Ebert si disperdono circa il 55% più tardi rispetto alle nubi omogenee. Questo spostamento sistematico nasce perché le strutture BE distribuiscono la stessa massa su un raggio fisico maggiore, richiedendo al guscio di percorrere una distanza maggiore.
4. Dominanza Pre-Supernova: Nello spazio dei parametri testato, la dispersione o il collasso della nube è determinato dal feedback pre-supernova (venti, radiazione, $P_{\text{H II}}$) prima dell'arrivo della prima supernova a collasso del nucleo (tipicamente $>3$ Myr).
5. Fuga LyC: Il codice traccia la frazione di fuga del continuum di Lyman ($f_{\text{LyC}}^{\text{esc}}$) come un "orologio" separato dalla dispersione dinamica, distinguendo il momento in cui il feedback altera la dinamica del guscio dal momento in cui i fotoni iniziano a fuoriuscire nel mezzo interstellare diffuso.

Significato e Affermazioni
Gli autori affermano che Trinity fornisce un framework altamente efficiente e interpretabile per mappare la dominanza del feedback attraverso lo spazio dei parametri delle nubi. La conclusione principale è che la dominanza del feedback non può essere dedotta dalla sola popolazione stellare. L'esito dipende criticamente da:

• Come i canali di feedback si accoppiano al guscio (in particolare il ruolo di $P_{\text{H II}}$).
• La tempistica della transizione di raffreddamento della bolla calda.
• La struttura di densità della nube natale.

Il documento afferma che assumere una nube omogenea canonica o un singolo pendenza a legge di potenza introduce errori sistematici nei tempi di comparsa e nelle soglie di dispersione dedotti. Offrendo uno strumento che accoppia in modo coerente venti, radiazione e fotoionizzazione con profili di densità flessibili, Trinity permette lo studio di regioni H II risolte e fornisce prescrizioni sub-griglia per simulazioni su scala galattica riguardanti la fuga di fotoni ionizzanti e la dispersione delle nubi.

Limitazioni
Gli autori notano esplicitamente che la simmetria sferica 1D media sulle strutture angolari (ad es. flussi champagne, rotture a bolle) e non può risolvere instabilità come quelle di Kelvin–Helmholtz o Rayleigh–Taylor che potrebbero aumentare il raffreddamento. Il trattamento di $P_{\text{H II}}$ utilizza un singolo raggio per il sistema bolla-fronte di ionizzazione anziché evolvere $R_b$ e $R_{\text{if}}$ come variabili dinamiche separate. Inoltre, il modello attualmente trascura il supporto magnetico, i gradienti di pressione dei raggi cosmici e l'evoluzione della frazione di copertura (gusci permeabili). Questi sono identificati come obiettivi per sviluppi futuri.