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

Planteamiento del Problema
El origen y la evolución de los agujeros negros supermasivos (MBH) siguen siendo inciertos, particularmente en lo que respecta a sus fusiones en el Universo Local. Los Pulsar Timing Arrays (PTA) han detectado recientemente una señal de fondo de ondas gravitacionales (GWB) consistente con una superposición estocástica de binarias de agujeros negros masivos (MBHB), pero existe una tensión persistente entre la amplitud observada y las predicciones teóricas, que a menudo favorecen amplitudes menores. La interpretación del espectro de la GWB se complica por la duración finita de los conjuntos de datos y la posible presencia de fuentes de onda continua (CW) individualmente resolubles que sobresalen por encima del fondo. Los métodos actuales para comparar modelos teóricos con datos de PTA suelen basarse en suposiciones simplificadas que pueden introducir sesgos. Este artículo tiene como objetivo abordar estos problemas utilizando la simulación cosmológica Horizon-AGN para estudiar las propiedades de las MBHB que contribuyen a la banda de los PTA, generar poblaciones realistas para las predicciones de GWB y CW, y evaluar críticamente las metodologías utilizadas para inferir las propiedades espectrales de la GWB a partir de los residuos de temporización.

Metodología
El estudio utiliza la simulación hidrodinámica cosmológica Horizon-AGN, que modela un volumen de $142 \, \text{cMpc}^3$ con una cosmología $\Lambda$CDM estándar. Los autores construyeron un catálogo de fusiones de MBH rastreando los ID de los MBH a través de la simulación e identificando fusiones numéricas. Para tener en cuenta el retraso físico entre la fusión numérica (a una separación de $\sim 4 \, \text{kpc}$) y la coalescencia real, aplicaron un flujo de trabajo de post-procesamiento que incorpora:

1. Fricción Dinámica: Calculada basada en la fricción dinámica estelar dentro de una aproximación de esfera isotérmica.
2. Endurecimiento de la Binaria (Binary Hardening): Modelado mediante endurecimiento estelar, torques viscosos en discos circumbinarios y emisión de ondas gravitacionales (GW), siguiendo prescripciones establecidas (p. ej., Sesana & Khan 2015; Dotti et al. 2015).
3. Evolución de la Excentricidad: Rastreada mediante la dispersión del loss-cone, arrastre viscoso y emisión de GW, asignando excentricidades iniciales basadas en resultados numéricos previos.

Para estimar la GWB, los autores emplearon dos enfoques distintos:

• Enfoque de Ensamble Gaussiano (GE): Un cálculo determinista del espectro de la deformación característica asumiendo una distribución continua de fuentes, a menudo utilizado en la literatura.
• Enfoque de Población: Un método estocástico que genera 2,000 "realizaciones del Universo" mediante el muestreo de binarias en fase de inspiración a partir de una distribución de Poisson basada en el catálogo de fusiones. Esto tiene en cuenta la varianza cósmica y la naturaleza discreta de las fuentes.

Los autores simularon los residuos de temporización de los PTA para estas poblaciones utilizando dos métodos:

1. Ensamble Gaussiano: Convirtiendo la deformación característica en la densidad de potencia espectral (PSD) de los residuos de temporización mediante una respuesta de púlsar promediada.
2. Enfoque de Población: Sumando las contribuciones individuales de la deformación de GW del término terrestre de todas las binarias en una realización, modelando explícitamente la duración de observación finita ($T_{\text{obs}}$) y los efectos de fuga espectral mediante funciones sinc.

Compararon estas predicciones teóricas contra conjuntos de datos de "juguete" de PTA que representan las sensibilidades actuales de EPTA (25 púlsares, 10.33 años) y futuras de SKA (80 púlsares, 10.33 años). También identificaron posibles candidatos a CW (fuentes individuales con Relación Señal-Ruido $> 3$) y analizaron sus propiedades electromagnéticas (ópticas, NIR, radio, rayos X) modelando la emisión de AGN y las propiedades de las galaxias anfitrionas.

Contribuciones Clave y Resultados

• Inferencia Espectral de la GWB y Sesgo: El estudio destaca una advertencia metodológica significativa al comparar modelos teóricos con datos de PTA. El enfoque estándar de Ensamble Gaussiano, que asume una duración de observación infinita, no tiene en cuenta la fuga espectral causada por una $T_{\text{obs}}$ finita. Cuando se utiliza el enfoque de población con un análisis de PTA realista (incluyendo la marginalización del modelo de temporización), el espectro de la GWB inferido está sesgado hacia pendientes más tenues y amplitudes más altas en $f = 1/\text{yr}$ en comparación con el espectro subyacente real. Este sesgo, impulsado por la fuga espectral de fuentes de baja frecuencia, puede explicar parcialmente la tensión entre las predicciones teóricas y las amplitudes observadas por los PTA.
• Estocasticidad y Recuento de Fuentes: A bajas frecuencias ($\sim 3 \, \text{nHz}$), la GWB está dominada por cientos a miles de binarias, confirmando la validez de la aproximación de fondo estocástico. Sin embargo, a frecuencias más altas ($\sim 10\text{--}50 \, \text{nHz}$), el número de fuentes contribuyentes disminuye significamente (a unas pocas decenas), lo que conduce a una pérdida de estocasticidad y un aumento de la varianza cósmica. Las poblaciones excéntricas mantienen la estocasticidad a frecuencias más altas que las circulares debido a la emisión en armónicos de orden superior.
• Propiedades de las Fuentes (Fondo vs. CW):
  • Fondo: La GWB está dominada por binarias con masas de chirp $M_c \sim 10^{8.5} - 10^{9.5} \, M_\odot$ en desplazamientos al rojo $z \in 0.05 - 1$. Estas se encuentran principalmente alojadas en galaxias masivas ($M_\star > 10^{11} \, M_\odot$) dentro de halos de masa $\sim 10^{13} \, M_\odot$.
  • Candidatos a CW: Las fuentes individuales resolubles (candidatos a CW) tienden a ser más masivas y estar localizadas en desplazamientos al rojo más bajos ($z < 0.4$) en comparación con el fondo general. Son predominantemente alojadas en las galaxias centrales de grupos y cúmulos masivos ($M_{\text{halo}} > 10^{13} \, M_\odot$).
• Contrapartidas Electromagnéticas:
  • Naturaleza: La mayoría de los candidatos a CW aparecen como Núcleos Galácticos Activos (AGN) en lugar de cuásares, caracterizados por bajas razones de Eddington ($f_{\text{Edd}} < 10^{-2}$).
  • Detectabilidad: En las bandas óptica y del infrarrojo cercano (NIR), la mayoría de los candidatos a CW son eclipsados por sus galaxias anfitrionas, lo que hace que las búsquedas de variabilidad en curvas de luz sean desafiantes y estén sesgadas hacia sistemas raros y luminosos.
  • Mejores Longitudes de Onda: Las bandas de radio y rayos X ofrecen las mejores perspectivas de identificación. Para candidatos con $f_{\text{Edd}} > 10^{-3}$, la emisión del AGN puede eclipsar a la galaxia anfitriona en estas bandas. Incluso para sistemas de baja acreción, estimaciones optimistas sugieren que hasta un 80% (radio) y un 30% (rayos X) podrían ser detectables por encima de la anfitriona.
• Probabilidades de Detección: Asumiendo una línea de base de 10 años, la probabilidad de detectar al menos un candidato a CW con $S/N > 3$ se estima en un 4% para sensibilidades tipo EPTA y un 20% para sensibilidades tipo SKA.

Significancia
El artículo afirma que su principal significancia radica en proporcionar un marco unificado de multi-mensajero para comprender la señal de los PTA utilizando una simulación cosmológica de vanguardia. Crucialmente, desafía la metodología estándar para comparar modelos teóricos de GWB con datos observacionales. Al demostrar que la duración finita de la observación y la fuga espectral pueden inflar artificialmente la amplitud inferida de la GWB y suavizar la pendiente espectral, los autores sugieren que la aparente tensión entre los modelos astrofísicos y las observaciones de los PTA puede ser parcialmente un artefacto de las técnicas de análisis y no un fallo de la astrofísica subyacente. Además, el estudio proporciona una guía específica para el seguimiento electromagnético, sugiriendo que las búsquedas deben priorizar las bandas de radio y rayos X y enfocarse en grupos y cúmulos de galaxias masivas, en lugar de depender únicamente de la variabilidad óptica en sistemas tipo cuásar.