Mirage Sources and Large TeV Halo-Pulsar Offsets: Exploring the Parameter Space

Resumen Técnico: Fuentes Espejismo y Desplazamientos de Halo de Pulsar de Gran Escala en TeV

Planteamiento del Problema
Observaciones recientes realizadas por el Observatorio de Chorrores de Aire de Gran Altitud LHAASO (Large High Altitude Air Shower Observatory) y HAWC han identificado emisiones extendidas de rayos gamma en el rango de TeV ("halos de púlsar") alrededor de púlsares de edad media. Un desafío significativo para interpretar estas fuentes es la observación frecuente de grandes desplazamientos espaciales entre la posición del púlsar y el centroide de la emisión de TeV. Los modelos tradicionales que asumen una difusión isotrópica de pares electrón-positrón ($e^\pm$) en el medio interestelar (ISM) no pueden explicar naturalmente estos desplazamientos, ya que la difusión isotrópica solo es una aproximación válida cuando la escala física supera múltiples longitudes de coherencia magnética. En consecuencia, para electrones de ultra alta energía ($>100$ TeV) con tasas de enfriamiento rápidas, la escala de propagación puede ser menor que la longitud de coherencia del campo magnético turbulento del ISM. Por consiguiente, el transporte de partículas puede ser altamente asimétrico, siguiendo las líneas del campo magnético local en lugar de difundirse uniformemente. Este artículo investiga si dicha propagación asimétrica puede generar fuentes "espejismo" (mirage sources)—picos de emisión aparentes desplazados del acelerador—y desplazamientos espaciales significativos.

Metodología
Los autores emplean una simulación de partículas de prueba acelerada por GPU, basada en primeros principios, para rastrear las trayectorias de electrones de 100 TeV inyectados en un campo magnético turbulento.

• Generación del Campo Magnético: El campo magnético turbulento se modela como una superposición de cuatro campos independientes de espectro Kolmogorov ($n_{in} = -11/3$) generados mediante un método de espacio de Fourier. Este enfoque multiescala (utilizando espaciamientos de rejilla desde $4 \times 10^{-9}$ hasta $5 \times 10^{-4}$ kpc) asegura un rango inercial realista y campos sin divergencia ($\nabla \cdot B = 0$). El campo total incluye componentes tanto regulares ($B_r$) como turbulentos ($B_t$).
• Propagación de Partículas: Los electrones se propagan utilizando un algoritmo de empuje de Boris para resolver la ecuación de movimiento relativista bajo la fuerza de Lorentz. La simulación contabiliza las pérdidas de energía radiativa mediante radiación de sincrotrón y dispersión de Compton inverso (IC) sobre la Radiación de Fondo de Microondas Cósmicas (CMB) y el Campo de Radiación Interestelar (ISRF).
• Muestreo Estadístico: Para muestrear eficientemente el espectro de energía inicial ($dN/dE \propto E^{-\alpha} \exp(-E/E_{cut})$), los autores utilizan un enfoque de macropartículas ponderadas. Se inyectan 8,192 macropartículas, cada una representando un intervalo de energía logarítmica con un peso proporcional al conteo físico de electrones en dicho intervalo.
• Síntesis Observacional: Las distribuciones de electrones simuladas se convierten en mapas sintéticos de rayos gamma (banda de 25–40 TeV) mediante una matriz de dispersión de Compton inverso (IC). Estos mapas se conllevan con una Función de Respuesta de Punto (PSF) gaussiana de $\sigma_{PSF} = 0.43^\circ$ (imitando a LHAASO-KM2A) y se añaden a un modelo de ruido de fondo de Poisson.
• Identificación de Fuentes: Se aplica un marco de comparación de modelos bayesianos mediante Cadenas de Markov de Monte Carlo (MCMC). Los modelos con un número creciente de fuentes gaussianas se ajustan a los datos sintéticos. Una fuente se considera "identificada" si el Estadístico de Prueba (TS) aumenta más de 25 ($\Delta \log L > 12.5$) al añadir un componente.

Contribuciones Clave y Resultados
El estudio explora sistemáticamente el espacio de parámetros, variando la distancia del púlsar ($d$), el índice espectral de inyección de electrones ($\alpha$), la fuerza del campo magnético total ($B$), la longitud de coherencia de la turbulencia ($L_c$) y la relación entre la fuerza del campo regular y el turbulento ($B_r/B_t$).

1. Formación de Fuentes Espejismo: Las simulaciones confirman que la propagación asimétrica puede crear múltiples fuentes de rayos gamma distintas ("fuentes espejismo") a partir de un solo acelerador. Esto ocurre cuando los electrones se propagan a lo largo de las líneas del campo magnético que, debido a efectos de proyección, parecen originarse en una ubicación distante del punto de inyección.
2. Multiplicidad de Fuentes ($N$):
   • El número de fuentes identificadas está influenciado más fuertemente por la fuerza del campo magnético ($B$). Un $B$ más alto conduce a un enfriamiento más rápido, confinando a los electrones a escalas menores y reduciendo la probabilidad de formar múltiples fuentes distintas (menor $N$).
   • El número de fuentes tiene una correlación negativa con la distancia ($d$) y la longitud de coherencia ($L_c$), y una correlación positiva con $B^2$.
3. Separación de Fuentes y Relación de Flujo:
   • Separación: La separación máxima entre las fuentes identificadas está positivamente correlacionada con $L_c$ y negativamente correlacionada con $B^2$.
   • Relación de Flujo ($\chi$): La relación entre el flujo de la fuente principal y el de la fuente espejismo disminuye con el aumento de $B^2$. En campos magnéticos fuertes, la fuente espejismo puede ser a veces más brillante que la fuente principal.
   • Campos Regulares vs. Turbulentos: La distancia de separación exhibe una dependencia no monotónica de $B_r/B_t$. La separación máxima ocurre en $B_r/B_t \approx 1$, donde el campo regular proporciona guía direccional mientras que la turbulencia permite suficiente dispersión para crear desplazamientos. Los campos puramente turbulentos o los campos regulares dominantes resultan en separaciones menores.
4. Desplazamientos Espaciales:
   • Los desplazamientos entre el púlsar y la fuente identificada más cercana se maximizan cuando $B_r/B_t \approx 1$ y $L_c$ es grande.
   • En regiones con grandes longitudes de coherencia (por ejemplo, regiones entre brazos o el halo galáctico), los desplazamientos pueden alcanzar hasta 50 pc.

Significado y Reivindicaciones
El artículo afirma proporcionar una explicación cuantitativa para las fuentes desconcertantes de LHAASO donde los halos de TeV aparecen espacialmente desconectados de sus púlsares asociados. Específicamente, los autores sugieren que fuentes tales como 1LHAASO J0216+4232u, J0212+4254u y J0206+4302u podrían ser halos espejismo alimentados por un solo púlsar de milisegundo (PSR J0218+4232), en lugar de aceleradores distintos.

El estudio enfatiza que la aproximación de difusión isotrópica es insuficiente para modelar el transporte de partículas cerca de los aceleradores en el ISM, particularmente para leptones de alta energía donde los tiempos de enfriamiento son más cortos que el tiempo de coherencia magnética. Al demostrar que los efectos de proyección en campos magnéticos turbulentos pueden producir naturalmente grandes desplazamientos y múltiples morfologías de fuentes, este trabajo ofrece un mecanismo para reconciliar el cielo de TeV observado con los modelos estándar de nebulosas de viento de púlsar. Los autores señalan que, aunque su turbulencia sintética asume isotropía (a diferencia de la turbulencia de Alfvén anisotrópica en MHD), captura cualitativamente la naturaleza filamentosa a gran escala de la turbulencia del ISM responsable de estas emisiones asimétricas.