Prospects of resolving and localising individual supermassive black hole binaries with pulsar timing arrays: the host ranking challenge

Resumo Técnico: Perspectivas de Resolução e Localização de Binárias de Buracos Negros Supermassivos Individuais com Matrizes de Temporização de Púlsares

Definição do Problema
As Matrizes de Temporização de Púlsares (PTAs) detectaram recentemente um fundo estocástico de ondas gravitacionais (GW), provavelmente originário de uma população de binárias de buracos negros supermassivos (SMBH). O próximo grande marco é a resolução de binárias individuais dentro deste fundo, o que permitiria a astronomia multi-mensageiro. No entanto, um desafio primário é a capacidade limitada de localização no céu das PTAs. A área de credibilidade de 90% ($A_{90}$) para a primeira binária resolvida individualmente espera-se que abranja centenas de graus quadrados, potencialmente contendo centenas de milhares de galáxias hospedeiras candidatas. Além disso, os catálogos de galáxias eletromagnéticas (EM) existentes são incompletos em regiões de alta contaminação estelar ou de poeira (ex: o Plano Galáctico), criando hospedeiros candidatos "ausentes". Atualmente, não existe um método padronizado para classificar eficientamente esse vasto número de candidatos com base nos posteriors de GW e propriedades EM para priorizar observações de acompanhamento.

Metodologia
Os autores simulam populações realistas de binárias de SMBH consistentes com o fundo de ondas gravitacionais medido pelo Conjunto de Dados de 15 anos do NANOGrav (NG15) usando o software holodeck. Eles geram 1.000 realizações de populações de binárias e selecionam a binária mais ruidosa de cada população que seja detectável em futuras configurações de PTA, mas indetectável no NG15.

Três configurações de PTA são simuladas para representar a sensibilidade futura:

• IPTA_20: Uma linha de base de 20 anos com 116 púlsares (semelhante à próxima liberação de dados 3 da IPTA).
• IPTA_25 (Fiducial): Uma linha de base de 25 anos com 158 púlsares (adicionando 42 novos púlsares à IPTA_20).
• IPTA_30: Uma linha de base de 30 anos com 200 púlsares.

Os autores injetam as binárias selecionadas em dados de temporização simulados contendo um fundo estocástico de ondas gravitacionais e utilizam o pipeline QuickCW para restringir os parâmetros da binária (localização no céu, massa de chirp, amplitude, frequência). Eles definem fontes "localizáveis" como aquelas com uma relação sinal-ruído (S/N) $\ge 8$, um limiar tipicamente requerido para posições no céu bem restringidas.

Para identificar hospedeiros, os autores cruzam os mapas de localização $A_{90}$ resultantes com dois catálogos EM de todo o céu abrangentes:

1. Catálogo ETG: Contém $\approx 19,4$ milhões de galáxias de tipo precoce (ETGs) quiescentes com redshifts fotométricos e massas estelares (derivados de 2MPZ e WISExSCOS).
2. Catálogo AGN: Contém $\approx 4,5$ milhões de candidatos a Núcleos Galácticos Ativos (AGN) (catálogo R90 de Assef et al. 2018) com magnitudes aparentes na banda W1.

Um sistema de classificação de hospedeiros é desenvolvido para atribuir uma pontuação ($H$) a cada galáxia candidata. A pontuação é o produto da probabilidade de posição no céu ($p_{sky}$) e da probabilidade de as propriedades EM da galáxia corresponderem aos posteriors de GW ($p_{cat}$).

• Para ETGs, $p_{cat}$ é derivado do redshift da galáxia (convertido em distância de luminosidade) e de sua massa de SMBH estimada (derivada da massa estelar via relações de escala).
• Para AGN, $p_{cat}$ é derivado da magnitude aparente W1, que é uma função da massa total e da distância de luminosidade.

Os autores calculam explicitamente o número de hospedeiros "ausentes" estimando a fração da área $A_{90}$ não coberta pelos catálogos EM.

Principais Resultados

• Detectabilidade e Localização:

  • Nos próximos 0–5 anos (IPTA_20), há uma probabilidade de $\approx 21\%$ de resolver uma binária ($S/N \ge 4$), mas apenas uma chance de $\approx 0,3\%$ de obter uma localização bem restringida ($S/N \ge 8$).
  • Nos próximos 5–10 anos (IPTA_25), a probabilidade de resolver uma binária sobe para $\approx 38\%$, com uma chance de $\approx 3,8\%$ de localização.
  • Em 10 anos (IPTA_30), a probabilidade de resolução atinge $\approx 51\%$, com uma chance de $\approx 14,1\%$ de localização.
  • Espera-se que as fontes localizáveis sejam relativamente próximas (mediana $d_{lum} \approx 400$ Mpc) e muito massivas ($M_{tot} \gtrsim 10^{9,8} M_{\odot}$).

• Áreas de Localização e Contagem de Hospedeiros:

  • Para a configuração fiducial (IPTA_25), 30 de 38 binárias potencialmente localizáveis possuem $A_{90} < 2.000$ deg$^2$. A mediana do $A_{90}$ é $408$ deg$^2$.
  • Estas áreas contêm uma mediana de $\approx 190.000$ ETGs e $\approx 40.000$ AGN.
  • Em IPTA_30, a mediana do $A_{90}$ encolhe para $212$ deg$^2$, reduzindo o número mediano de hospedeiros em $\approx 40\%$.

• Cobertura Incompleta:

  • A cobertura EM incompleta é uma limitação significativa. Na configuração fiducial, o número mediano de ETG hospedeiros ausentes é $\approx 20.000$ e de AGN ausentes é $\approx 6.000$.
  • Para alguns casos, o $A_{90}$ está parcialmente ou completamente fora da pegada dos catálogos (ex: sobrepondo o Plano Galáctico), tornando o verdadeiro hospedeiro inidentificável sem novos levantamentos.

• Eficiência de Classificação:

  • O método de classificação exclui com sucesso galáxias com propriedades inconsistentes com os posteriors de GW.
  • ETGs: Aproximadamente 47% (IPTA_25) a 61% (IPTA_30) das ETGs são excluídas devido a redshifts ou massas de SMBH incompatíveis. A eficiência de classificação ($E_{95-50}$, a razão entre o percentil 95 da pontuação e a mediana) é de $\approx 12,4$ para IPTA_25 e $\approx 17,2$ para IPTA_30.
  • AGN: Apenas $\approx 0,3\%$ (IPTA_25) a $\approx 3,8\%$ (IPTA_30) dos AGN são excluídos. A eficiência de classificação é significativamente menor ($\approx 3,7$ para IPTA_25), pois a magnitude aparente é uma combinação degenerada de massa e distância, e a massa de chirp é mal restringida por PTAs.

Significância e Alegações
O artigo afirma fornecer a primeira quantificação realista do número de candidatos hospedeiros e do impacto da cobertura incompleta do céu para a primeira binária de SMBH resolvida por PTA. Os autores demonstram que, embora o número de candidatos seja grande, um sistema de classificação baseado em posteriors de GW e propriedades EM pode reduzir efetivamente o espaço de busca, particularmente para ETGs onde a massa estelar e o redshift estão disponíveis.

O estudo destaca que:

1. A descoberta multi-mensageiro é plausível, mas desafiadora: A probabilidade de localizar uma fonte dentro de 10 anos é não-negligível ($\approx 14\%$), mas as áreas de localização permanecem grandes.
2. A completude do catálogo é crítica: A existência de milhares de hospedeiros ausentes devido a lacunas de levantamento (ex: o Plano Galáctico) apresenta um risco de falsos negativos na identificação do hospedeiro.
3. A qualidade dos dados dita o sucesso da classificação: A eficiência do método de classificação depende fortemente da disponibilidade de estimativas independentes de redshift e massa de SMBH. A incapacidade atual de classificar eficientemente AGN é atribuída à falta de dados de redshift/massa em catálogos de AGN e à dificuldade de restringir a massa de chirp por PTAs.
4. Melhorias futuras: Os autores sugerem que melhorias futuras nas medições de distância de púlsares (para utilizar o termo do púlsar em sinais de GW) e o desenvolvimento de catálogos de AGN com estimativas de redshift e massa aumentariam significativamente as capacidades de localização e classificação.

O artigo conclui que, embora a primeira detecção multi-mensageiro esteja ao alcance, o sucesso depende não apenas da sensibilidade da PTA, mas também da completude e uniformidade dos catálogos EM e do desenvolvimento de pipelines de classificação eficientes para priorizar observações de acompanhamento.