Efficient evaluation of the dark-matter two-loop power spectrum in the EFT of LSS

1. Il Problema

Le moderne indagini cosmologiche sulla Struttura su Grande Scala (LSS), come quelle condotte da DESI ed Euclid, richiedono previsioni teoriche estremamente precise per estrarre informazioni cosmologiche dai dati. La Teoria Efficace della Struttura su Grande Scala (EFTofLSS) è lo strumento standard per modellare le fluttuazioni di densità nella regime non lineare debole.

Tuttavia, l'uso pratico dell'EFTofLSS si scontra con due ostacoli principali quando si tenta di andare oltre il calcolo a un loop:

1. Complessità Computazionale: Il calcolo del potere spettrale a due loop (2-loop) per la materia oscura richiede l'integrazione numerica di diagrammi con due momenti interni. Questi integrali sono computazionalmente costosi e lenti da valutare, rendendo impossibile l'esplorazione efficiente dello spazio dei parametri cosmologici necessaria per l'analisi dei dati ("full-shape analysis").
2. Singularità UV e IR: Gli integrali di loop presentano divergenze ultraviolette (UV) dovute alla fisica non perturbativa a piccole scale e singolarità infrarosse (IR) che, sebbene si cancellino nella somma totale dei diagrammi, causano instabilità numerica se trattate localmente nell'integrando.

2. Metodologia Proposta

Gli autori sviluppano un metodo ibrido che combina tecniche analitiche di rinormalizzazione con un'integrazione numerica altamente ottimizzata. La strategia si articola in quattro pilastri fondamentali:

A. Rinormalizzazione UV e Sottrazioni Locali

Per gestire le divergenze UV, gli autori non si limitano a regolarizzare gli integrali, ma separano l'integrando in una parte "regolare" e una parte "UV" che cattura il comportamento asintotico.

• Sottrazione UV: Sottraggono approssimazioni asintotiche degli integrandi (basate su espansioni di Taylor per momenti grandi) direttamente a livello dell'integrando.
• Controtermini EFT: Dimostrano che le parti sottratte sono esattamente degenerate con i contotermini dell'EFT (8 controtermini necessari a due loop, inclusi effetti di vorticità generata). Questo permette di assorbire le divergenze nei parametri fisici misurabili, lasciando integrali finiti e regolari.

B. Cancellazioni IR Locali

Le singolarità IR (quando i momenti di loop tendono a zero) si cancellano solo globalmente tra diagrammi diversi, rendendo l'integrazione stocastica instabile.

• Partizioni dell'Unità: Utilizzano una tecnica basata su "partizioni dell'unità" (invece di funzioni gradino) per riscrivere gli integrandi in modo che tutte le singolarità IR siano mappate nella stessa regione dello spazio dei momenti (l'origine).
• Risultato: Questo rende l'integrando localmente liscio in tutto il dominio di integrazione, permettendo una convergenza rapida degli algoritmi Monte Carlo.

C. Espansione Intorno a un Cosmologia di Riferimento

Per decouplare la dipendenza dai parametri cosmologici dal calcolo degli integrali:

• Si sceglie una cosmologia di riferimento (es. Planck) con uno spettro di potenza lineare $P^{[0]}_{lin}(k)$.
• Qualsiasi altro modello cosmologico $j$ viene espresso come una piccola differenza rispetto al riferimento: $\Delta P^{[j]}_{lin}(k) = P^{[j]}_{lin}(k) - N^{[j]} P^{[0]}_{lin}(k)$.
• Poiché le differenze tra modelli realistici sono piccole, si espandono le correzioni a loop in serie di potenze di $\Delta P^{[j]}_{lin}$.

D. Pre-calcolo degli Integrali Universali

Sfruttando l'espansione sopra descritta, gli integrali a loop vengono decomposti in una combinazione lineare di tensori di integrali indipendenti dalla cosmologia.

• Gli integrali universali (che dipendono solo dalla forma delle funzioni di base usate per approssimare $\Delta P_{lin}$) vengono pre-calcolati una sola volta con alta precisione.
• Per qualsiasi nuova cosmologia, la previsione si ottiene semplicemente contraendo questi tensori pre-calcolati con i coefficienti di espansione specifici del modello.

3. Contributi Chiave

• Prima valutazione numerica efficiente a due loop: Il lavoro presenta il primo metodo capace di valutare lo spettro di potenza della materia oscura a due loop nell'EFTofLSS in modo sufficientemente veloce da essere utilizzato in analisi di dati reali.
• Identificazione di 8 Controtermini: Viene dimostrato che sono necessari otto controtermini EFT (inclusi quelli che catturano la vorticità generata e la non-località temporale) per rinormalizzare completamente la teoria a due loop. Questo supera le stime precedenti che ne richiedevano solo 4 o meno.
• Libreria di Integrali Pre-calcolati: Gli autori rilasciano pubblicamente i tensori di integrali pre-calcolati, permettendo alla comunità di calcolare rapidamente le previsioni per qualsiasi cosmologia senza dover rieseguire costose integrazioni Monte Carlo.
• Gestione rigorosa delle singolarità: La combinazione di sottrazioni UV analitiche e partizioni dell'unità per le cancellazioni IR risolve i problemi di stabilità numerica che avevano finora limitato i calcoli a due loop.

4. Risultati

• Precisione: Il metodo raggiunge un'accuratezza teorica inferiore all'1% (spesso al livello di per-mille) rispetto all'incertezza sperimentale prevista per survey come DESI.
• Validazione: I risultati sono stati validati confrontando il metodo di approssimazione (basato su funzioni di base) con l'integrazione numerica diretta degli spettri lineari esatti per un modello cosmologico di test. Le differenze sono risultate trascurabili.
• Scalabilità: Il metodo permette di esplorare lo spazio dei parametri cosmologici in tempi ridotti, rendendo fattibile l'analisi "full-shape" a due loop.
• Convergenza: Viene mostrato che l'espansione in $\Delta P_{lin}$ converge rapidamente; i termini di ordine superiore (cubici) sono trascurabili per la maggior parte dei modelli cosmologici realistici.

5. Significato

Questo lavoro rappresenta un passo fondamentale per la cosmologia di precisione di prossima generazione.

1. Abilitazione dell'Analisi a Due Loop: Permette di sfruttare la precisione aggiuntiva offerta dai calcoli a due loop, estendendo il regime di validità delle previsioni teoriche verso scale di lunghezza più piccole (k più alti), dove la quantità di informazione è maggiore.
2. Ottimizzazione delle Risorse: Decouplando la dipendenza cosmologica dagli integrali, il metodo riduce drasticamente il costo computazionale per l'inferenza dei parametri, un requisito essenziale per missioni come Euclid e DESI.
3. Completezza Teorica: La corretta identificazione e inclusione di tutti gli otto controtermini necessari assicura che le previsioni siano teoricamente consistenti e prive di dipendenze arbitrarie dalla fisica UV non controllata.

In sintesi, il paper fornisce sia il formalismo teorico completo sia gli strumenti pratici (codice e dati) per implementare l'EFTofLSS a due loop, colmando il divario tra la teoria avanzata e l'analisi dei dati osservativi su larga scala.