Accurate Galaxy Cluster Shear and Mass Calibration for LSST with AnaCal

Sintesi Tecnica: Calibrazione Accurata dello Shear e della Massa dei Cluster di Galassie per LSST con AnaCal

Problematica
La cosmologia dei cluster di galassie si basa sull'instaurazione di una relazione accurata tra gli osservabili dei cluster e le loro masse, un processo noto come calibrazione della massa. Il lensing gravitazionale debole è lo strumento standard per questa calibrazione perché sonda la massa totale (incluso la materia oscura) senza fare affidamento su assunzioni di equilibrio dinamico o termico. Tuttavia, una sfida significativa sorge nell'ambiente dei cluster: il segnale di shear vicino ai centri dei cluster può raggiungere il regime non lineare (ad esempio, $|g| \sim 0,15$), dove molti stimatori di shear si basano su assunzioni perturbative che richiedono una validazione esplicita. Inoltre, il lensing dei cluster introduce complessità uniche non presenti nelle misurazioni standard dello shear cosmico, tra cui:

• Regime di Shear Elevato: Lo shear tangenziale nelle regioni interne ($\sim 0,2$ Mpc) è diverse volte superiore ai valori tipicamente utilizzati per i test di calibrazione basati sulla perturbazione.
• Deflessione Inomogenea: Il campo di lensing del cluster sposta le posizioni apparenti delle galassie sullo sfondo, alterando le loro configurazioni di sovrapposizione e blending, il che impatta il rilevamento e la misurazione della forma.
• Dipendenza Radiale: Il campo di convergenza varia radialmente, potendo guidare una risposta media di shear radialmente decrescente.
• Effetti di Ordine Superiore: Vicino al centro del cluster, distorsioni di ordine superiore (ad esempio, la flessione) o la luce anisotropa dagli aloni intra-cluster possono introdurre segnali spuri.

Sebbene la LSST Dark Energy Science Collaboration (DESC) abbia selezionato AnaCal e metadetection come metodi fiduciali di stima dello shear, le loro prestazioni non erano state rigorosamente testate sotto le specifiche condizioni di alto shear e deflessione inomogenea caratteristiche dei lenti dei cluster.

Metodologia
Gli autori utilizzano simulazioni d'immagine rappresentative dei 10 anni di dati LSST per testare lo stimatore di shear AnaCal. La metodologia prevede:

1. Framework di Simulazione: Utilizzando descwl-shear-sims e xlens, il team ha generato simulazioni d'immagine realistiche del lensing dei cluster. Hanno modellato 9 masse di aloni da $10^{14} M_\odot$ a $10^{15} M_\odot$ come aloni NFW a un redshift del lens $z_l = 0,25$, con galassie sorgente a $z_s = 1,0$.

   • Le simulazioni includevano 2178 realizzazioni per massa, coprendo circa 1700 deg$^2$ totali.
   • Per isolare gli effetti di stima dello shear, le simulazioni hanno escluso le galassie membri dei cluster, la luce intra-cluster e gli errori di calibrazione del redshift, concentrandosi strettamente sull'impatto dell'alto shear e del campo di deflessione inomogeneo.
   • È stato utilizzato un PSF di tipo Moffat circolare costante ($FWHM = 0,8''$) per rappresentare la qualità dell'immagine attesa dell'Osservatorio Rubin.

2. Implementazione di AnaCal: Lo studio impiega il framework AnaCal, che deriva analiticamente la risposta di shear dell'intero processo di misurazione (osservabili, rilevamento e selezione) utilizzando l'auto-differenziazione e una base di shapelet.

   • Osservabili: Le immagini sono proiettate su kernel di rilevamento e funzioni di base di shapelet polari.
   • Calcolo della Risposta: La risposta di shear è calcolata differenziando le osservabili rispetto a perturbazioni di shear controllate.
   • Correzione del Rumore: Viene applicato uno schema di "renoising" per correggere il bias del rumore, aggiungendo realizzazioni di rumore sintetico alle immagini.
   • Pesatura: Vengono applicati pesi di selezione e rilevamento basati sul rapporto segnale-rumore (SNR) e sulla dimensione della galassia, progettati per essere invarianti sotto rotazioni di 90 gradi.

3. Analisi e Fitting:

   • Stima dello Shear: Le ellitticità tangenziali e trasversali sono calcolate, e la matrice di risposta viene ruotata nel frame tangenziale. Fondamentalmente, gli autori utilizzano la risposta di shear in ogni bin angolare separatamente, piuttosto che mediarla su tutti i bin, per evitare il bias derivante dalla dipendenza radiale del campo di convergenza.
   • Fitting della Massa: Un approccio di forward-modeling adatta i modelli NFW osservati al profilo di ridotto shear tangenziale ($g_T$) utilizzando una funzione di verosimiglianza che incorpora una matrice di covarianza analitica derivata dalle previsioni a 10 anni di LSST.
   • Quantificazione del Bias: Il bias moltiplicativo dello shear ($m$) viene adattato alla forma $m = a + b(g_{true})^2$, testando per bias lineari e non lineari di terzo ordine.

Risultati Chiave

• Recupero dello Shear: AnaCal recupera lo shear di input con un bias minimo anche a livelli di shear moderatamente elevati ($|g| \sim 0,15$).
• Trend della Risposta Radiale: Lo studio conferma una risposta di shear media radialmente decrescente, guidata dalla dipendenza radiale della convergenza. Gli autori dimostrano che se questo effetto non viene modellato (ovvero, se viene utilizzata una risposta media globale), esso introduce un bias nell'inferenza dello shear.
• Bias di Terzo Ordine: Un bias di stima dello shear positivo viene rilevato al terzo ordine nel ridotto shear vicino al centro del cluster.
• Accuratezza della Calibrazione della Massa: Nonostante la presenza di regimi di alto shear e bias di terzo ordine, l'impatto sulla stima finale della massa è trascurabile. Per le masse dei cluster di lente nell'intervallo $[10^{14} M_\odot, 10^{15} M_\odot]$ (adottando un taglio di scala di $\sim 0,2$ Mpc a $z=0,25$), il bias medio della massa del cluster risultante è $0,24 \pm 0,26\%$ in condizioni ideali. Questo piccolo bias è attribuito al fatto che solo una piccola frazione di galassie si trova nel regime di alto shear, e tali misurazioni sono pesate meno dalla matrice di covarianza.

Significato e Rivendicazioni
L'articolo sostiene che questi risultati dimostrano che AnaCal è uno strumento robusto per l'accurata calibrazione della massa dei cluster nell'era di LSST. Nello specifico:

• Valida che la stima dello shear basata su perturbazioni (tramite AnaCal) rimane accurata anche nel regime di alto shear tipico dei centri dei cluster, a condizione che la dipendenza radiale della risposta sia correttamente modellata.
• Conferma che le non linearità e gli effetti di alto shear non degradano significativamente la precisione delle stime di massa per campioni di cluster cosmologicamente rilevanti.
• Il lavoro fornisce un passaggio critico di validazione per la LSST DESC, assicurando che i metodi di stima dello shear selezionati possano soddisfare i rigorosi requisiti per la cosmologia dei cluster senza introdurre errori sistematici significativi.

Gli autori notano esplicitamente che le loro simulazioni non includevano effetti da membri del cluster, luce intra-cluster o distorsioni di ordine superiore come la flessione, isolando le prestazioni dello stimatore di shear contro la sfida primaria dell'alto shear e dei campi di deflessione inomogenei.