The Impact of Galaxy Formation on Galaxy Biasing, and Implications for Primordial non-Gaussianity Constraints

Sintesi Tecnica: L'impatto della formazione delle galassie sul biasing galattico e le implicazioni per i vincoli sulla non-gaussianità primordiale

Enunciato del problema
Il parametro $f_{\text{NL}}$ quantifica la non-gaussianità primordiale (PNG) locale nelle fluttuazioni di energia iniziali dell'Universo. Una rilevazione non nulla distinguerebbe tra modelli di inflazione a campo singolo e a più campi. Sebbene i prossimi sondaggi galattici mirino a vincolare $f_{\text{NL}}$ con errori statistici di $\sigma_{f_{\text{NL}}} \sim 1$, gli errori sistematici dominano attualmente questi vincoli. Nello specifico, la misurazione di $f_{\text{NL}}$ tramite l'clustering galattico è degenerata con i parametri di bias galattico $b_1$ (risposta alla densità di massa) e $b_\phi$ (risposta alle perturbazioni del potenziale primordiale).

Le attuali stime si basano spesso sulla "relazione di universalità", che assume una scala fissa tra $b_\phi$ e $b_1$ (ad esempio, $b_\phi = 2\delta_c(b_1 - 1)$) derivata dagli aloni di materia oscura. Tuttavia, lavori recenti (ad esempio, Barreira et al. 2020) indicano che questa relazione fallisce per le galassie, poiché $b_\phi$ dipende dalla selezione delle galassie e dalla specifica fisica della formazione galattica. La complessità della modellizzazione della formazione delle galassie introduce un'incertezza significativa in $b_\phi$, minacciando la precisione dei futuri vincoli su $f_{\text{NL}}$.

Metodologia
Per affrontare l'incertezza in $b_\phi$ derivante dalla modellizzazione della formazione delle galassie, gli autori utilizzano il Modello Semi-Analitico Santa Cruz (SC-SAM) integrato in simulazioni N-body a universi separati (separate-universe).

• Simulazioni: Sono state eseguite tre simulazioni N-body AREPO con dimensione della box $L = 205 h^{-1}$ cMpc e $N = 1280^3$ particelle. Queste simulazioni condividono identiche condizioni iniziali gaussiane, ma differiscono per l'ampiezza dello spettro di potenza primordiale ($\sigma_8$), utilizzando la tecnica dell'universo separato per misurare efficientemente la risposta della formazione galattica alle perturbazioni a lunga lunghezza d'onda.
• Variazioni del modello: Gli autori hanno eseguito 56 istanze uniche di SC-SAM. Queste variazioni comportano la modifica di tre parametri di feedback chiave:
  1. $\epsilon_{\text{SN}}$ (ampiezza dell'uscita di massa guidata dalle stelle, variata tramite il fattore $A_{\text{SN1}}$).
  2. $\alpha_{\text{rh}}$ (esponente della relazione di uscita di massa, variato tramite il fattore $A_{\text{SN2}}$).
  3. $\kappa_{\text{radio}}$ (ampiezza del feedback radio-mode dell'AGN, variata tramite il fattore $A_{\text{AGN}}$).
• Misurazioni:
  • $b_\phi$ è calcolato tramite differenziazione finita del conteggio numerico delle galassie tra simulazioni ad alta-$\sigma_8$, fiduciali e a bassa-$\sigma_8$.
  • $b_1$ è derivato dallo spettro di cross-potenza galassia-materia.
  • Le galassie sono state selezionate in base alla Massa Stellare ($M_*$), al Tasso di Formazione Stellare (SFR) e allo sSFR (sSFR specifico) a redshift $z=0, 1, 2$.
  • Un metric di "realismo" ($\chi^2$) è stato assegnato a ciascun modello in base al suo accordo con le funzioni di massa stellare osservate a $z=0$ (Bernardi et al. 2018).
  • Lo studio ha inoltre valutato l'impatto del bias di assemblaggio galattico mescolando le galassie in aloni di massa simile ma con storie di formazione differenti.

Contributi Chiave

1. Esplorazione sistematica del feedback: A differenza di studi precedenti che spesso utilizzavano una singola simulazione idrodinamica o un SAM fisso, questo lavoro varia sistematicamente i parametri di feedback stellare e AGN per mappare il risultante panorama di $b_\phi$ e $b_1$.
2. Applicazione del SAM a universi separati: Gli autori applicano la tecnica dell'universo separato a un modello semi-analitico, permettendo la misurazione efficiente di $b_\phi$ attraverso un ampio spazio di parametri che sarebbe computazionalmente proibitivo con simulazioni idrodinamiche complete.
3. Analisi della dipendenza dalla selezione: Il documento confronta esplicitamente come diversi criteri di selezione galattica ($M_*$, SFR, sSFR) rispondano alle variazioni della fisica della formazione galattica, identificando quali selezioni siano più robuste alle incertezze di modellizzazione.
4. Validazione tra modelli: I risultati sono confrontati con i precedenti risultati della simulazione idrodinamica IllustrisTNG, del Munich SAM e del Galacticus SAM per determinare se lo spazio dei parametri SC-SAM comprenda il comportamento di altri modelli di formazione galattica.

Risultati

• Deviazione dall'universalità: I modelli SC-SAM confermano che la relazione di universalità ($b_\phi = 2\delta_c(b_1 - 1)$) non è accurata per le galassie. La relazione $b_\phi(b_1)$ varia significativamente a seconda della selezione delle galassie e dei specifici parametri di feedback.
• Selezioni di Massa Stellare: $b_\phi$ aumenta generalmente con la massa stellare. Tuttavia, la relazione $b_\phi(b_1)$ è più ripida di quella trovata in IllustrisTNG (Barreira et al. 2020), particolarmente per le selezioni di bassa massa. La relazione è sensibile ai parametri di feedback stellare ($A_{\text{SN1}}, A_{\text{SN2}}$).
• Selezioni SFR: I valori di $b_\phi$ per le galassie selezionate tramite SFR sono caotici e spesso vicini allo zero. La relazione è altamente sensibile ai parametri di feedback perché diversi parametri spostano la relazione SFR–massa dell'alone in modi complessi.
• Selezioni sSFR (Risultato Chiave): Le selezioni basate sul tasso di formazione stellare specifico (sSFR) sono notevolmente robuste alle variazioni della modellizzazione della formazione galattica. In tutti i 56 modelli, le relazioni $b_\phi$ vs sSFR e $b_\phi(b_1)$ seguono un percorso stretto e consistente. Ciò è attribuito alla natura autoregolata dello SC-SAM, che mantiene una relazione sSFR–massa dell'alone coerente indipendentemente dalle variazioni dei parametri di feedback.
• Bias di Assemblaggio: Lo studio trova che il bias di assemblaggio galattico influenza le misurazioni di $b_1$ (riducendo il clustering per le selezioni di massa stellare, aumentandolo per le selezioni SFR/sSFR quando la storia di formazione viene rimescolata). Tuttavia, questi segnali di bias di assemblaggio hanno un impatto minimo sulle relazioni $b_\phi(b_1)$ misurate.
• Confronto con altri modelli: Le variazioni SC-SAM comprendono i comportamenti di $b_\phi$ misurati nel Munich SAM e nel Galacticus. Sebbene vi siano tensioni con i risultati di IllustrisTNG (probabilmente dovute ad effetti barionici assenti nei SAM o a differenze nella costruzione degli alberi di fusione), le variazioni SC-SAM coprono un ampio intervallo di valori di $b_\phi$ coerenti con altri modelli.

Significato e Rivendicazioni
Il documento afferma che le incertezze nella modellizzazione della formazione delle galassie rappresentano un sistema critico per i futuri vincoli su $f_{\text{NL}}$. Dimostrando che $b_\phi$ non è universale e dipende fortemente dalla selezione e dalla modellizzazione, gli autori sostengono che assumere una relazione di universalità fissa rischi di produrre vincoli distorti.

Il significato primario di questo lavoro risiede nell'identificare l'sSFR come una selezione galattica robusta per gli studi su $f_{\text{NL}}$. Poiché le selezioni sSFR producono valori di $b_\phi$ consistenti attraverso una vasta gamma di fisiche di formazione galattica, esse offrono una via per ridurre gli errori sistematici nei futuri sondaggi. Gli autori concludono che, sebbene nessun singolo modello possa descrivere perfettamente l'universo, la selezione delle galassie tramite sSFR minimizza la sensibilità dei vincoli su $f_{\text{NL}}$ ai dettagli specifici della fisica della formazione galattica. Inoltre, il lavoro fornisce un quadro per costruire prior informativi per $b_\phi$ basati su dati osservativi, muovendosi verso vincoli fisicamente motivati sulla non-gaussianità primordiale.