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

Enunciado do Problema
Com o lançamento do quarto catálogo de ondas gravitacionais LIGO–Virgo–Kagra (GWTC-4.0), o número de fusões de binárias de buracos negros (BBH) detectadas com confiança ultrapassou 150. Embora as expectativas iniciais sugerissem que um pequeno número de detecções poderia ser suficiente para identificar um único canal de formação primário, os dados agora revelam uma população complexa que não pode ser facilmente explicada por um único mecanismo. Tentativas anteriores de identificar subpopulações usando inclinações de spin (spin tilts) enfrentaram desafios devido a dificuldades de medição e incertezas teóricas. Além disso, embora existam correlações globais entre parâmetros, sua interpretação astrofísica permanece não trivial. O problema central abordado é se existem subpopulações distintas de BBHs em fusão, caracterizadas por fronteiras de transição claras na massa primária ($m_1$) e distribuições divergentes de razões de massa ($q$) e magnitudes de spin ($\chi$).

Metodologia
Os autores adotam uma abordagem fenomenológica para modelar a população de BBH usando os dados cumulativos do GWTC-4.0 (excluindo GW190814 devido à ambiguidade da secundária e GW231123 devido a preocupações com sistemáticas de forma de onda e spins extremos). Eles empregam inferência Bayesiana hierárquica para comparar modelos com diferentes níveis de complexidade.

Principais características metodológicas incluem:

• Parametrização: A distribuição de massa primária ($p(m_1)$) é modelada como uma lei de potência quebrada com dois picos Gaussianos. As magnitudes de spin ($p(\chi)$) e inclinações ($p(\cos \theta)$) são modeladas usando normais truncadas e misturas de componentes isotrópicos/alinhados, respectivamente. O redshift ($p(z)$) segue uma lei de potência.
• Massas de Transição: A inovação central é permitir que as distribuições de $q$, $\chi$ e $\cos \theta$ mudem de forma em até três massas de transição, criando efetivamente até quatro subpopulações. A forma funcional da distribuição pode permanecer a mesma através de uma transição, mas os hiperparâmetros que controlam a forma são permitidos a diferir.
• Seleção de Modelo: Os autores utilizam fatores de Bayes para realizar a seleção de modelos, aplicando uma penalidade de Occam para evitar o sobreajuste (overfitting). Eles comparam um modelo de "múltiplas populações" (multipop) contra um modelo fortemente parametrizado fiduciário (gwtc4-fiducial) e variantes onde características específicas são removidas.

Principais Contribuições e Resultados
A análise fornece evidências fortes para pelo menos três subpopulações distintas de binárias de buracos negros em fusão, separadas por duas massas de transição: $m_t^{AB} = 27,7^{+4,1}_{-3,4} M_\odot$ e $m_t^{BC} = 40,2^{+4,7}_{-3,2} M_\odot$. O modelo "multipop" é preferido em relação ao modelo fiduciário com um fator de Bayes de logaritmo natural de 12,3.

As três subpopulações são caracterizadas da seguinte forma:

1. Subpopulação A ($m_1 \lesssim 27,7 M_\odot$):

   • Razão de Massa: Exibe uma distribuição quase plana ($\beta_0 \approx -0,8$).
   • Spin: Caracterizada por pequenas magnitudes de spin, com suporte primário para $\chi \leq 0,5^{+0,1}_{-0,1}$.
   • Interpretação: Provavelmente inclui contribuições de múltiplos canais de formação, potencialmente incluindo binárias isoladas com transporte eficiente de momento angular.

2. Subpopulação B ($27,7 M_\odot \lesssim m_1 \lesssim 40,2 M_\odot$):

   • Razão de Massa: Mostra uma preferência muito mais nítida pelo emparelhamento de massas iguais ($q \approx 1$), com $\beta_1 \approx 3,6$. Esta distribuição de emparelhamento distinta resolve um pico anteriormente relatado na distribuição de massa primária em $\sim 35 M_\odot$, sugerindo que esse pico era um artefato de erro de especificação do modelo.
   • Spin: As distribuições de magnitude e inclinação de spin são estatisticamente indistinguíveis da Subpopulação A.
   • Interpretação: Os autores sugerem a Evolução Homogênea Química (CHE) ou estrelas População III como potenciais origens, embora notem desafios em reconciliar a falta de diferenças de spin com as previsões padrão da CHE. A formação dinâmica é considerada menos provável dada a falta de distribuições de inclinação de spin distintas.

3. Subpopulação C ($m_1 \gtrsim 40,2 M_\odot$):

   • Spin: Caracterizada por uma distribuição de magnitude de spin significativamente mais ampla com suporte substancial para $\chi > 0,5$ (com pico em torno de $\chi \sim 0,7$).
   • Razão de Massa: Mostra suporte tentativo para uma preferência em direção a $q \approx 0,5$.
   • Interpretação: As propriedades são amplamente consistentes com fusões hierárquicas (fusões de buracos negros de gerações anteriores), particularmente a ampla distribuição de spin e a preferência por $q \approx 0,5$. No entanto, a presença de spins pequenos dentro deste grupo e a falta de uma transição nas inclinações de spin sugerem que fusões hierárquicas em aglomerados densos podem não explicar totalmente esta subpopulação.

Significância e Alegações
O artigo afirma que os dados preferem fortemente a interpretação de que estas três subpopulações existem em conjunto, em vez de serem uma única população contínua. A identificação de fronteiras de transição nítidas nas distribuições de razão de massa e magnitude de spin implica que a suposição de separabilidade na distribuição populacional dos parâmetros de BBH quebra acima de $\sim 25 M_\odot$.

Os autores enfatizam que, embora seus resultados dependam de modelos fortemente parametrizados e estejam sujeitos a possíveis erros de especificação de modelo, a força estatística do sinal e a confirmação independente dessas massas de transição em outras análises da literatura sugerem que a existência destas subpopulações é real. Eles concluem que diferentes vias de formação provavelmente dominam diferentes faixas de massa, com a Subpopulação B potencialmente ligada à CHE ou estrelas População III, e a Subpopulação C ligada a fusões hierárquicas, embora cenários alternativos (como acreção em discos de AGN) permaneçam possíveis. O estudo afirma explicitamente que não encontra evidências para transições em torno do pico de $10 M_\odot$, deixando essa subestrutura para investigações futuras.