A Boundary-Consistent Two-Zone Electron Kernel for Distant Pulsar Contributions to Positron Flux and Anisotropy

Resumen Técnico: Un núcleo de electrones de dos zonas con consistencia de frontera para las contribuciones de púlsares distantes al flujo de positrones y la anisotropía

Planteamiento del problema
Los espectros locales de electrones y positrones de rayos cósmicos sirven como sondas de precisión para aceleradores cercanos de alta energía, particularmente los púlsares. Púlsares maduros como Geminga y Monogem son candidatos naturales para el flujo de leptones en la escala de los TeV, pero una cuestión crítica permanece: ¿deben los flujos de leptones locales estar dominados por los púlsares más cercanos, o pueden los púlsares más distantes contribuir colectivamente? Esta indagación se complica por dos factores: (1) el rápido enfriamiento radiativo (sincrotrón e Compton inverso) limita el horizonte de propagación de las partículas de alta energía, y (2) los límites de anisotropía de Fermi-LAT y AMS-02 restringen los modelos donde un único origen cercano domina de forma demasiado amplia. Estudios previos han sugerido que las fuentes distantes pueden reproducir los electrones de baja energía, pero faltaba un marco semi-analítico estable para cuantificar el factor de conteo de fuentes geométrico en presencia de una difusión lenta cerca de la fuente. Además, los pocos cientos de pársecs interiores de la Galaxia (la Burbuja Local) presentan una discontinuidad física donde una densidad continua y suave de púlsares es sospechosa debido a la baja densidad de gas molecular y los impulsos natales de los púlsares.

Metodología
Los autores desarrollan una solución de serie semi-analítica para la propagación de electrones y positrones en un modelo de difusión de dos zonas esférico.

• Núcleo de dos zonas: El modelo define una región de difusión lenta de radio $R_1$ que rodea a la fuente, donde el coeficiente de difusión es $D_1(E)$, y una región interestelar estándar ($R_1 < r < R_2$) con coeficiente $D_2(E)$. Los coeficientes comparten la misma dependencia energética ($D \propto f_D(E)$), lo que permite soluciones de autofunción analíticas.
• Condiciones de frontera: La solución impone la continuidad de la densidad y el flujo difusivo en $R_1$ y una frontera absorbente en $R_2$ (establecida en 20 kpc). La formulación evita las cancelaciones oscilatorias inherentes a las evaluaciones de integrales directas de modelos de dos zonas, preservando los límites nítidos de enfriamiento radiativo.
• Enfriamiento: Las pérdidas radiativas (sincrotrón y Compton inverso de Klein–Nishina sobre campos de microondas de fondo cósmico, infrarrojo y óptico) se tratan mediante características de energía, vinculando la energía de inyección con la edad de la fuente.
• Modelado de población: El núcleo se aplica a una población continua de púlsares en el disco. Los autores imponen una cavidad vacía para $R < 0.1$ kpc para dar cuenta de la Burbuja Local y los movimientos propios de los púlsares.
• Análisis estadístico: Se realiza un análisis de Monte Carlo por cadenas de Markov (MCMC) con un "solucionador en el bucle" para ajustar el modelo a los datos de flujo de positrones de AMS-02, las mediciones de amplitud de dipolo de electrones totales de CALET y los límites superiores de anisotropía de Fermi-LAT. El ajuste incluye historias de inyección de pérdida de rotación continuas, un espectro de inyección de pares de ley de potencia rota y perfiles sobre amplitudes lineales no negativas para una base de positrones de tipo secundario y un fondo de ley de potencia suave.

Contribuciones clave

1. Solución de serie estable: El artículo proporciona una solución de serie semi-analítica numéricamente estable para el problema de difusión de dos zonas con enfriamiento radiativo. Este método es computacionalmente eficiente y evita las inestabilidades numéricas de los métodos de volumen finito de fuerza bruta o las integrales oscilatorias, particularmente en el régimen de distancia media.
2. Cuantificación geométrica: El estudio cuantifica cómo el "factor de concha" (conteo de fuentes frente a distancia) interactúa con los núcleos de propagación. Demuestra que mientras una distribución 3D uniforme está dominada por conchas a escala de kiloparsecs a 10–100 GeV, una geometría de disco ($h=0.2$ kpc) desplaza la contribución hacia adentro; sin embargo, las fuentes más allá de 1 kpc aún suministran el 37–47% del flujo de 10–100 GeV.
3. Papel del halo lento: Los autores aclaran que la inclusión de un halo de difusión lenta no mejora el flujo de las fuentes distantes en relación con un modelo de una sola zona; más bien, introduce una penalización de enfriamiento dependiente de la energía que reduce el flujo absoluto de todas las distancias, con una supresión más fuerte a energías más altas.
4. Restricciones de anisotropía: El trabajo integra las restricciones de anisotropía directamente en el ajuste de la población, utilizando las amplitudes de dipolo de CALET y los límites superiores de Fermi-LAT para probar la viabilidad de los parámetros de halo común.

Resultos

• Descomposición de flujo: En la población de púlsares continua ajustada (excluyendo $R < 0.1$ kpc), la contribución de las fuentes más allá de 0.3 kpc domina el anillo local de 0.1–0.3 kpc en todas las energías desde 10 GeV hasta la escala de los TeV. A 10 GeV, la contribución de $R > 0.3$ kpc es aproximadamente un orden de magnitud mayor que la del anillo local.
• Parámetros del halo: El análisis MCMC arroja una distribución posterior para el radio del halo lento $R_1$ centrada en $0.0709^{+0.0510}_{-0.0407}$ kpc ($\approx 71$ pc) y un coeficiente de difusión $D_1(100 \text{ TeV}) \approx 2.08 \times 10^{28}$ cm$^2$ s$^{-1}$. Este resultado es consistente con los estudios de morfología de halos de TeV de Geminga y Monogem.
• Compatibilidad de anisotropía: El modelo ajustado, con los parámetros de halo derivados, sigue siendo compatible con los límites actuales de anisotropía (CALET y Fermi-LAT). Las predicciones de la amplitud del dipolo caen dentro de los límites superiores del 95%.
• Dependencia de la energía: Mientras que las fuentes distantes dominan el flujo sub-TeV, el enfriamiento radiativo causa que la contribución se localice cerca de la escala de los TeV. A 1 TeV, aproximadamente el 90% del flujo en el modelo de disco proviene de dentro de 1 kpc.

Significado y afirmaciones
El artículo concluye que una contribución colectiva de púlsares distantes (más allá de la vecindad inmediata de la Burbuja Local) es un contribuyente plausible y dominante para el flujo de positrones y electrones por debajo de unos pocos cientos de GeV. El componente de púlsares ajustado no está dominado por las fuentes más cercanas, sino por una población que se extiende más allá de 0.3 kpc.

Sin embargo, los autores mantienen una postura modesta respecto a la unicidad de sus hallazgos:

• Los datos de flujo y de anisotropía por sí solos no determinan de forma única el tamaño del halo ($R_1$) o el coeficiente de difusión ($D_1$); existen degeneraciones con las formas de inyección, el modelado del fondo, las historias de liberación y la incompletitud de los catálogos.
• El tamaño del halo derivado "a escala de Geminga" es una restricción condicional basada en el modelo de población y el conjunto de datos específicos utilizados, no una medición definitiva de los halos de sub-difusión físicos.
• Para pasar de una prueba de consistencia a una determinación de la región de difusión lenta, se requieren observables adicionales, específicamente perfiles de brillo superficial de Compton inverso en el TeV de halos individuales, datos detallados de la población de púlsares (movimientos propios, sitios de nacimiento) y verosimilitudes completas de anisotropía vectorial.

El trabajo sirve como una demostración controlada por la propagación que permite probar el espacio de parámetros del halo común, proporcionando un núcleo estable para futuras inferencias de población, al tiempo que destaca la necesidad de entradas multi-mensajero para una solución completa.