Empirical formula for total inelastic cross-section of proton-nucleus scattering

Dit artikel presenteert een universele empirische formule voor de totale inelastische proton-kernverstrooiingscross-sectie over een breed energiebereik van 15 MeV tot 1 TeV, die is gecalibreerd op experimentele data en gevalideerd tegen bestaande modellen en GEANT4-simulaties voor een scala aan elementen.

De kern

Stel je voor dat je een enorme, chaotische danszaal binnenloopt. In deze zaal zijn er twee soorten gasten: protonen (kleine, snelle dansers) en atoomkernen (grote, zware dansers van verschillende maten, van een kleine muis tot een olifant).

Wanneer een kleine proton-danser tegen een grote kern-danser botst, kan er van alles gebeuren. Soms stoten ze zachtjes aan elkaar en stuiteren ze terug (dit noemen we elastisch botsen). Maar vaak is de botsing zo heftig dat er nieuwe deeltjes ontstaan, de kern trilt, of er stukjes van de kern afvliegen. Dit noemen we inelastisch botsen.

De wetenschappers in dit artikel (Hemant Kumar, Tanmay Maji en collega's) hebben een nieuw recept bedacht om te voorspellen hoe vaak deze "ruwe" botsingen gebeuren. Ze noemen dit de "totale inelastische doorsnede". In de praktijk is dit een maatstaf voor de kans dat er iets interessants gebeurt tijdens een botsing.

Hier is hoe hun nieuwe recept werkt, vertaald naar alledaagse taal:

1. Het probleem met oude recepten

Vroeger hadden wetenschappers al enkele recepten om deze botsingen te voorspellen.

• De oude simpele regels: Sommige regels werkten goed voor lichte kernen (zoals koolstof), maar faalden volledig voor zware kernen (zoals uranium).
• De complexe regels: Andere regels waren zo ingewikkeld als een recept voor een 10-gangenmenu met 50 ingrediënten. Ze werkten soms goed, maar waren onhandig om te gebruiken en faalden bij lichte kernen.

De auteurs zeggen: "Waarom niet één recept dat voor iedereen werkt, van de kleinste muis tot de zwaarste olifant, en dat bovendien makkelijk te berekenen is?"

2. Het nieuwe recept: Drie ingrediënten

Hun nieuwe formule is als een slimme cocktail die bestaat uit drie delen, afhankelijk van hoe hard de protonen dansen (hun energie):

• De Basis (De Hoge Energie):
  Als de protonen extreem snel zijn (zoals een Formule 1-auto), is het niet meer belangrijk hoe de kern eruitziet, maar alleen hoe groot hij is. De formule kijkt dan alleen naar het gewicht van de kern (het aantal deeltjes erin) en berekent een basiswaarde. Dit is als zeggen: "Hoe groter de dansvloer, hoe meer kans er is dat iemand erop stapt."
• De Trilling (De Lage Energie):
  Bij lagere snelheden (zoals een wandelaar) is het gedrag chaotischer. De kans op een botsing gaat omhoog en omlaag, alsof er een onzichtbare trilling in de zaal zit. De formule voegt hier een "golf" aan toe (een wiskundige sinus-golf) om deze pieken en dalen na te bootsen.
• De Magnetische Muur (De Coulomb-barrière):
  Dit is het slimste deel. Protonen zijn positief geladen, en de kern is ook positief. Ze stoten elkaar af, net als twee noordpolen van magneetjes. Bij lage snelheden is deze afstoting heel sterk. De formule heeft een speciale "rem" toegevoegd die rekening houdt met hoe groot de kern is. Voor kleine kernen werkt dit anders dan voor zware kernen. Het zorgt ervoor dat het recept precies de juiste piek voorspelt bij lage snelheden.

3. De Test: Van Deuterium tot Uranium

De auteurs hebben hun nieuwe recept getest tegen de echte dansfeesten die in laboratoria zijn opgenomen. Ze hebben gekeken naar 33 verschillende soorten kernen, variërend van heel licht (zoals Deuterium, een lichte water-isotoop) tot heel zwaar (zoals Uranium).

• Het resultaat: Hun nieuwe formule bleek over het algemeen nauwkeuriger te zijn dan de oude recepten.
• De vergelijking: Ze hebben hun formule ook vergeleken met GEANT4. Stel je GEANT4 voor als een superkrachtige computersimulatie die alles tot in detail berekent (alsof je elke danser individueel scant). Hun nieuwe, simpele formule kwam verrassend dicht in de buurt van deze zware computer, maar dan veel sneller en eenvoudiger.

4. Waarom is dit nuttig?

Waarom moeten we hierover weten? Omdat dit soort botsingen overal om ons heen gebeurt en cruciaal is voor:

• Kankerbehandeling: Bij protonentherapie moet je precies weten hoe protonen door weefsel (kernen) reizen om kankercellen te doden zonder gezond weefsel te beschadigen.
• Ruimtevaart: Straling in de ruimte bestaat uit protonen die tegen de wanden van een ruimteschip (en de astronauten erin) botsen. Een goede formule helpt bij het ontwerpen van betere schilden.
• Cosmische straling: Het helpt ons te begrijpen hoe lang kosmische deeltjes door het heelal reizen voordat ze ergens tegenaan botsen.

Conclusie

Kortom, deze wetenschappers hebben een universele, simpele en accurate "recept" bedacht om te voorspellen wat er gebeurt als protonen tegen atoomkernen botsen. Het werkt goed voor lichte én zware kernen, voor lage én hoge snelheden, en het is veel makkelijker te gebruiken dan de oude, ingewikkelde modellen. Het is alsof ze een universele vertaler hebben gevonden die zowel de fluisterende muis als het brullende olifant perfect begrijpt.

Technische samenvatting

Titel: Empirische formule voor de totale inelastische doorsnede van proton-kernverstrooiing

Auteurs: Hemant Kumar, Tanmay Maji, Deepa Gupta en Ashavani Kumar.
Publicatiedatum: 25 februari 2026 (voorgesteld).

---

1. Probleemstelling

De interactie tussen protonen en atoomkernen is fundamenteel voor het begrijpen van hadronische reacties in een nucleair medium. De totale inelastische doorsnede (het totaal van alle niet-elastische processen zoals deeltjesproductie en nucleaire excitatie) is een cruciale grootheid voor toepassingen in medische stralingstherapie, stralingsbescherming in de ruimte, detectorontwikkeling en kosmische straling.

Hoewel er diverse empirische modellen bestaan (zoals die van Letaw et al., Shen, Tripathi et al. en Nakano et al.), vertonen deze beperkingen:

• Ze zijn vaak niet uniform nauwkeurig over het volledige energiebereik (van 15 MeV tot 1 TeV).
• Ze presteren vaak slecht bij lichte kernen (zoals Deuterium en Helium) of bij zware elementen (zoals Uranium).
• Sommige modellen vereisen specifieke aanpassingen per doelwit of hebben complexe wiskundige formuleringen met veel parameters.
• Er is een gebrek aan betrouwbare data voor niet-gemeten kernen, wat voorspellingen noodzakelijk maakt.

Het doel van dit onderzoek is het ontwikkelen van een eenvoudige, universele en nauwkeurige empirische formule die geldig is voor een breed scala aan doelwitten (van Deuterium tot Uranium) en energieën (15 MeV tot 1 TeV).

2. Methodologie

De auteurs stellen een nieuw model voor dat de totale inelastische doorsnede ($\sigma_{inel}$) factoriseert in een hoog-energetisch en een laag-energetisch regime. Het model is geïnspireerd op het werk van Letaw et al., maar bevat geactualiseerde experimentele data en fysieke correcties.

A. Hoog-energetisch regime (E > 5 GeV)

In dit regime wordt de interactie gedomineerd door geometrische overwegingen en is de doorsnede bijna energie-onafhankelijk (verzadiging). De formule voor de hoog-energetische bijdrage $\sigma_h(A)$ is:
$$\sigma_h(A) = a A^b [1 - c \cos(d \sqrt{A}) + e A^f]$$

• De term $A^b$ vertegenwoordigt de geometrische benadering.
• De cosinus-term houdt rekening met nucleaire structuur-effecten.
• De term $A^f$ is specifiek toegevoegd om trends bij zware elementen te vangen.
• De parameters ($a, b, c, d, e, f$) zijn bepaald via een $\chi^2$-fitting op experimentele data van Bobchenko et al. (5-9 GeV).

B. Laag-energetisch regime (E < 2 GeV)

Bij lagere energieën is de doorsnede sterk energie-afhankelijk en vertoont het een oscillatief gedrag met een scherpe piek onder de 200 MeV. De totale doorsnede wordt gegeven door:
$$\sigma_{inel}(Z, A, E) = \sigma_h(A) \cdot X(E) \cdot Y_{coul}(Z, E)$$

• $X(E)$: Een energieafhankelijke correctiefunctie die een exponentieel afnemende term en een gemoduleerde sinus-term bevat om de piek en het verval bij lagere energieën te modelleren.
• $Y_{coul}(Z, E)$: Een Coulomb-correctiefunctie die rekening houdt met de afstotende Coulomb-barrière tussen het proton en de kern bij lage energieën. Deze term hangt af van het atoomnummer $Z$.
• De parameters voor deze laag-energetische termen zijn bepaald door fitting op uitgebreide experimentele data voor Aluminium (Al) en Koolstof (C) met behulp van machine learning-methoden.

3. Belangrijkste Resultaten

Nauwkeurigheid en Vergelijking

Het model is getest tegen experimentele data voor 33 verschillende kernen, variërend van lichte elementen (Deuterium, Helium) tot zware elementen (Uranium, Lood).

• Universele Toepasbaarheid: Het model presteert consistent goed voor zowel lichte als zware kernen, in tegenstelling tot bestaande modellen die vaak falen bij de ene of de andere categorie.
• Vergelijking met Bestaande Modellen: De auteurs vergelijken hun resultaten met de modellen van Letaw, Shen, Tripathi, Nakano en GEANT4-simulaties (met de FTFP_BERT physics list).
  • Bestaande modellen (zoals Letaw) vertonen vaak grote afwijkingen bij lage energieën voor zware kernen of voorspellen onfysische waarden voor lichte kernen.
  • Het voorgestelde model levert de kleinste relatieve fouten (RE) voor de meeste van de 33 kernen, zoals weergegeven in Tabel V van het artikel.
• Lage Energie: Het model slaagt erin de scherpe piek bij energieën < 200 MeV nauwkeurig te reproduceren, wat een zwak punt is bij veel concurrenten.

Toepassingen

1. Levensduur van Kosmische Straling:
   Het model werd gebruikt om de interactie-snelheid en de gemiddelde vrije weg van kosmische straling in het interstellaire medium (ISM) te berekenen. De berekende levensduur toont een verbetering van ongeveer 20% ten opzichte van de voorspellingen van de modellen van Shen en Letaw, voornamelijk door de betere modellering van de lage-energiepiek bij lichte kernen (H en He).

2. Antiproton-productie (Mu2e Experiment):
   Voor het Mu2e-experiment bij Fermilab werd de productie van antiprotonen geschat. Het model voorspelt een inelastische doorsnede van 1689 mb voor Aluminium, wat resulteert in een antiproton-productie die zeer dicht bij de Letaw-waarde ligt (1% afwijking) maar ongeveer 11% afwijkt van de MCNP-simulatie. Dit onderstreept de bruikbaarheid voor het schatten van achtergronden in gevoelige experimenten.

4. Bijdragen en Significantie

• Eenvoud en Efficiëntie: Het model is computatievriendelijk en vereist minder parameters dan complexe modellen zoals dat van Nakano et al., terwijl het vergelijkbare of betere nauwkeurigheid biedt.
• Breed Energiebereik: Het is het eerste empirische model dat expliciet is ontworpen om het volledige bereik van 15 MeV tot 1 TeV te dekken met één samenhangende formulering.
• Voorspellend Vermogen: Omdat experimentele data ontbreekt voor veel kernen en energieën, biedt deze formule een betrouwbaar hulpmiddel om niet-gemeten nucleaire doorsnedes te voorspellen. Dit is essentieel voor het ontwerp van nieuwe versnellers, stralingsbescherming en astrofysische modellen.
• Fysische Inzicht: Door de factorisatie in geometrische, energie-afhankelijke en Coulomb-termen, behoudt het model een duidelijk fysiek verband met de onderliggende interactiemechanismen, in plaats van puur een "black box" fit te zijn.

Conclusie

De auteurs hebben een robuuste en universele empirische formule ontwikkeld voor de totale inelastische proton-kerndoorsnede. Door de factorisatie van hoog- en laag-energetische regimes en de integratie van een specifieke Coulomb-correctie, overtreft het model bestaande benaderingen in nauwkeurigheid en consistentie over een breed spectrum van doelwitten en energieën. Dit maakt het een waardevol instrument voor zowel fundamenteel nucleair onderzoek als toegepaste toepassingen in stralingsfysica en kosmische straling.