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

Resumen Técnico: Perspectivas de Resolución y Localización de Binarias de Agujeros Negros Supermasivos Individuales con Arreglos de Pulsares de Tiempo (PTA)

Planteamiento del Problema
Los Arreglos de Pulsares de Tiempo (PTA) han detectado recientemente un fondo estocástico de ondas gravitacionales (GW), probablemente originado por una población de binarias de agujeros negros supermasivos (SMBH). El próximo gran hito es la resolución de binarias individuales dentro de este fondo, lo que permitiría la astronomía multi-mensajero. Sin embargo, un desafío principal es la limitada capacidad de localización en el cielo de los PTA. El área de localización de credibilidad del 90% ($A_{90}$) para la primera binaria resuelta individualmente se espera que abarque cientos de grados cuadrados, conteniendo potencialmente cientos de miles de galaxias candidatas. Además, los catálogos de galaxias electromagnéticas (EM) existentes están incompletos en regiones de alta contaminación estelar o de polvo (por ejemplo, el plano galáctico), lo que crea candidatos a anfitriones "faltantes". Actualmente no existe un método estandarizado para clasificar eficientemente este vasto número de candidatos basándose en los posteriors de GW y las propiedades EM para priorizar las observaciones de seguimiento.

Metodología
Los autores simulan poblaciones realistas de binarias SMBH consistentes con el fondo de ondas gravitacionales medido por el conjunto de datos de 15 años de NANOGrav (NG15) utilizando el software holodeck. Generan 1,000 realizaciones de poblaciones de binarias y seleccionan la binaria más ruidosa de cada población que sea detectable en futuras configuraciones de PTA pero no detectable en NG15.

Se simulan tres configuraciones de PTA para representar la sensibilidad futura:

• IPTA_20: Una línea de base de 20 años con 116 púlsares (similar a la próxima entrega de datos 3 de IPTA).
• IPTA_25 (Fiducial): Una línea de base de 25 años con 158 púlsares (añadiendo 42 nuevos púlsares a IPTA_20).
• IPTA_30: Una línea de base de 30 años con 200 púlsares.

Los autores inyectan las binarias seleccionadas en datos de temporización simulados que contienen un fondo estocástico de ondas gravitacionales y utilizan el pipeline QuickCW para restringir los parámetros de la binaria (ubicación en el cielo, masa de chirp, amplitud, frecuencia). Definen las fuentes "localizables" como aquellas con una relación señal-ruido (S/N) $\ge 8$, un umbral típicamente requerido para posiciones en el cielo bien restringidas.

Para identificar anfitriones, los autores cruzan los mapas de localización $A_{90}$ resultantes con dos catálogos EM de todo el cielo exhaustivos:

1. Catálogo ETG: Contiene $\approx 19.4$ millones de galaxias tempranas (ETGs) quiescentes con desplazamientos al rojo fotométricos y masas estelares (derivadas de 2MPZ y WISExSCOS).
2. Catálogo AGN: Contiene $\approx 4.5$ millones de candidatos a Núcleos Galácticos Activos (AGN) (catálogo R90 de Assef et al. 2018) con magnitudes aparentes en la banda W1.

Se desarrolla un sistema de clasificación de anfitriones para asignar una puntuación ($H$) a cada galaxia candidata. La puntuación es el producto de la probabilidad de posición en el cielo ($p_{sky}$) y la probabilidad de que las propiedades EM de la galaxia coincidan con los posteriors de GW ($p_{cat}$).

• Para las ETGs, $p_{cat}$ se deriva del desplazamiento al rojo de la galaxia (convertido a distancia de luminosidad) y su masa de SMBH estimada (derivada de la masa estelar mediante relaciones de escala).
• Para los AGN, $p_{cat}$ se deriva de la magnitud aparente W1, que es una función de la masa total y la distancia de luminosidad.
  Los autores calculan explícitamente el número de anfitriones "faltantes" estimando la fracción del área $A_{90}$ no cubierta por los catálogos EM.

Resultados Clave

• Detectabilidad y Localización:

  • En los próximos 0–5 años (IPTA_20), hay una probabilidad de $\approx 21\%$ de resolver una binaria ($S/N \ge 4$), pero solo una probabilidad de $\approx 0.3\%$ de lograr una localización bien restringida ($S/N \ge 8$).
  • En los próximos 5–10 años (IPTA_25), la probabilidad de resolver una binaria aumenta a $\approx 38\%$, con un $\approx 3.8\%$ de probabilidad de localización.
  • En los próximos 10 años (IPTA_30), la probabilidad de resolución alcanza $\approx 51\%$, con un $\approx 14.1\%$ de probabilidad de localización.
  • Se espera que las fuentes localizables sean relativamente cercanas (mediana $d_{lum} \approx 400$ Mpc) y muy masivas ($M_{tot} \gtrsim 10^{9.8} M_{\odot}$).

• Áreas de Localización y Conteo de Anfitriones:

  • Para la configuración fiducial (IPTA_25), 30 de las 38 binarias potencialmente localizables tienen $A_{90} < 2,000$ deg$^2$. El área $A_{90}$ mediana es de $408$ deg$^2$.
  • Estas áreas contienen una mediana de $\approx 190,000$ ETGs y $\approx 40,000$ AGN.
  • En IPTA_30, el $A_{90}$ mediano se reduce a $212$ deg$^2$, reduciendo el número mediano de anfitriones en $\approx 40\%$.

• Cobertura Incompleta:

  • La cobertura EM incompleta es una limitación significativa. En la configuración fiducial, el número mediano de ETGs anfitriones faltantes es de $\approx 20,000$ y de AGN faltantes es de $\approx 6,000$.
  • Para algunas fuentes, el $A_{90}$ está parcial o completamente fuera de la huella de los catálogos (por ejemplo, solapándose con el plano galáctico), lo que hace que el verdadero anfitrión sea inidentificable sin nuevos sondeos.

• Eficiencia de Clasificación:

  • El método de clasificación excluye con éxito las galaxias con propiedades inconsistentes con los posteriors de GW.
  • ETGs: Aproximadamente el 47% (IPTA_25) al 61% (IPTA_30) de las ETGs son excluidas debido a desplazamientos al rojo o masas de SMBH incompatibles. La eficiencia de clasificación ($E_{95-50}$, la razón entre el percentil 95 de la puntuación y la mediana) es de $\approx 12.4$ para IPTA_25 y $\approx 17.2$ para IPTA_30.
  • AGN: Solo $\approx 0.3\%$ (IPTA_25) al $3.8\%$ (IPTA_30) de los AGN son excluidos. La eficiencia de clasificación es significativamente menor ($\approx 3.7$ para IPTA_25), ya que la magnitud aparente es una combinación degenerada de masa y distancia, y la masa de chirp está pobremente restringida por los PTA.

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera cuantificación realista del número de candidatos anfitriones y el impacto de la cobertura incompleta del cielo para la primera binaria de SMBH resuelta por PTA. Los autores demuestran que, si bien el número de candidatos es grande, un sistema de clasificación basado en los posteriors de GW y las propiedades EM puede reducir eficazmente el espacio de búsqueda, particularmente para las ETGs donde la masa estelar y el desplazamiento al rojo están disponibles.

El estudio destaca que:

1. El descubrimiento multi-mensajero es plausible pero desafiante: La probabilidad de localizar una fuente en 10 años es no despreciable ($\approx 14\%$), pero las áreas de localización siguen siendo grandes.
2. La completitud del catálogo es crítica: La existencia de miles de anfitriones faltantes debido a los vacíos de los sondeos (por ejemplo, el plano galáctico) plantea un riesgo de falsos negativos en la identificación del anfitrión.
3. La calidad de los datos dicta el éxito de la clasificación: La eficiencia del método de clasificación depende fuertemente de la disponibilidad de estimaciones independientes de desplazamiento al rojo y masa de SMBH. La incapacidad actual para clasificar eficientemente los AGN se atribuye a la falta de datos de desplazamiento al rojo y masa en los catálogos de AGN y a la pobre restricción de la masa de chirp por parte de los PTA.
4. Mejoras futuras: Los autores sugieren que las mejoras futuras en las mediciones de la distancia de los púlsares (para utilizar el término del púlsar en las señales de GW) y el desarrollo de catálogos de AGN con estimaciones de desplazamiento al rojo y masa mejorarían significativamente las capacidades de localización y clasificación.

El artículo concluye que, si bien la primera detección multi-mensajero está al alcance, el éxito depende no solo de la sensibilidad de los PTA, sino también de la completitud y uniformidad de los catálogos EM y del desarrollo de pipelines de clasificación eficientes para priorizar las observaciones de seguimiento.