ECHO21: a tool for modelling global 21-cm signal from dark ages to reionization

Technische Samenvatting van "ECHO21: Een instrument voor het modelleren van het globale 21-cm signaal van de dark ages tot reïonisatie"

Probleem en Motivatie
De dark ages en de cosmic dawn zijn slecht begrepen perioden in de kosmische tijdlijn vanwege een gebrek aan directe observationele sondes. De 21-cm hyperfine overgangslijn van neutraal waterstof biedt een potentiële sonde om de fysica van de eerste sterren, het intergalactische medium (IGM) en kosmologische parameters te ontrafelen. Hoewel numerieke en semi-numerieke codes zoals 21cmFAST en BEORN nauwkeurige, fysisch rijke voorspellingen bieden, zijn ze computationeel traag, wat efficiënte astrofysische of kosmologische inferentie uit grote datasets belemmert. Daarentegen bieden bestaande analytische codes zoals ARES en Zeus21 weliswaar snelheid, maar missen zij bepaalde fysische componenten, specifiek de opwarmende effecten van Lyα-fotonen. De auteurs identificeren dat het weglaten van Lyα-verwarming kan leiden tot sterke biases in signaalmodellering. Om te voldoen aan de behoefte aan een snelle, voldoende realistische en flexibele tool die in staat is om grote parameterruimten te behandelen voor inferentie, introduceren de auteurs ECHO21 (Exploring the Cosmos with Hydrogen Observation).

Methodologie
ECHO21 is een Python-package die gebruikmaakt van een analytisch kader om het globale 21-cm signaal te modelleren van roodverschuiving $z \approx 1500$ tot de huidige dag. De code werkt met een twee-zonesmodel van het IGM, waarbij onderscheid wordt gemaakt tussen het bulk-IGM en volledig geïoniseerde H II-regio's.

• Fysisch Kader: De code lost gekoppelde differentiaalvergelijkingen op voor de kinetische gastemperatuur ($T_k$), de elektronfractie van het bulk-IGM ($x_e$) en de volumefullingsfactor van H II-regio's ($Q$).

  • Thermische Geschiedenis: De evolutie van $T_k$ houdt rekening met adiabatische koeling, Compton-verwarming, X-ray verhitting en, cruciaal, Lyα-verhitting. De auteurs merken op dat waar voorgaande analytische codes Lyα-verhitting negeerden, ECHO21 dit wel bevat door de verwarmingssnelheid te berekenen op basis van de Lyα-specifieke intensiteit ($J_{Ly}$) afgeleid van de stervormingssnelheidsdichtheid (SFRD).
  • Ionisatiegeschiedenis: De code volgt de recombinatie van het bulk-IGM en de ionisatie gedreven door X-straling en UV-fotonen. Het modelleert de volumefullingsfactor van H II-regio's door de emissiesnelheid van ioniserende fotonen te vergelijken met radiatieve recombinaties.
  • Stervormingssnelheidsdichtheid (SFRD): Een nieuw kenmerk van ECHO21 is de inclusie van drie verschillende SFRD-modellen:
    1. Fysisch-gemotiveerd: Gebaseerd op de instorting van baryonen in donkere materie halo's (met gebruik van Press-Schechter, Sheth-Tormen of Tinker halo-massenfuncties) met een constante ster vormings-efficiëntie.
    2. Semi-empirisch: Een variant van het fysische model die een karakteristieke ster vormings-tijdschaal en een scherpe afkap bij een minimale viriale temperatuur incorporeert.
    3. Empirisch-gemotiveerd: Een model dat beperkt wordt door UV-luminositeitsfunctie-data (Madau-Dickinson vorm bij lage $z$, exponentiële daling bij hoge $z$).

• Parameters: Het model wordt beheerd door zes vrije astrofysische parameters: Lyα-normalisatie ($f_{Ly}$), Lyman-serie spectrale index ($s$), X-ray normalisatie ($f_X$), X-ray spectrale index ($w$), ontsnappingsfractie van ioniserende fotonen ($f_{esc}$), en de minimale viriale temperatuur ($min(T_{vir})$). Daarnaast staat de code variatie van standaard kosmologische parameters toe ($H_0, \Omega_m, \Omega_b, \sigma_8, n_s, T_{\gamma0}, Y_p$).

• Prestaties: De analytische benadering stelt ECHO21 in staat om een enkele globale 21-cm signaal instantie te genereren in ongeveer $O(1)$ seconde op een enkele CPU. De code is ontworpen voor MPI-parallellisatie, waardoor het genereren van duizenden modellen binnen minuten mogelijk is.

Belangrijkste Bijdragen

1. Inclusie van Lyα-verhitting: ECHO21 is de eerste analytische code die expliciet Lyα-verhitting insluit in de thermische evolutie van het IGM. De auteurs demonstreren dat het negeren van deze term kan leiden tot significante biases, met name in scenario's met lage X-ray emissiviteit en sterke Lyα-emissie.
2. Diversiteit in SFRD-modellen: De code biedt uniek drie verschillende benaderingen voor het modelleren van de SFRD (fysisch, semi-empirisch en empirisch), waardoor gebruikers kunnen beoordelen hoe de keuze van een ster vormings-voorschrift de resulterende 21-cm signaal beïnvloedt.
3. Flexibiliteit in Kosmologie: In tegen tegenover veel gespecialiseerde codes behandelt ECHO21 kosmologische parameters met dezelfde gemak als astrofysische parameters, wat de studie van kosmologische onzekerheden (bijv. de Hubble-spanning) op het 21-cm signaal faciliteert.
4. Snelheid en Schaalbaarheid: De code is geoptimaliseerd voor snelle inferentie, in staat om 15.625 modellen te genereren over een zes-dimensionale astrofysische parameterruimte in ~17 minuten met 50 CPU's.

Resultaten

• Fiducieel Model: Met gebruik van een fysisch-gemotiveerde SFRD en een Press-Schechter halo-massenfunctie, voorspelt het fiduciaele model een reïonisatiegeschiedenis die voltooid is tegen $z \approx 5.3$ met een 50% ionisatie roodverschuiving van $z_{re} = 7.42$. De elektron-verstrooiings optische diepte is berekend als $\tau_e \approx 0.0592$, consistent met de Planck 2018 restricties.
• Signaal Karakteristieken: Het fiduciaal globale 21-cm signaal vertoont een eerste absorptiepiek van $-41.1$ mK bij $z=87.5$ (dark ages) en een tweede, diepere absorptiepiek van $-142.4$ mK bij $z=19.2$ (cosmic dawn), gevolgd door een emissiepiek van $22.2$ mK bij $z=11.2$.
• Impact van Lyα-verhitting: Vergelijkingen met modellen die Lyα-verhitting uitsluiten, laten zien dat het weglaten hiervan leidt tot een sterker (dieper) 21-cm absorptiesignaal. De auteurs benadrukken dat Lyα-verhitting een cruciale component is voor nauwkeurige modellering, aangezien het koelingseffecten kan tegenwerken en de signaalsterkte niet-monotoon kan veranderen afhankelijk van de X-ray parameters.
• Impact van SFRD Modellen: De keuze van het SFRD-model verandert het signaal aanzienlijk. Empirisch-gemotiveerde modellen, die een latere start van stervorming voorspellen vergeleken met fysisch-gemotiveerde modellen, resulteren in signalen die met $\sim 8.8$ in roodverschuiving zijn vertraagd en met $\sim 39.4$ mK minder diep zijn. Ook de keuze van de Halo Mass Function (HMF) beïnvloedt de timing en diepte van de signaalfeatures.
• Parametergevoeligheid: Het signaal is het meest gevoelig voor $f_{Ly}$, $f_X$, $w$, en $min(T_{vir})$. Variaties in kosmologische parameters binnen de huidige observationele onzekerheden resulteren in relatief kleine variaties in de signaalamplitude (bijv. een bereik van $\sim 18$ mK in de piekabsorptiediepte).

Betekenis en Claims
De auteurs positioneren ECHO21 als een flexibele en uitbreidbare open-source tool, ontworpen voor het maken van "snelle maar voldoende realistische astrofysische inferenties" uit toekomstige 21-cm datasets (bijv. van REACH, SARAS, of toekomstige ruimte-gebaseerde experimenten). Zij claimen dat de code bijzonder geschikt is voor:

• Het uitvoeren van robuuste astrofysische en kosmologische inferentie via grootschalige parametereksploratie.
• Het testen van data-analyse pipelines (specifiek voor de REACH-collaboratie) met fysisch-gemotiveerde signalen in plaats van louter empirische fits.
• Het bouwen van synergieën tussen 21-cm observaties en andere kosmologische probes, zoals BBN-voorspellingen, CMB-anisotropieën en de Hubble-spanning.

Het paper erkent bescheiden dat hoewel ECHO21 minder nauwkeurig is dan semi-numerieke modellen vanwege de analytische benaderingen (zoals het verwaarlozen van meervoudige Lyα-verstrooiingen), de snelheid en modulariteit ervan het een ideaal instrument maken voor initiële inferenties en voor het integreren van nieuwe fysica in toekomstige versies. De code is publiekelijk beschikbaar via PyPI en GitHub.