Primordial Physics in the Nonlinear Universe: Towards particle constraints using the Weak lensing, Thermal SZ, and X-ray fields

Sintesi Tecnica: Fisica Primordiale nell'Universo Non Lineare

Problematica
Le non-gaussianità primordiali (PNG) nel campo di densità iniziale fungono da sonda critica per la fisica delle particelle nell'Universo primordiale, potenzialmente rivelando scale di energia inaccessibili agli esperimenti terrestri. Sebbene la Radiazione Cosmica di Fondo (CMB) abbia fornito i vincoli più stringenti ad oggi, essa è limitata alle scale lineari. Al contrario, la formazione delle strutture nel tempo tardivo (ad esempio, lensing debole, ammassi di galassie) sonda il regime altamente non lineare dove le firme delle PNG sono amplificate, in particolare nelle code della distribuzione di densità (picchi e vuoti). Tuttavia, modellare questi regimi non lineari è complesso a causa della complessa interazione tra l'evoluzione gravitazionale e i processi astrofisici (feedback barionico, allineamenti intrinseci, foreground). Le analisi esistenti spesso si affidano a modelli analitici che non possono catturare appieno la dinamica non lineare o trattare gli osservabili in isolamento, fallendo nel sfruttare l'informazione sinergica disponibile attraverso molteplici lunghezze d'onda. Questo lavoro affronta la necessità di un framework basato su simulazioni e multi-lunghezza d'onda per vincolare le PNG utilizzando campi di lensing debole (WL), effetto Sunyaev-Zeldovich termico (tSZ) e raggi X, modellando coerentemente i parametri di disturbo (nuisance) astrofisici.

Metodologia
Gli autori impiegano un approccio di modellazione diretta basato su simulazioni utilizzando la suite di simulazioni N-body Ulagam, che genera coni di luce a cielo intero con $512^3$ particelle. Per gestire il costo computazionale della simulazione di una vasta gamma di modelli di PNG, gli autori utilizzano algoritmi del loro lavoro precedente (Paper I e II) che approssimano qualsiasi bispettro tramite un insieme finito di termini fattorizzabili, consentendo la generazione efficiente delle condizioni iniziali per sei distinti modelli inflazionistici: Locale, Equilatero, Quasi-Single Field (QSF), Scalar I (SI), Scalar II (SII) e Non-Bunch-Davies (NBD2).

L'innovazione principale risiede nella pipeline Vaanam, che genera cieli sintetici multi-lunghezza d'onda applicando modelli barionici semi-analitici (SAM) agli output delle simulazioni dark-matter-only.

• Modellazione Astrofisica: Gli autori utilizzano il framework BaryonForge per modellare gli effetti barionici (profili del gas, frazioni stellari, gas espulso) in modo coerente attraverso tutti gli osservabili. Ciò consente l'auto-calibrazione dei parametri di disturbo rompendo le degenerazioni tra effetti cosmologici e astrofisici.
• Osservabili:
  • Lensing Debole (WL): Modellato come campo di convergenza ($\kappa$), includendo allineamenti intrinseci (modello NLA), clustering delle sorgenti, ingrandimento (magnification) e rumore di forma.
  • SZ Termico (tSZ): Modellato tramite l'integrale lungo la linea di vista della pressione elettronica termica, incorporando componenti di pressione non termica e foreground (CIB, sorgenti radio) che sono esplicitamente correlati al campo di densità simulato.
  • Raggi X: Modellato come tassi di conteggio dei fotoni da ammassi di galassie, AGN e binarie di raggi X. Crucialmente, gli autori introducono un profilo di metallicità radialmente variabile per gli ammassi e modellano le popolazioni di AGN/binarie di raggi X tramite campionamento di Poisson dal campo di densità, garantendo le correlazioni spaziali con il segnale cosmologico.
• Statistiche: L'analisi comprime le mappe in statistiche riassuntive utilizzando il secondo e il terzo momento (e un sottoinsieme sparso di quarti momenti) dei campi, mediati su 10 scale angolari ($3.2' < \theta < 200'$).
• Inferenza: I vincoli sono derivati utilizzando la matrice di informazione di Fisher, marginalizzando su un set completo di parametri di disturbo che includono la cosmologia ($\Omega_m, \sigma_8$), gli allineamenti intrinseci, i parametri SAM barionici e le ampiezze dei foreground.

Contributi Chiave

1. Framework Unificato Multi-Lunghezza d'Onda: Il documento presenta la prima analisi con modellazione diretta congiunta di campi WL, tSZ e raggi X per i vincoli sulle PNG, correlando esplicitamente tutti i foreground (CIB, radio, AGN) con il sottostante campo di densità.
2. Modellazione Barionica Semi-Analitica: Estendendo il modello BaryonForge per includere osservabili a raggi X e profili di metallicità, gli autori dimostrano un metodo per auto-calibrare i parametri di disturbo astrofisici attraverso diversi probe, anziché marginalizzare su di essi indipendentemente.
3. Vincoli sulle PNG Basati su Simulazioni: Il lavoro va oltre le approssimazioni analitiche utilizzando simulazioni N-body per i vincoli sulle PNG, consentendo l'uso di scale non lineari e statistiche di ordine superiore (momenti) che sono sensibili ai picchi del campo di densità dove i segnali PNG sono più forti.
4. Rilascio del Codice Pubblico: Gli autori rilasciano il codebase Vaanam per la generazione di cieli sintetici e il codebase Aarambam per la costruzione delle condizioni iniziali delle PNG, facilitando future inferenze basate su simulazioni.

Risultati

• Rottura delle Degenerazioni: L'analisi congiunta di WL, tSZ e raggi X rompe significativamente le degenerazioni dei parametri. Nello specifico, la combinazione dei dati tSZ e raggi X vincola i parametri barionici (ad esempio, scala di espulsione del gas $\theta_{ej}$ e massa dell'alone $M_c$) molto più strettamente rispetto al solo WL. Questa auto-calibrazione dei parametri di disturbo previene il degrado dei vincoli sulle PNG che tipicamente si verifica marginalizzando sulle incertezze astrofisiche.
• Potere di Vincolo:
  • WL da Solo: Marginalizzando su tutti i parametri di disturbo (cosmologia, IA, barioni), i vincoli di WL su $f_{NL}$ degradano di un fattore da 2 a 3 rispetto a uno scenario con parametri fissi.
  • Sinergia Multi-Probe: L'inclusione dei campi tSZ e raggi X recupera il potere di vincolo perduto. Un'analisi congiunta di tutti e tre i campi, pur marginalizzando completamente sui parametri di disturbo, produce vincoli su $f_{NL}$ che sono circa 2 volte più stretti rispetto ai vincoli WL-only con parametri di disturbo fissi.
  • Dipendenza dalla Scala: I campi tSZ e raggi X sono particolarmente sensibili agli aloni massicci, fornendo informazioni complementari al WL. L'inclusione di statistiche non gaussiane (terzi momenti) fornisce un miglioramento di un fattore 2 rispetto ai soli secondi momenti.
• Specifiche del Modello: I vincoli variano per modello di PNG. Il WL è più competitivo per i template Equilatero e NBD2 (che imprimono caratteristiche sulle piccole scale), mentre il template Locale rimane impegnativo a causa della sua dominanza sulle grandi scale dove il WL è meno sensibile.

Significatività e Rivendicazioni
Il documento sostiene che le analisi multi-lunghezza d'onda siano essenziali per elevare i vincoli sulle PNG allo stesso status multi-probe goduto dalle analisi $\Lambda$CDM. Sfruttando la sinergia tra WL, tSZ e raggi X, gli autori dimostrano che è possibile vincolare simultaneamente la fisica dell'universo primordiale e marginalizzare sui complessi sistemi astrofisici. I risultati suggeriscono che i prossimi survey (LSST, Simons Observatory, eROSITA) possano ottenere vincoli sulle PNG comparabili o migliori dei limiti attuali della CMB, in particolare per i modelli con firme su piccola scala. Il lavoro stabilisce una tabella di marcia per utilizzare l'inferenza basata su simulazioni per estrarre la massima informazione dall'universo non lineare, andando oltre i limiti dei modelli analitici e delle analisi a singolo probe. Gli autori notano che le loro previsioni sono conservative, poiché non includono ancora tutte le scale disponibili o l'intero suite di quarti momenti, suggerendo un potenziale ancora maggiore per i futuri vincoli.