Assembly bias and local Primordial non-Gaussianity from DESI DR1 Quasars

Titolo e Contesto

Il lavoro si concentra sull'analisi del clustering su larga scala dei quasar osservati dal Dark Energy Spectroscopic Instrument (DESI) per vincolare la non-gaussianità primordiale locale (PNG), parametrizzata dal parametro $f_{\text{NL}}$. Il documento affronta un problema fondamentale nell'interpretazione dei dati di DESI: la degenerazione tra $f_{\text{NL}}$ e il parametro di bias non-gaussianità $b_\phi$, che dipende a sua volta dal parametro di "assembly bias" $p$.

Il Problema Scientifico

La PNG locale induce un bias dipendente dalla scala nel clustering degli aloni di materia oscura e delle galassie. La grandezza misurabile è il prodotto $f_{\text{NL}} b_\phi$.
Il parametro $b_\phi$ è legato al bias lineare $b_1$ e al parametro $p$ dalla relazione:
$$b_\phi = 2\delta_c(b_1 - p)$$
dove $\delta_c \approx 1.686$ è la soglia di collasso sferico.

• Relazione di universalità: Se il campione di traccianti è selezionato solo in base alla massa dell'alone, si assume $p=1$.
• Bias di assemblaggio: Se la selezione dipende da altre proprietà (come la storia di fusione degli aloni), $p$ può deviare da 1. Per gli aloni che hanno subito fusioni recenti, studi precedenti suggerivano $p \approx 1.6$.
• La situazione DESI: Nell'analisi dei dati della prima release (DR1) di DESI, a causa dell'incertezza sul valore di $p$, sono stati forniti due vincoli separati su $f_{\text{NL}}$ (uno per $p=1$ e uno per $p=1.6$), rendendo i risultati parzialmente degeneri e meno informativi.

Metodologia

Gli autori hanno sviluppato un approccio basato su simulazioni idrodinamiche di stato dell'arte per costruire un prior fisico e robusto per il parametro $p$, specifico per i quasar selezionati da DESI.

1. Selezione del Campione in IllustrisTNG:

   • È stata utilizzata la simulazione idrodinamica TNG300-1 della suite IllustrisTNG.
   • È stato selezionato un campione di quasar che replica i criteri di selezione di DESI (densità numerica e massa dell'alone media a $z \approx 1.5$).
   • La selezione dei quasar è stata basata sul rapporto di Eddington ($\chi$), selezionando gli oggetti con accrescimento più efficiente (modalità "quasar"), che sono associati a fusioni recenti.
   • È stato confrontato con un campione di controllo (selezione casuale di aloni di massa simile) per isolare l'effetto della storia di fusione.

2. Calcolo del Parametro $p$:

   • Utilizzando il formalismo Extended Press-Schechter (ePS), gli autori hanno calcolato il parametro $p$ per ogni alone ospite, basandosi sulla sua storia di accrescimento di massa e sulla redshift di formazione ($z_f$).
   • La relazione analitica collega $p$ alla variabile $\omega_f$, che dipende dalla differenza tra la redshift di osservazione e quella di formazione.

3. Validazione con CAMELS:

   • Per quantificare la dipendenza dalle assunzioni astrofisiche e dalla varianza cosmica, è stata utilizzata la suite di simulazioni CAMELS (che include variazioni dei parametri cosmologici e della fisica sub-griglia, in particolare i feedback degli AGN).
   • Sono stati analizzati 26 volumi con condizioni iniziali diverse (CV set) e 36 realizzazioni con parametri astrofisici variati (1P set).

4. Applicazione ai Dati DESI DR1:

   • Il prior empirico $\Pi(p)$ ottenuto dalle simulazioni è stato propagato nell'analisi dei dati reali di DESI DR1.
   • È stata eseguita una marginalizzazione sulla distribuzione di $p$ per ottenere un singolo vincolo su $f_{\text{NL}}$, rompendo la degenerazione con $b_\phi$.

Risultati Chiave

• Distribuzione di $p$:

  • Il campione di quasar DESI-like in IllustrisTNG mostra una distribuzione di $p$ fortemente asimmetrica verso valori maggiori di 1, con una media $\bar{p} \simeq 1.4$ (specificamente 1.41).
  • Questo conferma che i quasar risiedono in aloni che hanno subito fusioni recenti, ma il valore medio è intermedio tra l'universalità ($p=1$) e il caso estremo di fusione recente ($p=1.6$).
  • La distribuzione è stabile rispetto al rumore di selezione: anche introducendo una variabilità del 34% nella selezione degli oggetti, il valore medio di $p$ rimane invariato.

• Robustezza:

  • L'analisi con CAMELS dimostra che il risultato è robusto rispetto alla varianza cosmica e alle incertezze nella fisica degli AGN (feedback, crescita dei buchi neri).
  • Variando i parametri astrofisici, la media di $p$ si sposta leggermente a $\bar{p} \simeq 1.38$, confermando la solidità del prior.
  • La dipendenza dalla redshift è debole nell'intervallo rilevante per DR1.

• Nuovi Vincoli su $f_{\text{NL}}$:

  • Applicando il prior $\Pi(p)$ ai dati DESI DR1, gli autori ottengono un vincolo aggiornato e unico:
    $$f_{\text{NL}} = -3.3 \pm 9.2$$
  • Questo risultato è confrontato con le analisi precedenti di DESI (che davano intervalli diversi per $p=1$ e $p=1.6$). L'uso del prior fisico porta a vincoli leggermente più stringenti e, soprattutto, fisicamente motivati.
  • È stato dimostrato che fissare $p$ al suo valore medio ($1.4$) durante l'inferenza produce risultati statisticamente indistinguibili dalla marginalizzazione completa sul prior, semplificando le procedure future.

Significato e Contributi

1. Risoluzione della Degenerazione: Il lavoro risolve il problema della degenerazione tra $f_{\text{NL}}$ e $b_\phi$ fornendo una stima fisica del parametro di bias di assemblaggio $p$ specifico per i quasar, invece di affidarsi a scelte arbitrarie ($p=1$ o $p=1.6$).
2. Validazione Astrofisica: Dimostra che la selezione dei quasar da parte di DESI è intrinsecamente legata alla storia di fusione degli aloni, validando l'uso di simulazioni idrodinamiche per calibrare i bias cosmologici.
3. Preparazione per il Futuro: Sebbene l'impatto sui dati DR1 sia moderato (a causa dell'attuale precisione statistica), questo approccio è cruciale per le future release di dati di DESI e per survey di prossima generazione (come Euclid o SKA), dove la precisione richiederà una modellazione accurata dei bias di assemblaggio.
4. Metodologia Riproducibile: Fornisce un framework robusto per derivare prior fisici per i parametri di bias non-gaussianità basati su simulazioni, riducendo l'incertezza sistematica nelle misure di cosmologia primordiale.

In conclusione, questo studio rappresenta un passo avanti significativo verso la misurazione precisa della non-gaussianità primordiale, trasformando un'incertezza sistematica (il valore di $p$) in un parametro vincolato dalla fisica della formazione delle strutture cosmiche.