Evidence for Three Subpopulations of Merging Binary Black Holes at Different Primary Masses

Problematica
Con il rilascio del quarto catalogo di onde gravitazionali LIGO–Virgo–KAGRA (GWTC-4.0), il numero di fusioni di buchi neri binari (BBH) confermate rilevate ha superato le 150. Sebbene le aspettative iniziali suggerissero che un piccolo numero di rilevazioni potesse bastare per identificare un singolo canale di formazione primario, i dati rivelano ora una popolazione complessa che non può essere spiegata facilmente da un unico meccanismo. I tentativi precedenti di identificare sottopopolazioni utilizzando le inclinazioni dello spin hanno affrontato sfide dovute alle difficoltà di misurazione e alle incertezze teoriche. Inoltre, sebbene esistano correlazioni globali tra i parametri, la loro interpretazione astrofisica rimane non banale. Il problema centrale affrontato è se esistano sottopopolazioni distinte di BBH in fusione, caratterizzate da chiari confini di transizione nelle masse primarie ($m_1$) e distribuzioni differenziate di rapporti di massa ($q$) e magnitudo dello spin ($\chi$).

Metodologia
Gli autori adottano un approccio fenomenologico per modellare la popolazione di BBH utilizzando i dati cumulativi GWTC-4.0 (escludendo GW190814 per ambiguità della componente secondaria e GW231123 per preoccupazioni riguardanti la sistematica delle forme d'onda ed estremi valori di spin). Impiegano l'inferenza bayesiana gerarchica per confrontare modelli con diversi livelli di complessità.

Caratteristiche metodologiche chiave includono:

• Parametrizzazione: La distribuzione della massa primaria ($p(m_1)$) è modellata come una legge di potenza spezzata con due picchi gaussiani. Le magnitudo dello spin ($p(\chi)$) e le inclinazioni ($p(\cos \theta)$) sono modellate utilizzando normali troncate e miscele di componenti isotrope/allineate, rispettivamente. Il redshift ($p(z)$) segue una legge di potenza.
• Masse di transizione: L'innovazione centrale consiste nel permettere alle distribuzioni di $q$, $\chi$ e $\cos \theta$ di cambiare forma in fino a tre masse di transizione, creando efficacemente fino a quattro sottopopolazioni. La forma funzionale della distribuzione può rimanere la stessa attraverso una transizione, ma gli iperparametri che controllano la forma possono differire.
• Selezione del modello: Gli autori utilizzano i fattori di Bayes per eseguire la selezione del modello, applicando una penalità di Occam per evitare l'overfitting. Confrontano un modello "multipop" (che permette transizioni) contro un modello fiduciale fortemente parametrizzato (gwtc4-fiducial) e varianti in cui specifiche caratteristiche vengono rimosse.

Contributi Chiave e Risultati
L'analisi fornisce prove forti di almeno tre distinte sottopopolazioni di buchi neri binari in fusione, separate da due masse di transizione: $m_t^{AB} = 27.7^{+4.1}_{-3.4} M_\odot$ e $m_t^{BC} = 40.2^{+4.7}_{-3.2} M_\odot$. Il modello "multipop" è preferito rispetto al modello fiduciale con un logaritmo naturale del fattore di Bayes di 12.3.

Le tre sottopopolazioni sono caratterizzate come segue:

1. Sottopopolazione A ($m_1 \lesssim 27.7 M_\odot$):

   • Rapporto di massa: Esibisce una distribuzione quasi piatta ($\beta_0 \approx -0.8$).
   • Spin: Caratterizzata da piccole magnitudo di spin, con supporto principalmente per $\chi \leq 0.5^{+0.1}_{-0.1}$.
   • Interpretazione: Include probabilmente contributi da molteplici canali di formazione, potenzialmente inclusi binari isolati con un efficiente trasporto del momento angolare.

2. Sottopopolazione B ($27.7 M_\odot \lesssim m_1 \lesssim 40.2 M_\odot$):

   • Rapporto di massa: Mostra una preferenza molto più netta per l'accoppiamento a massa uguale ($q \approx 1$), con $\beta_1 \approx 3.6$. Questa distinta distribuzione di accoppiamento risolve un picco precedentemente riportato nella distribuzione della massa primaria a $\sim 35 M_\odot$, suggerendo che tale picco fosse un artefatto di una errata specificazione del modello.
   • Spin: Le distribuzioni della magnitudo e dell'inclinazione dello spin sono statisticamente indistinguibili dalla Sottopopolazione A.
   • Interpretazione: Gli autori suggeriscono l'Evoluzione Chimicamente Omogenea (CHE) o le stelle di Popolazione III come potenziali origini, pur notando le sfide nel conciliare la mancanza di differenze di spin con le previsioni standard della CHE. La formazione dinamica è considerata meno probabile dato il mancato riconoscimento di distribuzioni di inclinazione dello spin distinte.

3. Sottopopolazione C ($m_1 \gtrsim 40.2 M_\odot$):

   • Spin: Caratterizzata da una distribuzione della magnitudo dello spin significativamente più ampia con un sostanziale supporto per $\chi > 0.5$ (con un picco intorno a $\chi \sim 0.7$).
   • Rapporto di massa: Mostra un supporto tenue per una preferenza verso $q \approx 0.5$.
   • Interpretazione: Le proprietà sono ampiamente coerenti con le fusioni gerarchiche (fusioni di buchi neri di generazione precedente), in particolare la distribuzione ampia dello spin e la preferenza per $q \approx 0.5$. Tuttavia, la presenza di piccoli spin all'interno di questo gruppo e la mancanza di una transizione nelle inclinazioni dello spin suggeriscono che le fusioni gerarchiche in ammassi densi potrebbero non spiegare completamente questa sottopopolazione.

Significatività e Rivendicazioni
L'articolo sostiene che i dati preferiscono fortemente l'interpretazione che queste tre sottopopolazioni esistano in congiunzione, piuttosto che come una singola popolazione continua. L'identificazione di netti confini di transizione nelle distribuzioni del rapporto di massa e della magnitudo dello spin implica che l'assunto di separabilità nella distribuzione della popolazione dei parametri BBH si interrompe sopra i $\sim 25 M_\odot$.

Gli autori sottolineano che, sebbene i loro risultati dipendano da modelli fortemente parametrizzati e siano soggetti a una potenziale errata specificazione del modello, la forza statistica del segnale e la conferma indipendente di queste masse di transizione in altre analisi della letteratura suggeriscono che l'esistenza di queste sottopopolazioni sia reale. Concludono che diversi percorsi di formazione probabilmente dominano diverse gamme di massa, con la Sottopopolazione B potenzialmente legata alla CHE o alle stelle di Popolazione III, e la Sottopopolazione C legata alle fusioni gerarchiche, sebbene scenari alternativi (come l'accrezione in dischi di AGN) rimangano possibili. Lo studio dichiara esplicitamente di non trovare evidenza di transizioni attorno al picco di $10 M_\odot$, lasciando questa sottostruttura per indagini future.