Pulsar searches of Fermi-LAT gamma-ray sources with the MWA

Resumen Técnico: Búsquedas de púlsares de fuentes de rayos gamma de Fermi-LAT con el MWA

Problema
Los púlsares de milisegundo (MSP) son herramientas críticas para la astrofísica, incluyendo la detección de ondas gravitacionales mediante arreglos de temporización de púlsares. Aproximadamente una quinta parte de los MSP galácticos conocidos han sido descubiertos mediante búsquedas de radio dirigidas a fuentes de rayos gamma no asociadas identificadas por el Telescopio de Gran Área de Fermi (Fermi-LAT). Sin embargo, estas búsquedas se han realizado históricamente a frecuencias de radio superiores a 300 MHz para mitigar la dispersión y la dispersión interestelar. Esto deja un vacío en el espacio de parámetros por debajo de los 300 MHz, donde podrían ser detectables los MSP con espectros de radio empinados. Las búsquedas a baja frecuencia enfrentan desafíos significativos debido a los efectos cromáticos en el medio interestelar ionizado frío (IISM), específicamente el costo computacional de corregir grandes retrasos dispersivos y el ensanchamiento intracanal (intrachannel smearing), lo que limita la sensibilidad a púlsares de período corto. Mientras que un estudio con el Low Frequency Array (LOFAR) utilizó exitosamente una estrategia de dedispersión semicoherente para descubrir tres MSP a 135 MHz, no se había realizado un estudio comparable utilizando el Murchison Widefield Array (MWA).

Metodología
Los autores desarrollaron un nuevo flujo de trabajo de búsqueda adaptado para el MWA para apuntar a fuentes de Fermi-LAT no asociadas utilizando datos de archivo del estudio SMART (Southern-sky MWA Rapid Two-metre) del cielo austral. La metodología involucra:

• Fuente de Datos: El estudio utilizó observaciones de formación de haces (beamformed) de 20 minutos del estudio SMART a 154 MHz, apuntando a 308 fuentes de Fermi-LAT no asociadas (seleccionadas del catálogo 4FGL-DR4 con declinaciones < +25°).
• Arquitectura del Pipeline: El pipeline emplea un esquema de dedispersión semicoherente para equilibrar el costo computacional y el ensanchamiento dispersivo.
  • Dedispersión de Fase Coherente: Utilizando el software CDMT acelerado por GPU, los voltajes complejos se dedispersan coherentemente a 27 pruebas de Medida de Dispersión (DM) que van desde 1.5 a 79.5 cm⁻³ pc en pasos de 3 cm⁻³ pc.
  • Dedispersión Incoherente: El resultado es luego dedispersado incoherentemente utilizando la biblioteca DEDISP para llenar los huecos entre los pasos coherentes, resultando en un total de 27,000 pruebas de DM.
  • Búsqueda de Aceleración: La rutina ACCELSEARCH de PRESTO realiza una búsqueda de aceleración en el dominio de Fourier (FDAS) para detectar púlsares binarios, buscando aceleraciones orbitales que corresponden a valores de $z$ entre -200 y +200.
  • Filtrado de Candidatos: Los candidatos se filtran según umbrales de relación señal-ruido, consistencia armónica y validez de DM. Los candidatos finales se pliegan y optimizan usando PREPFOLD.
• Validación: El pipeline fue validado mediante la detección ciega de cinco MSP conocidos (cuatro en binarias de órbita corta) dentro de los datos de SMART, confirmando su capacidad para detectar sistemas tipo spider a pesar de la baja frecuencia.

Contribuciones Clave

1. El mayor estudio de baja frecuencia: Este trabajo representa el mayor estudio de radio de fuentes de Fermi-LAT no asociadas hasta la fecha y el segundo estudio de este tipo realizado por debajo de 300 MHz.
2. Implementación del Pipeline: Los autores adaptaron e implementaron exitosamente un pipeline de dedispersión semicoherente para los datos del Sistema de Captura de Voltaje (VCS) del MWA, permitiendo búsquedas de MSP más profundas de lo que sería posible anteriormente con métodos totalmente incoherentes a estas frecuencias.
3. Análisis de Sensibilidad: El estudio proporciona estimaciones detalladas de la sensibilidad del MWA a 154 MHz, calculando densidades de flujo mínimas detectables ($S_{min}$) que oscilan entre aproximadamente 30 y 220 mJy (equivalente a 0.7–5.2 mJy a 1.4 GHz) para un púlsar de período de 2 ms con un ciclo de trabajo del 28%.

Resultados

• Rendimiento de Descubrimiento: No se identificaron nuevos púlsares en el estudio de las 308 fuentes. Los autores atribuyen este resultado nulo a una sensibilidad insuficiente.
• Detecciones Ciegas: El pipeline detectó exitosamente dos púlsares no reciclados conocidos (PSR J1752−2806 y PSR J2048−1616) en los lóbulos secundarios del haz de haz unido (tied-array beam), demostrando la funcionalidad del pipeline.
• Límites de Sensibilidad: El estudio fue sensible entre el 8% y el 44% de la población conocida de púlsares de rayos gamma con emisión de radio (por debajo del límite de DM impuesto por la dispersión). En contraste, el estudio de LOFAR fue sensible al 86–95% de la población y descubrió tres MSP. Los autores señalan que incluso si el estudio del MWA hubiera igualado el rendimiento de LOFAR, el número esperado de descubrimientos habría sido menor a uno, lo cual es consistente con el resultado nulo.
• Prueba Ciega: De los cinco MSP de prueba utilizados para validar el pipeline, cuatro fueron detectados ciegamente. Uno (PSR J0952−0607) no fue detectado, probablemente debido a una baja relación señal-ruido y un ensanchamiento fraccional significativo.

Significado y Perspectivas Futuras
El artículo concluye que, si bien el estudio piloto no produjo nuevos descubrimientos, establece un flujo de trabajo validado para búsquedas dirigidas de MSP con el MWA. Los autores argumentan que el resultado nulo es consistente con los límites de sensibilidad del estudio más que con un fallo del método.

La importancia de este trabajo radica en su preparación para futuras actualizaciones:

• MWA Fase III: La próxima actualización de la Fase III, que duplica el número de baldosas (tiles) a 256, se proyecta que mejorará la sensibilidad instantánea en un factor de dos. Esto aumentaría el rendimiento esperado de púlsares de rayos gamma detectables en aproximadamente un 30% para el mismo tiempo de integración (20 min).
• Formador de Haz en Tiempo Real (RTB): El desarrollo de un RTB permitirá observaciones dirigidas con el haz primario centrado más cerca de la fuente, reduciendo las pérdidas de sensibilidad por desplazamientos del haz y permitiendo tiempos de integración más largos.
• Integración Extendida: Aumentar el tiempo de integración a 60 minutos podría incrementar el rendimiento esperado en un ~60%, aunque esto reduce la sensibilidad a órbitas binarias muy cortas.
• Aplicaciones Más Amplias: El pipeline semicoherente desarrollado es aplicable a búsquedas de remanentes de supernova, cúmulos globulares y fuentes de radio de espectro empinado, proporcionando una base para futuros estudios de SKA-Low.

Los autores enfatizan que observar a frecuencias más bajas sigue siendo ventajoso para púlsares con espectros empinados, siempre que los límites de dispersión y los tiempos de integración se gestionen eficazmente mediante futuras actualizaciones de hardware y software.