From cosmological simulations to binary black hole mergers: The impact of using analytical star formation history models on gravitational-wave source populations

Sintesi Tecnica: Dalle Simulazioni Cosmologiche alle Fusioni di Buchi Neri Binari

Definizione del Problema
Le osservazioni delle fusioni di buchi neri binari (BBH) forniscono approfondimenti critici sull'evoluzione delle stelle massicce attraverso i tempi cosmici. Tuttavia, connettere le proprietà osservate delle fusioni ai loro astri progenitori richiede modelli accurati della storia della formazione stellare cosmica (SFRH) e dell'arricchimento chimico. Attualmente, il campo si affida frequentemente a fit analitici per la densità del tasso di formazione stellare dipendente dalla metallicità, $S(Z, z)$, come proxy per storie complesse sottostanti derivate da simulazioni cosmologiche o osservazioni. Queste prescrizioni analitiche sono spesso non vincolate e il loro impatto sulle popolazioni di BBH inferite rimane insufficientemente compreso. Inoltre, gli effetti della risoluzione e della dimensione della box delle simulazioni cosmologiche sulla derivazione di $S(Z, z)$ e sui successivi tassi di fusione BBH non sono stati quantificati sistematicamente.

Metodologia
Gli autori utilizzano la suite di simulazioni cosmologiche IllustrisTNG (specificamente TNG50-1, TNG100-1, TNG300-1), che variano in risoluzione e dimensione della box, per derivare direttamente dai dati delle celle di gas il valore $S(Z, z)$ basato sulla simulazione. Confrontano questi risultati diretti della simulazione con un modello di fit analitico ampiamente utilizzato (basato su L. A. C. van Son et al. 2023), che combina una forma funzionale di Madau & Dickinson (2014) per la densità totale del SFR, $S(z)$, con una distribuzione di metallicità log-normale asimmetrica.

Per valutare l'impatto di questi diversi input di $S(Z, z)$ sulle popolazioni BBH, gli autori accoppiano sia i modelli $S(Z, z)$ derivati dalla simulazione che quelli analitici al codice di sintesi di popolazione binaria rapida COMPAS (Compact Object Mergers: Population Astrophysics and Statistics). Utilizzano un framework di "integrazione cosmica" per calcolare i tassi di fusione e le distribuzioni di massa attraverso il redshift. Lo studio investiga specificamente:

1. La fedeltà del fit analitico nel riprodurre la forma e la normalizzazione di $S(Z, z)$ dalle diverse risoluzioni TNG.
2. Le conseguenti deviazioni nei tassi di fusione BBH e nelle distribuzioni di massa primaria ($M_{BH,1}$) utilizzando il fit analitico rispetto ai dati completi della simulazione.
3. L'influenza della risoluzione della simulazione (che spazia dall'alta risoluzione di TNG50-1 alla bassa risoluzione di TNG300-1) su tali risultati.

Risultati Chiave
Lo studio rileva che le prescrizioni analitiche per $S(Z, z)$ introducono errori sistematici significativi quando si modellano le popolazioni BBH:

• Discrepanze in $S(Z, z)$: Il fit analitico non riesce a catturare caratteristiche strutturali chiave della distribuzione di metallicità basata sulla simulazione. Nello specifico, manca di un "bump" secondario nel regime ad alta metallicità ($Z \gtrsim 3 Z_\odot$) e di una coda a bassa metallicità appiattita presenti nelle simulazioni. Il fit tende a sovrastimare il tasso di formazione stellare ad alti redshift ($z \gtrsim 6$) e alte metallicità, mentre sottostima i contributi a basse metallicità e redshift molto elevati, a seconda della risoluzione della simulazione.
• Impatto sui Tassi di Fusione: Ad alti redshift ($z \gtrsim 6$), il fit analitico sovrastima i tassi di fusione BBH con fattori che vanno da 10 a $10^4$ rispetto ai tassi basati sulla simulazione, con l'entità dell'errore che aumenta al diminuire della risoluzione della simulazione. A basso redshift ($z \approx 0$), i tassi di fusione locali derivati dal fit concordano generalmente con i tassi basati sulla simulazione entro un fattore di $\sim 1.5$, sebbene TNG50-1 mostri una sovrastima di $\sim 15$ Gpc$^{-3}$ yr$^{-1}$.
• Impatto sulle Distribuzioni di Massa: L'uso di fit analitici introduce caratteristiche spurie nella distribuzione di massa primaria BBH. In particolare, a $z \lesssim 2$, il fit produce un eccesso artificiale (bump) vicino a $8 M_\odot$ nella distribuzione di massa primaria, che è assente quando si utilizza il pieno $S(Z, z)$ della simulazione. Questa caratteristica è attribuita al canale di trasferimento di massa stabile ed è causata dall'incapacità del fit di modellare correttamente il picco di formazione stellare ad alta metallicità. Al contrario, il fit sopprime le caratteristiche ad alta massa (vicino a $33 M_\odot$ e $45 M_\odot$) per fattori di 3–5 a $M_{BH,1} \gtrsim 35 M_\odot$.
• Dipendenza dalla Risoluzione: La risoluzione della simulazione influenza significativamente la normalizzazione e la posizione del picco di $S(z)$. Simulazioni a risoluzione più alta (TNG50-1) producono tassi di formazione stellare più elevati e hanno un picco a redshift leggermente inferiori ($z \approx 2$) rispetto a quelle a risoluzione inferiore ($z \approx 2.5$). L'accuratezza del fit analitico degrada con la minore risoluzione della simulazione, particolarmente nel riprodurre la pendenza ad alto redshift.

Significato e Rivendicazioni
Gli autori concludono che, sebbene i fit analitici offrano una comodità computazionale, essi sono proxy inaffidabili per la complessa $S(Z, z)$ sottostante derivata dalle simulazioni cosmologiche, specialmente quando si interpretano le osservazioni di onde gravitazionali. Lo studio dimostra che i modelli analitici semplici non riescono a catturare caratteristiche fisiche critiche — come i bump ad alta metallicità e specifiche code a bassa metallicità — che influenzano direttamente i tassi di fusione BBH e le distribuzioni di massa.

L'articolo sottolinea che le discrepanze derivano dal fatto che le forme funzionali degli attuali fit analitici non possono accomodare la natura multimodale delle distribuzioni di metallicità trovate nelle simulazioni. Di conseguenza, fare affidamento su questi fit può portare a inferenze errate sui canali di formazione e sulla storia evolutiva dei BBH. Gli autori sostengono la necessità di una cauta integrazione di simulazioni cosmologiche, fit analitici e sintesi di popolazione per evitare l'introduzione di caratteristiche spurie nell'interpretazione dei dati di onde gravitazionali. Essi non affermano che i fit analitici siano interamente inutili, ma ne evidenziano i limiti e la necessità di validarli contro dati di simulazione ad alta fedeltà, specialmente man mano che i detector di prossima generazione esploreranno redshift più elevati.