CoLoRe-2LPT: Lyman-$\alpha$ mock catalogues for the validation of DESI cosmological analyses

Dichiarazione del Problema
La foresta Lyman-$\alpha$ (Ly$\alpha$) è una sonda critica per lo studio della struttura su larga scala dell'universo ad alti redshift ($z > 2$), offrendo vincoli sulle Oscillazioni Acustiche dei Barioni (BAO) e sul clustering di forma completa (FS). Tuttavia, poiché il Dark Energy Spectroscopic Instrument (DESI) sta entrando in un'era di precisione senza precedenti con il suo Data Release 2 (DR2), la sensibilità alle assunzioni di modellazione, ai sistemi astrofisici e agli effetti osservativi è aumentata. I cataloghi mock realistici sono essenziali per validare le pipeline di analisi, quantificare i bias sistematici e verificare il recupero dei parametri. I precedenti metodi di generazione di mock per la foresta Ly$\alpha$ si basavano su realizzazioni log-normali del campo di densità. Sebbene efficaci per studi BAO su larga scala, questi approcci sono inaffidabili per scale moderatamente non lineari ($k \sim 0.1 - 0.5 \, h/\text{Mpc}$) e non riescono a riprodurre accuratamente l'allargamento non lineare del picco BAO o il clustering su piccola scala richiesto per analisi FS robuste. Inoltre, i precedenti mock spesso mancavano di una modellazione realistica dei quasar (QSO) e dei contaminanti astrofisici come i sistemi ad Alta Colonna di Densità (HCD) e le righe metalliche.

Metodologia
Gli autori presentano "CoLoRe-2LPT", una nuova suite di cataloghi mock rapidi della foresta Ly$\alpha$ cosmologica progettata per affrontare i limiti degli approcci log-normali. La metodologia prevede una pipeline multi-stadio:

1. Campi di Densità e Velocità 2LPT: L'innovazione principale è l'uso della Teoria delle Perturbazioni Lagrangiane del Secondo Ordine (2LPT) per evolvere il campo di densità della materia, sostituendo la trasformazione log-normale. La pipeline utilizza il codice CoLoRe per generare coni di luce cosmologici. Inizia con un campo casuale gaussiano a $z=0$, evolve le particelle di test tramite 2LPT per ottenere spostamenti e velocità, e le riscala per generare coni di luce che coprano l'intervallo di redshift dei QSO di DESI ($0 < z < 3.8$).
2. Calibrazione del Clustering dei QSO: Per popolare i coni di luce con i QSO, gli autori impiegano un modello di bias soglia dipendente dalla densità. I parametri di questo modello sono calibrati empiricamente contro simulazioni N-body ad alta risoluzione Abacus. Ciò garantisce che il bias lineare dei QSO e le distorsioni nello spazio di redshift (RSD) corrispondano ai dati osservativi. Il modello include l'aggiunta di effetti "Fingers-of-God" (FoG) tramite velocità casuali per riprodurre lo smorzamento osservato nelle simulazioni N-body.
3. Generazione della Trasmissione Ly$\alpha$: Il modulo LyaCoLoRe trasforma gli skewer di densità e velocità 2LPT in skewer di trasmissione Ly$\alpha$. Questo processo prevede:
   • L'aggiunta di potenza su piccola scala ai campi di densità 2LPT utilizzando un potere spettrale 1D per far corrispondere il clustering caratteristico della foresta fino a $\sim 100 \, \text{kpc}/h$.
   • La conversione della densità in profondità ottica utilizzando l'Approssimazione di Gunn-Peterson Fluttuante (FGPA).
   • L'applicazione delle RSD per mappare la profondità ottica nello spazio reale verso lo spazio di redshift.
   • Il calcolo del flusso trasmesso.
   • I parametri per la potenza su piccola scala, la profondità ottica e la scala di velocità sono sintonizzati per riprodurre le misurazioni DESI DR1 del potere spettrale 1D ($P_{1D}$), il flusso medio, il bias Ly$\alpha$ ($b_{\text{Ly}\alpha}$) e il parametro RSD ($\beta_{\text{Ly}\alpha}$).
4. Contaminanti Astrofisici e Rumore: Per creare mock simili a quelli di DESI, la pipeline incorpora:
   • HCD: Modellati popolando le celle sopra una soglia di densità, con parametri sintonizzati per corrispondere al clustering e ai tassi di incidenza degli HCD.
   • Righe Metalliche: Aggiunte dopo le RSD utilizzando lo script quickquasars, modellando le transizioni (es. Si II, Si III) che creano picchi di correlazione spurii.
   • Continuo e Rumore: I continui dei QSO sono generati usando SIMQSO, e il rumore strumentale è aggiunto tramite specsim per corrispondere all'impronta digitale e alle caratteristiche di esposizione di DESI DR2.

Contributi Chiave

• Prima Applicazione di 2LPT a Mock Completi di Cono di Luce Ly$\alpha$: Mentre i metodi LPT sono stati utilizzati per survey galattiche e mock Ly$\alpha$ basati su box, questo lavoro costituisce la prima applicazione di 2LPT a mock completi di cono di luce per una survey della foresta Ly$\alpha$.
• Modellazione Non Lineare Migliorata: L'approccio 2LPT cattura naturalmente le non-linearità moderate e l'allargamento non lineare del picco BAO, che nei precedenti mock log-normali venivano modellati in modo ad-hoc.
• Validazione Comprensiva: I mock sono validati rispetto a molteplici statistiche DESI, inclusi il potere spettrale 1D, il flusso medio, i bias Ly$\alpha$ e QSO, e il parametro $\beta$.
• Modellazione Realistica dei Contaminanti: L'inclusione di HCD e righe metalliche consente l'esplorazione dei sistemi indotti da questi contaminanti astrofisici, il che è cruciale per le analisi FS.

Risultati

• Statistiche di Clustering: I nuovi mock riproducono le statistiche Ly$\alpha$ entro il 10% delle ultime misurazioni DESI. Nello specifico, il bias Ly$\alpha$ e il parametro $\beta$, il flusso trasmesso medio e il potere spettrale 1D sono ben riprodotti.
• Clustering dei QSO: Il bias lineare dei QSO è recuperato entro il 5% dei valori osservativi. L'inclusione delle velocità 2LPT e degli effetti FoG migliora significativamente il recupero delle distorsioni nello spazio di redshift rispetto agli approcci precedenti.
• Validazione BAO e FS: Quando adattati con il modello BAO completo di DESI DR2, le correlazioni dei mock mostrano un eccellente accordo con i dati osservativi. I parametri cosmologici derivati dai mock deviano dall'unità al massimo di 0.003, uno scostamento trascurabile rispetto alle incertezze statistiche. I mock riproducono con successo lo smorzamento anisotropo del picco BAO naturalmente, senza allargamento ad-hoc.
• Validazione dei Contaminanti: L'inclusione di metalli e HCD risulta in correlazioni dei mock che corrispondono qualitativamente e quantitativamente ai picchi e alle distorsioni spurie osservati nei dati DESI DR2. I bias dei metalli concordano con i parametri DESI DR2 entro circa il 20%.

Significatività
L'articolo sostiene che questo nuovo insieme di mock sia uno strumento "chiave" per potenziare la modellazione e la validazione dell'analisi cosmologica full-shape della foresta Ly$\alpha$ di DESI DR2. Fornendo uno strumento computazionalmente efficiente e con motivazione fisica che cattura accuratamente sia il clustering su larga scala che quello moderatamente non lineare, i mock CoLoRe-2LPT consentono una validazione robusta delle analisi FS e BAO per le attuali e le prossime survey Ly$\alpha$. Il lavoro dimostra che i mock basati su 2LPT possono sostituire gli approcci log-normali per fornire caratteristiche su piccola scala e l'allargamento del picco BAO più realistici, elementi essenziali man mano che le analisi cosmologiche si spostano oltre le semplici misurazioni BAO verso i vincoli di forma completa (full-shape).