Chasing Gamma-Ray Signals from Binary Neutron Star Coalescences with the Cherenkov Telescope Array: Prospects and Observing Strategies

Resumen Técnico: Perspectivas y Estrategias de Observación de Señales de Rayos Gamma de Coalescencias de Estrellas de Neutrones Binarias con el CTA

1. El Problema
La detección de la onda gravitacional GW170817 y su contraparte electromagnética (GRB 170817A) marcó el inicio de la astronomía multimensajero. Sin embargo, la detección de emisión en el rango de muy alta energía (VHE, >20 GeV) proveniente de estallidos de rayos gamma cortos (sGRB) asociados a fusiones de estrellas de neutrones binarias (BNS) sigue siendo un desafío. Aunque se han detectado fotones TeV en GRBs largos, la evidencia en sGRBs es limitada (solo un indicio en GRB 160821B por MAGIC).

El problema central abordado en este trabajo es determinar la viabilidad y las mejores estrategias para que el Observatorio del Telescopio de Cherenkov (CTAO) detecte estos contrapartes VHE durante el quinto ciclo de observación (O5) de la red LIGO-Virgo-KAGRA (LVK). Los desafíos incluyen:

• La gran incertidumbre en la localización del cielo de las alertas de ondas gravitacionales (que puede abarcar cientos de grados cuadrados).
• La naturaleza "desviada" (off-axis) de la mayoría de los eventos detectados, lo que reduce drásticamente el flujo de rayos gamma recibido.
• La necesidad de optimizar el tiempo de exposición y la estrategia de barrido (tiling) para maximizar la probabilidad de detección antes de que la señal decaiga.

2. Metodología
Los autores implementaron una cadena de simulación completa y realista que integra tres componentes principales:

• Simulación de Eventos de Ondas Gravitacionales (GW): Se utilizó una muestra de 2307 fusiones de BNS simuladas, basadas en distribuciones astrofísicas realistas de masas, espines y distancias, considerando la sensibilidad esperada de la red LVK en el ciclo O5. Para cada evento, se estimó la localización en el cielo (mapa de probabilidad del 90% de región creíble) utilizando el código BAYESTAR.
• Modelado de la Emisión Electromagnética (sGRB): Se asumió que las fusiones lanzan chorros relativistas estructurados (con perfiles gaussianos de energía y factor de Lorentz). La emisión en VHE se estimó mediante prescripciones fenomenológicas basadas en observaciones de GRBs largos y cortos, asumiendo que la luminosidad en TeV es comparable a la de rayos X blandos y decae a una tasa similar. Se calcularon curvas de luz y espectros tanto para observadores en el eje (on-axis) como fuera del eje (off-axis), considerando la atenuación por la luz de fondo extragaláctica (EBL).
• Simulación de la Respuesta del CTAO: Se utilizó la configuración "Alpha" del CTAO (dos arrays, Norte y Sur) y sus Funciones de Respuesta Instrumental (IRF) basadas en simulaciones Prod5. Se desarrolló un pipeline que simula la respuesta del telescopio a cada evento, calculando el tiempo de exposición necesario para alcanzar una significancia de 5$\sigma$ en función del retraso (latencia) desde la fusión.
• Estrategias de Observación y Programación: Se utilizó el software tilepy para planificar las observaciones. Se evaluaron cinco estrategias de seguimiento:
  1. Exposiciones fijas (1, 5 y 20 minutos).
  2. Exposiciones variables basadas en la evolución temporal de un GRB promedio.
  3. Exposiciones adaptativas que ajustan el tiempo según la curva de luz específica de cada evento (estrategia ideal para referencia).
     Además, se incorporó el sistema de Generación de Alertas Científicas (SAG) y el análisis en tiempo real (RTA), que permite detener el barrido del cielo si se detecta un "punto caliente" (hotspot) de 5$\sigma$ y reasignar recursos para observar esa fuente específica.

3. Contribuciones Clave

• Cadena de Simulación Integral: Es el primer estudio que conecta simulaciones de GW, modelos de emisión de sGRB estructurados y la respuesta instrumental detallada del CTAO con estrategias de programación realistas (incluyendo visibilidad, ciclo día/noche y condiciones lunares).
• Evaluación de Estrategias de "Tiling": Cuantificación comparativa de diferentes duraciones de exposición y estrategias de barrido para cubrir regiones de incertidumbre grandes.
• Análisis del Rol del RTA: Demostración cuantitativa de cómo el análisis en tiempo real transforma la eficiencia de la detección al permitir el "auto-disparo" de observaciones focalizadas.
• Identificación de Parámetros Críticos: Establecimiento de que el ángulo de visión efectivo ($\theta_{eff}$) y la latencia de la alerta son los factores determinantes más importantes para la detectabilidad.

4. Resultados Principales

• Tasa de Detección: Se estima que aproximadamente el 5% de los sGRB asociados a GW detectables por LVK en el ciclo O5 producirán radiación GeV-TeV detectable por el CTAO. Esta tasa aumenta significativamente si se considera solo la fracción de eventos que son teóricamente detectables (sin restricciones de programación), llegando a un ~18% en la muestra simulada.
• Impacto de la Latencia y el Ángulo de Visión:
  • La detección es extremadamente sensible a la latencia. Para latencias menores a ~1 hora, la mayoría de los eventos detectables pueden observarse con exposiciones cortas (1-10 min).
  • Después de 1-4 horas, la probabilidad de detección cae drásticamente (de ~20% a casi 0% en 24 horas).
  • Los eventos "on-axis" se detectan en minutos, mientras que los "off-axis" (que constituyen el 87% de la población) requieren ~10 minutos o más, y su detección es mucho más difícil.
• Estrategia Óptima:
  • Las campañas de seguimiento de 4 horas parecen ser el punto óptimo de equilibrio entre cobertura del cielo y tiempo de telescopio. Extender la campaña más allá de este tiempo ofrece retornos decrecientes.
  • La estrategia de exposición fija de 5 minutos resultó ser la más efectiva, logrando una tasa de detección del 5.1% con una campaña de 4 horas, superando ligeramente a las estrategias de 1 minuto (más rápidas pero menos sensibles) y 20 minutos (demasiado lentas para cubrir grandes áreas).
• Importancia del RTA: La implementación del análisis en tiempo real (RTA) es crucial. Si se detecta un hotspot de 5$\sigma$, detener el barrido y observar la fuente aumenta la significancia de la detección en un factor de 2 a 10, permitiendo estudios espectroscópicos detallados.
• Influencia de la Localización: Aunque las estrategias de barrido cubren grandes áreas, la detección real depende de si el chorro está orientado hacia la Tierra. Se sugiere que incluso estimaciones aproximadas del ángulo de visión en las alertas de GW mejorarían drásticamente la priorización de objetivos.

5. Significado e Implicaciones
Este trabajo proporciona un marco fundamental para la planificación de operaciones del CTAO en la era de la astronomía multimensajero.

• Optimización de Recursos: Define estrategias operativas concretas (ej. ventanas de 5 minutos, duración de campaña de 4 horas) para maximizar las oportunidades de detección en un escenario de alta incertidumbre.
• Validación de Modelos Físicos: La detección (o los límites superiores) de estos eventos permitirá restringir los modelos de formación de chorros, la estructura de los GRBs y la física de la aceleración de partículas en entornos de fusión de estrellas de neutrones.
• Preparación para el Futuro: El marco desarrollado es escalable para futuros detectores de ondas gravitacionales (como el Einstein Telescope) y para otros tipos de eventos (fusiones de estrellas de neutrones y agujeros negros), sentando las bases para la próxima generación de observaciones VHE.

En conclusión, el estudio demuestra que, aunque la detección conjunta GW-GRB en VHE es un evento raro (~5%), el CTAO, con estrategias de observación optimizadas y análisis en tiempo real, tiene el potencial único de capturar estos eventos, proporcionando información física inigualable sobre los remanentes de las fusiones de estrellas de neutrones.