Angular bispectrum of matter number counts in cosmic structures

Technische Samenvatting: Angulaire Bispectrum van Materie-aantalstellingen in Kosmische Structuren

Probleemstelling
De twee-punts correlatiefunctie (krachtspectrum) blijft het primaire instrument voor het beperken van het standaard $\Lambda$CDM-model. Echter, naarmate de grootschalige structuur (LSS) evolueert naar het zwak niet-lineaire regime, worden hogere-orde statistieken, specifiek de drie-punts correlatiefunctie en de bijbehorende Fourier-tegenhanger, het bispectrum, essentieel voor het onderzoeken van Fourier-fasecorrelaties en moduskoppelingen veroorzaakt door niet-lineaire gravitationele evolutie en primordiale niet-Gaussianiteit (PNG).

Huidige analyses van de bispectra van sterrenstelsels vertrouwen zwaar op de lokale vlakke-hemel benadering (flat-sky approximation), de Limber-benadering en Newtoniaanse dynamica. Deze benaderingen breken af voor eerstvolgende generatie surveys (bijv. Euclid, DESI, SPHEREx) die grote gebieden bestrijken en gebruikmaken van eindige roodverschuivingsbakken (tomografie). Bovendien, hoewel relativistische effecten (projectie, dynamische Algemene Relativiteitstheorie en straling) theoretisch bekend zijn, is hun exacte impact op het angulaire bispectrum in een setup met eindige roodverschuivingsbakken zonder de Limber-benadering nog niet volledig berekend. Daarnaast is het extraheren van deze subtiele relativistische signalen uit simulaties uitdagend vanwege numerieke artefacten en de dominantie van Newtoniaanse termen.

Methodologie
De auteurs presenteren een full-sky berekening van het angulaire bispectrum voor materie-aantalstellingen met behulp van tweede-orde perturbatietheorie. De methodologie wordt gekenmerkt door de volgende technische keuzes:

1. Raamwerk: De berekening maakt gebruik van "tomografische sferische harmonieken" (TSH), waarbij vensterfuncties worden gebruikt om eindige roodverschuivingsbakken te definiëren ($1.75 \le z \le 2.25$ en $0.55 \le z \le 0.65$). Deze aanpak vermijdt de Limber-benadering en de vlakke-hemel aanname.
2. Perturbatie-expansie: De aantalstelling-perturbatie $\Delta$ wordt uitgebreid tot de tweede orde. De auteurs organiseren de termen volgens de zwakke-veld expansieparameter $H/k$:
   • Newtoniaanse Termen ($O((H/k)^0)$): Bevat lineaire dichtheid, redshift-space distorties (RSD) en kwadratische termen (dichtheids-dichtheid en dichtheid-snelheid koppelingen).
   • Relativistische Projectie-effecten ($O(H/k)$): Leidende projectie-effecten (Doppler, Sachs-Wolfe, etc.) tot de eerste orde in de zwakke-veld expansie.
   • Dynamische GR en Stralingseffecten ($O((H/k)^2)$ en hoger): Bevat niet-lineaire dynamische GR-correcties en stralingscorrecties (volgens Tram et al. 2016).
   • Koppelings-termen: Bevat koppelingen tussen projectie en dynamische effecten tot $O((H/k)^3)$.
3. Uitsluitingen: Voor de eenvoud en om directe consistentiecontroles mogelijk te maken, verwaarlozen de auteurs tracer-bias, langs-de-gezichtslijn geïntegreerde effecten (lensing, Integrated Sachs-Wolfe) en waarnemer-positie termen. Ze laten ook de volgende-orde projectie-effecten op de tweede orde weg vanwege de complexiteit en het gebrek aan consensus in de literatuur.
4. Numerieke Implementatie: De auteurs breiden het efficiënte wiskundige raamwerk van Assassi et al. (2017) uit, waarbij gebruik wordt gemaakt van FFTLog-decomposities en hypergeometrische functie-evaluaties om het probleem te reduceren tot één-dimensionale integralen. Ze leveren een publieke code, ang_bispec, om de 14 verschillende tree-level bispectrum bijdragen te berekenen.
5. Validatie: Theoretische voorspellingen worden vergeleken met relativistische $N$-body simulaties (Montandon et al. 2023). Om een eerlijke vergelijking te garanderen, verwijderen de auteurs de geïntegreerde projectie-effecten uit de simulatiegegevens en passen ze binning en Gaussische smoothing toe op de theoretische voorspellingen om de simulatie-estimator (Bucher et al. 2016) te matchen.

Belangrijkste Bijdragen

• Eerste Full-Sky Berekening: Dit werk biedt de eerste full-sky berekening van het angulaire bispectrum in tweede-orde perturbatietheorie zonder de Limber-benadering, geformuleerd voor eindige roodverschuivingsbakken.
• Nieuwe Newtoniaanse Berekening: Zelfs het Newtoniaanse deel van het bispectrum binnen deze specifieke setup (eindige bakken, geen Limber) wordt als nieuw gerapporteerd.
• Omvattende Relativistische Inclusie: De auteurs bevatten, tot de tweede orde, alle dynamische GR-effecten, stralingscorrecties en leidende relativistische projectie-effecten binnen een verenigd raamwerk.
• Code-release: De ang_bispec code wordt publiekelijk beschikbaar gesteld, waarbij eerdere berekeningen worden uitgebreid naar alle dominante Newtoniaanse termen en relativistische effecten.

Resultaten
De auteurs evalueerden het bispectrum voor twee roodverschuivingsbakken en vergeleken de bijdragen van de verschillende termen:

1. Dominantie van Newtoniaanse Termen: Zoals verwacht, zijn de Newtoniaanse bijdragen (dichtheid, RSD, kwadratisch) doorgaans één of meerdere ordes van grootte groter dan de relativistische signalen over het gehele spectrum voor beide roodverschuivingen.
2. Relativistische Hiërarchie bij $z \approx 2$:
   • Projectie en dynamische GR-effecten hebben vergelijkbare amplitudes op grote schalen.
   • Onverwachte Bevinding: In de squeezed limit domineren stralingseffecten het relativistische signaal, waarbij ze de andere relativistische termen met een factor tien overtreffen.
3. Relativistische Hiërsarchie bij $z \approx 0.6$: De verwachte hiërarchie wordt hersteld waarbij projectie-effecten domineren. Echter, dynamische correcties (door niet-lineaire evolutie) blijven niet verwaarloosbaar en zijn slechts een factor 2–3 kleiner dan de projectie-effecten.
4. Simulatie Vergelijking:
   • Totaal Bispectrum: De theoretische resultaten komen redelijk goed overeen met de simulatie-metingen voor het totale bispectrum, behalve op kleine schalen waar de resolutie van de simulatie de kracht beperkt.
   • Relativistisch Signaal: Bij het vergelijken van het "zuivere" relativistische signaal (som van GR en stralingstermen) met de simulatie-metingen (afgeleid door het aftrekken van Newtonian van Relativistic simulaties), vinden de auteurs een discrepantie. De theoretische som is een factor $\sim 5$ kleiner dan de simulatie-meting bij een hoge roodverschuiving.
   • Interpretatie van de Discrepantie: De auteurs schrijven deze kloof toe aan numerieke artefacten in de simulaties, specifiek residuele gauge-effecten (mismatch tussen Poisson en N-body gauges) en backscaling-effecten (hoe straling wordt opgenomen in Newtoniaanse runs). Ze concluderen dat deze numerieke fouten vergelijkbaar zijn in amplitude met de fysieke dynamische relativistische effecten.
   • Smoothing Effecten: De analyse benadrukt dat smoothing (nodig voor het detecteren van minuscule relativistische signalen) systematische biases kan introduceren door kleine-schaal numerieke fouten naar grote schalen te laten lekken, wat mogelijk de geobserveerde offsets in lage-roodverschuiving vergelijkingen verklaart.

Betekenis en Claims
Het artikel claimt een robuust theoretisch proces te hebben vastgesteld voor de analyse van het angulaire bispectrum van materie-aantalstellingen in fotometrische surveys, waarmee de kloof tussen lensing en sterrenstelsel-clustering observabelen wordt overbrugd.

• Validatie van Amplitudes: Het werk bevestigt de amplitude-schattingen van de relativistische effecten afgeleid van de simulaties van Montandon et al. (2023), terwijl het een term-voor-term decompositie biedt die simulaties niet gemakkelijk kunnen isoleren.
• Relevantie van Dynamische Effecten: In tegengens van de aanname dat projectie-effecten de enige significante relativistische bijdrage zijn, tonen de auteurs aan dat dynamische GR en stralingseffecten van dezelfde orde grootte zijn als projectie-effecten bij een hoge roodverschuiving en nog steeds significant blijven (binnen een factor 2) bij lagere roodverschuivingen.
• Waarschuwing bij Simulatie Extractie: Een belangrijke bevinding is de identificatie van numerieke artefacten (gauge en backscaling) in huidige relativistische simulaties die even groot zijn als de fysieke signalen die men zoekt. Dit suggereert dat het extraheren van "zuivere" relativistische signalen uit huidige simulatie-aftrekmethoden extreme voorzichtigheid vereist.
• Toekomstige Bruikbaarheid: De auteurs stellen dat hun raamwerk essentieel is voor toekomstige eerstvolgende generatie surveys (zoals Euclid) om primordiale niet-Gaussianiteit te onderscheiden van intrinsieke niet-Gaussianiteit en om theorieën van zwaartekracht in het niet-lineaire regime te testen, mits simplificaties zoals sterrenstelsel-bias en langs-de-gezichtslijn integralen uiteindelijk worden opgenomen.