Westcott $g$ Factors Extended to Arbitrary Neutron Energy Spectra

Dit artikel introduceert een methode en open-source software om Westcott $g$-factoren te berekenen voor willekeurige neutronenspectra, in plaats van alleen voor Maxwelliaanse verdelingen, en past dit toe op gids- en koude-neutronenbundels van specifieke onderzoekreactoren.

De kern

De Westcott-g-factor: Een Reis door de Neutronenwereld

Stel je voor dat je een enorme, onzichtbare regen van deeltjes hebt die je wilt gebruiken om de samenstelling van een steen of metaal te analyseren. Dit zijn neutronen. Wetenschappers gebruiken deze neutronenregen om atomen in een monster te "vissen" (ze activeren), zodat ze kunnen meten wat erin zit. Dit heet Neutronenactivatie-analyse.

Maar hier is het probleem: niet alle neutronen zijn hetzelfde. Sommige rennen hard, sommigen lopen langzaam. En niet alle atomen reageren even goed op deze verschillende snelheden.

Deze paper is eigenlijk een update van de "vertaalwoordenlijst" die wetenschappers gebruiken om van de neutronenregen naar een nauwkeurig resultaat te komen. Hier is hoe het werkt, vertaald naar alledaagse taal:

1. Het Probleem: De "Regelmatige" vs. de "Eigenzinnige" Atomen

In de ideale wereld gedragen atomen zich als een trein: hoe langzamer de trein rijdt (de neutronen), hoe makkelijker het is om erin te stappen. Dit noemen ze de 1/v-wet. Als je dit weet, is het makkelijk om te berekenen hoeveel atomen je hebt "gevangen".

Maar sommige atomen (zoals Gadolinium of Lutetium) zijn eigenzinnige dansers. Ze hebben geen vaste regel. Soms dansen ze heel snel, soms heel langzaam, en soms hebben ze plotseling een "springtje" (een resonantie) waarbij ze ineens heel goed neutronen vangen, zelfs als die neutronen niet op de juiste snelheid zitten.

De Westcott-g-factor is de correctiefactor die wetenschappers gebruiken om deze eigenzinnige dansers te corrigeren. Het is een getal dat zegt: "Houd rekening met het feit dat deze atomen niet doen wat je verwacht."

2. Het Oude Mantra: Alles is een "Maxwelliaanse" Soep

Vroeger dachten wetenschappers dat ze de neutronenregen als een perfecte, warme soep konden beschouwen (een Maxwelliaanse verdeling). Ze dachten: "Als we de temperatuur van de soep weten, kunnen we de g-factor uit een boekje halen."

Het probleem? De echte neutronenregen in laboratoria (zoals in Hongarije en Duitsland) is geen soep. Het is meer als een geleidde stroom van water door een slang. De wanden van de slang reflecteren de langzamere deeltjes en veranderen de stroom. De oude "soep-formule" werkte niet goed voor deze specifieke situaties, vooral niet voor die eigenzinnige dansende atomen.

3. De Oplossing: Een Nieuwe Kaart voor Elke Reis

De auteurs van dit paper zeggen: "Stop met het raden van de temperatuur van de soep. Laten we gewoon kijken naar de echte stroom van neutronen die we hebben."

Ze hebben twee nieuwe, open-source softwaretools ontwikkeld (noem ze DeCE en WestcottFactors).

• DeCE is voor de oude garde die nog werkt met de traditionele data-indeling.
• WestcottFactors is de moderne, Python-gebaseerde tool die makkelijk te gebruiken is.

De analogie:
Stel je voor dat je een tol moet betalen om een brug over te steken.

• Vroeger: Je keek in een boekje: "Als het 20 graden is, kost het €5." (Onnauwkeurig als de brug anders is).
• Nu: Je gebruikt een app die de echte verkeersdrukte en de echte brug meet en zegt: "Voor deze brug, op dit moment, met deze auto's, kost het €4,73."

4. Waarom is dit belangrijk?

Als je de verkeerde correctiefactor gebruikt, kun je de hoeveelheid van een stof in een monster verkeerd inschatten. Voor sommige atomen kan het verschil 20% of meer zijn! Dat is alsof je denkt dat je 100 gram goud hebt, maar er eigenlijk maar 80 gram in zit.

De auteurs hebben getoond dat:

1. De oude methode (alleen kijken naar de temperatuur) vaak fouten maakt voor die eigenzinnige atomen.
2. De nieuwe methode (kijken naar de echte neutronenstroom) veel nauwkeuriger is.
3. Ze hun software gratis beschikbaar hebben gesteld, zodat iedereen (van studenten tot professionele laboranten) de juiste berekening kan doen zonder zelf de complexe wiskunde te hoeven begrijpen.

Samenvatting

Deze paper is een handboek-update. Het zegt tegen de wereld van de nucleaire wetenschap: "De wereld is complexer dan we dachten. Gebruik niet langer de oude, simpele formules voor alles. Gebruik onze nieuwe, slimme tools om de echte situatie te meten, zodat je resultaten kloppen."

Het is een eerlijke en transparante manier om ervoor te zorgen dat wanneer we de samenstelling van materialen meten, we de waarheid zien en niet een verzonnen versie ervan.

Technische samenvatting

Titel: Westcott 𝑔-factoren uitgebreid tot willekeurige neutronenspectra

Auteurs: D.A. Matters, A.M. Hurst en T. Kawano
Publicatie: Voorprent (Preprint) ingediend bij Elsevier (arXiv:2602.05995v2)

---

1. Het Probleem

In de Neutronenactivatie-analyse (NAA) en Prompt Gamma-ray Activation Analysis (PGAA) wordt de Westcott 𝑔-factor gebruikt om de impact van afwijkingen van de $1/v$-wet op de opbrengst van activatieproducten te kwantificeren. Voor de meeste kernen is het neutronenvangstwerkingsvermogen ($\sigma$) omgekeerd evenredig met de snelheid ($v$) van het invallende neutron ($1/v$-gedrag). Echter, bepaalde kernen (irreguliere kernen) vertonen neutronresonanties in het thermische en koude energiebereik (< 5 eV), wat leidt tot niet-$1/v$-gedrag.

Traditioneel worden Westcott 𝑔-factoren berekend onder de aanname van een Maxwelliaans neutronenfluxdistributie met een specifieke temperatuur ($T$). Dit is echter een vereenvoudiging. In moderne faciliteiten, zoals de Budapest Research Reactor (BRR) en de Forschungsreaktor München II (FRM II), worden gebruikgemaakt van geleide neutronenbundels (guided beams) en koude neutronenbronnen. De spectra van deze bundels wijken aanzienlijk af van een pure Maxwelliaanse verdeling. Het gebruik van tabellair vastgelegde 𝑔-factoren gebaseerd op Maxwelliaanse spectra kan leiden tot significante fouten bij de bepaling van effectieve werkingsvermogens voor irreguliere kernen, omdat de 𝑔-factor sterk afhankelijk is van de exacte vorm van het neutronenspectrum.

2. Methodologie

De auteurs hebben een nieuwe aanpak ontwikkeld om Westcott 𝑔-factoren te berekenen voor willekeurige neutronenspectra, niet beperkt tot Maxwelliaanse distributies.

• Definitie: De 𝑔-factor wordt gedefinieerd als de verhouding tussen het effectieve werkingsvermogen ($\sigma_{eff}$) en het werkingsvermogen bij thermische snelheid ($\sigma_0$), waarbij het effectieve werkingsvermogen wordt berekend door integratie over de werkelijke neutronenflux $\phi(v)$ en het werkingsvermogen $\sigma(v)$.
• Data Bronnen: Er zijn twee methoden gebruikt om de werkingsvermogensdata te verkrijgen:
  1. Irregulariteitsfunctie-methode: Een benadering die gebruikmaakt van de laagst-energetische resonantie (Lorentz-vorm) om het niet-$1/v$-gedrag te modelleren. Dit is de traditionele methode die vaak wordt gebruikt in handboeken.
  2. Directe Integratie-methode: Een nauwkeurigere methode waarbij de volledige neutronenvangstwerkingsvermogenskromme uit de Evaluated Nuclear Data File (ENDF/B-VIII.1) wordt gebruikt. Dit omvat alle resonanties en niet-resonantie bijdragen.
• Software Tools: Twee open-source tools zijn ontwikkeld of aangepast om deze berekeningen uit te voeren:
  • DeCE (Descriptive Correction of ENDF-6): Een C++-code die werkt met het legacy ENDF-6 formaat.
  • WestcottFactors: Een Python-toolkit die werkt met het modernere GNDS (Generalized Nuclear Database Structure) formaat en CSV-bestanden voor experimentele spectra. Deze tool maakt het mogelijk voor gebruikers om hun eigen gemeten neutronenspectra in te voeren.
• Validatie: De methoden zijn getest op een reeks irreguliere kernen (bijv. $^{149}$Sm, $^{151}$Eu, $^{157}$Gd, $^{176}$Lu) en vergeleken met bestaande tabellen (IAEA, PGAA Handbook) en met elkaar.

3. Belangrijkste Bijdragen

• Uitbreiding naar Willekeurige Spectra: De paper presenteert een robuuste methodologie om 𝑔-factoren te berekenen voor elk willekeurig neutronenspectrum, inclusief experimenteel gemeten spectra van geleide bundels.
• Open-Source Software: De ontwikkeling van de WestcottFactors pakket en de aangepaste DeCE fork, die beschikbaar zijn via GitHub en PyPI. Dit democratiseert de toegang tot nauwkeurige 𝑔-factor berekeningen voor NAA/PGAA-praktijk.
• Nieuwe Tabellen: Er zijn uitgebreide tabellen gegenereerd voor 𝑔-factoren bij verschillende Maxwelliaanse temperaturen en voor ideale geleide bundels, gebaseerd op de meest recente ENDF/B-VIII.1 data.
• Kritische Analyse van Bestaande Data: De auteurs tonen aan dat de traditionele "irregulariteitsfunctie"-methode (die slechts één resonantie gebruikt) onnauwkeurig kan zijn voor kernen met meerdere lage-energie resonanties.

4. Resultaten

• Vergelijking Methodes: Voor kernen met bijna-$1/v$-gedrag (zoals $^{83}$Kr) geven de irregulariteitsfunctie-methode en de directe integratie-methode vergelijkbare resultaten. Voor sterk irreguliere kernen (zoals $^{151}$Eu en $^{176}$Lu) vertoont de irregulariteitsfunctie-methode echter significante afwijkingen (tot wel 20% verschil bij hogere temperaturen) ten opzichte van de directe integratie over het volledige ENDF-spectrum.
• Invloed van het Spectrum: De berekeningen met de experimentele spectra van de BRR en FRM II tonen aan dat 𝑔-factoren voor irreguliere kernen sterk variëren afhankelijk van het gebruikte spectrum.
  • Bijvoorbeeld, voor $^{149}$Sm en $^{176}$Lu kunnen de verschillen tussen een benadering met een Maxwelliaans spectrum en het gebruik van het werkelijke gemeten spectrum 20% of meer bedragen.
  • De "ideale geleide bundel" benadering is vaak een betere schatting dan een Maxwelliaanse verdeling, maar het gebruik van het daadwerkelijk gemeten spectrum levert de meest accurate resultaten op.
• Software Consistentie: De resultaten van de C++ (DeCE) en Python (WestcottFactors) implementaties tonen een uitstekende overeenkomst (verschillen < 1%), wat de betrouwbaarheid van beide tools bevestigt.

5. Betekenis en Conclusie

Deze studie benadrukt dat het gebruik van standaard, tabellair vastgelegde Westcott 𝑔-factoren (die vaak gebaseerd zijn op ideale Maxwelliaanse spectra) kan leiden tot aanzienlijke systematische fouten in NAA en PGAA experimenten, vooral wanneer gebruik wordt gemaakt van geleide bundels of koude neutronenbronnen.

De belangrijkste conclusie is dat voor nauwkeurige analyses van irreguliere kernen, de 𝑔-factor moet worden berekend met het daadwerkelijk gemeten neutronenspectrum van de gebruikte faciliteit. De auteurs bevelen aan dat onderzoekers de open-source tools WestcottFactors of DeCE gebruiken om deze berekeningen uit te voeren in plaats van te vertrouwen op gepubliceerde tabellen. Dit verbetert de nauwkeurigheid van elementaire analyse en nucleaire data-evaluatie aanzienlijk.