A Multiwavelength Interpretation of HESS J1857+026 Emission Using the Fermi-LAT, VERITAS, and HAWC Observatories

Resumen Técnico: Una Interpretación Multilongitud de Onda de la Emisión de HESS J1857+026

Planteamiento del Problema
HESS J1857+026 es una fuente de rayos $\gamma$ de muy alta energía (VHE) extendida e no identificada, coincidente con el púlsar energético PSR J1856+0245. Aunque la fuente es una candidata a nebulosa de viento de púlsar (PWN), su naturaleza multi-longitud de onda sigue siendo compleja. Observaciones previas han revelado un desajuste espacial entre las extensiones de la emisión en GeV y TeV, así como una falta de contrapartes de rayos X confirmadas, lo que plantea interrogantes sobre el mecanismo de emisión dominante. Específicamente, no está claro si la emisión es puramente leptónica (dispersión Compton inversa de una PWN), puramente hadrónica (decaimiento de piones de una remanente de supernova interactuando con gas molecular), o una combinación de ambas. Además, las propiedades de difusión de partículas en el entorno inmediato de las PWN evolucionadas permanecen insuficientemente restringidas en comparación con el medio interestelar (ISM) promedio.

Metodología
Este estudio presenta un análisis exhaustivo multi-instrumental de HESS J1857+026 utilizando 15 años de datos de Fermi-LAT (300 MeV–2 TeV), datos de VERITAS (0.3–10 TeV) y nuevos datos de HAWC (1–316 TeV). El análisis procede a través de varias etapas:

1. Análisis de Datos:

   • Fermi-LAT: Se realizó un análisis de verosimilitud por bins (binned likelihood analysis) sobre datos de la clase SOURCE de Pass 8. La región de interés (ROI) se modeló utilizando el catálogo 4FGL-DR3 y plantillas difusas. Las pruebas de extensión espacial compararon plantillas de punto, disco radial y Gaussiana radial.
   • VERITAS: Se llevó a cabo un análisis de máxima verosimilitud en 3D sobre 30 horas de exposición utilizando el marco de trabajo Gammapy, modelando las fuentes como leyes de potencia con formas espaciales gaussianas simétricas.
   • HAWC: Se aplicó un análisis de máxima verosimilitud de plegado hacia adelante (forward-folding) al conjunto de datos Pass 5 (2860 días). El estudio probó múltiples plantillas espaciales (Gaussianas vs. basadas en difusión) y formas espectrales (ley de potencia, log-parábola, ley de potencia con corte exponencial) utilizando criterios de información Bayesiano y de Akaike (BIC/AIC) para seleccionar los mejores modelos de ajuste.

2. Modelado Radiativo:

   • Independiente del Tiempo: Las distribuciones espectrales de energía de banda ancha (SEDs) se ajustaron utilizando el paquete NAIMA para probar escenarios puramente hadrónicos (SNR) y lepto-hadrónicos (PWN/SNR).
   • Dependiente del Tiempo: Se empleó un modelo evolutivo (similar a Gelfand et al. 2009) para simular la expansión y el enfriamiento de la PWN dentro de una cáscara de SNR. Un procedimiento de Monte Carlo por cadenas de Markov (MCMC) restringió los parámetros físicos tales como la edad del sistema, la fuerza del campo magnético y los espectros de inyección de partículas.

3. Análisis de Difusión:

   • Se extrajeron perfiles de brillo superficial radiales de los tres instrumentos. Estos perfiles se ajustaron con un modelo analítico basado en la difusión (Ecuación 16) para restringir el coeficiente de difusión de las partículas, $D(E)$, y la supresión de la difusión con respecto al ISM.

Resultados Clave

• Características Espectrales y Espaciales:

  • Los datos de Fermi-LAT revelan una fuente extendida que se ajusta mejor mediante una Gaussiana radial con un radio de $0.38^\circ \pm 0.04^\circ$ y un índice espectral $\Gamma \approx 2.07$. Se observa una ruptura espectral, con un índice más suave ($\Gamma \approx 2.63$) por debajo de 5 GeV y uno más duro ($\Gamma \approx 1.99$) por encima de 5 GeV.
  • VERITAS detecta la fuente con una significancia de $15.2\sigma$, midiendo un índice espectral de $\Gamma = 2.48 \pm 0.06$ y un ancho Gaussiano de $0.236^\circ \pm 0.017^\circ$.
  • Los datos de HAWC muestran que la fuente se extiende hasta $\sim 37$ TeV. El mejor modelo espacial de ajuste favorece una plantilla de difusión sobre una simple Gaussiana para un campo magnético de $1 \mu$G, con una energía de corte $E_c \approx 14$ TeV.
  • El tamaño de la fuente disminuye con el aumento de la energía (de $0.38^\circ$ en energías de GeV a $0.26^\circ$ en energías de HAWC), lo que indica una morfología dependiente de la energía consistente con una PWN evolucionada.

• Origen Radiativo:

  • Lepto-hadrónico vs. Hadrónico: Aunque un modelo puramente hadrónico no puede descartarse por completo para el componente de baja energía ($<10$ GeV), el modelo lepto-hadrónico proporciona un ajuste superior a los datos de banda ancha. La falta de una detección de la cáscara de la SNR y la baja densidad de gas ambiental ($n_0 \lesssim 0.13$ cm$^{-3}$) derivada del modelo evolutivo desfavorecen un origen hadrónico dominante para la emisión de TeV.
  • Evolución de la PWN: El modelo evolutivo sugiere que el sistema es una PWN evolucionada con una edad de $\tau = [16, 21]$ kyr y una fuerza de campo magnético baja de $B = [0.4, 1.6] \mu$G. Estos parámetros explican la tenue emisión de rayos X (debido al enfriamiento por sincrotrón) y la morfología dependiente de la energía. El bajo campo magnético también respalda el predominio de la dispersión Compton inversa sobre los campos de fotones locales para la emisión de TeV.

• Difusión de Partículas:

  • Se encontró que el coeficiente de difusión alrededor de HESS J1857+026 está suprimido en relación con el promedio del ISM galáctico. A 50 TeV, el coeficiente de difusión es $D_{PWN} \sim 2.14 \times 10^{28}$ cm$^2$ s$^{-1}$ (para $B=1 \mu$G), que es aproximadamente dos órdenes de magnitud inferior al valor del ISM ($D_{ISM} \sim 1.2 \times 10^{30}$ cm$^2$ s$^{-1}$).
  • La supresión se encuentra principalmente en la normalización ($D_0$) más que en la dependencia de la rigidez (índice espectral $\delta$), que permanece consistente con los modelos de turbulencia de Kolmogorov ($\delta=1/3$) o Kraichnan ($\delta=1/2$). Esto respalda el escenario de turbulencia magnética localmente aumentada cerca de la fuente.

Significancia y Reivindicaciones
El artículo afirma que los datos combinados de múltiples longitudes de onda proporcionan una herramienta poderosa para restringir las propiedades de difusión de partículas en las proximidades de las PWN evolucionadas. La conclusión principal es que HESS J1857+026 es muy probablemente una PWN impulsada por PSR J1856+0245, donde la emisión de rayos $\gamma$ observada está dominada por electrones relativistas que experimentan dispersión Compton inversa.

El estudio destaca que el componente espectral de baja energía ($E < 10$ GeV) aún puede albergar una contribución hadrónica de una SNR huésped o cercana, o alternativamente, una población distinta de electrones dentro de la misma PWN. Sin embargo, la emisión de alta energía es firmemente leptónica. El factor de supresión de difusión derivado es consistente con otras PWN de TeV (por ejemplo, Geminga, Monogem), reforzando la hipótesis de que el transporte de partículas cerca de estas fuentes es significativamente dificultado por la turbulencia magnética local en comparación con el ISM general. Los autores señalan que, aunque se ha reportado un exceso de neutrinos con baja significancia, el gran límite superior de la región de emisión de neutrinos en comparación con la región de rayos $\gamma$ sugiere un origen complejo o no relacionado para cualquier señal potencial de neutrinos.