FolpsD: combining EFT and phenomenological approaches for joint power spectrum and bispectrum analyses

1. Il Problema

L'estrazione di informazioni cosmologiche dai dati di clustering delle galassie (come quelli del Dark Energy Spectroscopic Instrument, DESI) richiede modelli teorici accurati su un'ampia gamma di scale. Le analisi "full-shape" (che utilizzano l'intera forma dello spettro di potenza, non solo le oscillazioni acustiche barioniche - BAO) sono limitate dall'evoluzione gravitazionale non lineare e dagli effetti di distorsione dello spazio redshift (RSD).

• Limiti dell'EFT standard: L'approccio basato sulla Teoria dei Campi Effettivi (EFT) dello spazio delle strutture su larga scala è robusto ma viene tipicamente limitato a scale quasi lineari ($k \lesssim 0.2 \, h\text{Mpc}^{-1}$) per evitare sistematiche di modellazione.
• Necessità di estensione: Per massimizzare il potere statistico di survey come DESI, è necessario estendere l'analisi a scale più piccole (più non lineari) e includere statistiche di ordine superiore (come il bispettro) per rompere le degenerazioni tra parametri di bias e cosmologici.
• Sfida: Come combinare in modo controllato l'estensione delle scale nello spettro di potenza (tramite modelli fenomenologici) con l'inclusione del bispettro, mantenendo stime dei parametri robuste e non biasate.

2. Metodologia

Gli autori presentano FolpsD, un'estensione del codice esistente Folps (Fast Optimized Large-scale structure Perturbation Solver). Il modello combina l'EFT a un loop con approcci fenomenologici e statistiche di ordine superiore.

• Fattore di Smorzamento (Line-of-Sight Damping):

  • Viene introdotto un fattore di smorzamento fenomenologico (di tipo Lorentziano) lungo la linea di vista (LoS) per modellare l'effetto "Finger-of-God" (FoG) causato dalle dispersioni di velocità su piccola scala.
  • Questo permette di modellare lo spettro di potenza a valori di $k$ più elevati rispetto all'EFT standard, estendendo il range di validità fino a $k \sim 0.3 \, h\text{Mpc}^{-1}$.
  • Il modello è ispirato al modello TNS (Taruya–Nishimichi–Saito) e VDG (Velocity-Difference Generator), ma mantiene la struttura completa dei correlatori a un loop dell'EFT, garantendo il rispetto delle simmetrie di conservazione della quantità di moto e dell'invarianza galileiana (a differenza di alcune implementazioni TNS pure che possono violare queste simmetrie su scale intermedie).

• Modellazione del Bispettro:

  • Viene implementato un modello di bispettro delle galassie a livello "tree-level" (albero).
  • Il bispettro è proiettato sulla base degli armonici sferici tripolari (base Sugiyama), che è particolarmente adatta per analizzare la dipendenza dall'orientamento rispetto alla linea di vista.
  • Viene anche considerata una versione del bispettro nella base Scoccimarro per confronto.
  • Il fattore di smorzamento viene applicato anche al bispettro, proiettato su ciascun vettore d'onda del triangolo di configurazione.

• Analisi dei Dati e Simulazioni:

  • Mock: Il modello è stato testato su 25 realizzazioni indipendenti di cataloghi mock di seconda generazione (AbacusSummit), simulando Luminous Red Galaxies (LRG) e Quasar (QSO) simili a quelli di DESI.
  • Dati Reali: La componente di spettro di potenza è stata applicata ai dati reali di DESI Data Release 1 (DR1).
  • Metodi Statistici: Utilizzo di catene di Markov Monte Carlo (MCMC) con il sampler Cobaya. Sono state analizzate diverse combinazioni: solo spettro di potenza (EFT vs FolpsD) e analisi congiunta (Spettro + Bispettro).

3. Contributi Chiave

1. Sviluppo di FolpsD: Un codice pubblico che unisce l'EFT a un loop, un fattore di smorzamento fenomenologico LoS e un bispettro tree-level nella base Sugiyama.
2. Valutazione della Complementarità: Dimostrazione che l'analisi congiunta di spettro di potenza e bispettro riduce significativamente le degenerazioni tra parametri di bias di secondo ordine ($b_2, b_s$) e parametri cosmologici.
3. Analisi dei Limiti dello Smorzamento: Identificazione critica che, sebbene lo smorzamento migliori i vincoli per tracciatori ad alto rapporto segnale/rumore (come le LRG), può introdurre bias nei parametri cosmologici per tracciatori dominati dal rumore (come i QSO) o in scenari fisici dove la forma dello spettro è degenerata con lo smorzamento (es. massa dei neutrini).
4. Estensione delle Scale: Dimostrazione che, grazie alla correlazione incrociata tra monopolo dello spettro di potenza e monopolo del bispettro, il bispettro può essere utilizzato in modo affidabile fino a scale più piccole ($k \sim 0.24 \, h\text{Mpc}^{-1}$) senza introdurre bias significativi, superando i tagli conservativi precedenti.

4. Risultati Principali

• Su Mock LRG (Galassie Rosse Luminose):

  • L'analisi congiunta FolpsD-Pk+Bk (Spettro + Bispettro con smorzamento) fornisce vincoli fino al 30% più stretti su $A_s$ (ampiezza delle fluttuazioni primordiali) e $\omega_{cdm}$ (densità di materia oscura fredda) rispetto all'EFT standard con solo spettro di potenza.
  • L'estensione del bispettro fino a $k_{max} = 0.24 \, h\text{Mpc}^{-1}$ (invece di 0.12) aumenta ulteriormente il "Figure of Merit" (FoM) del 14% e riduce il "Figure of Bias" (FoB) del 30%.
  • I parametri di bias di secondo ordine sono vincolati molto meglio, il che si traduce in vincoli cosmologici più precisi.

• Su Mock QSO (Quasar):

  • Per tracciatori con basso rapporto segnale/rumore, il fattore di smorzamento è poco vincolato e tende ad assorbire le fluttuazioni di rumore, portando a bias significativi nei parametri cosmologici (es. $\omega_{cdm}$ spostato di oltre $1\sigma$).
  • Questo suggerisce cautela nell'uso di modelli di smorzamento per tracciatori rumorosi.

• Su Dati DESI DR1:

  • Applicando il modello di spettro di potenza FolpsD ai dati reali DR1 (con tagli conservativi $k < 0.3 \, h\text{Mpc}^{-1}$), si ottiene un miglioramento dei vincoli fino al 20% su $\sigma_8$ e $\Omega_m$ rispetto all'analisi EFT standard.
  • Modello $w_0 w_a \text{CDM}$: Utilizzando il metodo "ShapeFit" su FolpsD, si ottengono vincoli del 15% su $w_0$ e del 21% su $w_a$ più stretti rispetto all'EFT.
  • Risultato Sorprendente: L'analisi full-shape su dati DESI DR1 con FolpsD mostra una deviazione dall'energia oscura costante ($w = -1$) a un livello leggermente superiore a $1\sigma$, un risultato preliminare che merita ulteriori conferme.

• Massa dei Neutrini:

  • L'analisi mostra che l'uso di fattori di smorzamento può degradare i vincoli sulla massa dei neutrini a causa della degenerazione tra lo smorzamento fenomenologico e la soppressione dello spettro di potenza indotta dai neutrini. L'inclusione del bispettro aiuta a mitigare parzialmente questo problema.

5. Significato e Implicazioni

Il lavoro di FolpsD rappresenta un passo avanti significativo nell'analisi dei dati di survey di galassie di prossima generazione come DESI.

• Ottimizzazione delle Risorse: Dimostra che combinare statistiche di ordine superiore (bispettro) con modelli fenomenologici controllati (smorzamento) permette di estrarre più informazioni cosmologiche senza dover ricorrere a simulazioni N-body estremamente costose per ogni punto di dati.
• Robustezza: Fornisce una mappa chiara dei limiti di questi approcci: funzionano bene per tracciatori ad alto S/N (LRG) e per modelli cosmologici standard o di energia oscura, ma richiedono cautela per tracciatori rumorosi (QSO) o per fisica che introduce nuove scale (neutrini, gravità modificata).
• Impatto Futuro: I risultati suggeriscono che le analisi future su DESI DR2 e oltre dovrebbero adottare strategie congiunte di spettro e bispettro per massimizzare la precisione sui parametri cosmologici, potenzialmente portando a scoperte sulla natura dell'energia oscura o sulla fisica dei neutrini.

Il codice è reso pubblico per facilitare l'adozione da parte della comunità scientifica.