On the Contribution of Local Sources to the Galactic Cosmic-Ray Spectrum: An Exact Series Solution for Two-Zone Diffusion

Resumen Técnico: Contribución de Fuentes Locales al Espectro de Rayos Cósmicos Galácticos mediante Difusión de Dos Zonas

Planteamiento del Problema
Las mediciones de los espectros de protones y helio de los rayos cósmicos (RC) galácticos por debajo del "rodilla" (knee) revelan desviaciones de las leyes de potencia simples, incluyendo estructuras en el rango de los multi-TeV. Aunque una explicación potencial involucra la contribución dominante de una o unas pocas fuentes cercanas, estudios previos que utilizan modelos de difusión homogénea (HDM) han concluido que tal escenario es estadísticamente improbable. Específicamente, las simulaciones de Monte Carlo sugieren que la probabilidad de que una fuente cercana altere significamente el espectro local es muy baja (~0.4%). Sin embargo, estas conclusiones dependen de la suposición de que los RC se difunden eficientemente inmediatamente después de dejar sus fuentes. Observaciones recientes de emisión de rayos gamma extendida alrededor de aceleradores galácticos (por ejemplo, Cygnus X, halos de púlsares) sugieren que el transporte cerca de las fuentes puede estar inhibido, caracterizado por un coeficiente de difusión suprimido. Este artículo investiga cómo cambia la probabilidad estadística de que una fuente local dominante sea relevante si los RC experimentan un transporte ineficiente en una región localizada que rodea sus sitios de aceleración.

Metodología
Para abordar los desafíos computacionales de modelar la propagación de RC con una región de difusión lenta localizada, los autores derivaron una solución analítica exacta en serie (función de Green) para un problema de difusión esférica de dos zonas.

• Geometría del Modelo: El modelo consiste en una región interna ($0 < r < r_1$) con un coeficiente de difusión $D_1$ y una región externa ($r_1 < r < r_2$) con un coeficiente mayor $D_2$. Se impone la continuidad de la densidad y del flujo difusivo en la interfaz $r_1$.
• Solución Analítica: Los autores derivaron expansiones de modos propios tanto para escenarios de inyección tipo estallido (burst-like) como para inyección continua. Este núcleo analítico permite el cálculo rápido de las distribuciones de partículas, haciendo factibles las simulaciones de Monte Carlo a gran escala sin la ineficiencia de resolver simultáneamente escalas espaciales pequeñas y grandes en códigos numéricos.
• Simulaciones de Monte Carlo: Utilizando la función de Green derivada, los autores generaron $10^5$ realizaciones de una población de remanentes de supernova (SNR) galácticos siguiendo el modelo "Jelly" (adoptado de Evoli et al.). Compararon la contribución de flujo de la fuente local más fuerte (definida como $d < 0.5$ kpc y edad $t < 1$ Myr) contra el flujo de fondo de todas las demás fuentes a 10 TeV.
• Análisis de Candidatos: El estudio también examinó candidatos cercanos catalogados específicos (púlsares y Loop I) para determinar si fuentes conocidas podrían reproducir las características espectrales observadas bajo diferentes suposiciones de transporte (HDM, Modelo de Difusión de Dos Zonas [2ZDM], y Propagación Dependiente del Espacio [SDP]).

Contribuciones Clave

1. Solución de Serie Exacta: La derivación de una función de Green analítica y rápida para un problema de difusión de dos zonas, permitiendo estudios de Monte Carlo eficientes de poblaciones de fuentes con inhibición de transporte local.
2. Cambio Cuantitativo de Probabilidad: Una reevaluación rigurosa de la probabilidad estadística de que una fuente local domine el espectro de RC cuando el transporte cercano a la fuente está inhibido.
3. Demostración de la Dependencia del Modelo: Una comparación sistemática que muestra cómo la contribución predicha de fuentes cercanas específicas (por ejemplo, Vela, Geminga, Monogem) varía drásticamente dependiendo del modelo de transporte asumido y los parámetros del espectro de inyección.

Resultados

• Aumento de la Probabilidad: En el HDM estándar, la probabilidad de que la fuente local más fuerte contribuya con un flujo comparable al de fondo a 10 TeV es aproximadamente del 0.4%. Bajo el 2ZDM, esta probabilidad aumenta al 1.7% para la inyección tipo estallido y al 2.2% para la inyección continua. Aunque esto representa un aumento de cuatro a cinco veces, la probabilidad absoluta sigue siendo baja (del orden de porcentajes).
• Mecanismo Físico: La región interna de difusión lenta actúa como un reservorio temporal, retrasando la escape de las partículas y redistribuyéndolas en el tiempo y la energía. Este efecto es más pronunciado para observadores fuera de la región interna ($\xi > 1$) en escalas de tiempo comparables al tiempo de difusión a través de la zona interna ($\tau \sim 1$).
• Análisis de Fuentes Candidatas: Al aplicar el 2ZDM a candidatos cercanos conocidos (usando púlsares como trazadores de SNR), una característica espectral de 10 TeV solo puede reproducirse bajo suposiciones adicionales específicas:
  • Espectro de Inyección Más Duro: Las fuentes locales deben poseer un índice de inyección significativamente más duro ($\gamma \approx 2.08$) comparado con la población de fondo ($\gamma \approx 2.34$). Si el índice de inyección local coincide con el de fondo, la contribución se vuelve insignificante.
  • Sensibilidad del Modelo de Transporte: El flujo predicho de una fuente dada varía fuertemente entre los modelos HDM, 2ZDM y SDP. Por ejemplo, bajo el SDP, Geminga y Monogem muestran contribuciones relativas diferentes en comparación con el 2ZDM, resaltando que la "importancia" de una fuente no es intrínseca, sino que depende del entorno magnético asumido y de las propiedades de transporte.

Significancia y Reivindicaciones
El artículo concluye que, si bien el transporte inhibido cerca de la fuente debilita significativamente el argumento estadístico contra la dominancia de una fuente local, no elimina la dificultad. La probabilidad de una fuente local dominante sigue siendo baja (~1–2%). Por consiguiente, la interpretación de las estructuras espectrales de multi-TeV como originadas en fuentes locales específicas es altamente dependiente del modelo.

Los autores sostienen que una fuente cercana puede explicar una característica espectral solo si se cumplen simultáneamente condiciones favorables: una edad/distancia de la fuente específica, un espectro de inyección local más duro y un entorno de transporte específico. Por lo tanto, identificar una fuente local como el origen de una característica espectral basándose únicamente en el ajuste espectral no es robusto. El artículo argumenta que establecer un origen local definitivo requiere restricciones independientes sobre la historia de inyección de la fuente, los mecanismos de transporte de partículas y la turbulencia interestelar local, idealmente mediante el modelado conjunto de espectros, composición, anisotropía y datos de rayos gamma.