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

Technische Samenvatting: Bijdrage van lokale bronnen aan het spectrum van de Galactische kosmische straling via twee-zones diffusie

Probleemstelling
Metingen van Galactische kosmische-stralingsspectra (protonen en helium) onder de "knee" onthullen afwijkingen van eenvoudige machtswetten, inclusen structuren in het multi-TeV bereik. Hoewel een potentiële verklaring een dominante bijdrage van één of enkele nabijgelegen bronnen betreft, concludeerden eerdere studies die gebruikmaakten van homogene diffusiemodellen (HDM) dat een dergelijk scenario statistisch gezien onwaarschijnlijk is. Specifiek suggereren Monte Carlo-simulaties dat de waarschijnlijkheid dat een nabijgelegen bron het lokale spectrum significant verandert zeer laag is (~0,4%). Deze conclusies rusten echter op de aanname dat kosmische straling (CR) efficiënt diffundeert onmiddellijk nadat deze de bronnen verlaat. Recente observaties van uitgebreide gammastraling rond Galactische versnellers (bijv. Cygnus X, pulsar-halo's) suggereren dat transport nabij bronnen kan worden geremd, gekenmerkt door een onderdrukte diffusiecoëfficiënt. Dit artikel onderzoekt hoe de statistische waarschijnlijkheid van een dominante lokale bron verandert als CR een inefficiënt transport ervaart in een gelokaliseerd gebied rondom hun acceleratieplaatsen.

Methodologie
Om de computationele uitdagingen aan te pakken bij het modelleren van CR-propagatie met een gelokaliseerd gebied met trage diffusie, hebben de auteurs een exacte analytische reeksoplossing (Green's functie) afgeleid voor een sferisch twee-zones diffusieprobleem.

• Modelgeometrie: Het model bestaat uit een binnenregio ($0 < r < r_1$) met een diffusiecoëfficiënt $D_1$ en een buitenregio ($r_1 < r < r_2$) met een grotere coëfficiënt $D_2$. Continuïteit van dichtheid en diffusieve flux wordt afgedwongen op de interface $r_1$.
• Analytische oplossing: De auteurs hebben eigenmodi-expansies afgeleid voor zowel burst-achtige als continue injectiescenario's. Deze analytische kern maakt snelle berekeningen van deeltjesdistributies mogelijk, wat grootschalige Monte Carlo-simulaties haalbaar maakt zonder de inefficiëntie van het gelijktijdig resolveren van kleine en grote ruimtelijke schalen in numerieke codes.
• Monte Carlo-simulaties: Met behulp van de afgeleide Green's functie genereerden de auteurs $10^5$ realisaties van een populatie van Galactische supernova-restanten (SNR) volgens het "Jelly"-model (overgenomen van Evoli et al.). Ze vergeleken de fluxbijdrage van de sterkste lokale bron (gedefinieerd als $d < 0,5$ kpc en leeftijd $t < 1$ Myr) tegen de achtergrondflux van alle andere bronnen bij 10 TeV.
• Kandidaatanalyse: De studie onderzocht ook specifieke gecatalogiseerde nabije kandidaten (pulsars en Loop I) om te bepalen of bekende bronnen de geobserveerde spectrale kenmerken kunnen reproduceren onder verschillende transportaannames (HDM, Twee-zones Diffusiemodel [2ZDM], en Ruimtelijk Afhankelijke Propagatie [SDP]).

Belangrijkste Bijdragen

1. Exacte Reeksoplossing: De afleiding van een snelle, analytische Green's functie voor een twee-zones diffusieprobleem, wat efficiënte Monte Carlo-studies van bronpopulaties met lokale transportremming mogelijk maakt.
2. Kwantitatieve Waarschijnlijkheidsverschuiving: Een rigoureuze herwaardering van de statistische waarschijnlijkheid dat een lokale bron het CR-spectrum domineert wanneer transport nabij de bron wordt geremd.
3. Demonstratie van Modelafhankelijkheid: Een systematische vergelijking die aantoont hoe de voorspelde bijdrage van specifieke nabije bronnen (bijv. Vela, Geminga, Monogem) drastisch varieert afhankelijk van het aangenomen transportmodel en de injectiespectrumparameters.

Resultaten

• Waarschijnlijkheidsverhoging: In de standaard HDM is de waarschijnlijkheid dat de sterkste lokale bron een flux bijdragen die vergelijkbaar is met de achtergrond bij 10 TeV ongeveer 0,4%. Onder de 2ZDM stijgt deze waarschijnlijkheid naar 1,7% voor burst-achtige injectie en 2,2% voor continue injectie. Hoewel dit een vier- tot vijfmaal verhoging vertegenwoordigt, blijft de absolute waarschijnlijkheid laag (orde van grootte procenten).
• Fysiek Mechanisme: De binnenste regio met trage diffusie fungeert als een tijdelijk reservoir, dat het ontsnappen van deeltjes vertraagt en ze in tijd en energie herverdeelt. Dit effect is het meest uitgesproken voor waarnemers buiten de binnenregio ($\xi > 1$) op tijdschalen die vergelijkbaar zijn met de diffusietijd over de binnenregio ($\tau \sim 1$).
• Analyse van Kandidaatbronnen: Bij het toepassen van de 2ZDM op bekende nabije kandidaten (waarbij pulsars als SNR-tracers worden gebruikt), kan een 10 TeV spectrale kenmerk alleen worden gereproduceerd onder specifieke aanvullende aannames:
  • Hardere Injectiespectrum: De lokale bronnen moeten een aanzienlijk harder injectie-index bezitten ($\gamma \approx 2,08$) vergeleken met de achtergrondpopulatie ($\gamma \approx 2,34$). Als de lokale injectie-index overeenkomt met de achtergrond, wordt de bijdrage verwaarloosbaar.
  • Gevoeligheid van het Transportmodel: De voorspelde flux van een gegeven bron varieert sterk tussen HDM, 2ZDM en SDP-modellen. Bijvoorbeeld, onder SDP vertonen Geminga en Monogem andere relatieve bijdragen vergeleken met 2ZDM, wat benadrukt dat de "belangrijkheid" van een bron niet intrinsiek is, maar afhankelijk is van de aangenomen magnetische omgeving en transporteigenschappen.

Betekenis en Claims
Het artikel concludeert dat hoewel geremde transport nabij de bron de statistische argumentatie tegen lokale bron-dominantie aanzienlijk verzwakt, het de moeilijkheid niet wegneemt. De waarschijnlijkheid van een dominante lokale bron blijft laag (~1–2%). Bijgevolg is de interpretatie van multi-TeV spectrale structuren als afkomstig van specifieke lokale bronnen sterk modelafhankelijk.

De auteurs stellen dat een nabije bron een spectrale kenmerk alleen kan verklaren als gunstige voorwaarden simultaan worden vervuld: een specifieke bronleeftijd/afstand, een hard lokaal injectiespectrum en een specifieke transportomgeving. Daarom is het identificeren van een lokale bron als de oorsprong van een spectrale kenmerk op basis van enkel spectrale fitting niet robuust. Het artikel betoogt dat het vaststellen van een definitieve lokale oorsprong een onafhankelijke beperking vereist van de broninjectiegeschiedenis, deeltjes-transportmechanismen en lokale interstellaire turbulentie, bij voorkeur door gezamenlijke modellering van spectra, compositie, anisotropie en gamma-gegevens.