POLAR-II: modeling star formation history of galaxies on the 21-cm signal from Epoch of Reionization

Dit artikel toont aan dat variaties in de sterformatiegeschiedenis van sterrenstelsels, in plaats van een constante vormingsrate, leiden tot grotere en warmere geïoniseerde gebieden en een andere topologie en timing van de 21-cm-signalen tijdens het tijdperk van reionisatie.

De kern

De Sterrenkijkers van de Oertijd: Hoe de "Levensloop" van Sterren de Kosmische Verlichting Beïnvloedt

Stel je het heelal voor als een gigantisch, donker zwembad. Duizenden jaren geleden was dit zwembad vol met koud, onzichtbaar water (het neutrale gas). Toen begonnen de eerste sterren en sterrenstelsels te ontstaan. Ze fungeerden als enorme lichten die het water begonnen te verwarmen en te "ontleden" (ioniseren), waardoor het zwembad helder en warm werd. Dit proces heet de Epoch of Reionization (de Tijdperk van Herionisatie).

Astronomen proberen dit proces te begrijpen door naar het heelal te kijken met speciale radio-antennes. Ze zoeken naar een heel zwak signaal, de "21-cm-lijn", dat ons vertelt hoe koud of warm het water is en hoe ver het licht al is doorgedrongen.

Het Probleem: Te Simpele Verhalen
In dit artikel, getiteld POLAR-II, kijken onderzoekers naar een foutje in de verhalen die we tot nu toe hebben verteld over hoe sterrenstelsels zich gedragen.

Stel je voor dat je een film maakt over een sterrenstelsel. Tot nu toe hebben wetenschappers vaak een heel simpel verhaal gebruikt: "Dit sterrenstelsel heeft gedurende zijn hele leven precies evenveel sterren per seconde gemaakt." Alsof een fabriek 24 uur per dag, 7 dagen per week, met exact hetzelfde tempo werkt.

Maar in het echte leven is dat niet zo. Sterrenstelsels hebben een levensloop (in het Engels: Star Formation History). Soms hebben ze een enorme uitbarsting van sterrengeboorte (een "starburst"), net als een feestje waar iedereen tegelijkertijd begint te dansen. Soms gaan ze "op slot" en stoppen ze met het maken van sterren (een "quenched" periode), alsof de fabriek tijdelijk sluit voor onderhoud.

De Nieuwe Aanpak: Een Realistischer Film
De onderzoekers van dit papier hebben een nieuwe, realistischere manier bedacht om deze sterrenstelsels te modelleren. Ze gebruiken een supercomputer-simulatie (Jiutian-300) en een slim rekenmodel (L-Galaxies 2020) om te kijken hoe sterrenstelsels écht groeien en verouderen.

Ze kijken specifiek naar één eigenschap: de gemiddelde leeftijd van de sterren in een sterrenstelsel.

• Jonge sterrenstelsels: Hebben veel jonge sterren. Ze stralen veel UV-licht uit en verwarmen het gas snel.
• Oude sterrenstelsels: Hebben veel oude sterren. Ze stralen minder UV-licht uit, maar hun "restwarmte" (röntgenstraling) verwarmt het gas op een andere manier.

Wat Vonden Ze? De "Grote Verlichting"
Toen ze deze realistische levenslopen in hun simulaties stopten, zagen ze interessante verschillen vergeleken met de oude, simpele modellen:

1. Grotere "Bellen" van Licht: Sterrenstelsels met een realistische levensloop creëren iets grotere gebieden waar het gas al is "ontleden" (ioniseerd) dan die met een constant tempo. Het is alsof een sterrenstelsel met een echte geschiedenis een groter lichtgebied kan verlichten dan een fabriek die altijd even hard werkt.
2. Warmteverschillen: De manier waarop het gas wordt verwarmd, hangt af van hoe oud de sterren zijn. Als een sterrenstelsel zijn sterren vroeg heeft gemaakt (een hoge "gemiddelde leeftijd"), verwarmt het het omringende gas minder goed dan een stelsel dat recentelijk veel nieuwe sterren heeft gemaakt.
3. Het Signaal Verandert: Omdat het gas op een andere manier wordt verwarmd en verlicht, verandert het 21-cm-signaal dat we met onze radio-antennes zouden moeten zien. Het signaal van een "realistisch" sterrenstelsel ziet er anders uit dan dat van een "simpele" sterrenstelsel.

De Grootste Les
De belangrijkste conclusie is dat we niet kunnen doen alsof alle sterrenstelsels hetzelfde werken. Net zoals mensen verschillende levenslopen hebben (sommige beginnen vroeg met werken, anderen later; sommigen werken hard, anderen nemen pauze), hebben ook sterrenstelsels hun eigen unieke geschiedenis.

Als we deze geschiedenis niet goed meenemen in onze berekeningen, krijgen we een verkeerd beeld van hoe het heelal is opgewarmd en verlicht. Het is alsof je probeert het weer te voorspellen door alleen te kijken naar de gemiddelde temperatuur, zonder rekening te houden met stormen, hittegolven of koude dagen.

Kortom:
Dit artikel laat zien dat als we de complexe "levensgeschiedenis" van sterrenstelsels serieus nemen, we een scherpere foto krijgen van hoe het jonge heelal van donker naar licht is gegaan. Dit helpt ons om de data van toekomstige telescopen (zoals het SKA) beter te begrijpen en te weten wat we precies moeten zoeken in de diepe duisternis van de kosmos.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "POLAR-II: modeling star formation history of galaxies on the 21-cm signal from Epoch of Reionization" in het Nederlands.

Probleemstelling

De Epoch of Reionization (EoR) is een cruciale fase in de geschiedenis van het heelal waarbij het neutrale intergalactische medium (IGM) wordt geïoniseerd en opgewarmd door de eerste sterren en sterrenstelsels. Hoewel semi-numerieke simulaties vaak worden gebruikt om de 21-cm-signaal van dit proces te modelleren, hebben ze een belangrijke beperking: ze nemen doorgaans een constante sterformatiesnelheid (SFR) aan voor sterrenstelsels.

In werkelijkheid ondergaan sterrenstelsels complexe evoluties met periodes van sterburst en uitdoving (quenching) als gevolg van fusies en feedback. Deze variaties in de Sterformatiegeschiedenis (SFH) beïnvloeden hoe ioniserende fotonen en X-stralen over tijd worden vrijgegeven. Hoewel de totale ioniserende fotonenbudget over de levensduur van een sterrenstelsel met dezelfde sterrenmassa vergelijkbaar kan zijn, kan het tijdsverloop van deze emissie de recombinatie en afkoeling van het IGM beïnvloeden. Bestaande modellen negeren deze dynamiek vaak, wat kan leiden tot onnauwkeurigheden in de voorspelde 21-cm-signatuur.

Methodologie

De auteurs hebben de semi-numerieke code POLAR geüpdatet om de SFH van sterrenstelsels expliciet te incorporeren. De methodologie omvat de volgende stappen:

1. Simulatie van Donkere Materie: Gebruik van de Jiutian-300 N-body simulatie (boxgrootte 300 $h^{-1}$ cMpc, $6144^3$ deeltjes) om de dichtheidsvelden van donkere materie te genereren.
2. Zelfvorming van Sterrenstelsels: Toepassing van het semi-analytische model L-Galaxies 2020 (LG20) op de halo's uit de Jiutian-300 simulatie. Dit model berekent de evolutie van sterrenstelsels, inclusief gaskoeling, sterformatie, fusies en feedback, en levert gedetailleerde SFH's op.
3. Bronnen van Straling:
   • Sterren: Berekening van de geïntegreerde Spectrale Energieverdeling (iSED) op basis van de leeftijd en metaliciteit van de sterrenpopulatie (met BPASS).
   • X-ray Binaries (XRB): Modellering van HMXB en LMXB emissie gebaseerd op SFR, sterrenmassa en metaliciteit.
   • Shock-verwarmd ISM: Berekening van zachte X-stralen veroorzaakt door supernova-explosies.
4. Radiative Transfer (RT): Gebruik van de 1-D RT-code Grizzly om de ionisatie en verwarming van het IGM te modelleren. In plaats van een constante SFR, wordt nu de specifieke SFH uit LG20 gebruikt. De auteurs hebben zes simulaties uitgevoerd:
   • Drie met een constante sterrenleeftijd ($\tau_{age}$ = 0.02, 0.1, 0.2 Gyr).
   • Eén met de werkelijke SFH uit LG20 (fiduciale model).
   • Twee vergelijkingsmodellen waarbij de SFR constant wordt gehouden ($SFR = M_*/\tau_{age}$) voor de LG20-leeftijden.

Belangrijkste Bijdragen

• Update van POLAR: De eerste implementatie van een realistische, evoluerende SFH (gebaseerd op semi-analytische modellen) in een semi-numerieke reionisatiecode, in plaats van een vereenvoudigde constante SFR.
• Koppeling van SFH aan 21-cm: Het kwantificeren van hoe variaties in de sterrenleeftijd ($\tau_{age}$) en de vorm van de SFH de ionisatiebellen en de thermische geschiedenis van het IGM beïnvloeden.
• Validatie van Aannames: Het tonen aan dat hoewel de totale ioniserende budget gelijk kan zijn, de tijdsverdeling van de emissie significant van invloed is op de topologie van de reionisatie.

Resultaten

De simulaties leverden de volgende kernresultaten op:

1. Invloed van $\tau_{age}$ op Ionisatiebellen:

   • Voor een vaste sterrenmassa ($M_*$) leiden sterrenstelsels met een grotere gemiddelde sterrenleeftijd ($\tau_{age}$) tot kleinere ionisatiebellen. Dit komt doordat bij een langere leeftijd de recombinatie van geïoniseerd waterstof en vrije elektronen over een langere periode plaatsvindt, wat de effectieve ionisatie vermindert.
   • Sterrenstelsels met een kleinere $\tau_{age}$ produceren grotere ionisatiebellen en warmen het IGM sterker op.

2. Vergelijking: LG20 SFH vs. Constante SFR:

   • Sterrenstelsels met een SFH uit LG20 produceren grotere en warmer ionisatiebellen dan die met een equivalente constante SFR.
   • Dit verschil is het grootst bij jonge sterrenstelsels (kleine $\tau_{age}$) en neemt af naarmate $\tau_{age}$ toeneemt.

3. Invloed op het 21-cm Signaal:

   • Temperatuur en Ionisatie: De verschillen in SFH leiden tot verschillende temperatuur-ionisatie relaties ($T_k$ vs. $x_{HII}$). Simulaties met LG20-SFH tonen hogere gas temperaturen voor dezelfde ionisatiegraad.
   • Timing van Reionisatie: Simulaties met een SFH uit LG20 (of een kleine constante $\tau_{age}$) sluiten de EoR eerder af (ongeveer $\Delta z \approx 0.5$ eerder) dan simulaties met een grote $\tau_{age}$ of constante SFR.
   • Krachtenspectrum ($\Delta^2_{21cm}$):
     • Bij hoge roodverschuivingen ($z > 9$), waar X-stralenverwarming domineert, leidt een zwakkere verwarming (grote $\tau_{age}$) tot een hoger machtsspectrum.
     • Bij lagere roodverschuivingen ($z < 8$), waar ionisatiebellen domineren, verschuiven de pieken in het machtsspectrum naar lagere $z$ voor modellen met grotere $\tau_{age}$.
     • Modellen met LG20-SFH vertonen over het algemeen een lager machtsspectrum dan modellen met constante SFR, vanwege de efficiëntere verwarming en grotere ionisatiebellen.

Significantie

De studie benadrukt dat de aannames over de sterformatiegeschiedenis van sterrenstelsels kritiek zijn voor het interpreteren van toekomstige 21-cm waarnemingen (bijv. van SKA, HERA, LOFAR).

• Het negeren van de evoluerende SFH kan leiden tot fouten in de geschatte timing en topologie van de EoR.
• De resultaten tonen aan dat de 21-cm-signatuur niet alleen afhangt van de totale hoeveelheid ioniserende straling, maar ook van hoe en wanneer deze straling wordt vrijgegeven.
• Voor een nauwkeurige interpretatie van de data van de James Webb Space Telescope (JWST) en toekomstige radiotelescopen is het noodzakelijk om realistische SFH-modellen te gebruiken in simulaties van de EoR.