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

Titolo: L'approssimazione del cielo piatto per il clustering delle galassie fotometriche di Euclid

1. Il Problema

Il satellite Euclid genererà un catalogo di circa $1,5 \times 10^9$ galassie con redshift fotometrici. Per sfruttare al meglio questi dati, è necessario calcolare gli spettri di potenza angolari ($C_\ell$) in modo sia veloce che accurato per l'inferenza dei parametri cosmologici (ad esempio tramite analisi MCMC).
Attualmente, l'approssimazione di Limber è lo standard per le survey fotometriche perché offre un enorme vantaggio computazionale (circa 100 volte più veloce del calcolo completo). Tuttavia, è noto che l'approssimazione di Limber fallisce quando le galassie sono distribuite in bin di redshift stretti (narrow bins), sottostimando o sovrastimando il segnale per multipoli bassi ($\ell \lesssim 100$) con errori che possono superare un fattore 2.
Il problema centrale di questo studio è determinare se l'approssimazione di Limber è sufficientemente accurata per le diverse configurazioni di rilascio dati di Euclid (DR1 e DR3) o se è necessario adottare un approccio alternativo, come l'approssimazione del cielo piatto (flat-sky), per evitare bias nelle stime cosmologiche.

2. Metodologia

Gli autori hanno confrontato tre metodi di calcolo degli spettri di potenza angolare per il numero di galassie:

1. Calcolo Completo (Full Calculation): Utilizza il codice Boltzmann CLASS per calcolare gli spettri esatti integrando su tutte le scale, inclusi effetti relativistici (distorsioni spaziali, lente gravitazionale, magnificazione).
2. Approssimazione di Limber: Sostituisce l'integrale sul numero d'onda $k$ con una funzione delta di Dirac, riducendo l'integrale doppio su $z, z'$ a uno singolo. È molto veloce ma trascura la componente parallela di $k$.
3. Approssimazione del Cielo Piatto (Flat-Sky): Approssima la sfera celeste come un piano, sostituendo gli armonici sferici con trasformate di Fourier 2D. Mantiene la componente parallela di $k$ e utilizza una formulazione che permette di pre-calcolare e memorizzare l'integrale su $k$, migliorando l'accuratezza su larga scala rispetto a Limber.

Configurazioni di Studio:

• DR1 (First Data Release): Configurazione ottimistica con 6 bin di redshift equi-populati e un taglio di magnitudine $I_E < 23.5$.
• DR3 (Third/Final Data Release): Configurazione finale con 13 bin di redshift equi-populati e un taglio di magnitudine più profondo $I_E < 24.5$.

Validazione:
È stato utilizzato il catalogo simulato "Flagship" di Euclid per definire le distribuzioni di redshift, i bias delle galassie e i bias di magnificazione. Per valutare le prestazioni, sono state generate 1000 realizzazioni di rumore gaussiano basate sulla matrice di covarianza teorica. È stato calcolato il $\chi^2$ per ciascuna approssimazione confrontandolo con il risultato "vero" (calcolo completo) per determinare la probabilità che le approssimazioni siano compatibili con i dati osservati.

3. Contributi Chiave

• Estensione dell'approssimazione Flat-Sky: Gli autori applicano una versione "ricalibrata" dell'approssimazione del cielo piatto (basata su Gao et al. 2024) specificamente per i dati fotometrici di Euclid, includendo termini di distorsione spaziale (RSD) e effetti di lente gravitazionale.
• Analisi Comparativa su Scenari Reali: Prima valutazione dettagliata delle prestazioni di queste approssimazioni su scenari realistici di Euclid (DR1 vs DR3), evidenziando come la precisione richiesta cambi drasticamente con la risoluzione in redshift.
• Quantificazione del Bias: Dimostrazione quantitativa che l'approssimazione di Limber introduce errori sistematici significativi nei bin stretti della DR3, che superano gli errori statistici attesi ($1\sigma$).

4. Risultati

• Caso DR1 (6 bin): L'approssimazione di Limber è sufficientemente accurata. Gli errori relativi rispetto al calcolo completo sono all'interno della banda di errore osservativa ($1\sigma$). Il test $\chi^2$ mostra che l'approssimazione di Limber non può essere rifiutata con alta confidenza ($p \approx 0.038$, vicino ma accettabile nel contesto delle incertezze), e non è necessario usare il calcolo flat-sky per questa fase.
• Caso DR3 (13 bin): L'approssimazione di Limber fallisce. I bin più stretti rendono l'approssimazione troppo grezza. Gli errori relativi superano il 5-10% e spesso escono dalla banda $1\sigma$, specialmente per i bin a redshift più alti e per le correlazioni cross-bin. Il test $\chi^2$ rifiuta l'ipotesi nulla per Limber con una p-value $< 0.001$ (99.9% di confidenza), indicando un bias inaccettabile.
• Prestazioni Flat-Sky: L'approssimazione del cielo piatto è altamente accurata in entrambi i casi, recuperando gli spettri di potenza con un errore inferiore al 5% rispetto al calcolo completo. Per la DR3, il test $\chi^2$ non rifiuta l'ipotesi nulla ($p \approx 0.131$), confermando che l'approssimazione non introduce bias significativi.
• Efficienza Computazionale:
  • Limber: ~200 volte più veloce del calcolo completo.
  • Flat-Sky (implementazione attuale): ~2 volte più veloce del calcolo completo (un miglioramento modesto rispetto a Limber, ma sufficiente per evitare bias).
  • Nota: L'implementazione flat-sky richiede ancora un'integrazione tripla, ma può essere ottimizzata (es. pre-calcolo dell'integrale $k$) per ridurre i tempi in analisi MCMC.

5. Significato e Implicazioni

Questo studio ha implicazioni cruciali per l'analisi dei dati di Euclid:

1. Scelta del Metodo per DR3: Per il rilascio finale dei dati (DR3) con 13 bin, l'uso dell'approssimazione di Limber porterebbe a risultati cosmologici distorti. È necessario adottare l'approssimazione del cielo piatto (o calcoli completi ottimizzati) per garantire l'accuratezza richiesta.
2. Strategia Ibrida: Poiché le deviazioni di Limber sono maggiori per $\ell \lesssim 20$, si potrebbe considerare una strategia ibrida: usare Limber per le scale piccole e il calcolo flat-sky (o completo) per le scale grandi. Inoltre, per i termini di lente gravitazionale (L-L) e densità-lente (D-L), Limber rimane accurato anche nei bin stretti.
3. Ottimizzazione del Codice: Sebbene il calcolo flat-sky sia attualmente più lento di Limber, la sua accuratezza lo rende indispensabile. Gli autori suggeriscono che ulteriori ottimizzazioni (simili a quelle del codice BLAST) sono necessarie per rendere il calcolo flat-sky competitivo in termini di velocità per le analisi MCMC su larga scala.
4. Impatto sulla Survey Spettroscopica: I risultati suggeriscono che anche per l'analisi combinata (6x2pt) che include i dati spettroscopici, se si utilizzano bin di redshift stretti, l'approssimazione di Limber sarà insufficiente, rendendo il flat-sky un'alternativa promettente anche per quel settore.

In sintesi, il lavoro stabilisce un confine chiaro: Limber va bene per i primi dati (bin larghi), ma il cielo piatto è la soluzione obbligata per l'analisi finale di Euclid (bin stretti) per evitare errori sistematici che comprometterebbero la precisione cosmologica della missione.