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

Sintesi Tecnica: Un Nucleo di Elettroni a Due Zone con Coerenza di Confine per i Contributi di Pulsar Distanti al Flusso di Positroni e all'Anisotropia

Problema
Gli spettri locali di elettroni e positroni nei raggi cosmici servono come sonde di precisione per acceleratori di alta energia nelle vicinanze, in particolare per le pulsar. Mentre pulsar mature come Geminga e Monogem sono candidati naturali per il flusso leptonico su scala TeV, una questione critica rimane: il flusso leptonico locale deve essere dominato dalle pulsar più vicine, o possono le pulsar più distanti contribuire collettivamente? Questa indagine è complicata da due fattori: (1) il rapido raffreddamento radiativo (sincrotrone e Compton inverso) limita l'orizzonte di propagazione delle particelle ad alta energia, e (2) i limiti di anisotropia da Fermi-LAT e AMS-02 vincolano i modelli in cui una singola sorgente vicina domina troppo ampiamente. Studi precedenti hanno suggerito che le sorgenti distanti possono riprodurre gli elettroni a bassa energia, ma mancava un quadro semi-analitico stabile per quantificare il fattore geometrico del conteggio delle sorgenti in presenza di una diffusione lenta vicino alla sorgente. Inoltre, i primi centinaia di parsec della Galassia (la Bolla Locale) presentano una discontinuità fisica dove una densità continua di pulsar è sospetta a causa della bassa densità di gas molecolare e dei calci natali delle pulsar.

Metodologia
Gli autori sviluppano una soluzione in serie semi-analitica per la propagazione di elettroni e positroni in un modello di diffusione a due zone sferiche.

• Nucleo a due zone: Il modello definisce una regione a diffusione lenta di raggio $R_1$ che circonda la sorgente, dove il coefficiente di diffusione è $D_1(E)$, e una regione interstellare standard ($R_1 < r < R_2$) con coefficiente $D_2(E)$. I coefficienti condividono la stessa dipendenza energetica ($D \propto f_D(E)$), permettendo soluzioni analitiche degli autovettori.
• Condizioni al contorno: La soluzione impone la continuità della densità e del flusso diffusivo a $R_1$ e un confine assorbente a $R_2$ (impostato a 20 kpc). La formulazione evita le cancellazioni oscillatorie inerenti alle valutazioni integrali dirette dei modelli a due zone, preservando i confini netti del raffreddamento radiativo.
• Raffreddamento: Le perdite radiative (sincrotrone e Compton inverso Klein–Nishina sui campi CMB, infrarossi e ottici) sono trattate tramite caratteristiche energetiche, collegando l'energia di iniezione all'età della sorgente.
• Modellazione della Popolazione: Il nucleo è applicato a una popolazione continua di pulsar nel disco. Gli autori impongono una cavità vuota per $R < 0.1$ kpc per tenere conto della Bolla Locale e dei moti propri delle pulsar.
• Analisi Statistica: Viene eseguita un'analisi Monte Carlo a catena di Markov (MCMC) con "solver-in-the-loop" per adattare il modello ai dati del flusso di positroni di AMS-02, alle misure di anisotropia del dipolo degli elettroni totali di CALET e ai limiti superiori di anisotropia di Fermi-LAT. Il fit include storie di iniezione di spin-down continue, uno spettro di iniezione di coppie a legge di potenza spezzata e profili su ampiezze lineari non negative per una base di positroni simile a quella secondaria e un background a legge di potenza liscio.

Contributi Chiave

1. Soluzione in Serie Stabile: Il documento fornisce una soluzione in serie semi-analitica numericamente stabile per il problema della diffusione a due zone con raffreddamento radiativo. Questo metodo è computazionalmente efficiente ed evita le instabilità numeriche dei metodi a volume finito brute-force o degli integrali oscillanti, particolarmente nel regime di media distanza.
2. Quantificazione Geometrica: Lo studio quantifica come il "fattore di shell" (conteggio delle sorgenti rispetto alla distanza) interagisca con i kernel di propagazione. Dimostra che mentre una distribuzione 3D uniforme è dominata da gusci (shell) su scala kiloparsec a 10–100 GeV, una geometria a disco ($h=0.2$ kpc) sposta il contributo verso l'interno, eppure le sorgenti oltre 1 kpc forniscono ancora il 37–47% del flusso a 10–100 GeV.
3. Ruolo dell'Alone Lento: Gli autori chiariscono che l'inclusione di un alone a diffusione lenta non aumenta il flusso dalle sorgenti distanti rispetto a un modello a una zona; piuttosto, introduce una penalità di raffreddamento dipendente dall'energia che riduce il flusso assoluto da tutte le distanze, con una soppressione più forte alle energie più elevate.
4. Vincoli di Anisotropia: Il lavoro integra direttamente i vincoli di anisotropia nella popolazione del fit, utilizzando le ampiezze del dipolo di CALET e i limiti superiori di Fermi-LAT per testare la viabilità dei parametri dell'alone comune.

Risultati

• Decomposizione del Flusso: Nella popolazione di pulsar continua (escludendo $R < 0.1$ kpc), il contributo delle sorgenti oltre 0.3 kpc domina l'anello locale di 0.1–0.3 kpc a tutte le energie da 10 GeV al TeV. A 10 GeV, il contributo di $R > 0.3$ kpc è circa un ordine di grandezza superiore rispetto all'anello locale.
• Parametri dell'Alone: L'analisi MCMC fornisce una distribuzione a posteriori per il raggio dell'alone lento $R_1$ centrata a $0.0709^{+0.0510}_{-0.0407}$ kpc ($\approx 71$ pc) e un coefficiente di diffusione $D_1(100 \text{ TeV}) \approx 2.08 \times 10^{28}$ cm$^2$ s$^{-1}$. Questo risultato è coerente con gli studi di morfologia dell'alone TeV di Geminga e Monogem.
• Compatibilità con l'Anisotropia: Il modello adattato, con i parametri dell'alone derivati, rimane compatibile con gli attuali limiti di anisotropia (CALET e Fermi-LAT). Le previsioni dell'ampiezza del dipolo rientrano nei limiti superiori del 95%.
• Dipendenza dall'Energia: Sebbene le sorgenti distanti dominino il flusso sub-TeV, il raffreddamento radiativo causa la localizzazione del contributo vicino alla scala TeV. A 1 TeV, circa il 90% del flusso nel modello a disco proviene da entro 1 kpc.

Significato e Rivendicazioni
Il documento conclude che un contributo collettivo da parte di pulsar distanti (oltre l'immediato vicinato della Bolla Locale) è un contributore plausibile e dominante per il flusso di positroni ed elettroni al di sotto di alcune centinaia di GeV. La componente pulsar adattata non è dominata dalle sorgenti più vicine, ma da una popolazione che si estende oltre 0.3 kpc.

Tuttovi, gli autori mantengono una posizione modesta riguardo alla unicità delle loro scoperte:

• I dati sul flusso e sull'anisotropia da soli non determinano univocamente la dimensione dell'alone ($R_1$) o il coefficiente di diffusione ($D_1$); esistono degenerazioni con le forme di iniezione, la modellazione del background, le storie di rilascio e l'incompletezza dei cataloghi.
• La dimensione dell'alone derivata "su scala Geminga" è un vincolo condizionale basato sul modello di popolazione e sul set di dati specifico utilizzati, non una misurazione definitiva di aloni di sub-diffusione fisici.
• Per passare da un test di coerenza alla determinazione della regione a diffusione lenta, sono necessari ulteriori osservabili, specificamente profili di luminosità superficiale Compton inverso TeV di singoli aloni, dati dettagliati sulla popolazione di pulsar (moti propri, siti di nascita) e likelihood vettoriali complete dell'anisotropia.

Il lavoro funge da dimostrazione controllata dalla propagazione per testare lo spazio dei parametri dell'alone comune, fornendo un nucleo stabile per future inferenze di popolazione, evidenziando al contempo la necessità di input multi-messaggero per una soluzione completa.