The multimessenger view of Pulsar Timing Array black holes with the Horizon-AGN simulation

Enunciado do Problema
A origem e a evolução dos buracos negros supermassivos (MBHs) permanecem incertas, particularmente em relação às suas fusões no Universo Local. Enquanto as Pulsar Timing Arrays (PTAs) detectaram recentemente um sinal de fundo de ondas gravitacionais (GWB) consistente com uma superposição estocástica de binários de buracos negros massivos (MBHBs), existe uma tensão persistente entre a amplitude observada e as previsões teóricas, que frequentemente favorecem amplitudes menores. A interpretação do espectro do GWB é complicada pela duração finita dos conjuntos de dados e pela potencial presença de fontes de ondas contínuas (CW) individualmente resolvíveis que se destacam acima do fundo. Os métodos atuais para comparar modelos teóricos com dados de PTA dependem de suposições simplificadas que podem introduzir vieses. Este artigo visa abordar essas questões utilizando a simulação cosmológica Horizon-AGN para estudar as propriedades dos MBHBs que contribuem para a banda de PTA, gerar populações realistas para previsões de GWB e CW, e avaliar criticamente as metodologias utilizadas para inferir as propriedades espectrais do GWB a partir dos resíduos de tempo.

Metodologia
O estudo utiliza a simulação hidrodinâmica cosmológica Horizon-AGN, que modela um volume de $142 \, \text{cMpc}^3$ com uma cosmologia $\Lambda$CDM padrão. Os autores construíram um catálogo de fusões de MBH rastreando IDs de MBH através da simulação e identificando fusões numéricas. Para contabilizar o atraso físico entre a fusão numérica (a uma separação de $\sim 4 \, \text{kpc}$) e a coalescência real, aplicaram um pipeline de pós-processamento incorporando:

1. Fricção Dinâmica: Calculada com base na fricção dinâmica estelar dentro de uma aproximação de esfera isotérmica.
2. Endurecimento de Binário (Binary Hardening): Modelado via endurecimento estelar, torques viscosos em discos circumbinários e emissão de ondas gravitacionais (GW), seguindo prescrições estabelecidas (ex: Sesana & Khan 2015; Dotti et al. 2015).
3. Evolução da Excentricidade: Rastreada via espalhamento de perda de cone (loss-cone scattering), arrasto viscoso e emissão de GW, atribuindo excentricidades iniciais baseadas em resultados numéricos anteriores.

Para estimar o GWB, os autores empregaram duas abordagens distintas:

• Abordagem de Ensamble Gaussiano (GE): Um cálculo determinístico do espectro de deformação característica assumindo uma distribuição contínua de fontes, frequentemente usado na literatura.
• Abordagem de População: Um método estocástico gerando 2.000 "realizações do Universo" ao desenhar binários em fase de inspiração a partir de uma distribuição de Poisson baseada no catálogo de fusões. Isso contabiliza a variância cósmica e a natureza discreta das fontes.

Os autores simularam resíduos de tempo de PTA para essas populações usando dois métodos:

1. Ensamble Gaussiano: Convertendo a deformação característica em densidade de potência espectral (PSD) de resíduos de tempo via uma resposta de pulsar média.
2. Abordagem de População: Somando as contribuições individuais de deformação de onda gravitacional do termo da Terra de todos os binários em uma realização, modelando explicitamente a duração de observação finita ($T_{\text{obs}}$) e efeitos de vazamento espectral via funções sinc.

Eles compararam essas previsões teóricas contra conjuntos de dados de "toy" PTA representando as sensibilidades atuais da EPTA (25 pulsares, 10,33 anos) e futuras da SKA (80 pulsares, 10,33 anos). Eles também identificaram potenciais candidatos a CW (fontes individuais com Relação Sinal-Ruído $> 3$) e analisaram suas propriedades eletromagnéticas (ópticas, NIR, rádio, raios-X) ao modelar a emissão de AGN e as propriedades das galáxias hospedeiras.

Principais Contribuições e Resultados

• Inferência Espectral do GWB e Viés: O estudo destaca uma importante ressalva metodológica ao comparar modelos teóricos com dados de PTA. A abordagem padrão de Ensamble Gaussiano, que assume uma duração de observação infinita, falha em contabilizar o vazamento espectral causado pela $T_{\text{obs}}$ finita. Quando a abordagem de população é usada com análise de PTA realista (incluindo marginalização do modelo de tempo), o espectro de GWB inferido é enviesado para inclinações mais rasas e amplitudes mais altas em $f = 1/\text{ano}$ comparado ao espectro subjacente real. Este viés, impulsionado pelo vazamento espectral de fontes de baixa frequência, pode explicar parcialmente a tensão entre as previsões teóricas e as amplitudes observadas pelas PTA.
• Estocasticidade e Contagem de Fontes: Em baixas frequências ($\sim 3 \, \text{nHz}$), o GWB é dominado por centenas a milhares de binários, confirmando a validade da aproximação de fundo estocástico. No entanto, em frequências mais altas ($\sim 10\text{--}50 \, \text{nHz}$), o número de fontes contribuintes cai significativamente (para algumas dezenas), levando à perda de estocasticidade e ao aumento da variância cósmica. Populações excêntricas mantêm a estocasticidade em frequências mais altas do que as circulares devido à emissão em harmônicos orbitais superiores.
• Propriedades das Fontes (Fundo vs. CW):
  • Fundo: O GWB é dominado por binários com massas de chirp $M_c \sim ^{10^{8.5} - 10^{9.5}} \, M_\odot$ em redshifts $z \in 0.05 - 1$. Estes são primordialmente hospedados em galáxias massivas ($M_\star > 10^{11} \, M_\odot$) dentro de halos de massa $\sim 10^{13} \, M_\odot$.
  • Candidatos a CW: Fontes individuais resolvíveis (candidatos a CW) tendem a ser mais massivas e localizadas em redshifts mais baixos ($z < 0.4$) comparadas ao fundo geral. Elas são predominantemente hospedadas nas galáxias centrais de grupos e aglomerados massivos ($M_{\text{halo}} > 10^{13} \, M_\odot$).
• Contrapartidas Eletromagnéticas:
  • Natureza: A maioria dos candidatos a CW aparece como Núcleos Ativos de Galáxias (AGN) em vez de quasares, caracterizados por baixas razões de Eddington ($f_{\text{Edd}} < 10^{-2}$).
  • Detectabilidade: Nas bandas óptica e infravermelho próximo (NIR), a maioria dos candidatos a CW é ofuscada por suas galáxias hospedeiras, tornando buscas de variabilidade de curva de luz desafiadoras e enviesadas para sistemas raros e luminosos.
  • Melhores Comprimentos de Onda: As bandas de rádio e raios-X oferecem as melhores perspectivas de identificação. Para candidatos com $f_{\text{Edd}} > 10^{-3}$, a emissão de AGN pode ofuscar a galáxia hospedeira nessas bandas. Mesmo para sistemas de baixa acreção, estimativas otimistas sugerem que até 80% (rádio) e 30% (raios-X) poderiam ser detectáveis acima da hospedeira.
• Probabilidades de Detecção: Assumindo uma linha de base de 10 anos, a probabilidade de detectar pelo menos um candidato a CW com $S/N > 3$ é estimada em 4% para sensibilidades do tipo EPTA e 20% para sensibilidades do tipo SKA.

Significância
O artigo afirma que sua principal significância reside em fornecer um framework multi-mensageiro unificado para entender o sinal de PTA usando uma simulação cosmológica de última geração. Crucialmente, ele desafia a metodologia padrão para comparar modelos teóricos de GWB com dados observacionais. Ao demonstrar que a duração finita de observação e o vazamento espectral podem inflar artificialmente a amplitude inferida do GWB e achatar a inclinação espectral, os autores sugerem que a aparente tensão entre modelos astrofísicos e observações de PTA pode ser parcialmente um artefato das técnicas de análise, e não uma falha da astrofísica subjacente. Adicionalmente, o estudo fornece orientação específica para o acompanhamento eletromagnético, sugerindo que as buscas devem priorizar as bandas de rádio e raios-X e focar em grupos e aglomerados de galáxias massivas, em vez de depender apenas da variabilidade óptica em sistemas do tipo quasar.