Measuring local primordial non-Gaussianity from the clustering of DESI DR1 LRGs and QSOs

Sintesi Tecnica: Misurazione della Non-Gaussianità Primordiale Locale da DESI DR1 Clustering

Problematica
La natura dell'epoca inflazionistica nell'Universo primordiale rimane una delle principali questioni aperte in cosmologia. Sebbene il Fondo Cosmico a Microonde (CMB) abbia fornito i vincoli più sensibili sulla non-gaussianità primordiale (PNG) fino ad oggi, i sondaggi della Grande Struttura su Larga Scala (LSS) offrono una sonda complementare e potenzialmente competitiva. La PNG, parametrizzata dal parametro di non-linearità locale $f_{\text{NL}}$, induce un bias dipendente dalla scala nel clustering dei traccianti. Tuttavia, misurare $f_{\text{NL}}$ nello spazio di configurazione (la funzione di correlazione a due punti, o 2pcf) è notoriamente difficile a causa della contaminazione da sistematica di imaging. Variazioni nella profondità del sondaggio, l'estinzione da polvere e la densità stellare possono indurre fluttuazioni della densità angolare che mimano il segnale di clustering su larga scala della PNG. Questo articolo affronta la sfida di estrarre una misura robusta di $f_{\text{NL}}$ dai dati DESI DR1, mitigando rigorosamente questi effetti sistematici.

Metodologia
Gli autori impiegano un framework statistico basato su simulazioni per modellare il monopolo della 2pcf di Galassie Rosse Luminose (LRG) e Quasar (QSO) dal campione DESI DR1.

• Elaborazione dei Dati: L'analisi utilizza l'algoritmo ConKer per stimare la 2pcf. I traccianti sono mappati su una griglia cartesiana 3D con pesi ($w_{\text{tot}}$) che tengono conto delle sistematiche di imaging ($w_{\text{sys}}$), dell'incompletezza del catalogo ($w_{\text{comp}}$), dei fallimenti di redshift ($w_{\text{zfail}}$) e dei pesi FKP ($w_{\text{FKP}}$). L'analisi copre i Cap di Nord e Sud (NGC/SGC) per le LRG ($0.4 < z < 1.1$) e i QSO ($0.8 < z < 3.1$).
• Modello Teorico: Il modello è interamente basato su simulazioni, ancorato alle simulazioni N-body ad alta risoluzione AbacusSummit per il caso fiduciale ($f_{\text{NL}}=0$) e sulle simulazioni "FastPM" per la dipendenza dalla PNG. La 2pc è modellata come un'interpolazione lineare tra le realizzazioni con $f_{\text{NL}}=0$ e $f_{\text{NL}}=100$.
• Bias e Universalità: Il modello incorpora il bias lineare delle galassie ($\tilde{b}_1$) e il bias della PNG ($b_\phi$). La relazione tra questi è governata dalla relazione di Universalità: $b_\phi = 2\delta_c(b_1 - p)$, dove $p$ è un parametro dipendente dal tracciante. L'analisi impone prior gaussiani su $p$: $p_{\text{LRG}} = 1.0 \pm 0.1$ e $p_{\text{QSO}} = 1.6 \pm 0.1$.
• Mitigazione delle Sistematiche: Una componente critica della metodologia è la modellazione esplicita delle sistematiche di imaging. Gli autori introducono parametri di disturbo ($K^R_{\text{sys}}$) che rappresentano l'applicazione frazionaria eccessiva o insufficiente dei pesi di imaging in specifiche regioni del sondaggio (MzLS, (N)DEC, (S)DEC). Questi parametri sono vincolati tramite un approccio basato su simulazioni in cui l'effetto della variazione dei pesi sulla 2pcf è pre-calcolato e adattato a una funzione quadratica.
• Inferenza Statistica: Una funzione di verosimiglianza viene costruita utilizzando una statistica $\chi^2$ che include la matrice di covarianza derivata da 1000 mock EZ (Zel'dovich effettivi). I parametri di interesse ($f_{\text{NL}}$ e bias lineare $b_1$) e sette parametri di disturbo vengono campionati tramite Markov Chain Monte Carlo (MCMC). Le scale delle Oscillazioni Acustiche dei Barioni (BAO) ($100 < s < 120 \, h^{-1}\text{Mpc}$) sono mascherate per evitare bias di modellazione non lineare.

Risultati Chiave
L'analisi produce la prima misura di $f_{\text{NL}}$ nello spazio di configurazione 2pcf utilizzando i dati DESI:

• Campione LRG: La misura fornisce $f_{\text{NL}} = -3^{+22}_{-21}$ al livello di confidenza (CL) del 68%.
• Campione QSO: La misura fornisce $f_{\text{NL}} = 0^{+17}_{-16}$ al CL del 68%.
• Vincolo Congiunto: Combinando entrambi i campioni di traccianti, si ottiene un vincolo congiunto di $f_{\text{NL}} = -3^{+12}_{-12}$ al CL del 68%.
• Parametri di Bias: I bias lineari sono vincolati a $b_1 = 2.01 \pm 0.06$ per le LRG e $b_1 = 2.57 \pm 0.10$ per i QSO.
• Sistematica: I parametri di disturbo che governano le sistematiche di imaging ($K_{\text{sys}}$) e il parametro di formazione $p$ rimangono ampiamente coerenti con i loro prior, indicando che i dati non prediligono fortemente deviazioni dai modelli sistematici fiduciali.

Significato e Rivendicazioni
L'articolo rivendica tre conclusioni primarie:

1. Robustezza dell'Analisi nello Spazio di Configurazione: Il lavoro dimostra che uno studio robusto della PNG richiede un'analisi sia nello spazio di Fourier che in quello di configurazione. Sebbene gli effetti della PNG differiscano in questi domini, la coerenza di questi risultati nello spazio di configurazione con le precedenti analisi nello spazio di Fourier (Chaussidon et al. 2025) valida l'approccio.
2. Criticità della Mitigazione delle Sistematiche: Lo studio sottolinea che il controllo delle sistematiche di imaging è fondamentale per le misure di PNG in LSS. L'accordo tra questo risultato nello spazio di configurazione e i precedenti risultati nello spazio di Fourier (che utilizzano diverse strategie di mitigazione delle sistematiche) serve come una forte validazione della capacità della pipeline di gestire l'incompletezza del sondaggio e le variazioni di densità angolare.
3. Potere di Vincolo della LSS: La rivendicazione più significativa è che i sondaggi LSS si stanno rapidamente avvicinando alla sensibilità delle sonde CMB. Con solo il primo anno di dati DESI, il vincolo congiunto ($\sigma_{f_{\text{NL}}} \approx 12$) supera significativamente i precedenti vincoli LSS e si avvicina alla sensibilità della misura CMB di Planck ($\sigma_{f_{\text{NL}}} \approx 5$). Gli autori postulano che il campo sia entrato in un'era in cui la caratterizzazione della struttura cosmologica tardiva sarà una sonda principale della fisica primordiale.

I risultati sono coerenti con il modello standard di inflazione "slow-roll" ($f_{\text{NL}} \approx 0$) e non mostrano evidenza di scenari di inflazione a campo multiplo entro i limiti di precisione attuali.