Is GW231123 a hierarchical merger?

Resumo Técnico: GW231123 como uma Fusão Hierárquica

1. O Problema

O evento de onda gravitacional GW231123 representa um desafio significativo para os modelos atuais de formação de binárias de buracos negros. Este evento destaca-se por duas características extremas:

• Massa: É o evento mais massivo observado até a data, com massas dos componentes $m_1 \approx 137 M_\odot$ e $m_2 \approx 103 M_\odot$. Essas massas situam-se no "gap de massa" (entre 50 e 130 $M_\odot$), onde a formação estelar isolada é teoricamente suprimida devido à instabilidade de pares (pair-instability).
• Spins (Rotação): Apresenta os spins componentes mais altos já medidos, com valores dimensionais $\chi_1 = 0.90^{+0.10}_{-0.19}$ e $\chi_2 = 0.80^{+0.20}_{-0.51}$ (intervalos de credibilidade de 90%).

A física estelar padrão (colapso de estrelas isoladas) prevê que buracos negros de massa estelar nasçam com spins muito baixos ($\chi < 0.01$) devido ao transporte eficiente de momento angular. Portanto, explicar a origem de buracos negros tão massivos e com rotação extrema exige cenários astrofísicos alternativos. Uma hipótese promissora é a montagem hierárquica, onde buracos negros resultantes de fusões anteriores (gerações 2G, 3G, etc.) se fundem novamente em ambientes densos, como aglomerados estelares globulares. Nesses cenários, o remanescente herda momento angular orbital, resultando em spins elevados ($\chi_f \approx 0.7$).

O objetivo deste trabalho é testar estatisticamente se o GW231123 pode ser explicado como a fusão de dois buracos negros de segunda geração (2G+2G) em um aglomerado estelar, considerando as incertezas nos modelos de waveform e ruído do detector.

2. Metodologia

Os autores empregaram uma abordagem estatística rigorosa baseada no formalismo desenvolvido por Passenger et al. (2024) para avaliar a consistência de eventos extremos com modelos populacionais.

• Modelo Populacional: Utilizaram a distribuição de spins de remanescentes ($\chi_f$) para buracos negros de segunda geração em aglomerados globulares, derivada de simulações do código Cluster Monte Carlo (CMC) por Borchers et al. (2025). O modelo assume spins de nascimento ($\chi_{birth}$) conservadores de 0.2 e marginaliza sobre a velocidade de escape do aglomerado e razões de massa.
• Estatística de Evidência Normalizada ($Z$): Definiram uma estatística $Z$ que compara a evidência de que o sinal é extraído do modelo populacional hierárquico ($\pi_{2G}$) contra uma prior uniforme ($\pi_U$).
  $$Z \equiv \frac{\int d\theta \, \mathcal{L}(\tilde{d}|\theta) \pi_{2G}(\chi_1, \chi_2) \pi(\eta)}{\int d\theta \, \mathcal{L}(\tilde{d}|\theta) \pi_U(\chi_1, \chi_2) \pi(\eta)}$$
  Onde $\mathcal{L}$ é a verossimilhança (likelihood) e $\theta$ são os parâmetros da binária.
• Simulação e Validação:
  1. Geraram injeções de sinais simulados baseados na distribuição $\pi_{2G}$ e injetaram em ruído gaussiano colorido (usando a densidade espectral de potência real do GW231123).
  2. Realizaram análises de parâmetros para essas injeções e para o evento real GW231123.
  3. Calcularam o valor $p$ como a probabilidade de obter um valor de $Z$ menor ou igual ao observado para o GW231123, dado o modelo populacional: $p = \int_{Z \le Z_{event}} dZ \, \pi_{pop}(Z)$.
• Análise de Waveforms: O estudo foi realizado com dois modelos de waveform distintos para abordar as incertezas sistemáticas:
  • IMRPhenomXO4a (XO4a): Modelo favorecido pelo fator de Bayes na análise original, mas com calibração limitada em spins extremos.
  • NRSur7dq4 (NRSur): Modelo baseado em Relatividade Numérica (NR), mais preciso fisicamente para spins altos, mas com limitações no intervalo de massas e razões de massa.

3. Resultados Principais

Os resultados demonstram uma forte dependência do modelo de waveform utilizado na conclusão:

• Análise com XO4a: O evento GW231123 é excluído da distribuição de fusões hierárquicas 2G+2G com um valor $p \ll 1\%$. O modelo XO4a infere spins tão altos ($\chi_1 \gtrsim 0.8$) que são inconsistentes com a distribuição prevista para remanescentes 2G+2G em aglomerados com velocidades de escape típicas.
• Análise com NRSur: O evento é compatível com o modelo hierárquico, com um valor $p \approx 5\%$. O modelo NRSur permite uma faixa de spins ligeiramente menor para o componente primário, tornando possível (embora pouco provável) que os spins observados sejam consistentes com a fusão de dois buracos negros de segunda geração.
• Incertezas Sistemáticas: A discrepância entre os resultados destaca a sensibilidade das inferências a:
  1. Sistemáticas de Waveform: Modelos não calibrados para spins extremos ($\chi > 0.8$) podem enviesar as estimativas.
  2. Ruído Não Estacionário: O GW231123 é um sinal curto (poucos ciclos) e sofreu interferência de "glitches" (falhas transitórias) nos detectores H1 e L1, o que pode contaminar a inferência de spin.

4. Contribuições Chave

• Teste Quantitativo de Cenários Hierárquicos: Fornecem um método estatístico robusto para testar se eventos extremos são consistentes com modelos de montagem hierárquica, indo além de comparações qualitativas.
• Identificação de Dependência de Modelo: Demonstram que a conclusão sobre a natureza hierárquica do GW231123 não é absoluta, mas depende criticamente da escolha do modelo de waveform e da qualidade dos dados.
• Análise de Limites de Calibração: Evidenciam a necessidade urgente de desenvolver e calibrar modelos de waveform em regimes de spin extremo e massa alta, onde os modelos atuais divergem significativamente.
• Avaliação de Cenários Alternativos: Discutem brevemente que fusões envolvendo gerações mais altas (3G+2G) ou mistas (2G+1G) também são possíveis, mas podem ter taxas de ocorrência muito baixas ou exigir ambientes com velocidades de escape ainda maiores (como núcleos de aglomerados estelares).

5. Significado e Conclusão

O estudo conclui que, embora seja possível (com $p \approx 5\%$ sob o modelo NRSur) que o GW231123 seja o resultado de uma fusão hierárquica de dois buracos negros de segunda geração em um aglomerado estelar, a evidência não é conclusiva devido às grandes incertezas sistemáticas.

A interpretação do evento como uma fusão hierárquica é plausível, mas não definitiva. A conclusão final depende de:

1. Melhorias nos modelos de waveform para spins extremos.
2. Modelos de ruído mais flexíveis que tratem adequadamente os glitches.
3. Futuras detecções de eventos similares que possam restringir melhor as distribuições populacionais.

O trabalho reforça que o GW231123 é um candidato forte para a física de montagem hierárquica, mas sua confirmação exige avanços na modelagem de ondas gravitacionais e na compreensão do ruído dos detectores.