The Cosmic Baryon Cycle in IllustrisTNG: flows of mass, energy, and metals

Dit artikel analyseert de stromen van massa, energie en metalen door het interstellaire en circumgalactische medium in de IllustrisTNG100-simulaties om de dominante terugkoppingsmechanismen te identificeren en te verklaren hoe deze de groei en uitdoving van sterrenstelsels reguleren.

De kern

De Kosmische Baryoncyclus: Hoe Sterrenstelsels Aderen, Ademen en Verouderen

Stel je voor dat het heelal een enorme, levende stad is. In deze stad zijn de sterrenstelsels de gebouwen. Maar hoe groeien deze gebouwen? Hoe krijgen ze nieuwe materialen, en waarom stoppen sommige met groeien terwijl andere blijven bloeien?

Dit artikel van Yossi Oren en zijn collega's is als een gedetailleerde inspectie van de afvoer- en aanvoersystemen van deze kosmische stad. Ze kijken naar de "TNG100"-simulatie, een superkrachtige computerwereld die nadoet hoe het heelal zich gedraagt. Hun doel? Begrijpen hoe massa (stof), energie (kracht) en metalen (zware elementen) binnenkomen en weer vertrekken bij sterrenstelsels.

Hier is de uitleg in simpele taal, met een paar creatieve vergelijkingen.

1. De twee grote problemen: Te veel of te weinig sterren

In de theorie zou er veel meer gas moeten zijn om sterren te maken dan er eigenlijk is. Als we alleen naar de zwaartekracht kijken, zouden sterrenstelsels gigantisch moeten zijn. Maar dat zijn ze niet. Er moet dus een rem zijn die voorkomt dat al dat gas in sterren verandert.

De auteurs noemen dit "feedback". Het is alsof een sterrenstelsel twee manieren heeft om te reageren op te veel gas:

• De "Uitwerpsel"-manier (Ejectief): Het stelsel gooit het gas eruit, zoals een vuilniswagen die zijn lading leegt.
• De "Sluit de Deur"-manier (Preventief): Het stelsel houdt het gas buiten de deur, zodat het nooit binnenkomt om sterren te maken.

2. De twee hoofdpersoonnen: Sterren en Zwarte Gaten

Er zijn twee krachten die deze feedback regelen:

1. Supernova's (Sterren): Wanneer zware sterren ontploffen, blazen ze als een gigantische vuilniswagen het gas weg. Dit werkt vooral goed in kleine sterrenstelsels.
2. Actieve Galactische Kernen (AGN / Zwarte Gaten): In het centrum van grote sterrenstelsels zitten superzware zwarte gaten. Als ze eten (gas opzuigen), kunnen ze twee dingen doen:
   • Thermische mode: Ze verwarmen het gas (alsof je de verwarming op 100 zet), zodat het niet kan afkoelen en sterren kan vormen.
   • Kinetische mode: Ze schieten krachtige stralen (jets) weg die het gas letterlijk uit het stelsel blazen. Dit is de "zware artillerie" die grote stelsels tot rust brengt.

3. De Reis van het Gas: Van de rand tot het centrum

De auteurs kijken naar twee plekken in het sterrenstelsel:

• De ISM-schaal (Het centrum): Dit is de "woonkamer" van het sterrenstelsel, waar de sterren en het koude gas zitten.
• De Halo-schaal (De rand): Dit is de "tuin" of het "omheining" ver weg van het centrum. Hier zit het warme gas dat nog niet helemaal binnen is gekomen.

Wat ontdekten ze?

• In het jonge heelal (hoge roodshift, $z > 2$):
  Het is een drukke bouwplaats. Er komt veel nieuw gas binnen (aanvoer), en er wordt ook wel wat weggeblazen, maar de aanvoer wint het. De sterrenstelsels groeien snel. De "tuin" (halo) zit vol met gas dat binnenkomt.

  • Vergelijking: Het is alsof er een enorme kraan openstaat die het huis vult met water. De afvoer is er ook, maar de kraan wint.

• In het oude heelal (lage roodshift, $z < 2$):
  Dan gebeurt er iets interessants. In de zware sterrenstelsels schakelt het zwarte gat over op de "kinetische mode". Het begint het gas hard weg te blazen.

  • Het evenwicht: De aanvoer en de afvoer komen bijna in evenwicht. Het huis groeit niet meer snel. Dit is het moment waarop sterrenstelsels "stopten met groeien" (quenching).
  • De "Fontein": In kleinere stelsels gebeurt er iets grappigs. Het gas wordt wel weggeblazen, maar het is niet zwaar genoeg om het heelal in te vliegen. Het koelt af en valt weer terug naar het stelsel. Dit noemen ze een "Galactische Fontein". Het gas wordt hergebruikt, net als regenwater dat in een emmer valt en weer verdampt.

4. De Belangrijke Verschillen tussen Massa, Energie en Metalen

• Massa (Het stof):
  In kleine stelsels is de "lading" (hoeveel gas eruit gaat per ster die erin gaat) groot. Naarmate stelsels zwaarder worden, wordt deze lading kleiner. Maar zodra het zwarte gat de macht overneemt, wordt de lading weer enorm groot, omdat het zwarte gat alles wegblaast.
• Energie (De kracht):
  Dit is verrassend: de hoeveelheid energie die het stelsel verliest of wint, is bijna altijd hetzelfde, ongeacht hoe groot het stelsel is of hoe oud het heelal is. Het is alsof de "motor" altijd even hard draait, maar de "wielomvang" (de massa) verandert.
• Metalen (De zware elementen):
  Metalen (zoals ijzer, koolstof, etc.) worden gemaakt in sterren. Als gas het stelsel verlaat, neemt het deze metalen mee.
  • In kleine stelsels is het uitgaande gas erg rijk aan metalen (het is puur stoffen uit het stelsel).
  • In grote stelsels, waar het zwarte gat blaast, wordt het gas "verdund" met armere gas uit de omgeving. Het is alsof je een glas sinaasappelsap (rijk aan metalen) opent en er een emmer water bij doet; de smaak wordt zwakker.

5. Waarom is dit belangrijk?

Deze studie is als een handleiding voor bouwers.
Vroeger maakten wetenschappers modellen (simulaties) van sterrenstelsels, maar ze wisten niet precies hoe ze de "feedback" (de remmen en versnellingen) moesten instellen. Ze deden het op gevoel.

Nu hebben Oren en zijn team de exacte cijfers geleverd:

• Hoeveel gas gaat eruit?
• Hoeveel energie gaat eruit?
• Hoe ver reist het?
• Wat gebeurt er als het zwarte gat wakker wordt?

Met deze gegevens kunnen andere wetenschappers hun modellen veel beter maken. Ze kunnen nu zeggen: "Oké, we moeten de afvoer zo instellen dat hij op dit moment in de tijd precies deze hoeveelheid gas wegblaast."

Conclusie in één zin

Sterrenstelsels zijn als levende organismen die ademen: ze zuigen gas in, bouwen sterren, en als ze te groot worden, blazen hun zwarte gaten het gas eruit om te voorkomen dat ze uit elkaar vallen, waardoor ze uiteindelijk rustig en "dood" (gequenched) worden, terwijl kleinere stelsels blijven groeien door een continue cyclus van uitblazen en terugvallen (de fontein).

Deze paper geeft ons de blauwdruk van die ademhaling, zodat we beter begrijpen waarom ons heelal er precies zo uitziet als het nu doet.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "The Cosmic Baryon Cycle in IllustrisTNG: Flows of Mass, Energy, and Metals" in het Nederlands.

Probleemstelling

De evolutie van sterrenstelsels binnen het $\Lambda$CDM-model vereist feedbackmechanismen om de waargenomen lage "baryon-conversie-efficiëntie" (het percentage baryonische massa dat in sterren wordt omgezet) te verklaren. Zonder feedback zouden zowel lage- als hoge-massa stelsels overproduktie van sterren vertonen. De centrale vraag is hoe feedback (voornamelijk door supernova's [SN] en actieve galactische kernen [AGN]) de instroom en uitstroom van massa, energie en metalen reguleert op verschillende schalen (ISM en halo).

Hoewel numerieke simulaties zoals IllustrisTNG deze processen modelleren via "subgrid"-recipes, blijft de gedetailleerde kwantificering van de baryoncyclus (in- en uitstroom van massa, energie en metalen) als functie van halo-massa en roodverschuiving onvoldoende onderzocht. Bestaande semi-analytische modellen (SAMs) en regulator-modellen missen vaak gedetailleerde data over energiestromen en de specifieke afhankelijkheid van feedbackmechanismen, wat nodig is om deze modellen te verbeteren.

Methodologie

De auteurs analyseren de TNG100-simulatie uit de IllustrisTNG-suite, een kosmologische hydrodynamische simulatie met een kubusgrootte van 110,7 Mpc en hoge resolutie.

• Steekproef: Er zijn 9.522 halo's geanalyseerd over 10 roodverschuivingen ($z = 0, 0.1, 0.2, 0.5, 1, 2, 4, 6, 8, 10$) en een bereik van viriale massa's ($M_{vir}$) van $10^{10}$tot$10^{14} M_\odot$.
• Schalen: De stromen worden gemeten op twee specifieke schalen:
  1. ISM-schaal: Een schaal van $0,05 R_{vir}$tot$0,15 R_{vir}$ (grens tussen interstellair medium en circumgalactisch medium).
  2. Halo-schaal: Een schaal van $0,95 R_{vir}$tot$1,05 R_{vir}$ (viriale straal).
• Feedback-dominantie: Om de dominante feedbackmechanismen te identificeren, berekenen de auteurs de geïntegreerde energie-injectie van SN, thermische AGN-feedback en kinetische AGN-feedback over tijdsintervallen die overeenkomen met de reistijd van uitstroom naar de viriale straal.
• Gedefinieerde Grootheden:
  • Massa-, Energie- en Metaalstroom: Berekend op basis van de radiale snelheid en massa van deeltjes.
  • Loading Factors: De verhouding tussen uitstroom (massa, energie, metalen) en de stervormingsrate (SFR) of SN-energie-injectie.
  • Metalen Verrijking: De verhouding tussen de metalenconcentratie in de stroming en die van het reservoir waaruit deze komt.

Belangrijkste Bijdragen

1. Gedetailleerde Kwantificering: Voor het eerst worden in TNG100 simultaan massa-, energie- en metaalstromen gemeten op zowel ISM- als halo-schalen, met een fijne binningsindeling in massa en roodverschuiving.
2. Scheiding van Feedback-mechanismen: De studie onderscheidt systematisch tussen halo's gedomineerd door SN-feedback, thermische AGN-feedback en kinetische AGN-feedback, in plaats van ze te aggregeren.
3. Energiebalans: In tegenstelling tot traditionele SAM's die zich vaak beperken tot massa en metalen, wordt hier expliciet de rol van energiestromen in de regulatie van koeling en accretie onderzocht.
4. Validatie voor SAM's: De resultaten dienen als een benchmark voor de ontwikkeling van nieuwe semi-analytische modellen die de baryoncyclus fysiek beter moeten kunnen nabootsen.

Resultaten

1. Dominantie van Feedback en Overgang:

• Bij hoge roodverschuiving ($z \gtrsim 2$) en lage massa's is SN-feedback dominant.
• Er is een duidelijke overgang naar AGN-dominantie bij $M_{vir} \approx 10^{12} M_\odot$ bij $z \approx 2$.
• Thermische AGN-feedback (hoge accretie) injecteert veel energie, maar drijft geen significante uitstroom aan en remt stervorming niet effectief.
• Kinetische AGN-feedback (lage accretie, "jet-mode") is de dominante drijvende kracht voor uitstroom en quenching in massieve halo's ($M_{vir} \gtrsim 10^{12} M_\odot$) bij lage roodverschuiving.

2. Massa- en Energiestromen:

• Hoge Roodverschuiving ($z \gtrsim 2$): Sterke netto instroom van massa op halo-schaal. Op ISM-schaal is de instroom ook positief, maar zwakker. Er is weinig bewijs voor preventieve feedback op halo-schaal.
• Lage Roodverschuiving ($z \lesssim 2$): In massieve halo's balanceren instroom en uitstroom elkaar bijna op ISM-schaal, wat leidt tot quenching. Op halo-schaal treden er significante netto uitstroom op bij de overgang naar kinetische AGN-feedback.
• Preventieve Feedback: Sterk aanwezig in halo's met $M_{vir} \approx 2 \times 10^{13} M_\odot$ (voornamelijk door kinetische AGN), waar verse IGM-gas de halo binnendringt maar niet het ISM bereikt.
• Recycling ("Galactic Fountains"): In halo's met $M_{vir} \sim 10^{11} M_\odot$ is er een snelle recycling van uitgestoten gas dat terugvalt op het stelsel, wat leidt tot een instroom op ISM-schaal die de instroom op halo-schaal overtreft.

3. Loading Factors:

• Massa Loading ($\eta_M$): Neemt af met toenemende halo-massa en toenemende roodverschuiving. Bij $z=0$ kan $\eta_M$ oplopen tot $\sim 100$ voor lage massa's. Bij de overgang naar kinetische AGN-feedback stijgt de gemeten loading factor dramatisch (tot >1000), wat aangeeft dat SN's geen rol meer spelen in het aandrijven van uitstroom.
• Energie Loading ($\eta_E$): Blijft bijna onafhankelijk van zowel massa als roodverschuiving (waarden rond 0,5–1,0), wat suggereert dat de specifieke energie van de wind constant blijft.
• Metaal Loading ($\eta_Z$): Is vrijwel onafhankelijk van de halo-massa maar neemt toe bij afnemende roodverschuiving.

4. Metalen en Verrijking:

• Uitstroom op ISM-schaal heeft een metalenverrijking die overeenkomt met het bronreservoir: $\sim 0,4 Z_{ISM}$ voor SN-gedreven winden en $\sim 1,0 Z_{ISM}$ voor kinetische AGN-gedreven uitstroom (die ISM-deeltjes direct wegduwen).
• Uitstroom op halo-schaal volgt de metalenconcentratie van het CGM, wat aangeeft dat er veel menging plaatsvindt.
• Instromend gas naar het ISM is verrijkt ten opzichte van het gemiddelde CGM, consistent met radiale metalen-gradiënten.

Significantie en Conclusies

De studie biedt een dieper inzicht in hoe sterren- en AGN-feedback de groei en quenching van sterrenstelsels reguleren in IllustrisTNG.

• Mechanisme: Kinetische AGN-feedback is cruciaal voor het stoppen van stervorming in massieve halo's, voornamelijk door preventieve feedback (het verhinderen van gasaccretie) en ejective feedback (het uitdrijven van gas).
• Subgrid-beperkingen: De resultaten tonen aan dat de emergente uitstroomeigenschappen (zoals massa-loading) sterk afwijken van de input-parameters van de subgrid-recipes, voornamelijk door interacties met het CGM (stalling bij hoge dichtheid/roodverschuiving en entrainment bij lage dichtheid).
• Toekomstige Modellen: De gedetailleerde data over energiestromen en loading factors als functie van massa en tijd biedt een essentiële richtlijn voor het ontwikkelen van de volgende generatie semi-analytische modellen (SAMs) en gasregulatormodellen, die nu fysiek onderbouwde schaalrelaties kunnen gebruiken in plaats van puur empirische aanpassingen.

Kortom, dit werk legt de kwantitatieve basis voor het begrijpen van de kosmische baryoncyclus in TNG100 en benadrukt de noodzaak om energiestromen en de dynamiek van het CGM expliciet mee te nemen in theoretische modellen van galactische evolutie.