Euclid preparation. The impact of redshift interlopers on the two-point correlation function analysis

Titolo: L'impatto degli interlopers di redshift sull'analisi della funzione di correlazione a due punti (Euclid)

1. Il Problema

La missione spaziale Euclid dell'ESA mira a mappare la struttura su larga scala dell'universo misurando il redshift spettroscopico di milioni di galassie a riga di emissione (ELG) tramite la linea H$\alpha$. Tuttavia, l'uso della spettroscopia senza slit (slitless) con risoluzione media-bassa ($R < 1000$) e un intervallo di banda limitato introduce ambiguità nell'identificazione delle linee spettrali.
Questo porta alla presenza di galassie "interlopers" nel campione spettroscopico, ovvero oggetti il cui redshift misurato ($z_{meas}$) è catastroficamente errato rispetto a quello vero ($z_{true}$). Esistono due tipi principali di interlopers:

• Interlopers di linea: Galassie con linee di emissione reali (es. [OIII] $\lambda 5008$ o [SIII] $\lambda 9531$) erroneamente identificate come H$\alpha$.
• Interlopers di rumore: Oggetti privi di linee di emissione significative o stelle, dove un picco di rumore nello spettro viene scambiato per una linea di emissione.

Questi errori di redshift causano uno spostamento radiale delle posizioni comoving delle galassie, diluendo e distorcendo la statistica di clustering (in particolare la Funzione di Correlazione a Due Punti, 2PCF) e potenzialmente introducendo bias nei parametri cosmologici inferiti, come il tasso di crescita delle strutture ($f\sigma_8$) e i parametri di Alcock-Paczynski ($\alpha_{\parallel}, \alpha_{\perp}$).

2. Metodologia

Lo studio si basa sull'analisi di 1000 cataloghi sintetici (EuclidLargeMocks) che mimano l'area e la funzione di selezione della prima release di dati (DR1) di Euclid. Questi mock incorporano un modello realistico di errori di redshift derivato da simulazioni end-to-end, includendo le frazioni attese di interlopers di linea e di rumore.

L'approccio metodologico include:

• Formalismo Teorico: Derivazione di un'espressione analitica per la 2PCF misurata ($\xi_m$) in presenza di contaminanti, decomponendola in termini di auto-correlazione delle galassie corrette, auto-correlazione degli interlopers (di linea e di rumore) e termini di cross-correlazione tra le diverse popolazioni (Eq. 17).
• Analisi dei Residui: Valutazione dell'ampiezza e della rilevanza di ciascun termine nell'equazione della 2PCF misurata, confrontando modelli teorici di complessità crescente con i dati simulati.
• Inferenza Cosmologica: Esecuzione di analisi MCMC (Monte Carlo Markov Chain) per stimare i parametri cosmologici ($f\sigma_8$, $b\sigma_8$, $\sigma_p$) e i parametri AP ($\alpha_{\parallel}, \alpha_{\perp}$). Vengono testati diversi modelli di fitting:
  1. Modello minimale: considera solo le galassie corrette con un fattore di attenuazione (damping) dell'ampiezza.
  2. Modello intermedio: include l'auto-correlazione degli interlopers di linea.
  3. Modello completo: include anche le cross-correlazioni (es. corrette-noise).
• Confronto: I risultati ottenuti su campioni contaminati vengono confrontati con un "caso di riferimento" (campione puro al 100%) per quantificare il bias sistematico introdotto da modelli incompleti.

3. Contributi Chiave

• Quantificazione delle Contribuzioni: Il lavoro dimostra che, sebbene la popolazione di galassie corrette sia dominante, i termini di auto-correlazione degli interlopers di linea e la cross-correlazione tra galassie corrette e interlopers di rumore sono i contributi secondari più significativi, con un'importanza che varia in funzione del redshift.
• Validazione di Modelli Semplificati: Lo studio identifica che un modello minimale, che tratta gli interlopers semplicemente come un'attenuazione dell'ampiezza del segnale di clustering (senza modellare esplicitamente la loro struttura interna o le cross-correlazioni), è sufficiente per recuperare i parametri cosmologici corretti nel contesto della DR1.
• Robustezza dei Parametri AP: Viene dimostrato che i parametri di Alcock-Paczynski sono estremamente robusti alla presenza di interlopers, poiché il termine di "broadband" polinomiale nei modelli BAO assorbe efficacemente gli effetti residui.

4. Risultati Principali

• Parametri di Crescita ($f\sigma_8$):
  • L'adozione di un modello minimale (solo attenuazione dell'ampiezza) induce un errore sistematico sulla stima di $f\sigma_8$ compreso tra 1% e 3%, a seconda del redshift.
  • Questo errore sistematico è inferiore all'errore statistico atteso per l'analisi della DR1 di Euclid.
  • L'aggiunta di termini di auto-correlazione degli interlopers o di cross-correlazione non migliora significativamente la precisione o l'accuratezza per la DR1, sebbene possa diventare rilevante per analisi future (es. DR3) con statistiche molto più elevate.
• Parametri di Alcock-Paczynski ($\alpha_{\parallel}, \alpha_{\perp}$):
  • L'analisi BAO (Oscillazioni Acustiche Barioniche) mostra una robustezza eccezionale.
  • L'uso di un modello che ignora completamente la presenza degli interlopers non introduce bias significativi nei valori di $\alpha_{\parallel}$ e $\alpha_{\perp}$.
  • L'incertezza statistica aumenta leggermente (fino al ~20% in alcuni casi di redshift più alti dove gli interlopers di linea sono abbondanti), ma il bias sistematico rimane trascurabile rispetto all'errore statistico atteso sia per la DR1 che per la DR3.
• Dipendenza dal Redshift: L'impatto varia con il redshift. A redshift più bassi ($z \approx 0.9-1.1$), gli interlopers di rumore dominano e la cross-correlazione è significativa; a redshift più alti ($z \approx 1.3-1.5$), gli interlopers di linea ([OIII], [SIII]) diventano più rilevanti, distorcendo la forma della 2PCF (spostamento e allargamento del picco BAO), ma senza compromettere l'estrazione dei parametri cosmologici con i modelli attuali.

5. Significato e Conclusioni

Questo lavoro fornisce una base fondamentale per la strategia di analisi dei dati della prima release (DR1) della missione Euclid.

• Efficienza Operativa: Dimostra che non è necessario sviluppare modelli estremamente complessi per gestire gli interlopers nella fase iniziale della missione. Un approccio "minimale" che corregga solo l'attenuazione dell'ampiezza del clustering è sufficiente per garantire risultati cosmologici affidabili.
• Pianificazione Futura: Sebbene i risultati siano ottimi per la DR1, lo studio avverte che per le release future (DR3), con errori statistici ridotti di un fattore $\sqrt{6}$, la precisione richiesta potrebbe necessitare di modelli più sofisticati che includano esplicitamente le cross-correlazioni e la struttura degli interlopers.
• Monitoraggio della Purezza: L'accuratezza finale dipenderà dalla capacità di stimare con precisione le frazioni di contaminazione (obiettivo di purezza < 20%). Se la contaminazione reale nella DR1 supera le aspettative (come suggerito da dati preliminari Q1), sarà necessario rivedere i cataloghi simulati e aggiornare le stime di bias.

In sintesi, lo studio rassicura la comunità scientifica sulla fattibilità di ottenere misure cosmologiche di alta precisione con Euclid, anche in presenza di un livello di contaminazione di redshift realistico, confermando la robustezza delle tecniche di analisi della struttura su larga scala.