AGN-driven BBH mergers: Black hole populations and hierarchical growth across the AGN parameter space

Sintesi Tecnica: Fusioni di BBH guidate da AGN: popolazioni di buchi neri e crescita gerarchica attraverso lo spazio dei parametri degli AGN

Problema e Motivazione
I dischi di nuclei galattici attivi (AGN) sono emersi come un canale dinamico promettente per la formazione e la fusione di sistemi binari di buchi neri (BBH). Sebbene i buchi neri (BH) di massa stellare immersi in questi dischi possano migrare, accoppiarsi e fondersi tramite coppie di torchi di gas ed emissione di onde gravitazionali (GW), le proprietà previste di questa popolazione rimangono altamente incerte. Studi precedenti hanno evidenziato una forte dipendenza da ingredienti fisici quali la struttura del disco, le coppie di migrazione, l'efficienza dell'indurimento del gas e la durata degli AGN, adottando spesso assunzioni differenti che ostacolano il confronto diretto. Vincoli osservativi recenti suggeriscono che il canale AGN non possa spiegare l'intera popolazione di fusioni BBH rilevate dagli osservatori di onde gravitazionali, ma il suo specifico contributo e le sue firme distintive richiedono un'esplorazione sistematica dello spazio dei parametri degli AGN per identificare caratteristiche osservative robuste e fonti di incertezza.

Metodologia
Gli autori presentano un aggiornato framework semi-analitico implementato all'interno dell'ambiente software fastcluster per modellare la formazione, l'evoluzione e la crescita gerarchica dei BBH nei dischi di AGN. L'approccio traccia il ciclo di vita dei buchi neri di origine stellare (prima generazione, 1g) attraverso la cattura, la migrazione radiale, la formazione di coppie binarie, l'indurimento guidato dal gas, le perturbazioni dinamiche e l'eventuale fusione.

Componenti metodologiche chiave includono:

• Modellazione del Disco: Il disco di accrezione è modellato utilizzando il modello analitico a stato stazionario di Sirko & Goodman (2003) (SG). Gli autori risolvono in modo autoconsistente le equazioni del disco monodimensionali per variabili di massa del buco nero supermassiccio ($M_\bullet$), rapporto di Eddington ($\lambda_\bullet$) e parametri di viscosità del disco ($\alpha$). Vengono testati due valori di viscosità: $\alpha = 0.01$ (bassa viscosità) e $\alpha = 0.1$ (alta viscosità).
• Condizioni Iniziali: Le simulazioni inizializzano le popolazioni di BH basandosi sulle relazioni di scala dei cluster stellari nucleari (NSC). I BH di prima generazione (1g) sono estratti da un catalogo di sintesi di popolazione (sevn) assumendo metallicità solare. Gli spin sono assegnati in base a una distribuzione Maxwelliana, con orientamenti allineati al momento angolare del disco.
• Fisica Evolutiva:
  • Migrazione: I BH migrano attraverso i regimi Type I e Type II, con prescrizioni derivate da Grishin et al. (2024) e Kanagawa et al. (2018).
  • Accoppiamento: La formazione di binari segue un criterio assistito dal gas, tenendo conto della frazione di orbite prograde ($f_{progr}$).
  • Indurimento ed Incontri: Dopo la formazione, i binari evolvono sotto l'effetto combinato di torchi di gas (interazione con il disco circumbinario) e incontri a tre corpi con singoli BH. Gli autori integrano le equazioni di evoluzione orbitale combinando l'indurimento guidato dal gas (usando Ishibashi & Gröbner 2024) e l'emissione di GW (Peters 1964). Gli incontri sono modellati utilizzando una griglia pre-calcolata di esperimenti di scattering post-newtoniani (tsunami).
  • Crescita Gerarchica: I resti di fusione vengono calcolati utilizzando prescrizioni di spin allineato e velocità di rinculo (kick). I resti che rimangono legati al disco partecipano a generazioni successive di fusioni ($N$-esima generazione), soggetti a limiti di durata dell'AGN e criteri di espulsione.
• Sintesi della Popolazione: Per costruire una popolazione intrinseca di BBH rappresentativa dell'Universo locale ($z \lesssim 0.3$), gli autori pesano gli output delle simulazioni secondo le distribuzioni osservate degli AGN di $M_\bullet$ e $\lambda_\bullet$ (Greene & Ho 2007).

Risultati Chiave
Lo studio esplora sistematicamente la dipendenza delle popolazioni di fusione BBH da $M_\bullet$, $\lambda_\bullet$ e $\alpha$, producendo i seguenti risultati:

1. Crescita Gerarchica e Distribuzione di Massa:

   • Gli dischi di AGN possono produrre fusioni ripetute, costruendo una coda ad alta massa che si estende oltre il gap della supernova a instabilità di coppia (PISN) ($50 \lesssim m_1/M_\odot \lesssim 130$) nell'intervallo dei buchi neri di massa intermedia (IMBH).
   • Dipendenza dalla Viscosità: La crescita gerarchica è significativamente più efficiente nei dischi a bassa viscosità ($\alpha = 0.01$), dove le densità superficiali del gas più elevate accorciano le scale temporali dinamiche, permettendo un numero maggiore di generazioni di fusione ($N_{gen} > 100$). I dischi ad alta viscosità ($\alpha = 0.1$) sopprimono la formazione di resti massicci, saturando tipicamente a generazioni inferiori ($N_{gen} \approx 20$).
   • Rapporti di Massa: La popolazione evolve verso rapporti di massa sempre più diseguali ad alte masse primarie. Mentre le masse primarie crescono attraverso fusioni successive, le masse secondarie si saturano a valori inferiori (tipicamente sotto il gap PISN), poiché le primarie di alta generazione si fondono con secondarie di bassa generazione.

2. Efficienza di Fusione e Spazio dei Parametri:

   • L'efficienza di fusione ($\eta$) aumenta generalmente con il rapporto di Eddington ($\lambda_\bullet$).
   • La dipendenza da $M_\bullet$ è non banale; l'efficienza e la massa massima del residuo tendono a raggiungere il picco a masse di SMBH intermedie e grandi, ma nessun comportamento caratteristico singolo copre l'intero spazio dei parametri. Alcune configurazioni ad alta massa non mostrano fusioni, indicando che un'elevata $M_\bullet$ da sola non garantisce un assemblaggio efficiente.

3. Firme di Spin ed Eccentricità:

   • Spin Effettivo ($\chi_{eff}$): La distribuzione di $\chi_{eff}$ dipende fortemente dalla frazione prograda iniziale ($f_{progr}$). Le configurazioni prograde producono $\chi_{eff}$ positivi, mentre quelle retrograde producono valori negativi. Le frazioni intermedie producono una distribuzione bimodale. Esiste una correlazione tra massa primaria e $|\chi_{eff}|$, poiché le fusioni ripetute rinforzano le orientazioni coerenti dello spin.
   • Eccentricità: La maggior parte delle fusioni entra nella banda dei detector a terra (a $f_{peak} = 10$ Hz) con basse eccentricità ($e \sim 10^{-3}–10^{-2}$). Tuttavia, una sottopopolazione mantiene un'alta eccentricità ($e \gtrsim 0.8$), essendosi fusa prima che la circolarizzazione fosse significativa.

4. Confronto con Eventi GW:

   • Il canale AGN può riprodurre sistemi coerenti con gli eventi massicci GW190521 e GW231123, sebbene questi occupino le code delle distribuzioni previste.
   • I sistemi corrispondenti ai parametri di massa di GW190521 si verificano a un tasso di $\sim 2 \times 10^{-3}$ nei modelli a bassa viscosità, mentre i sistemi simili a GW231123 sono più rari ($\sim 3 \times 10^{-4}$).
   • Il modello tende a sottostimare la massa secondaria ($m_2$) per questi specifici eventi rispetto ai parametri inferiti, suggerendo che, sebbene il canale AGN possa produrre la massa primaria e le caratteristiche di spin, far corrispondere l'intera distribuzione di massa dei componenti è una sfida.

5. Robustezza del Modello:

   • Le proprietà generali della popolazione (coda ad alta massa, correlazioni di spin) sono robuste rispetto alle variazioni nei formalismi di torque di migrazione (es. correzioni termiche) e nelle prescrizioni di indurimento del gas (IG20 vs. IG24 vs. C24).
   • Dipendenza Critica: La formazione di resti massicci e la coda ad alta massa sono completamente soppresse quando l'indurimento del gas viene disattivato, confermando che l'evoluzione binaria guidata dal gas è un ingrediente chiave per la crescita gerarchica nei dischi di AGN.

Significato e Rivendicazioni
L'articolo sostiene di fornire un framework completo e aggiornato che mappa sistematicamente lo spazio dei parametri degli AGN per valutare la robustezza delle previsioni di fusione BBH. Gli autori affermano che il canale AGN non produce una singola popolazione caratteristica, ma piuttosto una vasta diversità di risultati. Lo studio identifica firme distintive a livello di popolazione — specificamente la correlazione tra massa primaria e $|\chi_{eff}|$, la tendenza verso rapporti di massa diseguali ad alte masse e la conservazione della memoria riguardo alle configurazioni orbitali iniziali nella distribuzione dello spin.

Il lavoro sottolinea che, sebbene il canale AGN possa spiegare eventi BBH massicci come GW190521 e GW231123, questi eventi rappresentano probabilmente occorrenze rare all'interno del canale. Gli autori concludono che l'evoluzione binaria guidata dal gas è essenziale per un'efficiente crescita gerarchica, e che le specifiche assunzioni fisiche riguardanti la struttura del disco e l'indurimento influenzano significamente l'efficienza della formazione di resti massicci, sebbene l'esistenza qualitativa di una coda ad alta massa rimanga una caratteristica robusta del canale.