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

Resumo Técnico: Fusões de BBH impulsionadas por AGN: Populações de buracos negros e crescimento hierárquico através do espaço de parâmetros de AGN

Problema e Motivação
Os discos de núcleos galácticos ativos (AGN) emergiram como um canal promissor para a formação e fusão de binárias de buracos negros (BBHs). Embora buracos negros (BHs) de massa estelar inseridos nesses discos possam migrar, formar pares e fundir-se via torques gasosos e emissão de ondas gravitacionais (GW), as propriedades previstas para essa população permanecem altamente incertas. Estudos anteriores destacaram uma forte dependência de ingredientes físicos, como a estrutura do disco, torques de migração, eficiência de endurecimento gasoso e tempos de vida de AGN, frequentemente adotando diferentes suposições que dificultam a comparação direta. Restrições observacionais recentes sugerem que o canal de AGN não pode explicar a totalidade da população de fusões de BBH detectada por observatórios de ondas gravitacionais, mas sua contribuição específica e assinaturas distintas requerem uma exploração sistemática do espaço de parâmetros de AGN para identificar características observacionais robustas e fontes de incerteza.

Metodologia
Os autores apresentam um arcabouço semi-analítico atualizado, implementado no ambiente de software fastcluster, para modelar a formação, evolução e crescimento hierárquico de BBHs em discos de AGN. A abordagem rastreia o ciclo de vida de buracos negros de origem estelar (primeira geração, 1g) através de captura, migração radial, formação de binárias, endurecimento impulsionado por gás, perturbações dinâmicas e eventual fusão.

Principais componentes metodológicos incluem:

• Modelagem do Disco: O disco de acreção é modelado usando o modelo analítico de estado estacionário de Sirko & Goodman (2003) (SG). Os autores resolvem de forma autossuficiente as equações unidimensionais do disco para variando massas de buracos negros supermassivos ($M_\bullet$), razões de Eddington ($\lambda_\bullet$) e parâmetros de viscosidade do disco ($\alpha$). Dois valores de viscosidade são testados: $\alpha = 0,01$ (baixa viscosidade) e $\alpha = 0,1$ (alta viscosidade).
• Condições Iniciais: As simulações inicializam populações de BH com base em relações de escala de aglomerados estelares nucleares (NSC). BHs de 1g são extraídos de um catálogo de síntese populacional (sevn) assumindo metalicidade solar. Os spins são atribuídos com base em uma distribuição de Maxwell, com orientações alinhadas ao momento angular do disco.
• Física Evolutiva:
  • Migração: Os BHs migram através dos regimes de Tipo I e Tipo II, com prescrições de Grishin et al. (2024) e Kanagawa et al. (2018).
  • Formação de Pares: A formação de binárias segue um critério assistido por gás, contabilizando a fração de órbitas prógradas ($f_{progr}$).
  • Endurecimento e Encontros: Após a formação, as binárias evoluem sob torques gasosos combinados (interação com o disco circumbinário) e encontros de três corpos com BHs únicos. Os autores integram equações de evolução orbital combinando o endurecimento impulsionado por gás (usando Ishibashi & Gröbner 2024) e emissão de GW (Peters 1964). Os encontros são modelados usando uma grade pré-computada de experimentos de espalhamento pós-newtonianos (tsunami).
  • Crescimento Hierárquico: Os remanescentes de fusão são calculados usando prescrições de spin alinhado e velocidades de recuo (kick). Os remanescentes que permanecem ligados ao disco participam de gerações subsequentes de fusão ($N$-ésima geração), sujeitos a limites de tempo de vida de AGN e critérios de ejeção.
• Síntese Populacional: Para construir uma população intrínseca de BBH representativa do Universo local ($z \lesssim 0,3$), os autores ponderam os resultados das simulações de acordo com as distribuições observadas de AGN de $M_\bullet$ e $\lambda_\bullet$ (Greene & Ho 2007).

Principais Resultados
O estudo explora sistematicamente a dependência das populações de fusão de BBH em relação a $M_\bullet$, $\lambda_\bullet$ e $\alpha$, obtendo os seguintes achados:

1. Crescimento Hierárquico e Distribuição de Massa:

   • Discos de AGN podem produzir fusões repetidas, construindo uma cauda de alta massa que se estende além da lacuna de supernova de instabilidade de par (PISN) ($50 \lesssim m_1/M_\odot \lesssim 130$) para a faixa de buracos negros intermediários (IMBH).
   • Dependência da Viscosidade: O crescimento hierárquico é significativamente mais eficiente em discos de baixa viscosidade ($\alpha = 0,01$), onde densidades superficiais de gás mais altas encurtam as escalas de tempo dinâmicas, permitindo um maior número de gerações de fusão ($N_{gen} > 100$). Discos de alta viscosidade ($\alpha = 0,1$) suprimem a formação de remanescentes massivos, tipicamente saturando em gerações mais baixas ($N_{gen} \approx 20$).
   • Razões de Massa: A população evolui para razões de massa cada vez mais desiguais em altas massas primárias. Enquanto as massas primárias crescem através de fusões sucessivas, as massas secundárias saturam em valores mais baixos (tipicamente abaixo da lacuna PISN), conforme primárias de alta geração fundem-se com secundárias de baixa geração.

2. Eficiência de Fusão e Espaço de Parâmetros:

   • A eficiência de fusão ($\eta$) geralmente aumenta com a razão de Eddington ($\lambda_\bullet$).
   • A dependência de $M_\bullet$ é não trivial; a eficiência e a massa máxima do remanescente tendem a atingir o pico em massas de SMBH de intermediárias a grandes, mas nenhum comportamento característico único abrange todo o espaço de parâmetros. Algumas configurações de alta massa não mostram fusões, indicando que um grande $M_\bullet$ por si só não garante uma montagem eficiente.

3. Assinaturas de Spin e Excentricidade:

   • Spin Efetivo ($\chi_{eff}$): A distribuição de $\chi_{eff}$ depende fortemente da fração prógrada inicial ($f_{progr}$). Configurações prógradas geram $\chi_{eff}$ positivos, enquanto retrógradas geram valores negativos. Frações intermediárias produzem uma distribuição bimodal. Existe uma correlação entre a massa primária e $|\chi_{eff}|$, pois fusões repetidas reforçam as orientações de spin coerentes.
   • Excentricidade: A maioria das fusões entra na banda dos detectores baseados em solo (em $f_{peak} = 10$ Hz) com excentricidades baixas ($e \sim 10^{-3}–10^{-2}$). No entanto, uma subpopulação retém alta excentricidade ($e \gtrsim 0,8$), tendo fundido antes que a circularização significativa ocorresse.

4. Comparação com Eventos de GW:

   • O canal de AGN pode reproduzir sistemas consistentes com os eventos massivos GW190521 e GW231123, embora estes ocupem as caudas das distribuições previstas.
   • Sistemas que correspondem aos parâmetros de massa de GW190521 ocorrem a uma taxa de $\sim 2 \times 10^{-3}$ em modelos de baixa viscosidade, enquanto sistemas do tipo GW231123 são mais raros ($\sim 3 \times 10^{-4}$).
   • O modelo tende a subestimar a massa secundária ($m_2$) para esses eventos específicos em comparação com os parâmetros inferidos, sugerindo que, embora o canal de AGN possa produzir a massa primária e as características de spin, corresponder à distribuição completa de massa dos componentes é desafiador.

5. Robustez do Modelo:

   • As propriedades gerais da população (cauda de alta massa, correlações de spin) são robustas frente a variações nos formalismos de torque de migração (ex: correções térmicas) e prescrições de endurecimento gasoso (IG20 vs. IG24 vs. C24).
   • Dependência Crítica: A formação de remanescentes massivos e a cauda de alta massa são inteiramente suprimidas quando o endurecimento gasoso é desligado, confirmando que a evolução binária impulsionada por gás é um ingrediente chave para o crescimento hierárquico em discos de AGN.

Significância e Alegações
O artigo afirma fornecer um arcabouço abrangente e atualizado que mapeia sistematicamente o espaço de parâmetros de AGN para avaliar a robustez das previsões de fusão de BBH. Os autores afirmam que o canal de AGN não produz uma população característica única, mas sim uma ampla diversidade de resultados. O estudo identifica assinaturas distintivas ao nível populacional — especificamente a correlação entre massa primária e $|\chi_{eff}|$, a tendência para razões de massa desiguais em altas massas, e a retenção de memória em relação às configurações orbitais iniciais na distribuição de spin.

O trabalho enfatiza que, embora o canal de AGN possa explicar eventos de BBH massivos como GW190521 e GW231123, esses eventos provavelmente representam ocorrências raras dentro do canal. Os autores concluem que a evolução binária impulsionada por gás é essencial para um crescimento hierárquico eficiente, e que as suposições físicas específicas sobre a estrutura do disco e o endurecimento influenciam significativamente a eficiência da formação de remanescentes massivos, embora a existência qualitativa de uma cauda de alta massa permaneça uma característica robusta do canal.