Integrated Sachs-Wolfe maps from the Gower Street $w$CDM simulations

Sintesi Tecnica: Mappe dell'effetto Sachs–Wolfe integrato da simulazioni wCDM di Gower Street

Problema e Motivazione
L'effetto Sachs–Wolfe integrato (ISW) lineare a tempi tardivi funge da sonda diretta dell'accelerazione cosmica e della natura dell'energia oscura. Tuttavia, rilevare queste deboli anisotropie secondarie di temperatura nella Radiazione Cosmica di Fondo (CMB) richiede previsioni teoriche accurate su un ampio spettro di modelli cosmologici. Mentre le analisi precedenti si sono concentrate prevalentemente sul modello standard $\Lambda$CDM, la natura dell'energia oscura rimane sconosciuta, con alternative teoriche che spaziano dai modelli fantasma ($w < -1$) a quelli di quintessenza ($w > -1$). Gli strumenti esistenti, come il pacchetto pyGenISW, erano in precedenza limitati alle cosmologie $\Lambda$CDM. Per far fronte alla necessità di una modellazione teorica robusta su uno spazio dei parametri dell'energia oscura più ampio, questo lavoro mira a generare mappe ISW a tutto cielo per una serie di cosmologie $w$CDM, consentendo analisi che vanno oltre il modello standard.

Metodologia
Gli autori hanno esteso il pacchetto open-source pyGenISW per gestire i modelli $w$CDM e i neutrini massivi, superando la precedente limitazione al solo $\Lambda$CDM. La metodologia si basa sulla suite di simulazioni Gower Street (GS), composta da 791 simulazioni N-body a tutto cielo, contenenti solo materia oscura, eseguite con pkdgrav3. Queste simulazioni coprono una cosmologia $w$CDM piatta con parametri di equazione di stato dell'energia oscura $-1.79 \lesssim w \lesssim -0.34$, variando sette parametri cosmologici tra cui $\Omega_m$, $\sigma_8$, $w$ e la massa del neutrino $m_\nu$.

La pipeline di costruzione delle mappe ISW prevede i seguenti passaggi:

1. Dati in ingresso: La pipeline utilizza cataloghi di coni di luce dalle simulazioni GS, convertendo i conteggi grezzi di particelle in mappe di contrasto di densità per ogni strato fino al redshift $z \simeq 1.0$.
2. Quadro teorico: Per calcolare le necessarie storie di espansione e crescita, gli autori hanno implementato tre metodi complementari per calcolare il fattore di crescita lineare ($D(z)$) e il tasso di crescita ($f(z)$):
   • Un approccio basato su Equazioni Differenziali Ordinarie (ODE) (fiduciale), che risolve l'equazione di perturbazione standard.
   • Un approccio basato su CAMB, che estrae i fattori di crescita dallo spettro di potenza della materia lineare.
   • Un'approssimazione analitica che utilizza il formalismo dell'indice di crescita.
     L'approccio ODE è stato validato contro la libreria pyCCL ed è risultato in eccellente accordo.
3. Ricostruzione delle mappe: Utilizzando pyGenISW esteso, la pipeline esegue una trasformata sferica di Bessel discreta (SBT) sulle mappe di contrasto di densità. Ciò comporta la proiezione dei campi di densità su una base di funzioni di Bessel sferiche adattate all'intervallo radiale finito del cono di luce, applicando una condizione al contorno di Neumann. Le anisotropie di temperatura ISW vengono quindi ricostruite integrando l'evoluzione temporale del potenziale gravitazionale lungo la linea di vista.
4. Validazione: La pipeline è stata validata confrontando gli spettri di potenza angolare misurati (auto-correlazione ISW, auto-correlazione galassie e incrocio ISW-galassie) con le previsioni teoriche ottenute da TheoryCL e pyCCL.

Contributi Chiave

• Estensione del codice: Gli autori hanno esteso con successo pyGenISW per supportare le cosmologie $w$CDM e i neutrini massivi, rendendo il codice e i dati ISW risultanti disponibili al pubblico.
• Utilizzo della suite di simulazioni: Hanno generato mappe ISW a tutto cielo per 777 simulazioni valide su 791 disponibili nella suite Gower Street, coprendo un'ampia gamma di equazioni di stato dell'energia oscura.
• Validazione della modellazione della crescita: Lo studio ha confermato che l'approccio basato su ODE per il calcolo dei fattori di crescita produce risultati coerenti con pyCCL nell'intervallo di redshift rilevante ($0 \le z \le 1$) per tutti i modelli $w$CDM testati.

Risultati

• Accordo degli spettri di potenza: Gli spettri di potenza angolare ISW misurati ($C_{\ell}^{ISW-ISW}$) e le correlazioni incrociate ISW-densità ($C_{\ell}^{g-ISW}$) mostrano un eccellente accordo con le previsioni della teoria lineare (calcolate tramite pyCCL e TheoryCL) nell'intervallo di multipolo $2 \le \ell \le 200$. Le discrepanze a $\ell < 10$ sono attribuite alla dimensione finita della scatola di simulazione, mentre i deficit a $\ell \ge 100$ sono legati agli effetti di risoluzione e al collasso delle assunzioni della teoria lineare.
• Dipendenza dal modello: L'intensità del segnale ISW varia significativamente con l'equazione di stato dell'energia oscura. I modelli simili alla quintessenza ($w > -1$, ad es. Sim 401 con $w \approx -0.34$) mostrano ampiezze ISW più elevate e distribuzioni di fluttuazioni di temperatura più ampie rispetto ai modelli fantasma ($w < -1$, ad es. Sim 127 con $w \approx -1.79$). Ciò è coerente con il decadimento potenziato dei potenziali gravitazionali a tempi tardivi negli scenari di quintessenza.
• Evoluzione del redshift: Sebbene i modelli con diversi valori di $w$ siano quasi indistinguibili a redshift molto bassi ($z < 0.2$) a causa di tassi di decadimento potenziale attuali simili (controllati principalmente da $\Omega_m$), chiare distinzioni emergono a redshift intermedi ($0.3 \lesssim z \lesssim 1.0$) dove il kernel ISW è più sensibile alla storia della crescita.
• Estremi di densità: L'analisi degli estremi di densità (vuoti e superammassi) rivela che il segnale ISW è potenziato in queste regioni. Il modello di quintessenza produce i maggiori spostamenti di temperatura ($\Delta T_{ISW} \sim \pm 15 \, \mu\text{K}$), mentre il modello fantasma mostra segnali più deboli ($\Delta T_{ISW} \sim \pm 5 \, \mu\text{K}$).
• Rapporti di potenza: Calcolando il rapporto degli spettri di potenza ISW rispetto a un modello $\Lambda$CDM fiduciale ($A_{ISW}^{Tg}$), gli autori hanno scoperto che i modelli con $w > -1$ e $\Omega_m$ basso possono produrre ampiezze coerenti con alcune affermazioni osservative di segnali in eccesso (ad es. $A_{ISW} \approx 5$), sebbene nessun modello in questa suite spieghi il riportato inversione di segno del segnale ISW a redshift specifici.

Significato e Affermazioni
Il paper conclude che la coerenza tra le mappe ISW $w$CDM generate e le previsioni teoriche conferma la robustezza della metodologia estesa. Gli autori affermano che questo quadro stabilisce uno strumento affidabile per le analisi teoriche e osservative ISW-Struttura su Grande Scala (LSS). Nello specifico, le mappe sono destinate a supportare:

• Calcoli di covarianza per sondaggi di prossima generazione.
• Vari statistiche basate su mappe, incluse analisi di stacking e inferenza basata su simulazioni.
• Previsioni teoriche per l'impronta ISW dipendente dal redshift su un ampio spazio dei parametri che va dai domini fantasma a quelli di quintessenza.

Il lavoro non afferma di risolvere le tensioni osservative esistenti (come l'inversione di segno del segnale ISW), ma fornisce piuttosto l'infrastruttura teorica necessaria per modellare i segnali ISW attraverso cosmologie diverse, facilitando i futuri confronti con i dati di sondaggi come Euclid, DESI e Rubin-LSST.