Stochastic analysis of ultra-high energy cosmic ray interactions

Dit artikel introduceert een analytische probabilistische beschrijving van de stochastische interacties van ultra-hoog-energetische kosmische stralingen met fotonvelden, gekoppeld aan Markov-sprongprocessen, om de beperkingen van huidige deterministische benaderingen en Monte Carlo-simulaties te overwinnen en zo een nauwkeurigere analyse van de samenstelling en het spectrum van kosmische stralingen mogelijk te maken.

De kern

De Kosmische Reis van een Sterrenstofdeeltje: Een Verhaal over Wiskunde en Toeval

Stel je voor dat het heelal een enorme, onzichtbare oceaan is. In deze oceaan zwemmen er gigantische deeltjes, de Ultra-Hoge Energie Kosmische Straling (UHECR). Deze deeltjes zijn vaak zware atoomkernen, zoals ijzer, die met bijna de lichtsnelheid door de ruimte vliegen. Ze komen van verre sterrenstelsels en zijn de snelste deeltjes die we kennen.

Maar hier is het probleem: als we deze deeltjes op Aarde opvangen, zijn ze vaak veranderd. Ze zijn niet meer de "oorspronkelijke" deeltjes die ze waren toen ze vertrokken. Onderweg botsen ze tegen onzichtbare muur van lichtdeeltjes (fotonen) aan, zoals de restwarmte van de Oerknal (de CMB).

Het oude verhaal: Een voorspelbare trein
Vroeger dachten wetenschappers dat dit proces heel voorspelbaar was. Ze dachten: "Als een ijzerkern een botsing heeft, verliest hij precies één deeltje, en dan nog één, en dan nog één." Het was alsof een trein op een spoor rijdt: elke stop is vastgelegd. Ze gebruikten simpele formules om te berekenen hoe ver zo'n trein kon komen voordat hij helemaal uit elkaar viel.

Het nieuwe verhaal: Een chaos van dobbelstenen
Dit nieuwe paper van Leonel Morejon en Karl-Heinz Kampert zegt: "Nee, het is niet zo simpel. Het is veel chaotischer."

Stel je voor dat de kosmische straling niet een trein is, maar een zwerm mieren die door een woud loopt.

• Soms botst een mier tegen een steen en valt hij uit elkaar in twee stukken.
• Soms valt hij in drie stukken uit elkaar.
• Soms valt hij in tien stukken uit elkaar.
• Soms valt hij helemaal niet uit elkaar, maar verliest hij alleen wat energie.

Elke keer dat er een botsing is, is het alsof je een dobbelspel speelt. Je weet niet precies wat er gebeurt, je weet alleen de kans dat er iets gebeurt. Dit noemen ze stochastisch (wiskundig toeval).

De uitdaging: Hoe bereken je het toeval?
Tot nu toe hadden wetenschappers twee manieren om dit te bestuderen:

1. De "Gok-methode" (Monte Carlo): Ze lieten een computer duizenden keren het dobbelspel spelen om een gemiddeld beeld te krijgen. Dit werkt goed, maar het is extreem traag en je krijgt geen diep inzicht in waarom het zo gaat. Het is alsof je duizenden keren een munt opgoopt om te zien of hij eerlijk is, in plaats van de wiskunde erachter te begrijpen.
2. De "Gemiddelde-methode": Ze negeerden het toeval en keken alleen naar het gemiddelde. Dit is snel, maar onnauwkeurig, vooral voor zware atomen zoals ijzer.

De oplossing: De wiskundige "Sieradenkist"
De auteurs van dit paper hebben een nieuwe, elegante manier bedacht. Ze hebben bewezen dat je dit hele chaotische proces kunt beschrijven met een wiskundig model dat precies werkt, zonder dat je hoeft te gokken.

Ze vergelijken het met een Sieradenkist:

• Stel je hebt een grote kist met verschillende vakjes (de verschillende soorten atoomkernen).
• Elke keer dat een deeltje een botsing heeft, springt het van het ene vakje naar het andere.
• De auteurs hebben een formule bedacht die precies beschrijft: "Als je begint met een zwaar ijzerdeeltje, wat is de kans dat je na 100 miljoen lichtjaar nog een ijzerdeeltje hebt? Wat is de kans dat het nu een koolstofatoom is? En wat is de kans dat het helemaal is verdwenen?"

Waarom is dit belangrijk?

1. Het is een GPS voor sterrenstof:
   Omdat we nu precies weten hoe snel en op welke manier deze deeltjes uit elkaar vallen, kunnen we beter berekenen hoe ver ze hebben gereisd. Het is alsof we een spoor van deeltjes vinden en kunnen zeggen: "Deze deeltjes komen niet van de buren, maar van een sterrenstelsel dat 500 miljoen lichtjaar verder weg ligt."

2. Het verklaren van de "Kleur" van het heelal:
   Als een zwaar deeltje uit elkaar valt, ontstaan er lichtere deeltjes. Soms zijn dit lichte atomen, soms zijn het protonen of neutronen. De auteurs laten zien dat de "samenstelling" van de straling die we zien, een code is. Als we de code kunnen kraken met hun nieuwe wiskunde, weten we wat er in de bron (bijvoorbeeld een explosie van een ster) is gebeurd.

3. Het kompas van de sterren:
   Deeltjes worden ook afgebogen door magnetische velden in het heelal. Als een deeltje uit elkaar valt, verandert zijn gewicht en lading, en dus verandert hoe sterk hij door het magnetische veld wordt afgebogen. De auteurs laten zien dat als je dit toeval meerekent, de richting waaruit het deeltje komt, heel anders kan zijn dan wat we dachten. Het is alsof je denkt dat een bal rechtaf rolt, maar door de wind en hobbels in het gras eindigt hij ergens anders.

Conclusie
Kortom: Dit paper geeft ons een nieuwe, super-nauwkeurige kaart voor de reis van kosmische deeltjes. In plaats van te gokken of te middelen, gebruiken ze een slimme wiskundige methode om het toeval te temmen. Hierdoor kunnen we beter begrijpen waar deze deeltjes vandaan komen, wat ze onderweg hebben meegemaakt, en wat er in de verre hoeken van het heelal gebeurt.

Het is alsof we van een ruwe schets van een landschap zijn gegaan naar een gedetailleerde 3D-kaart, compleet met de waarschijnlijke routes van elke reiziger.

Technische samenvatting

Titel: Stochastische analyse van interacties van ultra-hoge-energie kosmische straling

Auteurs: Leonel Morejon en Karl-Heinz Kampert
Instituut: Bergische Universiteit Wuppertal, Duitsland
Datum: 23 februari 2026 (conceptversie)

---

1. Het Probleem

Ultra-hoge-energie kosmische straling (UHECR) bestaat uit atoomkernen die door het heelal reizen en interageren met omringende fotonvelden (zoals de kosmische microgolfachtergrond - CMB, en de kosmische infraroodachtergrond - IRB). Deze interacties leiden tot energieverlies en fotodisintegratie (het afbreken van kernen).

De huidige uitdagingen zijn:

• Stochastische aard: Interacties zoals fotodisintegratie zijn fundamenteel probabilistisch. Een ijzerkern kan op verschillende manieren uiteenvallen, afhankelijk van de inelasticiteit en het aantal vrijgemaakte nucleonen.
• Beperkingen van bestaande methoden:
  • Deterministische benaderingen (ODE's): Deze behandelen energieverlies als continu en negeren fluctuaties. Ze zijn snel maar onnauwkeurig voor zware kernen en complexe cascades.
  • Monte Carlo-simulaties: Deze simuleren de stochastische aard nauwkeurig maar zijn computergewijs zeer duur. Ze bieden beperkt theoretisch inzicht in de correlatie tussen invoerparameters (zoals kernkruisdoorsneden) en uitkomsten, en convergeren soms traag.
• Gebrek aan analytisch kader: Er ontbreekt een formeel theoretisch raamwerk dat de stochastische interacties van UHECR's analytisch beschrijft, waardoor het moeilijk is om de connectie tussen de bronnen en de waarnemingen op Aarde exact te modelleren.

2. Methodologie

De auteurs presenteren een analytisch probabilistisch raamwerk dat UHECR-interacties beschrijft als Markov-sprongprocessen (Markov jump processes).

• Kernconcept: In plaats van het volgen van de gemiddelde dichtheid van kernen (zoals bij ODE's), modelleren ze de evolutie van de waarschijnlijkheidsverdeling van de toestanden (kernsoorten) als een functie van de afgelegde afstand (of "target thickness").
• Mathematische Formulering:
  • Het proces wordt beschreven door Kolmogorov's differentiaalvergelijking: $\frac{d}{dt}P_t = G P_t$, waarbij $P_t$ de matrix is van overgangswaarschijnlijkheden en $G$ de infinitesimale generator is (vergelijkbaar met het interactiematrix $\Lambda$).
  • De oplossing voor homogene gevallen is $P_t = e^{Gt}$, wat leidt tot gesloten vormen voor waarschijnlijkheidsdichtheden.
• Classificatie van Cascades: De auteurs onderscheiden drie types cascades op basis van de complexiteit van de interactiekanalen:
  1. Serieel en Regelmatig (RSeC): Eén kernsoort per massa, één kanaal per stap (bijv. één nucleon verlies per interactie). Dit leidt tot hypo-exponentiële verdelingen.
  2. Serieel en Irregulier (ISeC): Meerdere kernen per massa of variërende kruisdoorsneden die afwijken van de regelmaat (bijv. door spontane verval of specifieke kernstructuren).
  3. Concurrente Cascades (CoC): Meerdere kanalen per stap, waardoor een netwerk van kruisende paden ontstaat (bijv. emissie van 1, 2 of 6+ nucleonen). Dit vereist een volledige interactiematrix $\Lambda$.
• Inhomogeniteiten: Het model behandelt ook "Coherente Inhomogeniteiten" (CI), waarbij de interactiegraad verandert door kosmologische expansie (roodverschuiving), en "Dispersieve Inhomogeniteiten" (DI), waarbij energieverlies afhankelijk is van de specifieke kernsoort (bijv. synchrotronverlies), wat de Markov-eigenschap verstoort.

3. Belangrijkste Bijdragen

• Analytische Gesloten Vormen: Voor het eerst worden gesloten uitdrukkingen afgeleid voor de waarschijnlijkheidsverdelingen van de afstand tot volledige disintegratie en de samenstelling van secundaire deeltjes.
• Verbinding met Markov-processen: Het artikel legt een formele link tussen UHECR-propagatie en de theorie van Markov-sprongprocessen, wat nieuwe wiskundige tools (zoals "transformaties via beloningen") beschikbaar maakt voor het berekenen van grootheden zoals hoekafwijkingen.
• Open Source Implementatie: De methoden zijn geïmplementeerd in een open-source Python-pakket genaamd CRISP (Cosmic Ray Stochastic Interactions for Propagation), wat de reproduceerbaarheid en toepasbaarheid voor de gemeenschap vergroot.
• Reverse Propagation: Een methode om de oorspronkelijke samenstelling bij de bron te reconstrueren op basis van waarnemingen op Aarde, door het Markov-proces om te keren (quasi-stationaire verdeling).

4. Resultaten

• Horizon en Disintegratie: De berekende "horizon" (de afstand waarbij de kans op volledige disintegratie 99% is) verschilt aanzienlijk van traditionele schattingen gebaseerd op gemiddeld energieverlies (LEL). Vooral bij lagere Lorentz-factoren (waar IRB-interacties domineren) zijn de stochastische effecten significant.
• Massa-evolutie: Er wordt een "disciplined disintegration" (DD) effect bevestigd, maar de auteurs tonen aan dat bij concurrente cascades (CoC) de relatie tussen massa en disintegratie-afstand lineair is in plaats van logaritmisch (zoals bij seriële cascades).
• Invloed van Roodverschuiving: Cosmologische effecten (roodverschuiving van fotonvelden) worden nauwkeurig gemodelleerd via een effectieve "dikte" ($\delta$) van het medium, wat leidt tot kortere horizons dan bij homogene benaderingen.
• Aankomstrichting en Magnetische Afwijking: De stochastische disintegratie beïnvloedt de magnetische afwijking van kosmische straling. Omdat de rigiditeit ($R=E/Z$) verandert tijdens de cascade, is de hoekafwijking geen vaste waarde maar een verdeling. De auteurs tonen aan dat het negeren van disintegratie leidt tot grote fouten in de geschatte bronpositie (tot enkele graden verschil).
• Bronmodellen: Het model kan worden toegepast op bronnen (zoals GRB's) met tijdsafhankelijke injectie en complexe ontsnappingsmechanismen (bijv. diffuus vs. advecctie), wat leidt tot nauwkeurigere voorspellingen van het spectraal en samenstellingsprofiel.

5. Betekenis en Conclusie

Dit werk markeert een paradigmaverschuiving in de UHECR-fysica door de afhankelijkheid van Monte Carlo-simulaties voor probabilistische vragen te verminderen ten gunste van een efficiënter en inzichtelijker analytisch kader.

• Wetenschappelijke Impact: Het biedt een manier om de onzekerheden in kernkruisdoorsneden en fotonvelden systematisch te bestuderen zonder kostbare parameter-scan.
• Toekomstige Toepassingen: Het raamwerk maakt het mogelijk om de samenstelling van UHECR's op Aarde direct te koppelen aan de emissie-eigenschappen van bronnen, rekening houdend met de volledige stochastische geschiedenis van de deeltjes.
• Praktisch Nut: De beschikbaarheid van CRISP stelt onderzoekers in staat om complexe scenario's (zoals de oorsprong van de "Amaterasu"-deeltjes) te analyseren met een precisie die voorheen alleen met enorme rekenkracht mogelijk was.

Samenvattend biedt deze paper een fundamenteel nieuw instrument om de reis van ultra-hoge-energie deeltjes door het heelal te begrijpen, waarbij de inherente willekeur van kerninteracties niet als een hinderpaal, maar als een kwantificeerbaar en voorspelbaar aspect wordt behandeld.