Euclid preparation: The flat-sky approximation for the clustering of Euclid's photometric galaxies

Probleemstelling

De Euclid-missie zal een enorme catalogus van ongeveer $1,5 \times 10^9$ fotometrische sterrenstelsels genereren. Een van de belangrijkste statistieken voor deze survey is de analyse van de hoekvermogensspectra ($C_\ell$) van het aantal sterrenstelsels (galaxy number counts) in verschillende roodverschuivingsbinnen. Om de data optimaal te benutten voor het afleiden van kosmologische parameters (bijvoorbeeld via MCMC-analyses), moeten deze theoretische spectra snel en nauwkeurig worden berekend.

De standaardbenadering in de huidige analysepakketten is de Limber-benadering. Deze benadering versnelt de berekening aanzienlijk (ongeveer 100 keer sneller dan de volledige berekening) door de integraal over de golflengte $k$ te vervangen door een Dirac-delta-functie. Hoewel de Limber-benadering uitstekend werkt voor weak lensing (waarbij de vensterfuncties breed zijn), staat bekend dat deze tekortschiet voor fotometrische surveys met smalle roodverschuivingsbinnen. In dergelijke gevallen kan de Limber-benadering afwijken met een factor 2 of meer van de ware waarde bij lage multipoles ($\ell \lesssim 100$). Het artikel onderzoekt of deze onnauwkeurigheid kritiek is voor de toekomstige Euclid-data releases (DR1 en DR3) en vergelijkt de Limber-benadering met de flat-sky-benadering.

Methodologie

De auteurs vergelijken drie methoden voor het berekenen van de hoekvermogensspectra:

1. Volledige berekening: De exacte oplossing (Eq. 7) berekend met de Boltzmann-solver CLASS, inclusief alle relativistische effecten (dichtheidsfluctuaties, roodverschuivingsruimtestorringen (RSD) en lensing-magnificatie).
2. Limber-benadering: De standaard snelle benadering (Eq. 9).
3. Flat-sky-benadering: Een alternatieve benadering (Eq. 18) die de hemel benadert als een vlak en de $(\ell, m)$-coördinaten vervangt door een 2D-vector $\boldsymbol{\ell}$. Deze methode behoudt de parallelle component van de golfvector $k_\parallel$ en is nauwkeuriger bij grote schalen, terwijl ze nog steeds sneller is dan de volledige berekening.

Simulatie-instellingen:
De studie gebruikt het Flagship-simulatie (een N-body simulatie) om twee scenarios te definiëren:

• DR1 (Data Release 1): Een optimistisch scenario met 6 equi-populaire roodverschuivingsbinnen en een helderheidscutoff van $I_E < 23.5$.
• DR3 (Data Release 3): Het definitieve scenario aan het einde van de missie met 13 equi-populaire binnen en een diepere cutoff ($I_E < 24.5$).

Evaluatiemethode:
Om de prestaties te beoordelen in de context van echte metingen, simuleerden de auteurs 1000 ruisrealisaties van de vermogensspectra gebaseerd op de theoretische covariantiematrix (Gaussisch). Ze berekenden de $\chi^2$-verdeling voor zowel de Limber- als de flat-sky-benadering ten opzichte van de "ware" volledige berekening. De afstand tussen de verdelingen ($\Delta\chi^2$) en de bijbehorende p-waarden geven aan of een benadering statistisch significant afwijkt van de ware oplossing binnen de verwachte meetfouten.

Belangrijkste Resultaten

1. DR1 Scenario (6 binnen):

• Voor het eerste data release met brede binnen is de Limber-benadering voldoende nauwkeurig.
• De fouten die door de Limber-benadering worden gemaakt, vallen binnen de $1\sigma$ meetfouten van de survey.
• De $\chi^2$-analyse toont aan dat de verdeling van de Limber-benadering ($p \approx 0.038$) net op de rand van significantie ligt, maar de flat-sky-benadering ($p \approx 0.39$) is niet significant verschillend van de ware oplossing.
• Conclusie: Voor DR1 is de extra complexiteit van de flat-sky-benadering niet nodig; Limber is acceptabel.

2. DR3 Scenario (13 binnen):

• Voor het definitieve scenario met smalle binnen faalt de Limber-benadering karakteristiek.
• De relatieve fouten van de Limber-benadering overschrijden de $1\sigma$-foutbanden, vooral bij gelijke roodverschuivingen in hogere binnen en bij lage $\ell$-waarden.
• De $\chi^2$-analyse toont een enorme afwijking: de Limber-benadering resulteert in een $\Delta\chi^2 \sim 410$ met een p-waarde $< 0.001$. Dit betekent dat de nullhypothese (dat de benadering consistent is met de ware data) wordt verworpen met meer dan 99% zekerheid.
• De flat-sky-benadering presteert uitstekend: de $\Delta\chi^2$ is gemiddeld slechts $\sim 5$ met een p-waarde van $0.131$. De verdeling is consistent met nul, wat aangeeft dat de benadering de spectra nauwkeurig herstelt zonder bias.

3. Fysieke Oorzaak van de Afwijking:
De onnauwkeurigheid van de Limber-benadering in de DR3-configuratie wordt voornamelijk veroorzaakt door het falen om de Roodverschuivingsruimtestorringen (RSD) correct te modelleren. Bij smalle binnen wordt de RSD-bijdrage vergelijkbaar met de dichtheidsbijdrage (D-D) op grote schalen, en omdat de Limber-benadering de parallelle component $k_\parallel$ negeert, worden deze termen verkeerd geschat.

Significantie en Conclusie

• Noodzaak van een nieuwe benadering: Voor de uiteindelijke analyse van Euclid (DR3) met 13 smalle binnen is de Limber-benadering te grof en zal deze leiden tot gefabriceerde resultaten en bias in de kosmologische parameters. De flat-sky-benadering is de vereiste vervanging.
• Nauwkeurigheid: De flat-sky-benadering kan de hoekvermogensspectra van het aantal sterrenstelsels herstellen met een nauwkeurigheid van beter dan 5% ten opzichte van de volledige berekening.
• Berekeningskosten: Hoewel de flat-sky-benadering nauwkeuriger is, is deze nog steeds langzamer dan de Limber-benadering (ongeveer 2x sneller dan de volledige berekening, terwijl Limber 200x sneller is). De auteurs stellen echter vast dat verdere optimalisaties mogelijk zijn (bijvoorbeeld door pre-computing van de $k$-integralen of het combineren van benaderingen per schaal).
• Toekomstperspectief: De flat-sky-benadering biedt een haalbare weg voor de data-analyse van Euclid, waardoor de berekening van de volledige vermogensspectra in harmonische ruimte mogelijk wordt zonder de noodzaak van de onnauwkeurige Limber-benadering voor smalle binnen. Dit is ook relevant voor toekomstige 6x2pt-analyses die spectroscopische en fotometrische data combineren.

Samenvattend concludeert het artikel dat de flat-sky-benadering een cruciale verbetering is voor de Euclid-missie, specifiek voor de analyse van de finale data release met hoge resolutie in roodverschuiving, waar de standaard Limber-benadering niet langer toereikend is.