The THESAN project: Lyman-alpha emitters as probes of ionized bubble sizes

Probleemstelling

Tijdens het Tijdperk van de Reionisatie (EoR, ongeveer $z \approx 5-20$) ioniseerden de eerste generaties sterrenstelsels het neutrale waterstof in het intergalactische medium (IGM), wat leidde tot de vorming van ionisatiebellen. Het begrijpen van de grootte en topologie van deze bellen is cruciaal voor het reconstrueren van de geschiedenis van de reionisatie.
De belangrijkste observatieproef hiervoor is de Lyman-alfa (Ly$\alpha$) emissielijn. Echter, de detecteerbaarheid van Ly$\alpha$-emitters (LAE's) wordt niet alleen bepaald door de grootte van de ionisatiebel, maar ook door complexe factoren zoals:

• Intrinsieke Ly$\alpha$-productie in sterrenstelsels.
• Absorptie door stof (dust) en resonante verstrooiing in het interstellair medium (ISM) en circumgalactisch medium (CGM).
• De ionisatietoestand van het omringende IGM (belgrootte en neutrale fractie).
• Peculiare snelheden die de lijnprofielen beïnvloeden.

Bestaande analytische modellen gebruiken vaak vereenvoudigde, bolvormige bel-geometrieën, wat kan leiden tot degeneraties en onnauwkeurigheden. Er is behoefte aan een robuust, op simulaties gebaseerd kader om de relatie tussen waargenomen Ly$\alpha$-eigenschappen en de werkelijke grootte van ionisatiebellen te kalibreren.

Methodologie

De auteurs gebruiken de thesan-suite van kosmologische stralingshydrodynamische (RHD) simulaties, die de vorming van sterrenstelsels koppelt aan on-the-fly stralingstransport (met de code arepo-rt).

1. Simulatiebasis: De flagship-run (thesan-1) heeft een kubusgrootte van 95,5 cMpc en een hoge resolutie ($m_{gas} \approx 5.82 \times 10^5 M_\odot$), wat het mogelijk maakt om atoomkoelingshalo's en de evolutie van ionisatiebellen rond sterrenstelsels op te lossen.
2. Intrinsieke Emissie: De simulaties berekenen de intrinsieke Ly$\alpha$-luminositeit gebaseerd op recombinatie, collisionele excitatie en sterrenemissie, gekoppeld aan de ionisatievelden van het IGM.
3. Transmissie door IGM: Voor ongeveer 90% van de zwaarste centrale sterrenstelsels wordt de transmissie van Ly$\alpha$-fotonen berekend langs 768 stralen (ray-tracing) door het IGM, rekening houdend met Doppler-verschuivingen en de neutrale waterstofkolom.
4. Empirische Kalibratie (Nieuwe Bijdrage): Omdat de simulaties sub-galactische processen (zoals stofabsorptie en winds) niet volledig oplossen, wordt een empirisch model toegepast om de waargenomen Ly$\alpha$-eigenschappen te voorspellen:
   • Lijnprofiel: Een analytische oplossing voor een asymmetrisch, dubbel-gepiekt profiel binnen een uitdijende wolk, gekoppeld aan de cirkelsnelheid van het sterrenstelsel.
   • Stofcorrectie: Een geüpdate model voor stofextinctie ($A_{1500}$) dat is gekalibreerd op ALMA-observaties, waarbij de Ly$\alpha$-ontsnappingsfractie beperkt wordt tot de continuüm-ontsnappingsfractie.
   • Optimalisatie: De parameters van het model (voor piekpositie, windsterkte en stof) worden geoptimaliseerd door de gesimuleerde Ly$\alpha$-luminositeitsfuncties (LF) af te stemmen op waarnemingen bij $z=5.5$ en $z=6.6$ (Ouchi et al., Konno et al.).

Belangrijkste Resultaten

1. Kalibratie en Luminositeitsfuncties
Het geüpdate model reproduceert succesvol de waargenomen Ly$\alpha$-luminositeitsfuncties. Het model toont aan dat stof en uitstromingen (winds) de waargenomen luminositeit van heldere sterrenstelsels significant onderdrukken (tot twee ordes van grootte ten opzichte van de intrinsieke waarde). Er is een lichte spanning bij $z=6.6$ voor zeer heldere bronnen, waarschijnlijk door de beperkte kubusgrootte van de simulatie die zeldzame, extreme bronnen niet volledig kan genereren.

2. Relatie tussen Eigenschappen en Omgeving

• UV-Brightness vs. Ly$\alpha$: Helderere sterrenstelsels (in UV) hebben doorgaans lagere Ly$\alpha$-ontsnappingsfracties en equivalente breedtes (EW) door sterkere stofabsorptie en winds.
• Evolutie: Naarmate de reionisatie vordert (lagere $z$), neemt de IGM-transmissie toe. Fijne sterrenstelsels tonen een sterkere toename in transmissie dan heldere sterrenstelsels, omdat ze minder beïnvloed worden door interne ISM-processen.

3. LAE-Fractie en Reionisatiegeschiedenis
De voorspelde fractie van sterrenstelsels die als LAE's worden gedetecteerd (met $L_{Ly\alpha} > 10^{41.5}$ erg s$^{-1}$ en $EW > 25$ Å) komt goed overeen met JWST-observaties bij hoge roodverschuivingen ($z \gtrsim 8$). Bij lagere roodverschuivingen ($z < 7$) overschat het model de LAE-fractie licht, wat suggereert dat de reionisatie in thesan iets te vroeg voltooid is vergeleken met de nieuwste waarnemingen.

4. Correlatie met Belgrootte (Kernbevinding)
De studie onthult een sterke statistische relatie tussen de eigenschappen van LAE's en de grootte van de ionisatiebel waarin ze zich bevinden:

• Vóór het middenpunt van reionisatie ($z \gtrsim 7$): Sterrenstelsels in grotere bellen ($\gtrsim 10$ cMpc) hebben een hogere waargenomen Ly$\alpha$-luminositeit en een hogere equivalente breedte. Grote bellen staan toe dat Ly$\alpha$-fotonen uit resonantie verschuiven voordat ze neutraal gas tegenkomen.
• Selectiecriteria: De correlatie is het sterkst voor selectie op Ly$\alpha$-luminositeit en equivalente breedte (EW), en zwakker voor selectie op UV-magnitude. Vooral bij $z \gtrsim 7$ is een hoge EW een betrouwbare indicator voor een grote ionisatiebel.
• Verdwijning van correlatie: Naarmate de bellen samensmelten (percolatie) en het IGM volledig ioniseert, zwakt deze correlatie af. Uiteindelijk kunnen ook sterrenstelsels in kleinere bellen (of waar de bellen al samengesmolten zijn) als LAE's worden waargenomen.

Significantie en Conclusie

Dit artikel biedt een cruciaal theoretisch kader om de complexe waarnemingen van LAE's tijdens de reionisatie te interpreteren:

1. Statistische Tracer: De evoluerende fractie van LAE's en hun eigenschappen (vooral EW en $L_{Ly\alpha}$) fungeren als praktische statistische tracers voor de grootte van ionisatiebellen, vooral in de vroege fasen van de reionisatie.
2. Correctie voor Bias: Het werk identificeert en corrigeert voor biases in idealistische bolmodellen door gebruik te maken van realistische 3D-simulaties die inhomogene dichtheidsstructuren en complexe stralingstransport meenemen.
3. Toekomstige Toepassingen: De gegenereerde catalogi met gekalibreerde eigenschappen zullen dienen als referentie voor het interpreteren van huidige en toekomstige waarnemingen (o.a. van JWST en narrow-band surveys). Dit helpt bij het vertalen van observaties naar de topologie van de reionisatie en het beperken van de reionisatiegeschiedenis van het heelal.

Kortom, de studie bevestigt dat Ly$\alpha$-emitters, mits correct gekalibreerd voor stof en winds, krachtige instrumenten zijn om de grootte en evolutie van ionisatiebellen in het vroege heelal in kaart te brengen.