Investigating the emission signatures of pulsar halo candidate HESS J1813-126

Resumen Técnico: Investigación de las firmas de emisión del candidato a halo de púlsar HESS J1813-126

1. Problema y Contexto
Las fuentes de rayos gamma extendidas que rodean a púlsares de edad media (típicamente >20 kyr) se interpretan como "halos de púlsar". En estos sistemas, pares relativistas $e^\pm$ inyectados por el púlsar difunden en el medio interestelar (ISM) y producen emisión gamma mediante dispersión Compton inversa (IC) sobre campos de radiación de fondo. Sin embargo, la naturaleza exacta de muchos candidatos sigue siendo incierta debido a la degeneración entre modelos de transporte de partículas y la falta de observaciones morfológicas detalladas a altas energías.

El objeto de estudio es HESS J1813-126, una fuente extendida de rayos gamma de teraelectronvoltios (TeV) asociada al púlsar PSR J1813-1246. Este púlsar es rápido, energético, radio-silencioso y de edad media ($\sim$43 kyr), ubicado fuera del plano galáctico. Aunque se ha sugerido que HESS J1813-126 es un halo de púlsar, su naturaleza no está confirmada, especialmente dado que no se ha detectado ninguna nebulosa de viento de púlsar (PWN) en radio o rayos X, lo que apoyaría un origen de halo en lugar de una nebulosa confinada.

2. Metodología
El estudio emplea un enfoque multifacético combinando análisis de datos observacionales y modelado teórico:

• Análisis de Datos Observacionales:
  • Fermi-LAT: Se realizó un análisis estándar y uno basado en fases (phased analysis) para el púlsar. Se buscó emisión residual fuera del pulso (off-pulse) para descartar la presencia de una PWN en GeV. Se determinó la localización, significancia y morfología de la contraparte en GeV.
  • Datos Multi-frecuencia: Se integraron datos de H.E.S.S., HAWC y LHAASO para cubrir el espectro desde GeV hasta >100 TeV.
• Modelado de la Emisión del Púlsar (Magnetosfera):
  • Se utilizó el modelo de radiación sincro-curvatura (SC) para explicar el espectro de energía (SED) combinado de rayos X y gamma del púlsar. Este modelo asume que los pares $e^\pm$ son acelerados por campos eléctricos paralelos en la magnetosfera y emiten mientras espiralizan y se deslizan a lo largo de líneas de campo curvas.
• Modelado del Halo de Púlsar (Transporte de Partículas):
  • Se resolvieron ecuaciones de transporte y pérdida dependientes del tiempo para los pares $e^\pm$ inyectados.
  • Se probaron tres modelos de transporte distintos para explicar la morfología y el SED de HESS J1813-126:
    1. Difusión Suprimida Isótropa de Dos Zonas (2ZISD): Un modelo con una zona interna de difusión lenta y una externa más rápida.
    2. Balístico a Difusivo (B2D): Donde las partículas viajan inicialmente de forma balística antes de entrar en régimen difusivo.
    3. Difusión Anisotrópica (AD): Considera la difusión a lo largo y perpendicular a las líneas del campo magnético ordenado.
  • Se calcularon los perfiles de brillo superficial (SBP) y la emisión dependiente de la apertura para cada modelo, comparándolos con los datos observados.
• Balance de Pares: Se comparó la multiplicidad de pares de alta energía en la magnetosfera (calculada a partir del modelo SC) con la requerida para alimentar el halo, normalizando ambas a la tasa de flujo de Goldreich-Julian.

3. Contribuciones Clave y Resultados

• Caracterización de la Contraparte en GeV:
  • La fuente en GeV (4FGL J1813.4-1246) es puntual y coincide con el púlsar, con un SED curvado bien descrito por una función de corte exponencial.
  • El análisis de fases reveló una fuente débil fuera del pulso (PS J1813.3-1246), pero su espectro (ley de potencia suave) y morfología puntual sugieren que es contaminación residual de la magnetosfera y no una PWN en GeV. Esto refuerza la hipótesis de halo.
• Éxito del Modelo Sincro-Curvatura:
  • El modelo SC reprodujo exitosamente el SED combinado de rayos X (componente sincrotrón) y rayos gamma (componente curvatura) de PSR J1813-1246.
  • Se identificó que la transición entre los regímenes de emisión ocurre a una distancia específica de la luz, explicando la emisión en diferentes bandas.
• Modelado del Halo y Comparación de Transportes:
  • Modelo 2ZISD: Proporcionó la descripción más consistente con los datos multi-instrumento (H.E.S.S., HAWC, LHAASO). Predice que el SED aumenta con la apertura de integración principalmente a bajas energías (debido a electrones de larga vida que se difunden más lejos), saturándose a altas energías. Los SBP muestran un pico central brillante que se estrecha con la energía.
  • Modelo B2D: Aunque ajusta los datos de H.E.S.S. dentro de una pequeña apertura, predice un flujo excesivo en HAWC y LHAASO al aumentar la apertura debido a perfiles de brillo superficial más planos y una difusión cuasi-balística. Además, requiere una eficiencia de conversión inusualmente alta ($\eta \gtrsim 10\%$).
  • Modelo AD: Captura las huellas geométricas de la difusión anisotrópica. La configuración alineada con el campo ($\zeta=0^\circ$) es consistente con los datos, mientras que la desalineada ($\zeta=10^\circ$) tiende a sobreestimar el flujo en aperturas grandes, sugiriendo que el halo es probablemente esférico o débilmente asimétrico.
• Consistencia Energética:
  • Se calculó la multiplicidad de pares de alta energía: $\kappa_{mag} \approx 5.9 \times 10^4$ (magnetosfera) y $\kappa_{halo} \approx 2.2 \times 10^2$ (halo).
  • Estos valores indican que la magnetosfera tiene un presupuesto energético suficiente para alimentar tanto la emisión SC como el halo sin necesidad de canales de creación de pares adicionales o aceleradores exóticos.

4. Significancia e Implicaciones

• Confirmación de Halo de Púlsar: El estudio proporciona evidencia sólida de que HESS J1813-126 es un halo de púlsar alimentado por PSR J1813-1246, basado en la ausencia de una PWN, la morfología extendida, la coincidencia posicional y la capacidad de los modelos de halo para reproducir los datos.
• Diagnósticos Observacionales Futuros: El trabajo destaca que los perfiles de brillo superficial (SBP) y la evolución dependiente de la apertura del SED son firmas observables críticas. Futuros telescopios de efecto Cherenkov de imagen (IACTs) de alta resolución, como el CTA, podrán distinguir entre los modelos de transporte (especialmente 2ZISD vs. B2D/AD) midiendo cómo cambia la extensión de la fuente con la energía y el tamaño de la apertura.
• PeVatrones Galácticos: La detección de emisión >100 TeV por LHAASO confirma que este sistema es un candidato a PeVatron galáctico, demostrando que los halos de púlsar pueden ser fuentes significativas de partículas de ultra-alta energía en la Vía Láctea.
• Unificación de Fenomenología: El estudio conecta exitosamente la física de la magnetosfera del púlsar (emisión SC) con la física del medio interestelar (halo), ofreciendo una imagen coherente y auto-consistente del sistema.

En conclusión, el papel establece un marco robusto para interpretar HESS J1813-126 como un halo de púlsar y proporciona herramientas diagnósticas clave para desentrañar los mecanismos de transporte de partículas en estos entornos extremos mediante observaciones futuras.