The 'Forgotten' Neutrons: Implications for the Propagation of High-Energy Cosmic Rays in Magnetized Astrophysical and Cosmological Structures

Probleemstelling

Cosmische straling (CR), voornamelijk bestaande uit geladen protonen, wordt doorgaans gemodelleerd als deeltjes die magnetisch worden opgesloten in grote astrofysische structuren zoals sterrenstelsels, sterrenstelselclusters en kosmische filamenten. Traditionele propagatiemodellen gaan ervan uit dat deze deeltjes gedurende hun hele reis geladen blijven en zich voortbewegen via diffusie of stroming, beïnvloed door magnetische velden.

Echter, bij ultra-hoge energieën (UHE, $>10^{18}$ eV) kunnen hadronische interacties (zoals $pp$- en $p\gamma$-botsingen) protonen omzetten in neutronen. Neutronen zijn elektrisch neutraal en worden daarom niet beïnvloed door magnetische velden; ze bewegen zich ballistisch voort. Omdat neutronen een eindige levensduur hebben (ongeveer 10 minuten in rust), maar door relativistische tijddilatatie bij hoge energieën over enorme afstanden kunnen reizen voordat ze vervallen, vormen ze een potentieel "ontsnappingskanaal".

Het probleem is dat deze "geladen-neutrale schakeling" (charged-neutral switching) in de meeste modellen voor het transport tussen magnetische reservoirs wordt genegeerd. Dit leidt tot een onderschatting van de ontsnapping van ultra-hoge-energie CR's uit structuren waar ze anders sterk zouden worden geabsorbeerd of opgesloten.

Methodologie

De auteurs ontwikkelen een nieuw propagatiemodel dat de tijdelijke conversie van protonen naar neutronen expliciet behandelt als een stochastisch proces.

1. Transportvergelijking met Poisson-sprongen:
   In plaats van een gekoppeld systeem van vergelijkingen voor protonen en neutronen op te lossen, gebruiken de auteurs een stochastische differentiaalvergelijking (SDE).

   • Deterministisch deel: Beschrijft de diffusie en stroming van geladen protonen binnen het magnetische veld.
   • Stochastisch deel: Neutronproductie en verval worden gemodelleerd als een Poisson-"jump"-proces. Wanneer een proton een hadronische interactie ondergaat, wordt er met een bepaalde waarschijnlijkheid een neutron gegenereerd dat ballistisch beweegt.
   • De vergelijking wordt numeriek opgelost met een hybride schema dat een stijve ODE-oplosser (RADAU5) combineert met Monte Carlo-sampling voor de stochastische sprongen.

2. Omgevingen en Parameters:
   Het model wordt toegepast op vier representatieve omgevingen met verschillende schalen en magnetische eigenschappen:

   • Een normaal sterrenstelsel (Milky Way-achtig).
   • Een sterrenburst-stelsel (hoge dichtheid, sterke magnetische velden).
   • Een sterrenstelselcluster.
   • Een kosmisch filament.
     De simulaties lopen over energieën van $10^{12}$ eV tot $10^{21}$ eV en nemen rekening met interacties met materie ($pp$) en stralingsvelden ($p\gamma$, inclusief de CMB en sterrenlicht).

3. Berekening van Ontsnapping:
   De auteurs definiëren een "ontsnappingsfractie" ($f_{esc}$) en een "ontsnappingsfactor" ($A$) voor het neutronkanaal. Ze analyseren of neutronen de structuur kunnen verlaten voordat ze vervallen tot protonen, en hoe dit de totale baryonflux beïnvloedt.

Belangrijkste Bijdragen

• Formulering van een hybride transportmodel: De eerste expliciete implementatie van een stochastisch "jump"-proces voor neutronen in de context van CR-transport tussen grote kosmische structuren.
• Kwantificering van het "neutron-bypass"-effect: Het aantonen dat neutronen een ballistische route bieden die magnetische opsluiting en stralingsabsorptie kan omzeilen, zelfs in omgevingen waar geladen deeltjes gevangen zouden zitten.
• Redshift-afhankelijkheid: Het analyseren van hoe dit effect verandert bij hogere roodverschuivingen, waar de CMB-dichtheid hoger is en $p\gamma$-interacties frequenter worden.

Resultaten

De simulaties leiden tot de volgende cruciale bevindingen:

1. Kwalitatieve verschuiving in ontsnapping:

   • In sterrenstelsels (vooral sterrenburs) is het neutronkanaal minder dominant bij lage energieën omdat neutronen hier te snel vervallen om de structuur te verlaten.
   • In grootschalige structuren (clusters en filamenten) is het effect dramatisch. Bij ultra-hoge energieën ($>10^{20}$ eV) kunnen neutronen de hele structuur doorkruisen voordat ze vervallen. Dit herstelt een ontsnappingskanaal dat in standaardmodellen (alleen geladen deeltjes) volledig zou worden onderdrukt door magnetische opsluiting en $p\gamma$-absorptie.
   • In clusters kan de ontsnappingsfractie bij $z=1$ met een factor twee toenemen door het neutronkanaal.

2. Richtingsafhankelijkheid en Geometrie:

   • Filamenten: In tegenstelling tot het idee dat filamenten als "snelwegen" fungeren voor CR's langs hun as, blijkt uit de simulaties dat de geometrie (smalle diameter vs. lange lengte) ervoor zorgt dat CR's voornamelijk via de zijwanden ontsnappen. Filamenten gedragen zich meer als "lekke pijpen" dan als gesloten snelwegen, vooral bij hoge energieën.
   • Sterrenstelsels: Ontsnapping vindt voornamelijk plaats via de boven- en onderkappen van het schijfje, tenzij de energie zeer hoog is.

3. Invloed van Roodverschuiving:
   Bij hogere roodverschuivingen (bijv. $z=1$) wordt de CMB-dichtheid hoger, wat de $p\gamma$-interactietijd verkort en de ontsnapping van geladen protonen verder onderdrukt. Het neutronkanaal compenseert dit echter gedeeltelijk, waardoor er toch een significante flux van ultra-hoge-energie deeltjes de intergalactische ruimte in kan stromen.

4. Implicaties voor Calorimetrie:
   In sterrenstelsels blijft het calorimetrisch gedrag (energieopname) grotendeels ongewijzigd, maar in clusters en filamenten verandert het beeld: clusters zijn niet langer noodzakelijk perfecte "calorimeters" die alle energie omzetten in secundaire straling (gamma/neutrino's), omdat een deel van de energie als neutronen kan ontsnappen.

Significantie en Conclusie

De studie concludeert dat het negeren van neutronen leidt tot een fundamenteel verkeerd beeld van de distributie van ultra-hoge-energie kosmische straling in het universum.

• Herschikking van CR-reservoirs: Neutronen fungeren als een transportmiddel dat baryonen uit magnetisch opgesloten reservoirs (clusters) naar de omringende filamenten en vacuümruimtes (voids) kan vervoeren. Dit betekent dat filamenten en voids een veel rijkere populatie van UHE-CR's kunnen bevatten dan eerder werd gedacht.
• Bronnen van UHE-CR's: Het model suggereert dat sterrenstelselclusters, die eerder werden gezien als te "dicht" voor het ontsnappen van UHE-CR's, toch een plausibele bron kunnen zijn voor de op aarde waargenomen ultra-hoge-energie deeltjes, dankzij dit neutron-bypass-mechanisme.
• Observatie-implicaties: Dit effect zou zichtbaar moeten zijn in de anisotropie (richtingsafhankelijkheid) van de aankomstrichtingen van UHE-CR's. Correlaties met nabije grote structuren (zoals clusters) zouden zwakker of anders van aard zijn dan voorspeld door modellen zonder neutronen, vooral bij de hoogste energieën.

Kortom, de "vergeten neutronen" zijn een cruciaal schakelstuk in de kosmische web-ecologie van kosmische straling, die de grenzen tussen magnetische reservoirs en de lege ruimte effectief doorbreekt.