Implications of the nanoHertz Gravitational-Wave Background for Galactic Feedback and Massive Black Hole Growth

Resumo Técnico: Implicações do Fundo de Ondas Gravitacionais em NanoHertz para o Feedback Galáctico e o Crescimento de Buracos Negros Massivos

Enunciado do Problema
Medidas recentes por Arrays de Cronometragem de Pulsares (PTAs) forneceram evidências robustas de um fundo estocástico de ondas gravitacionais em nanoHertz (GWB), provavelmente gerado pela superposição incoerente de binárias de buracos negros supermassivos (SMBHBs). A amplitude deste fundo é sensível à função de massa dos buracos negros supermassivos (SMBHs). No entanto, existe uma tensão: simulações hidrodinâmicas cosmológicas (como IllustrisTNG, MillenniumTNG e Simba) consistentemente subestimam a amplitude do GWB inferida pelos dados das PTAs. Esta discrepância sugere que os modelos atuais de crescimento, semente e mecanismos de feedback que os regulam (feedback de AGN e estelar) podem estar incompletos. Especificamente, não está claro como variações nas prescrições de feedback galáctico alteram o extremo de alta massa da Função de Massa de Buracos Negros (BHMF) e, consequentemente, o sinal do GWB previsto.

Metodologia
Os autores empregam um quadro semianalítico baseado em Casey-Clyde et al. (2022) para ligar as BHMFs previstas pelas simulações à amplitude do GWB. A metodologia procede da seguinte forma:

1. Suites de Simulação: O estudo analisa quatro suites principais de simulação: IllustrisTNG, MillenniumTNG, Simba e o projeto CAMELS (Cosmology and Astrophysics with MachinE Learning Simulations).
2. Tamanho da Caixa e Convergência: Para abordar a sensibilidade do GWB ao extremo de alta massa da BHMF, os autores realizam estudos de convergência usando múltiplos tamanhos de caixa (variando de 25 Mpc/h a 500 Mpc/h). Eles determinam que caixas menores que ~35 Mpc/h (para TNG) e ~50 Mpc/h (para Simba) sofrem de ruído de Poisson e amostragem insuficiente de halos massivos ($M_{BH} > 10^9 M_\odot$).
3. Variações de Feedback: O estudo explora o impacto do feedback galáctico por:
   • Comparar modelos fiduciais entre diferentes códigos de simulação.
   • Analisar variações "1P" (um parâmetro) nas suites CAMELS, onde 23 parâmetros astrofísicos que controlam o feedback estelar e de AGN são variados individualmente.
   • Examinar variantes específicas de feedback do Simba onde os modos de jato de AGN são removidos ou alterados (ex: "SW only", "SW + Rad").
4. Cálculo do GWB: As BHMFs das simulações são convertidas em uma população de binárias usando uma Função de Distribuição de Razão de Eddington (ERDF) para identificar quasares ativos. Esta população é escalada por um fator de normalização ($N_{cal} \approx 4$) para contabilizar a razão binária-quasar, e massas secundárias são atribuídas via uma distribuição de razão de massas. A população de binárias resultante é integrada para calcular o espectro de deformação característica e a amplitude do GWB em uma frequência de referência de $f = 1 \text{ yr}^{-1}$.

Principais Contribuições

• Quantificação do Impacto do Feedback: O artigo demonstra que a escolha dos modelos de feedback estelar e de AGN altera significativamente o extremo de alta massa da BHMF. Feedback de AGN ineficiente ou ausente leva a uma abundância excessiva de SMBHs massivos, aumentando a amplitude prevista do GWB por fatores de 2 (em variações fiduciais) a 10 (em variações extremas do CAMELS).
• Análise de Convergência: Os autores estabelecem os tamanhos mínimos de caixa necessários para previsões robustas do GWB, identificando que as estatísticas de $M_{BH} > 10^9 M_\odot$ não estão convergidas nas caixas padrão de 25 Mpc/h do CAMELS.
• Subestimação Sistemática: O estudo confirma que os modelos fiduciais de IllustrisTNG, MillenniumTNG e Simba subestimam a amplitude do GWB inferida pelas PTAs por um fator de aproximadamente 2. Para corresponder às observações, estes modelos exigiriam de 1,5 a 5,5 vezes mais SMBHs na faixa de $M_{BH} > 10^8 M_\odot$.

Resultados

• Modelos Fiduciais: Todos os modelos de simulação padrão (fiduciais) subestimam a amplitude do GWB em comparação com os dados de 15 anos do NANOGrav. Isso implica que as prescrições de feedback atuais suprimem o crescimento de SMBHs com eficiência excessiva ou que a semente ocorre muito tarde.
• Variantes de Feedback: Simulações com feedback de AGN reduzido ou ausente (especificamente execuções do Simba sem modos de jato) produzem amplitudes de GWB consistentes com os dados das PTAs. No entanto, estes modelos falham em reproduzir populações galácticas realistas (ex: falham em extinguir a formação estelar efetivamente), tornando-os fisicamente irreais.
• Sensibilidade de Parâmetros: Nos conjuntos 1P do CAMELS, parâmetros específicos impulsionam a amplitude do GWB:
  • Para CAMELS-TNG, o parâmetro $\beta$ (controlando a troca entre estados de AGN de baixa e alta acreção) tem o maior impacto. Um valor de $\beta=0$ corresponde aos dados das PTAs, mas subestima a função de massa estelar de baixa massa.
  • Para CAMELS-Simba, o parâmetro de razão de massa de semente ($A_{BHseed}$) é o mais crítico. Reduzir o limiar de semente permite um crescimento de SMBH mais precoce, correspondendo aos dados das PTAs, mas pode perturbar a formação estelar inicial.
• Feedback Estelar: Embora frequentemente assumido como afetando apenas buracos negros de baixa massa, o estudo encontra que certos parâmetros de feedback estelar (ex: normalização da velocidade do vento e fração térmica) também influenciam a BHMF de alta massa e a amplitude do GWB, regulando o combustível disponível para o crescimento do buraco negro.

Significado e Alegações
O artigo alega que as medições das PTAs fornecem uma sonda poderosa e independente da física do feedback galáctico, complementando restrições eletromagnéticas tradicionais. A incompatibilidade entre as previsões das simulações e os dados das PTAs sugere que o crescimento de SMBHs no universo primordial é mais eficiente ou ocorre mais cedo do que capturado pelos modelos atuais.

Os autores concluem que, embora modificar parâmetros de feedback ou prescrições de semente possa alinhar artificialmente as simulações com os dados das PTAs, tais modificações frequentemente quebram o acordo com estatísticas galácticas observadas (como a função de massa estelar). Portanto, a tensão entre a amplitude do GWB e as previsões das simulações indica uma necessidade fundamental de reavaliar as prescrições de semente e crescimento de SMBHs em simulações hidrodinâmicas cosmológicas. O estudo enfatiza que modelos futuros devem simultaneamente reproduzir populações galácticas realistas e a amplitude observada do GWB, sugerindo que os modelos de feedback atuais podem estar suprimindo excessivamente o crescimento de SMBHs. Os autores notam que seu quadro depende de suposições simplificadas (ex: negligenciar o endurecimento ambiental e detalhes do tempo de fusão), e que um tratamento autoconsistente da evolução binária é necessário para precisão futura.