Towards an application of fourth-order shear statistics I. The information content of $\langle M_\mathrm{ap}^4 \rangle$

1. Il Problema e il Contesto

Negli ultimi decenni, l'analisi del lensing debole (cosmic shear) ha fornito informazioni cruciali sulla struttura su larga scala dell'universo e sui parametri cosmologici. Tradizionalmente, le analisi si sono basate sulle statistiche del secondo ordine (funzioni di correlazione a due punti, 2PCF), che catturano le informazioni gaussiane del campo di materia. Tuttavia, con l'avvento di survey di "Stage III" (come DES, KiDS, HSC) e l'imminente arrivo di survey di "Stage IV" (come Euclid), la precisione osservativa è aumentata al punto da rendere rilevanti le statistiche di ordine superiore.

Le statistiche di terzo ordine (bispettri, 3PCF) sono state già utilizzate per catturare le non-gaussianità del campo di materia. Questo lavoro si concentra sul passo successivo: lo sviluppo e l'applicazione delle statistiche di quarto ordine (4PCF e aperture di quarto ordine). L'obiettivo è determinare se queste statistiche possano:

1. Rompere le degenerazioni tra i parametri cosmologici rispetto alle statistiche di ordine inferiore.
2. Fornire un nuovo test per i sistematici e i parametri di disturbo (nuisance parameters), come l'allineamento intrinseco o i modelli barionici, anche nel caso in cui siano degeneri con gli ordini precedenti.

2. Metodologia

Gli autori hanno sviluppato un quadro teorico e numerico completo per l'analisi delle statistiche di quarto ordine, strutturato come segue:

• Quadro Teorico:

  • Hanno generalizzato la relazione tra le statistiche di apertura di ordine $n$ ($\langle M_n^{ap} \rangle$) e le funzioni di correlazione a $n$ punti (nPCF) del campo di shear, nonché gli $n$-polispettri di convergenza.
  • Per il caso specifico del quarto ordine, hanno derivato la forma funzionale dei filtri necessari per trasformare le 4PCF nelle statistiche di apertura.
  • Hanno definito le componenti naturali della 4PCF (8 componenti complesse indipendenti) e le relative proiezioni (in particolare la proiezione "×") per gestire le coordinate e le singolarità.
  • Hanno derivato le espressioni per separare i modi $E$ (segnale di lensing puro) e $B$ (rumore o sistematici) a partire dalle 8 misure di aperture complesse.

• Pipeline Numerica e Integrazione:

  • Hanno sviluppato una pipeline di integrazione numerica per calcolare le statistiche di apertura integrando le 4PCF su configurazioni di quadrilateri.
  • Hanno testato la precisione di questa pipeline utilizzando un Campo Casuale Gaussiano (GRF). In un GRF, le statistiche di quarto ordine sono puramente "disconnesse" (possono essere espresse come combinazioni di statistiche del secondo ordine), permettendo un confronto diretto con i risultati attesi dalle statistiche 2PCF.
  • Hanno studiato gli effetti della discretizzazione (binning) e dell'intervallo di integrazione radiale, identificando problemi di mescolamento tra modi $E$ e $B$ dovuti a intervalli di integrazione finiti e binning insufficiente.

• Previsione Cosmologica (Fisher Forecast):

  • Hanno eseguito una previsione basata sulla matrice di Fisher per valutare il potere vincolante delle statistiche di quarto ordine in un setup simile a DES-Y3.
  • Hanno utilizzato simulazioni cosmo-SLICS per stimare le derivate delle osservabili rispetto ai parametri cosmologici ($\Omega_m, S_8, h, w_0$).
  • Hanno introdotto un metodo innovativo per stimare le derivate da un insieme di simulazioni con distribuzioni cosmologiche arbitrarie (latin hypercube), minimizzando una funzione di perdita per costruire un piano tangente robusto, superando l'instabilità delle differenze finite su dati rumorosi.

3. Risultati Chiave

• Precisione Numerica:

  • Analizzando e mitigando gli effetti numerici (come il mescolamento $E/B$ dovuto a intervalli di integrazione limitati e la discretizzazione), gli autori hanno raggiunto una precisione al 2% sulle statistiche di apertura di quarto ordine per un campo GRF.
  • Questa precisione è ben al di sotto del budget di rumore previsto per le survey di Stage IV, rendendo la pipeline affidabile per l'uso su dati reali.
  • È stato dimostrato che l'integrazione delle 4PCF all'interno di bin radiali (simulando la media reale dei dati) riduce la necessità di una risoluzione di binning estremamente fine, rendendo il calcolo fattibile anche per dati di Stage III.

• Contenuto Informativo (Fisher Forecast):

  • Nel contesto di un setup simile a DES-Y3 (senza tomografia e con statistiche di apertura a scala uguale), l'aggiunta delle statistiche di quarto ordine a un'analisi congiunta di secondo e terzo ordine ($\langle M_2^{ap} \rangle + \langle M_3^{ap} \rangle$) porta a un miglioramento minimo nel potere vincolante sui parametri cosmologici.
  • Gli autori attribuiscono questo risultato a diverse limitazioni delle assunzioni fatte:
    1. L'assunzione di una distribuzione gaussiana delle osservabili (non valida per ordini superiori).
    2. L'uso di statistiche a scala uguale, che catturano solo una frazione dell'informazione totale disponibile (le statistiche multi-scala potrebbero essere più informative).
    3. La marginalizzazione su parametri come $h$ e $w_0$ e la mancanza di tomografia, che riducono la capacità delle statistiche di ordine superiore di rompere le degenerazioni.
    4. Il fatto che le statistiche di ordine superiore decadono rapidamente con la scala, concentrando l'informazione su scale piccole difficili da modellare e non sempre accessibili nelle simulazioni usate per le derivate.

4. Contributi Principali

1. Formalismo Teorico: Derivazione completa dei filtri per le statistiche di apertura di quarto ordine in termini di 4PCF e trispettri.
2. Pipeline di Integrazione: Sviluppo e validazione di un metodo numerico robusto per calcolare le statistiche di apertura da 4PCF, con una stima accurata degli errori sistematici numerici.
3. Metodo di Derivazione: Proposta di una tecnica basata sulla minimizzazione della perdita per stimare le derivate cosmologiche da suite di simulazioni irregolarmente campionate, applicabile a spazi parametrici ad alta dimensionalità.
4. Valutazione dell'Impatto: Una valutazione realistica del valore aggiunto delle statistiche di quarto ordine per survey attuali (DES-Y3) rispetto a quelle future (Euclid).

5. Significato e Conclusioni

Il lavoro conclude che, sebbene le statistiche di quarto ordine contengano informazioni non-gaussiane preziose, il loro impatto immediato sulle survey di Stage III (come DES-Y3) è limitato a causa delle incertezze statistiche e delle limitazioni dei modelli teorici attuali.

Tuttavia, il quadro teorico e la pipeline numerica sviluppati sono fondamentali per le survey di Stage IV (come Euclid). In futuro, con una precisione osservativa molto superiore, le statistiche di quarto ordine potrebbero:

• Rompere le degenerazioni tra parametri cosmologici.
• Funzionare come un potente test di coerenza per i sistematici osservativi e i parametri di disturbo, anche se risultassero degeneri con gli ordini inferiori.

Questo studio getta le basi per l'applicazione pratica di queste statistiche avanzate, che verranno sfruttate nella parte II dello studio (Porth et al. 2025) per la rilevazione effettiva delle statistiche di quarto ordine sui dati reali del Dark Energy Survey.