Pressure-regulated feedback-modulated star formation as a subgrid model for galaxy formation simulations

Sintesi Tecnica: Modello Subgrid di Formazione Stellare Modulata dal Feedback e Regolata dalla Pressione

Problema
Le simulazioni cosmologiche a grande volume della formazione delle galassie (ad es. IllustrisTNG, EAGLE, SIMBA) riproducono con successo le proprietà statistiche delle popolazioni galattiche, ma non possono risolvere le scale fisiche delle nubi molecolari giganti (GMC) o del mezzo interstellare (ISM) multifase. Di conseguenza, queste simulazioni si affidano a prescrizioni subgrid empiricamente calibrate per i tassi di formazione stellare (SFR) e l'equazione di stato (EoS) effettiva. I modelli tradizionali assumono tipicamente che lo SFR specifico scali con la densità del gas come $\dot{m}_*/m_g \propto \rho^{1/2}$ (tempo di caduta libera inverso), con una normalizzazione impostata dalle osservazioni dell'universo locale.

Gli autori identificano due limitazioni primarie in questo approccio:

1. Mancanza di Connessione Fisica: La scalatura empirica manca di un legame diretto con i processi fisici che regolano la formazione stellare, specificamente l'equilibrio tra gravità e feedback stellare.
2. Sensibilità alla Risoluzione e all'Ambiente: Queste prescrizioni spesso falliscono in regimi in cui la densità del gas e la pressione sono significativamente più elevate rispetto all'universo locale (ad es. galassie ad alto redshift). Inoltre, i modelli basati sulla densità volumetrica locale sono sensibili alla risoluzione numerica; se lo spessore di scala del disco di gas non è risolto, la densità volumetrica misurata diventa inaccurata, portando a previsioni di formazione stellare inconsistenti tra diverse risoluzioni.

Metodologia
Il documento introduce un nuovo framework subgrid basato sulla teoria della formazione stellare Regolata dalla Pressione e Modulata dal Feedback (PRFM) (Ostriker & Kim 2022). In questo framework, l'SFR si adatta per mantenere l'equilibrio termico e dinamico verticale nell'ISM, dove la pressione totale del piano medio ($P_{tot}$) bilancia il peso dell'ISM ($W$) ed è proporzionale allo SFR per unità di area ($\Sigma_{SFR}$) tramite una resa di feedback ($\Upsilon_{tot}$).

Gli autori implementano due versioni complementari di questo modello nel codice a mesh mobile Arepo, calibrate utilizzando le simulazioni magnetoidrodinamiche ad alta risoluzione della suite TIGRESS:

1. PRFM-vol (Volumetrico): Progettato per simulazioni in cui lo spessore verticale del disco di gas è risolto numericamente (tipicamente risoluzione di massa $\sim 10^5 M_\odot$).

   • Calcola la pressione effettiva ($P_{eff}$) in ogni cella utilizzando una EoS calibrata da TIGRESS ($P_{tot} \propto n_H^{1.8}$).
   • Il tempo dinamico ($t_{dyn}$) è derivato dalle densità volumetriche locali di gas, stelle e materia oscura.
   • Il tempo di esaurimento ($t_{dep}$) è impostato sulla base della pressione locale e di un fattore di rinormalizzazione ($R_f \approx 2$) per tenere conto dell'integrazione verticale della legge volumetrica.

2. PRFM-int (Integrato): Progettato per risoluzioni più grossolane dove lo spessore di scala del disco è non risolto (tipicamente risoluzione di massa $\gtrsim 10^6 M_\odot$).

   • Invece di fare affidamento sulla densità volumetrica locale, ricostruisce la densità e la pressione del piano medio considerando l'equilibrio verticale tramite le densità superficiali ($\Sigma_g$, $\Sigma_*$) e il peso dell'ISM ($W$).
   • Risolve per lo spessore di scala del gas in equilibrio ($H_g$) e deriva $t_{dep}$ direttamente da queste quantità integrate, evitando l'uso di densità volumetriche sensibili alla risoluzione nella prescrizione dell'SFR.

Contributi Chiave

• Calibrazione Fondata sulla Fisica: Il modello deriva l'SFR e la EoS dalla condizione di equilibrio PRFM ($W = P_{tot} = \Upsilon_{tot}\Sigma_{SFR}$), calibrata contro le simulazioni TIGRESS piuttosto che su fit empirici di galassie locali.
• Robustezza della Risoluzione: Lo sviluppo dell'approccio PRFM-int affronta specificamente il fallimento dei modelli volumetrici a bassa risoluzione, utilizzando densità superficiali e vincoli di equilibrio, che sono meno sensibili all'incapacità di risolvere la struttura verticale del disco.
• Efficienza Variabile: A differenza dei modelli tradizionali che assumono un'efficienza costante per tempo di caduta libera, il modello PRFM prevede che l'efficienza della formazione stellare per tempo dinamico ($\epsilon_{dyn}$) aumenti con la pressione ($\epsilon_{dyn} \propto P_{tot}^{0.43}$).

Risultati
Gli autori hanno testato entrambe le implementazioni in simulazioni di dischi isolati simili alla Via Lattea attraverso risoluzioni di massa di $10^5$, $10^6$ e $10^7 M_\odot$, confrontandole con la prescrizione standard di IllustrisTNG.

• Accordo ad Alta Risoluzione: A $10^5 M_\odot$, sia PRFM-vol che PRFM-int concordano strettamente tra loro e con i risultati TIGRESS per gli SFR e i tempi di esaurimento. Entrambi recuperano con successo la relazione Ostriker-Kim calibrata da TIGRESS ($\Sigma_{SFR} \propto P_{tot}^{1.2}$).
• Sensibilità alla Risoluzione dei Modelli Volumetrici: Il PRFM-vol e i modelli TNG mostrano una significativa sensibilità alla risoluzione. Man mano che la risoluzione si fa più grossolana a $10^6$ e $10^7 M_\odot$, l'incapacità di risolvere lo spessore di scala del disco porta a misurazioni inaccurate della densità volumetrica. Ciò causa al modello PRFM-vol di sovrastimare i tempi di esaurimento e di non riuscire a corrispondere alla storia dell'SFR ad alta risoluzione.
• Robustezza di PRFM-int: L'implementazione PRFM-int rimane robusta attraverso tutte le risoluzioni testate. Essa segue con successo la storia globale della formazione stellare e le distribuzioni di densità superficiale della simulazione PRFM-vol ad alta risoluzione ($10^5 M_\odot$), anche a $10^7 M_\odot$.
• Confronto con TNG: Rispetto al modello TNG, i modelli PRFM prevedono tempi di esaurimento del gas più brevi, particolarmente nelle regioni ad alta pressione e alta densità. I modelli PRFM mostrano un intervallo dinamico molto più ampio nei tempi di esaurimento e una dipendenza più ripida dalla pressione, mentre il modello TNG mantiene un'efficienza costante per tempo di caduta libera.

Significatività
L'articolo sostiene che le prescrizioni subgrid derivate da PRFM forniscono un framework fisicamente fondato e numericamente stabile per la formazione stellare attraverso l'intervallo dinamico delle simulazioni di formazione galattica.

• Per Simulazioni ad Alta Risoluzione: Il modello PRFM-vol offre un trattamento più realistico della fisica dell'ISM rispetto alle leggi empiriche tradizionali, collegando l'SFR direttamente al bilancio dinamico della galassia.
• Per Simulazioni Cosmologiche su Grande Scala: Il modello PRFM-int risolve il problema critico della dipendenza dalla risoluzione. Permette alle simulazioni cosmologiche a box grandi (che tipicamente non possono risolvere gli spessori di scala dei dischi) di utilizzare leggi di formazione stellare motivate fisicamente che rimangono coerenti con la fisica ad alta risoluzione.
• Implicazioni per l'Alto Redshift: Gli autori notano che l'efficienza dipendente dalla pressione del modello PRFM suggerisce che la formazione stellare nelle galassie ad alto redshift (dove le pressioni sono più elevate) possa essere significativamente più efficiente rispetto a quanto previsto dagli attuali modelli tipo TNG, influenzando potenzialmente l'interpretazione delle galassie luminose osservate da JWST.

Il lavoro è presentato come un passo verso l'obiettivo della collaborazione "Learning the Universe" di sviluppare modelli subgrid predittivi, piuttosto che empiricamente tarati, che possano essere aggiornati man mano che le simulazioni ISM ad alta risoluzione (come TIGRESS) si espandono per coprire spazi di parametri più ampi (ad es. variazione della metallicità).