Searching For Fast Radio Transients And Radio Pulsars Using SPOTLIGHT

Resumen Técnico: Búsqueda de Transitorios de Radio y Púlsares de Radio Rápidos Utilizando SPOTLIGHT

Planteamiento del Problema
El campo de la radioastronomía está transitando del estudio de fuentes "estáticas" a la exploración de fenómenos transitorios que ocurren en escalas de tiempo de milisegundos a microsegundos, incluyendo púlsares, Transitorios de Radio en Rotación (RRATs) y Explosiones de Radio Rápidas (FRBs). Aunque el Telescopio de Radio Gigante de Metros (GMRT) posee alta sensibilidad, sus capacidades tradicionales de sondeo están limitadas por su modo de arreglo de fase, que restringe el número de haces simultáneos a un máximo de cuatro. Esta limitación obstaculiza la velocidad de sondeo necesaria para detectar eventos transitorios raros en todo el cielo. Además, la llegada de la computación de alto rendimiento (HPC) y las Unidades de Procesamiento Gráfico (GPUs) ha creado una avalancha de datos que el software existente de radioastronomía en el dominio del tiempo no puede manejar eficientemente. Existe una necesidad urgente de actualizar las tuberías de software para procesar volúmenes masivos de datos en tiempo real, automatizar la selección de candidatos y mitigar la Interferencia de Radiofrecuencia (RFI) sin intervención humana.

Metodología
La tesis detalla el desarrollo y la prueba de dos tuberías de software principales para el proyecto SPOTLIGHT (Sondeo de Explosiones de Radio Esporádicas mediante un Arreglo de Haces Múltiples Commensal Alimentado por GPU en HPC en el GMRT), que utiliza el GMRT en un nuevo modo de Arreglo de Fase Post-Correlación (PCPA) para cubrir el campo de visión primario con aproximadamente 2.000 haces.

1. Tubería de Búsqueda de Transitorios en Tiempo Real (Detección de FRB):

   • Arquitectura: La tubería se interconecta con un formador de haces de múltiples haces que se ejecuta en un clúster de 60 nodos de computación. Los datos se almacenan en búfer en estructuras de memoria compartida (SHM) (TEL_SHM y FRB_SHM) para manejar el alto rendimiento de ~2.000 haces con una resolución temporal de ~1,3 ms.
   • Flujo de Procesamiento:
     • Mitigación de RFI: Utiliza el software gptool para la corrección de la respuesta en banda y el marcado iterativo en el dominio del tiempo/frecuencia.
     • Des-dispersión: Procesa los datos de los haces en un rango de Medida de Dispersión (DM) de prueba de 0–2000 cm⁻³pc.
     • Búsqueda de Pulsos Únicos (SPS): Emplea filtrado de caja en el dominio del tiempo para identificar pulsos estadísticamente significativos.
     • Filtrado de Candidatos: Utiliza el agrupamiento DBSCAN para agrupar candidatos en los planos DM-tiempo y Ascensión Recta-Declinación (RA-Dec), reteniendo solo los candidatos con la mayor Relación Señal-Ruido (SNR).
     • Clasificación: Extrae características (espectros dinámicos des-dispersados y transformadas DM/diagramas de "lazo de corbata") utilizando candies y clasifica los candidatos como interferencia o señales reales utilizando la Red Neuronal Convolucional (CNN) FETCH.
   • Localización: Para los candidatos clasificados positivamente, la tubería genera imágenes de las visibilidades de ~1,3 ms para lograr una localización de segundos de arco y estima el flujo integrado.

2. Tubería de Búsqueda de Púlsares en Tiempo Casi Real:

   • Arquitectura: Diseñada para ciclos de mantenimiento semanales donde los datos se almacenan y procesan fuera de línea.
   • Flujo de Procesamiento:
     • Des-dispersión: Implementada como un algoritmo de fuerza bruta O(m × N) acelerado mediante CUDA en GPUs NVIDIA, reduciendo el tiempo de ejecución en un factor de ~6 en comparación con la ejecución en CPU de un solo hilo.
     • Des-enrojecimiento: Aplica la sustracción de la mediana móvil para eliminar el ruido de base de variación lenta (ruido rojo) y la RFI, preferido sobre el blanqueamiento espectral para púlsares de largo período.
     • Re-muestreo: Corrige las derivadas del período (hasta la segunda derivada) re-muestreando la serie temporal para tener en cuenta la deriva de frecuencia, permitiendo el plegado de señales de sistemas binarios o púlsares en desaceleración.
     • Plegado: Implementa un Algoritmo de Plegado Rápido (FFA) acelerado por GPU. A diferencia de los métodos basados en Fourier, el FFA pliega los datos en períodos de prueba estrechamente espaciados, acercándose a la sensibilidad óptima teórica para todas las señales periódicas independientemente del ciclo de trabajo.
   • Clasificación: Los perfiles plegados se pasan a un clasificador Árbol de Decisión Muy Rápido de Hellinger Gaussiano (GH-VFDT).

Contribuciones Clave
Las contribuciones específicas del autor a la colaboración SPOTLIGHT incluyen:

• Integración de la Tubería: Interconexión exitosa de la tubería de búsqueda de transitorios con el nuevo formador de haces de múltiples haces, gestionando estructuras complejas de memoria compartida para el almacenamiento en búfer de datos en tiempo real.
• Simulación y Pruebas: Desarrollo y depuración del software arachne para inyectar FRBs simulados (utilizando formatos de matriz dispersa para eficiencia) en datos digitalizados a 8 bits del GMRT. Esto permitió pruebas rigurosas de extremo a extremo de la completitud y funcionalidad de la tubería antes del despliegue completo.
• Desarrollo de Algoritmos: Creación de una implementación acelerada por GPU del Algoritmo de Plegado Rápido (FFA) para búsquedas de púlsares, incluyendo la lógica para re-muestrear series temporales para corregir derivadas de frecuencia.
• Investigación de Mitigación de RFI: Investigación de la integración de gptool en la interfaz SHM para la mitigación de RFI de múltiples haces, aunque esta integración específica no se completó dentro del cronograma del proyecto.

Resultados

• Preparación de la Tubería: La tubería de búsqueda de transitorios en tiempo real ha sido probada exhaustivamente utilizando inyecciones simuladas y se considera lista para operaciones commensales. También existe una tubería de respaldo fuera de línea.
• Rendimiento: El paso de des-dispersión acelerado por GPU logró una mejora en el factor de tiempo real de ~100 a ~600. El software de búsqueda de púlsares basado en FFA es funcional, con el plegado ejecutándose actualmente en CPUs en espera de la aceleración completa por GPU.
• Capacidad de Detección: El sistema está diseñado para detectar y localizar aproximadamente 300 FRBs durante tres años de operación en el rango de frecuencias de 300 MHz a 1460 MHz.
• Localización: La tubería demuestra exitosamente la capacidad de generar imágenes de pulsos únicos a partir de visibilidades, logrando una localización de segundos de arco (sujeta a errores de procesamiento interno notados en los casos de prueba).

Significado
El artículo posiciona a SPOTLIGHT como un paso crítico hacia la "astronomía basada en datos", aprovechando las tecnologías HPC y GPU para superar las limitaciones de velocidad de sondeo de los interferómetros actuales. Al permitir que el GMRT opere con ~2.000 haces simultáneos, el proyecto busca aumentar significativamente la tasa de descubrimiento de FRBs y púlsares. El desarrollo de tuberías automatizadas de alto rendimiento se presenta como esencial para la próxima generación de radiotelescopios (incluyendo el Cuadrado Kilómetro). La tesis concluye que estos avances en software facilitarán el estudio estadístico de las galaxias anfitrionas de FRB y las propiedades de los progenitores, lo que potencialmente llevará a una comprensión unificada de los transitorios de radio que van desde púlsares hasta transitorios de período ultra-largo. El trabajo es modesto en sus afirmaciones, centrándose en la preparación de la infraestructura de software para el próximo ciclo de observación 48 del GMRT en lugar de afirmar descubrimientos astronómicos inmediatos.