Euclid preparation. 3D reconstruction of the cosmic web with simulated Euclid Deep spectroscopic samples

Probleemstelling

De Euclid-missie heeft tot doel de aard van donkere energie te karakteriseren en de verdeling van donkere materie te begrijpen. Een cruciaal "bijproduct" van deze missie is de verzameling van spectroscopische roodverschuivingen voor ongeveer twee miljoen sterrenstelsels in de Euclid Deep Fields (EDF). Deze data zullen een diepte bereiken van $5 \times 10^{-17} \text{ erg cm}^{-2} \text{ s}^{-1}$ voor de H$\alpha$-emissielijn in het roodverschuivingsbereik $0.4 < z < 1.8$.

Het centrale wetenschappelijke probleem dat in dit paper wordt aangepakt, is de vraag in welke mate het kosmisch web (de grote schaalstructuur bestaande uit filamenten, wanden, knooppunten en holtes) nauwkeurig kan worden gereconstrueerd op basis van deze specifieke spectroscopische dataset. De uitdagingen zijn:

1. Selectie-effecten: De data zijn beperkt door een flux-limiet voor H$\alpha$ (sterrenstelsels met hoge sterrenvorming), in plaats van een massa-limiet. Dit leidt tot een onvolledige steekproef van de totale populatie sterrenstelsels.
2. Roodverschuivingsfouten: Meetonzekerheden in de roodverschuiving ($\sigma_z \approx 0.001(1+z)$) kunnen de 3D-positie van sterrenstelsels verstoren.
3. Redshift-space distortions (RSD): Vooral het "Fingers-of-God" (FoG) effect, veroorzaakt door de hoge snelheden van sterrenstelsels binnen virialiseerde groepen, vervormt de waargenomen structuur langs de lijn van zicht, waardoor filamenten kunstmatig verlengd lijken.

Methodologie

De auteurs evalueren de kwaliteit van de reconstructie door gebruik te maken van twee verschillende sets van gesimuleerde sterrenstelselcatalogi (mocks):

1. Euclid Flagship: Een N-body simulatie met $4 \times 10^{12}$ donkere materie-deeltjes, waarbij sterrenstelsels worden toegewezen aan halos via Halo Occupation Distribution (HOD) en Abundance Matching.
2. GAEA (GAlaxy Evolution and Assembly): Een semi-analytisch model dat draait op de Millennium-simulatie, met een focus op fysiek gemotiveerde voorschriften voor gas, metalen en energie-recycling.

De simulatie van de EDF-data:
Om de toekomstige Euclid-data na te bootsen, passen de auteurs de volgende beperkingen toe op de mocks:

• Flux-limiet: Selectie op H$\alpha$-flux ($6 \times 10^{-17} \text{ erg cm}^{-2} \text{ s}^{-1}$).
• Completie: Willekeurige verwijdering van 40% van de sterrenstelsels om een completie van 60% te simuleren.
• Roodverschuivingsruis: Toevoeging van Gaussische fouten aan de roodverschuivingen.
• Correctie voor FoG: Gebruik van een anisotrope "Friends-of-Friends" (FoF) algoritme om groepen te identificeren en vervolgens de groepen radiaal te comprimeren om de FoG-effecten te minimaliseren.

Reconstructie-algoritme:
Voor het extraheren van het kosmisch web wordt DisPerSE (Discrete Persistent Structures Extractor) gebruikt, gebaseerd op discrete Morse-theorie en Delaunay-tessellatie. De auteurs vergelijken twee benaderingen:

1. Niet-gewogen: Gebaseerd op de dichtheid van het aantal sterrenstelsels.
2. Gewogen: De Delaunay-tessellatie wordt gewogen met de sterrenmassa van de sterrenstelsels om de onderliggende dichtheidsveld beter te benaderen.

Evaluatiemetrics:
De kwaliteit wordt gemeten aan de hand van:

• Geometrische statistieken: Lengte van filamenten, oppervlakte van wanden, volume van holtes.
• Topologische statistieken: Connectiviteit (aantal filamenten verbonden met een knooppunt) en multipliciteit (lokaal aantal verbindingen).
• Transversale gradiënten: Het herstel van de trend dat zwaardere sterrenstelsels dichter bij filamenten voorkomen dan lichtere.

Belangrijkste Bijdragen

1. Validatie van reconstructiemethoden voor Euclid: Het paper biedt de eerste uitgebreide technische richtlijnen voor het reconstrueren van het 3D-kosmisch web met de specifieke beperkingen van de Euclid EDF-data.
2. Rol van Sterrenmassa-Weigting: Het demonstreert dat het wegen van de tessellatie met sterrenmassa essentieel is om de bias van een H$\alpha$-flux-gelimiteerde steekproef te corrigeren, hoewel dit nieuwe systematische fouten introduceert die rekening moeten houden bij vergelijking met theorie.
3. Robuustheid van Multipliciteit: Het identificeert "multipliciteit" als een robuustere topologische maatstaf dan "connectiviteit" in het licht van ruis en onvolledige steekproeven.
4. Beperkingen van Roodverschuivingsfouten: Het kwantificeert dat roodverschuivingsonzekerheden de grootste uitdaging vormen voor het herstel van gradiënten in sterrenstelsel-eigenschappen, zelfs na correctie voor FoG.

Resultaten

• Geometrie: Het FoG-effect heeft een sterke invloed op de lengte van filamenten (deze worden te lang gereconstrueerd). Na compressie van groepen komt de verdeling van de lengtes goed overeen met de referentie, vooral bij gebruik van sterrenmassa-weigting. De onvolledigheid van de steekproef (60%) verschuift de verdelingen naar hogere waarden.
• Connectiviteit: Roodverschuivingsruis en onvolledigheid verminderen de kwaliteit van de reconstructie aanzienlijk. Sterrenmassa-weigting helpt de correlatie tussen sterrenmassa en connectiviteit te herstellen, maar introduceert een bias waarbij de connectiviteit kunstmatig wordt verhoogd bij lagere massa's.
• Multipliciteit: In tegenstelling tot connectiviteit, is multipliciteit minder gevoelig voor het FoG-effect en roodverschuivingsfouten. Zelfs zonder weging kan de relatie tussen sterrenmassa en multipliciteit redelijk goed worden hersteld na FoG-correctie.
• Sterrenmassa-gradiënten: De verwachte trend (zwaardere sterrenstelsels dichter bij filamenten) is aanwezig in de referentie. Echter, bij de H$\alpha$-flux-gelimiteerde steekproef met roodverschuivingsfouten en onvolledigheid, wordt dit signaal sterk verzwakt. Sterrenmassa-weigting verbetert het signaal, maar introduceert een bias omdat zwaardere sterrenstelsels per definitie dichter bij het web worden geplaatst.

Betekenis en Conclusie

Dit onderzoek is cruciaal voor de voorbereiding van de Euclid-missie. Het concludeert dat:

• Een 3D-reconstructie van het kosmisch web met Euclid-data mogelijk is, maar vereist zorgvuldige behandeling van selectie-effecten en roodverschuivingsfouten.
• Sterrenmassa-weigting een noodzakelijke maar dubbelzinnige tool is: het verbetert de geometrische reconstructie en het herstel van gradiënten, maar vervormt topologische maten zoals connectiviteit.
• Multipliciteit de voorkeur verdient boven connectiviteit als topologische maatstaf voor deze specifieke dataset.
• De grootste beperking voor het bestuderen van de co-evolutie van sterrenstelsels en het kosmisch web op hoge roodverschuiving ($z > 1$) de roodverschuivingsonzekerheid blijft.

De resultaten bieden een blauwdruk voor toekomstige analyses van de echte Euclid-data, waarbij wetenschappers kunnen focussen op het bestuderen van de invloed van de grote schaalstructuur op eigenschappen zoals sterrenvorming, morfologie en spin-oriëntatie van sterrenstelsels tijdens de "kosmische middag" ($z \approx 1-2$).