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

Technische Samenvatting: Drukgereguleerde, Feedback-gemoduleerde Stervorming als Subgrid-model

Probleemstelling
Grootschalige kosmologische simulaties van de vorming van sterrenstelsels (bijv. IllustrisTNG, EAGLE, SIMBA) reproduceren succesvol de statistische eigenschappen van populaties sterrenstelsels, maar kunnen de fysieke schalen van reusachtige moleculaire wolken (GMCs) of het meerfasige interstellaire medium (ISM) niet oplossen. Bijgevolg vertrouwen deze simulaties op empirisch gekalibreerde subgrid-voorschriften voor sterformatiesnelheden (SFR) en de effectieve toestandsvergelijking (EoS). Traditionele modellen gaan er doorgaans van uit dat de specifieke SFR schaalt met de gasdichtheid als $\dot{m}_*/m_g \propto \rho^{1/2}$ (inverse vrije-valtijd), waarbij een normalisatie wordt ingesteld op basis van waarnemingen in het lokale universum.

De auteurs identificeren twee primaire beperkingen in deze aanpak:

1. Gebrek aan Fysieke Koppeling: De empirische schaling mist een directe link met de fysieke processen die de stervorming reguleren, specifiek de balans tussen zwaartekracht en stellaire feedback.
2. Resolutie en Omgevingsgevoeligheid: Deze voorschriften falen vaak in regimes waar de gasdichtheid en druk aanzienlijk hoger zijn dan in het lokale universum (bijv. sterrenstelsels met een hoge roodverschuiving). Bovendien zijn modellen die vertrouwen op lokale volumetrische dichtheid gevoelig voor numerieke resolutie; als de verticale dikte van de gasschijf niet wordt opgelost, wordt de gemeten volumetrische dichtheid onnauwkeurig, wat leidt tot inconsistente stervormingsvoorspellingen over verschillende resoluties heen.

Methodologie
Het artikel introduceert een nieuw subgrid-raamwerk gebaseerd op de Pressure-Regulated, Feedback-Modulated (PRFM) theorie van stervorming (Ostriker & Kim 2022). In dit raamwerk past de SFR zich aan om het thermische en verticale dynamische evenwicht in het ISM te handhaven, waarbij de totale midplanedruk ($P_{tot}$) het gewicht van het ISM ($W$) balanceert en proportioneel is aan de SFR per oppervlakte ($\Sigma_{SFR}$) via een feedback-opbrengst ($\Upsilon_{tot}$).

De auteurs implementeren twee complementaire versies van dit model in de moving-mesh code Arepo, gekalibreerd met hoog-resolutie magnetohydrodynamische simulaties uit de TIGRESS-suite:

1. PRFM-vol (Volumetrisch): Ontworpen voor simulaties waarbij de verticale dikte van de gasschijf numeriek wordt opgelost (typisch massaresolutie $\sim 10^5 M_\odot$).

   • Het berekent de effectieve druk ($P_{eff}$) in elke cel met behulp van een TIGRESS-gekalibreerde EoS ($P_{tot} \propto n_H^{1.8}$).
   • De dynamische tijd ($t_{dyn}$) wordt afgeleid van lokale volumetrische dichtheden van gas, sterren en donkere materie.
   • De uitputtingstijd ($t_{dep}$) wordt ingesteld op basis van de lokale druk en een renormalisatiefactor ($R_f \approx 2$) om de verticale integratie van de volumetrische wet te compenseren.

2. PRFM-int (Geïntegreerd): Ontworpen voor grovere resoluties waar de verticale schaalhoogte van de schijf niet wordt opgelost (typisch massaresolutie $\gtrsim 10^6 M_\odot$).

   • In plaats van te vertrouwen op lokale volumetrische dichtheid, reconstrueert het de midplanedichtheid en -druk vanuit verticale evenwichtsoverwegingen met behulp van oppervlaktedichtheden ($\Sigma_g$, $\Sigma_*$) en het ISM-gewicht ($W$).
   • Het lost de evenwichtsschaalhoogte van het gas ($H_g$) op en leidt $t_{dep}$ direct af uit deze geïntegreerde grootheden, waardoor het gebruik van resolutiegevoelige volumetrische dichtheden in het SFR-voorschrift wordt vermeden.

Belangrijkste Bijdragen

• Fysisch Onderbouwde Kalibratie: Het model leidt de SFR en EoS af uit de PRFM-evenwichtsconditie ($W = P_{tot} = \Upsilon_{tot}\Sigma_{SFR}$), gekalibreerd tegen TIGRESS-simulaties in plaats van empirische fits aan lokale sterrenstelsels.
• Resolutie-robuustheid: De ontwikkeling van de PRFM-int aanpak adresseert specifiek het falen van volumetrische modellen bij lage resolutie door gebruik te maken van oppervlaktedichtheden en evenwichtsrestricties, die minder gevoelig zijn voor het onvermogen om de verticale structuur van de schijf op te lossen.
• Variabele Efficiëntie: In tegen tegenstelling tot traditionele modellen die een constante efficiëntie per vrije-valtijd aannemen, voorspelt het PRFM-model dat de stervormingsefficiëntie per dynamische tijd ($\epsilon_{dyn}$) toeneemt met de druk ($\epsilon_{dyn} \propto P_{tot}^{0.43}$).

Resultaten
De auteurs hebben beide implementaties getest in geïsoleerde, op de Melkweg achtige schijven, over massaresoluties van $10^5$, $10^6$, en $10^7 M_\odot$, waarbij ze deze vergeleken met het standaard IllustrisTNG-voorschrift.

• Overeenstemming bij Hoge Resolutie: Bij $10^5 M_\odot$ komen zowel PRFM-vol als PRFM-int nauw overeen met elkaar en met de TIGRESS-resultaten voor SFR's en uitputtingstijden. Beiden reproduceren succesvol de TIGRESS-gekalibreerde Ostriker-Kim relatie ($\Sigma_{SFR} \propto P_{tot}^{1.2}$).
• Resolutiegevoeligheid van Volumetrische Modellen: De PRFM-vol en TNG-modellen vertonen een significante gevoeligheid voor resolutie. Naarmate de resolutie grover wordt naar $10^6$ en $10^7 M_\odot$, leidt het onvermogen om de verticale schaalhoogte van de schijf op te lossen tot onnauwkeurige metingen van de volumetrische dichtheid. Dit veroorzaakt dat het PRFM-vol model de uitputtingstijden overschat en de hoge-resolutie SFR-geschiedenis niet slaagt te reproduceren.
• Robuustheid van PRFM-int: De PRFM-int implementatie blijft robuust over alle geteste resoluties. Het volgt succesvol de globale stervormingsgeschiedenis en oppervlaktedichtheidsverdelingen van de hoge-resolutie ($10^5 M_\odot$) PRFM-vol simulatie, zelfs bij $10^7 M_\odot$.
• Vergelijking met TNG: Vergeleken met het TNG-model voorspellen de PRFM-modellen kortere gasuitputtingstijden, met name in regio's met hoge druk en hoge dichtheid. De PRFM-modellen vertonen een veel breder dynamisch bereik in uitputtingstijden en een steilere afhankelijkheid van druk, terwijl het TNG-model een constante efficiëntie per vrije-valtijd behoudt.

Betekenis
Het artikel stelt dat de door PRFM afgeleide subgrid-voorschriften een fysisch onderbouwd en numeriek stabiel raamwerk bieden voor stervorming over het dynamische bereik van de vorming van sterrenstelsels.

• Voor Simulaties met Hoge Resolutie: Het PRFM-vol model biedt een realistischerere behandeling van de ISM-fysica dan traditionele empirische wetten, door de SFR direct te koppelen aan de dynamische balans van het sterrenstelsel.
• Voor Grootschalige Kosmologische Simulaties: Het PRFM-int model lost het cruciale probleem van resolutieafhankelijkheid op. Het stelt grootschalige kosmologische simulaties (die doorgaans de verticale schaalhoogte van schijven niet kunnen oplossen) in staat om fysisch gemotiveerde stervormingswetten te gebruiken die consistent blijven met hoge-resolutie fysica.
• Implicaties voor Hoge Roodverschuiving: De auteurs merken op dat de drukafhankelijke efficiëntie van het PRFM-model suggereert dat stervorming in sterrenstelsels met een hoge roodverschuiving (waar de drukken hoger zijn) aanzienlijk efficiënter kan zijn dan voorspeld door huidige TNG-achtige modellen, wat potentieel de interpretatie van lumineuze sterrenstelsels geobserveerd door JWST kan beïnvloeden.

Dit werk wordt gepresenteerd als een stap naar het doel van de "Learning the Universe Collaboration" om voorspellende, in plaats van empirisch afgestemde subgrid-modellen te ontwikkelen die kunnen worden bijgewerkt naarmate hoog-resolutie ISM-simulaties (zoals TIGRESS) een breder scala aan parameters (bijv. variërende metalliciteit) dekken.