Shell effects and the neutron emission within the multi-dimensional Langevin model for fission

Dit artikel beschrijft een multidimensionaal Langevin-model dat de tijdsafhankelijke vormevolutie van een splijtingskern simuleert om de emissie van neutronen te modelleren en zo de pre-splijtingsneutrone multipliciteit, massaverdeling en energiedistributie te berekenen en te vergelijken met experimentele data.

De kern

De Grote Splitsing: Hoe een atoomkern barst en deeltjes spuugt

Stel je een atoomkern voor als een enorme, trillende bal van klei. Deze bal is zo volgepropt met deeltjes (neutronen en protonen) dat hij eigenlijk niet stabiel is. Hij wil liever in twee kleinere ballen uiteenvallen. Dit proces noemen we kernsplitsing (of fissie).

De auteurs van dit artikel, een team van fysici uit Oekraïne en Japan, hebben een heel slim computerprogramma geschreven om te simuleren hoe deze "kleibal" precies uit elkaar valt. Maar ze hebben iets toegevoegd wat in eerdere modellen vaak werd vergeten: het spugen van neutronen tijdens het proces.

Hier is hoe ze het hebben aangepakt, vertaald in alledaagse taal:

1. De Langevin-beweging: Een dronken wandelaar

Het hart van hun onderzoek is iets dat ze de "Langevin-vergelijking" noemen.

• De Analogie: Stel je voor dat de atoomkern een dronken wandelaar is die over een heuvelachtig landschap loopt.
  • Het landschap is de energie: er zijn diepe dalen (waar de kern veilig is) en hoge heuvels (de barrière die hij moet overwinnen om te splijten).
  • De wandelaar heeft een stok (de wrijving) die zijn beweging vertraagt.
  • Er waait een willekeurige wind (de random force) die hem soms naar links en soms naar rechts duwt.
• Het doel: De computer berekent stap voor stap hoe deze wandelaar over het landschap hobbelt tot hij de top van de heuvel bereikt en in het dal aan de andere kant valt. Dat moment van vallen is de splitsing.

2. Het nieuwe geheim: De "neutronen-schroef"

In het verleden keken wetenschappers alleen naar hoe de wandelaar over de heuvel liep. Maar in werkelijkheid gebeurt er iets belangrijks onderweg: de wandelaar is het heet en begint zweetdruppels (neutronen) af te geven.

• Het probleem: Als de kern neutronen kwijtraakt, koelt hij af. Hij verliest energie. Als hij te veel energie verliest, kan hij misschien de heuvel niet meer over en valt hij terug in het dal (geen splitsing).
• De oplossing: De auteurs hebben hun computerprogramma zo aangepast dat het bij elke stap van de wandeling checkt: "Gaat er nu een neutron weg?"
  • Als het antwoord ja is, trekken ze direct de energie van de wandelaar af.
  • Ze veranderen de kaart (het landschap) een klein beetje, omdat de kern nu één neutron minder heeft.
  • Dan laten ze de wandelaar verder lopen op de nieuwe kaart.

3. Waarom is dit belangrijk?

Dit klinkt als een klein detail, maar het heeft grote gevolgen, net als het afwerpen van gewicht tijdens een bergbeklimming.

• De vorm van de brokken: Als de kern neutronen kwijtraakt voordat hij uiteenvalt, verandert dit hoe de twee nieuwe stukken (de splijtingsfragmenten) eruitzien.
  • Vergelijking: Stel je voor dat je een taart in tweeën snijdt. Als je er eerst wat room van afkrabt (neutronen), is de taart lichter en valt hij anders uiteen. De auteurs ontdekten dat door rekening te houden met deze "afgekrabde room", hun berekeningen veel beter overeenkwamen met wat mensen in het echt in het lab zien.
• Wanneer gebeurt het? Ze ontdekten iets verrassends:
  • Bij lage energie (een "koudere" wandelaar) worden de neutronen pas gespuugd als de wandelaar al over de heuveltop is.
  • Bij hoge energie (een "heetere" wandelaar) beginnen ze al te zweten terwijl ze nog in het dal zitten, voordat ze de heuvel opgaan.

4. De resultaten: Een betere voorspelling

De auteurs hebben hun model getest op het isotoop Uranium-236.

• Ze keken naar hoeveel neutronen er vrijkwamen.
• Ze keken naar de energie van die neutronen.
• Ze keken naar de verdeling van de gewichten van de twee splijtingsstukken.

Het verdict: Hun nieuwe model, dat rekening houdt met het "tijdig" verliezen van neutronen, gaf resultaten die veel dichter bij de echte meetdata lagen dan oude modellen. Het verklaarde ook waarom bepaalde schalen (de "schaleffecten") in de kern belangrijk blijven, zelfs als de kern heet wordt.

Samenvatting in één zin

Deze wetenschappers hebben een simulatie van een atoomkern die uit elkaar valt verbeterd door er rekening mee te houden dat de kern onderweg "zwetend" neutronen kwijtraakt, wat zorgt voor een veel nauwkeurigere voorspelling van hoe de kern uiteenvalt en welke brokstukken er overblijven.

Het is alsof ze niet alleen de route van de wandelaar hebben getekend, maar ook precies hebben berekend hoe zijn gewicht verandert door het zweet dat hij onderweg verliest, en hoe dat zijn loopstijl beïnvloedt.

Technische samenvatting

Titel: Shell-effecten en neutronenemissie binnen het multi-dimensionale Langevin-model voor kernsplijting

Auteurs: F.A. Ivanyuk, S.V. Radionov, C. Ishizuka, S. Chiba
Instituut: Instituut voor Kernonderzoek (Kiev, Oekraïne) en Instituut voor Wetenschap Tokio (Japan)

---

1. Probleemstelling

De theoretische beschrijving van kernsplijting, met name het proces van de samengestelde kern tot het scheidingspunt (scission point), is complex. Traditionele Langevin-benaderingen hebben de evolutie van de kervorm succesvol beschreven en kunnen grootheden zoals de massa-verdeling van splijtingsfragmenten en kinetische energieën goed reproduceren.

Echter, een belangrijk aspect dat in veel eerdere modellen werd verwaarloosd of vereenvoudigd, is de tijdsafhankelijke emissie van neutronen tijdens het splijtingsproces zelf (pre-scission neutronen).

• Bij excitatie-energieën onder de 80 MeV is neutronverdamping de dominante de-excitatiemodus.
• Bestaande modellen berekenen vaak alleen de totale emissie of nemen emissie alleen voor het scheidingspunt in beschouwing zonder de dynamische koppeling met de vormverandering van de kern.
• Er is behoefte aan een model dat de simultane evolutie van de kervorm en de neutronenemissie beschrijft, zodat men kan bepalen wanneer (op welk stadium van vervorming) en hoeveel neutronen worden uitgezonden, en hoe dit de overgebleven excitatie-energie en de uiteindelijke splijtingskarakteristieken beïnvloedt.

2. Methodologie

De auteurs hebben een uitgebreid vijf-dimensionaal Langevin-model aangepast om neutronenemissie dynamisch te integreren.

• Langevin-vergelijkingen: De tijdsevolutie van de collectieve variabelen ($q_\mu$) die de vorm van de kern beschrijven (gebaseerd op het Two-Center Shell Model - TCSM), wordt opgelost. De variabelen omvatten de afstand tussen de twee zwaartepunten, de vervorming van de linker- en rechterdelen, de massasymmetrie en de nekparameter.
  • Vereenvoudiging: Om de rekentijd beheersbaar te houden bij de toevoeging van emissie, wordt de nekparameter ($\epsilon$) vastgehouden op 0,25 (i.p.v. de gebruikelijke 0,35), wat bleek te leiden tot betere overeenstemming met experimentele data voor $^{236}\text{U}$.
• Dynamische koppeling met neutronenemissie:
  • Op elke integratiestap wordt de kans op neutronenemissie berekend.
  • Als een neutron wordt uitgezonden, wordt de lokale excitatie-energie ($E^*$) verlaagd met de neutronen-scheidingsenergie ($S_n$) plus de gemiddelde kinetische energie van het uitgezonden neutron.
  • De kern "springt" vervolgens naar een potentiaalenergie-oppervlak met één neutron minder.
  • De trajecten worden verder geïntegreerd tot het scheidingspunt (nekstraal < 2 fm).
• Emissieformule ($dN/dt$):
  • In plaats van empirische modellen, hebben de auteurs een analytische benadering voor de emissiesnelheid afgeleid op basis van de continuïteitsvergelijking en het Fermi-gas model.
  • De snelheid wordt berekend als de flux van deeltjes die de kern verlaten, afhankelijk van de lokale vorm, excitatie-energie en collectieve snelheden.
  • De formule houdt rekening met de diepte van het potentiaalput ($U_0$) en de chemische potentiaal ($\mu$), die vervormingsafhankelijk is.
• Stochastische emissie: Op elke stap wordt een willekeurig getal gegenereerd. Als de berekende emissie $\Delta N$ groter is dan dit getal, wordt een neutron uitgezonden. De energie van het uitgezonden neutron wordt getrokken uit een verdeling rondom $\langle e_n \rangle = 3/2 T$.

3. Belangrijkste Bijdragen

1. Integratie van emissie in de dynamica: Voor het eerst wordt de neutronenemissie niet als een statisch proces behandeld, maar als een dynamisch proces dat de excitatie-energie en de baan van de kern tijdens de splijting direct beïnvloedt.
2. Analytische uitdrukking voor $dN/dt$: De ontwikkeling van een formule voor de emissiesnelheid die lokaal is (afhankelijk van de huidige vorm en energie) en niet afhankelijk is van globale eigenschappen zoals de hoogte van de splijtingsbarrière of de terugwaartse reactie.
3. Gedetailleerde analyse van emissiestadia: Het model maakt het mogelijk om te onderscheiden of neutronen worden uitgezonden vanuit de grondtoestand, bij de zadel-punt (saddle), of tussen het zadelpunt en het scheidingspunt.
4. Toepassing op $^{236}\text{U}$: Het model is getest voor de splijting van $^{236}\text{U}$ bij excitatie-energieën tussen 10 en 45 MeV.

4. Resultaten

De berekeningen zijn vergeleken met beschikbare experimentele data:

• Massa-verdeling van fragmenten:
  • De rekeninghouding met neutronenemissie verbetert de overeenstemming tussen theorie en experiment aanzienlijk.
  • De emissie versterkt het effect van schaleffecten (shell effects) op de opbrengst (yield).
  • De demping van schaleffecten bij toenemende excitatie-energie wordt correct gereproduceerd.
• Tijdstip van emissie (vervormingsafhankelijkheid):
  • Bij lage excitatie-energie ($E^* = 10$ MeV) worden geen neutronen uitgezonden vanuit het grondtoestandsgebied. Als er een neutron zou worden uitgezonden, zou de resterende energie onder de splijtingsbarrière dalen, waardoor de kern niet zou splijten (evaporatie-residu). Emissie vindt hier alleen plaats na het passeren van de barrière.
  • Bij hogere energieën ($E^* \ge 20$ MeV) wordt een klein deel van de neutronen wel vanuit het grondtoestandsgebied uitgezonden.
  • Het merendeel van de pre-scheidingsneutronen wordt uitgezonden bij grotere vervormingen, tussen het zadelpunt en het scheidingspunt.
• Neutronspectra:
  • De berekende spectra van pre-scheidingsneutronen komen redelijk overeen met de standaard Watt-verdeling die wordt gebruikt voor prompte splijtingsneutronen.
• Pre-scheidings neutronenmultipliciteit ($M_{pre}$):
  • De berekende multipliciteit (het gemiddelde aantal neutronen per splijting) komt goed overeen met de beperkte experimentele data die beschikbaar is voor $^{236}\text{U}$ en verwante reacties.

5. Betekenis en Conclusie

Dit werk biedt een meer realistisch en gedetailleerd beeld van het kernsplijtingsproces. Door de koppeling tussen vormevolutie en de-excitatie (neutronenemissie) te versterken, kunnen theoretici:

• Beter voorspellen hoe de excitatie-energie tijdens het proces afneemt.
• De invloed van schaleffecten op de uiteindelijke fragmentverdeling nauwkeuriger modelleren.
• Inzicht krijgen in de tijdschaal van neutronenemissie, wat experimenteel moeilijk te meten is maar cruciaal is voor het begrijpen van de dynamica van de kern.

De methode bewijst dat een Langevin-aanpak, zelfs met beperkte dimensies (4D in dit geval door vastzetten van de nekparameter), zeer nauwkeurige resultaten kan leveren voor zowel massa-verdelingen als neutronenemissie-eigenschappen, mits de emissie dynamisch wordt gekoppeld aan de trajecten. Dit is een belangrijke stap vooruit voor de modellering van splijting in zowel fundamenteel onderzoek als toepassingen (zoals kernenergie en nucleaire veiligheid).