Optimal Transport Reconstruction of Biased Tracers in Primordial Non-Gaussian Fields

Sintesi Tecnica: Ricostruzione tramite Trasporto Ottimale di Traccianti Biasati in Campi Primordiali Non Gaussiani

Definizione del Problema
Il Trasporto Ottimale (OT) si è affermato come un framework efficiente per ricostruire il campo di spostamento cosmico, mappando le posizioni iniziali nelle posizioni attuali (o viceversa). Lavori precedenti hanno dimostrato che l'OT recupera accuratamente gli spostamenti per i traccianti biasati sotto l'assunzione di condizioni iniziali gaussiane. Tuttavia, una questione aperta critica rimane: questo framework rimane valido quando le condizioni iniziali sono non gaussiane? Questa è una sfida tempestiva, dato il crescente sforzo per vincolare la non-gaussianità primordiale (PNG) utilizzando i prossimi sondaggi cosmologici. Una difficoltà specifica sorge perché le implementazioni dell'OT spesso assumono una distribuzione spaziale iniziale uniforme. Sebbene ciò sia valido per l'intero campo di materia oscura, non è il caso dei sottoinsiemi biasati (ad esempio, gli aloni). Negli scenari gaussiani, la massa non osservata ("polvere") può essere modellata efficacemente; tuttavia, negli scenari non gaussiani, la polvere stessa acquisisce specifici bias dipendenti dalla scala che devono essere presi in considerazione per garantire una ricostruzione accurata.

Metodologia
Gli autori investigano l'applicabilità della ricostruzione OT a campi non gaussiani primordiali "locali", dove il potenziale gravitazionale $\phi_{NG}$ è modellato come una serie di potenze di un campo gaussiano $\phi_G$ (troncata al secondo ordine, $f_{NL} \neq 0$). In questo modello, il contrasto di densità $\delta_{NG}$ contiene un termine proporzionale al prodotto di campi gaussiani, portando a un distintivo bias dipendente dalla scala nel clustering dei traccianti biasati alle grandi scale (piccolo $k$).

Lo studio utilizza il set di simulazioni Quijote-PNG, composto da 20 box periodici ($L = 1 h^{-1}$ Gpc) con $512^3$ particelle, evoluti da $z=127$ a $z=0$ sotto una cosmologia $\Lambda$CDM piatta. I traccianti biasati sono identificati come aloni con massa $> 20 m_p$.

L'innovazione metodologica principale consiste nel modellare la "polvere" non osservata (la frazione di massa $1-f_b$ non presente nel catalogo degli aloni) per soddisfare il vincolo che l'intero campo di materia debba esibire un bias indipendente dalla scala. Gli autori derivano che, se i traccianti biasati esibiscono un bias dipendente dalla scala caratterizzato dai parametri $b_\delta$ e $b_\phi$ (dove il bias scala come $b(k) \approx b_\delta + 2f_{NL}b_\phi \alpha(k)$), la polvere deve possedere un bias dipendente dalla scala compensativo con segno opposto.

Due approcci di ricostruzione vengono testati:

1. OT-K: Utilizza un modello di filtro di Wiener per stimare la distribuzione della polvere basandosi sui traccianti osservati. Ciò richiede di ipotizzare la combinazione $2f_{NL}b_\phi$ (che può essere vincolata dai dati dei sondaggi) e la frazione di massa $f_b$.
2. OT-W: Utilizza le vere particelle di polvere della simulazione (un campione casuale) combinata con gli aloni pesati in massa per eseguire la ricostruzione, servendo come limite superiore delle prestazioni.

La qualità della ricostruzione viene valutata confrontando le posizioni ricostruite con le effettive posizioni dei protoaloni e analizzando la dipendenza dalla scala della divergenza del campo di spostamento ($\theta_{OT}$), che funge da proxy per il campo di densità iniziale.

Contributi Chiave e Risultati

• Caratterizzazione della Polvere: Il documento dimostra che nei modelli di PNG locale, la "polvere" non è un background passivo ma esibisce il proprio caratteristico bias dipendente dalla scala. Nello specifico, l'ampiezza del bias della polvere è determinata dal requisito che la somma dei campi dei traccianti biasati e della polvere produca un campo totale indipendente dalla scala.
• Estensione dell'OT alla Non-Gaussianità: Gli autori mostrano che la metodologia OT si estende con successo alle condizioni iniziali non gaussiane. Incorporando un modello realistico per la polvere che tiene conto del bias dipendente da $f_{NL}$, la ricostruzione dei traccianti biasati rimane accurata.
• Accuratezza della Ricostruzione:
  • Posizioni dei Traccianti: La ricostruzione OT-K (utilizzando il modello di polvere del filtro di Wiener) recupera correttamente il bias dipendente dalla scala dei protoaloni. La ricostruzione OT-W (usando la vera polvere) mostra un accordo ancora migliore, indicando che il metodo è robusto ma sensibile alla fedeltità del modello di polvere.
  • Proxy del Campo Completo: La divergenza del campo di spostamento OT ($\theta_{OT}$) mostra quasi nessuna dipendenza dalla scala a piccolo $k$, rispecchiando il comportamento del vero campo iniziale. Il modello OT-W estende questa indipendenza dalla scala a valori di $k$ più grandi rispetto a OT-K, ottenendo un miglioramento di un fattore 2 nell'intervallo di $k$ dove $b(k) \approx 1$.
• Robustezza rispetto a $f_{NL}$: Il valore medio del quadrato del campo ricostruito $\langle \theta_{OT}^2 \rangle$ rimane stabile attraverso diversi valori di $f_{NL}$ ($+100, 0, -100$), con variazioni entro la deviazione standard, suggerendo che il metodo non è influenzato (biased) dalla presenza di non-gaussianità in sé.

Significatività e Rivendicazioni
Il documento afferma che l'Optimal Transport fornisce un framework utile per ricostruire le strutture cosmiche anche quando le condizioni iniziali sono non gaussiane, a condizione che il modello per la "polvere" non osservata sia sufficientemente realistico. Gli autori sottolineano che il successo del metodo risiede nel fatto che la distintiva firma dipendente dalla scala della PNG è presente sia nei traccianti biasati che nella polvere, ma con segni opposti, permettendo al campo completo di rimanere non biasato.

Gli autori dichiarano con modestia di non aver dimostrato un miglioramento diretto dei vincoli su $f_{NL}$ da questo lavoro da solo. Invece, forniscono un framework fondamentale che permette alle analisi basate sulla ricostruzione di essere applicate coerentemente nelle cosmologie PNG. Suggeriscono che questo lavoro è un passo necessario verso l'integrazione delle tecniche di ricostruzione nelle future analisi ad alta precisione della PNG, potenzialmente aiutando a connettere più strettamente i traccianti biasati osservati ai loro controparti lagrangiani primordiali e mitigando sottili contributi non lineari o stocastici rilevanti per i prossimi sondaggi. Il documento conclude che un migliore modellazione della polvere (ad esempio, tramite metodi iterativi o machine learning) potrebbe produrre ulteriori guadagni.