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

Die Studie löst mehrdimensionale Langevin-Gleichungen, um die zeitliche Entwicklung von Spaltprozessen unter Berücksichtigung von Schaleffekten und der Emission von Neutronen zu simulieren, wobei berechnete Größen wie die prä-szissionale Neutronenmultiplizität und die Fragmentmassenverteilung mit experimentellen Daten verglichen werden.

Das Wesentliche

🌌 Die Reise eines Atomkerns: Wenn ein riesiger Ballon platzt

Stellen Sie sich vor, ein Atomkern ist wie ein riesiger, überdehnter Wasserballon, der mit Wasser (den Neutronen und Protonen) gefüllt ist. Normalerweise ist dieser Ballon stabil. Aber wenn man ihm genug Energie zuführt (wie wenn man ihn stark aufbläst), wird er instabil. Irgendwann wird er so dünn, dass er in zwei Teile reißt. Das nennt man Kernspaltung.

Die Autoren dieses Papers haben sich gefragt: Was passiert genau in diesem Moment des Platzens? Und noch wichtiger: Wie viele kleine Wassertropfen (Neutronen) fliegen dabei weg, bevor der Ballon endgültig in zwei Hälften zerfällt?

1. Der Simulator: Ein chaotischer Tanz

Die Wissenschaftler nutzen ein mathematisches Werkzeug namens „Langevin-Gleichung". Das können Sie sich wie einen sehr detaillierten Simulator für einen chaotischen Tanz vorstellen.

• Die Tänzer: Der Atomkern ist der Tänzer. Er verändert ständig seine Form (er wird langgezogen, dünn am Hals, dann reißt er).
• Die Musik: Die Musik ist die Energie, die im Kern steckt.
• Der Zufall: Der Tanz ist nicht perfekt choreografiert. Es gibt kleine Stöße (wie Windböen), die den Tänzer unvorhersehbar bewegen. Das Modell berechnet Millionen von möglichen Tänzen (Trajektorien), um zu sehen, wie oft der Kern in zwei Stücke zerfällt und wie viele „Tropfen" dabei verloren gehen.

2. Das Problem: Die „verlorenen" Tropfen

In früheren Modellen haben die Wissenschaftler oft nur auf das Endergebnis geschaut: Wie sehen die zwei Hälften aus? Sie haben aber oft vergessen, dass der Kern während des Tanzens ständig kleine Neutronen „verliert" (verdampft), genau wie ein heißer Kaffee Dampftropfen verliert, bevor er abkühlt.

Diese Arbeit füllt diese Lücke. Sie fragen: Verliert der Kern Neutronen, während er noch im Tal der Stabilität ist? Oder erst, wenn er den Berg der Instabilität (die Barriere) überquert hat?

3. Die neue Methode: Ein Zähler in Echtzeit

Die Autoren haben ihren Simulator verbessert. Statt nur am Ende zu zählen, schauen sie bei jedem einzelnen Schritt des Tanzes:

• Ist gerade genug Energie da, um einen Neutronen-Tropfen abzuschnallen?
• Ja? Dann wird die Energie des Kerns sofort verringert (wie wenn der Ballon etwas Luft verliert und weniger stark aufbläht).
• Nein? Der Tanz geht weiter.

Dadurch können sie genau sehen, wann und wo die Neutronen wegfliegen.

4. Die wichtigsten Entdeckungen (Die Geschichte der 236-Uran-Kerne)

Die Forscher haben sich den Kern von Uran-236 angesehen und verschiedene Energiestufen getestet:

• Bei niedriger Energie (10 MeV):
  Stellen Sie sich vor, der Kern hat nur wenig Energie. Er versucht, den Berg der Instabilität zu überqueren. Aber sobald er einen Neutronen verliert, hat er nicht mehr genug Energie, um den Berg zu erklimmen. Er fällt zurück ins Tal und verdampft einfach (wird zu einem anderen Kern), ohne zu spalten.

  • Ergebnis: Bei niedriger Energie fliegen keine Neutronen weg, bevor der Kern den Berg überquert. Alle Neutronen kommen erst, wenn der Kern schon auf dem Weg zum Zerreißen ist.

• Bei hoher Energie (20–40 MeV):
  Hier hat der Kern genug Kraft. Er kann Neutronen verlieren und trotzdem noch genug Energie haben, um den Berg zu überqueren und zu spalten.

  • Ergebnis: Ein kleiner Teil der Neutronen fliegt schon weg, während der Kern noch im Tal ist. Aber die meisten Neutronen fliegen weg, wenn der Kern schon stark in die Länge gezogen ist, kurz bevor er reißt (zwischen dem Gipfel des Berges und dem Riss).

5. Warum ist das wichtig?

Stellen Sie sich vor, Sie wollen wissen, wie ein Kuchen backt. Wenn Sie nur das fertige Ergebnis sehen, wissen Sie nicht, ob der Teig während des Backens aufgegangen ist oder ob er zusammengefallen ist.

Durch dieses neue Modell können die Wissenschaftler:

1. Genau vorhersagen, wie viele Neutronen bei einer Spaltung entstehen (wichtig für Reaktorsicherheit).
2. Verstehen, wie die Form des Kerns (seine „Tanzbewegung") die Verteilung der Fragmente beeinflusst.
3. Die Theorie mit der Realität abgleichen: Die Ergebnisse des Simulators stimmen sehr gut mit echten Messdaten aus dem Labor überein. Das bedeutet, ihr „Tanz-Simulator" funktioniert wirklich gut!

Zusammenfassung in einem Satz

Die Autoren haben einen hochmodernen Simulator gebaut, der nicht nur das Platzen eines Atomkerns nachahmt, sondern auch genau verfolgt, wie viele kleine Neutronen dabei „wegfliegen" und wann genau das passiert – und haben damit bewiesen, dass die Form des Kerns während des Tanzes entscheidend dafür ist, wie die Spaltung abläuft.

Technische Zusammenfassung

Titel: Schaleffekte und Neutronenemission im mehrdimensionalen Langevin-Modell für die Kernspaltung

Autoren: F. A. Ivanyuk, S. V. Radionov, C. Ishizuka, S. Chiba

---

1. Problemstellung und Motivation

Die theoretische Beschreibung der Kernspaltung, insbesondere bei Anregungsenergien oberhalb der Spaltbarriere, stößt oft an Grenzen, wenn die Emission von leichten Teilchen (vorwiegend Neutronen) während des Spaltprozesses vernachlässigt wird. Bisherige Ansätze basierten häufig auf der Annahme der "ersten Chance-Spaltung" (First Chance Fission), bei der die De-Exzitation durch Teilchenemission vor dem Erreichen des Spaltpunkts (Scission Point) ignoriert wurde.

Das Hauptproblem besteht darin, dass die Emission von Neutronen die Anregungsenergie des Kerns reduziert, was wiederum die Wahrscheinlichkeit für das Erreichen des Spaltpunkts und die Eigenschaften der Spaltfragmente (z. B. Massendistribution) beeinflusst. Es fehlte an einem konsistenten Modell, das die zeitliche Entwicklung der Kernform (beschrieben durch Langevin-Gleichungen) und die stochastische Emission von Neutronen in Echtzeit koppelt, um zu bestimmen, wann und wie viele Neutronen in welchem Stadium der Spaltung emittiert werden.

2. Methodik

Die Autoren erweitern den bestehenden mehrdimensionalen Langevin-Ansatz um eine Kopplung an die Neutronenemission.

• Langevin-Gleichungen: Die zeitliche Entwicklung der kollektiven Variablen $q_\mu$ (Formparameter) und der konjugierten Impulse $p_\mu$ wird durch stochastische Differentialgleichungen gelöst.
  • Formparametrisierung: Es wird das Zwei-Zentren-Schalenmodell (TCSM) verwendet. Um die Rechenzeit zu begrenzen, wird die Anzahl der Variablen auf vier reduziert, indem der Halsparameter $\epsilon$ fest auf $0.25$ gesetzt wird (leicht unter dem Standardwert von $0.35$).
  • Potenzial: Die freie Energie $F(q, T)$ wird im makroskopisch-mikroskopischen Modell berechnet, wobei Schalenkorrekturen (Strutinsky-Methode) und Paarungsenergie berücksichtigt werden. Die Dämpfung der Schaleffekte bei steigender Temperatur wird durch eine spezifische Funktion $\Phi(E^*)$ modelliert.
• Kopplung der Neutronenemission:
  • Emissionsrate ($dN/dt$): Anstelle von globalen Modellen wird eine analytische Näherung für die momentane Emissionsrate entwickelt. Diese basiert auf der Kontinuitätsgleichung und dem Fermi-Gas-Modell. Die Rate hängt von lokalen Größen ab: der Kernform, der lokalen Anregungsenergie und der Temperatur.
  • Algorithmus: In jedem Integrationsschritt wird die Wahrscheinlichkeit der Neutronenemission berechnet.
    • Falls ein Neutron emittiert wird, wird die lokale Anregungsenergie $E^*$ reduziert um die Neutronen-Separationsenergie ($S_n = -\mu$, wobei $\mu$ das chemische Potential ist) plus die kinetische Energie des emittierten Neutrons.
    • Der Kern wechselt auf eine Potentialenergiefläche mit einer geringeren Neutronenzahl ($A-1$).
    • Die Trajektorie wird mit den neuen Parametern weiter integriert.
  • Stochastik: Die kinetische Energie der emittierten Neutronen wird als Gauß-Verteilung um den Mittelwert $\langle e_n \rangle = 3/2 T$ mit einer Varianz simuliert.

3. Wichtige Beiträge und Innovationen

• Echtzeit-Kopplung: Die Studie führt eine dynamische Kopplung zwischen der Formevolution und der Teilchenemission ein, anstatt diese als nachgelagerten Prozess zu betrachten.
• Analytische Emissionsrate: Entwicklung einer Formel für $dN/dt$, die nur von lokalen Größen abhängt und somit in den stochastischen Langevin-Prozess integrierbar ist, ohne auf globale Barrierenhöhen oder inverse Reaktionsquerschnitte zurückzugreifen.
• Untersuchung der Emissionsphasen: Das Modell erlaubt erstmals eine detaillierte Aufschlüsselung, in welchem Stadium der Spaltung (Grundzustandsbereich, Sattelpunkt, Über-Sattelpunkt) Neutronen emittiert werden.
• Anwendung auf $^{236}\text{U}$: Die Methode wird erfolgreich auf die Spaltung von $^{236}\text{U}$ bei Anregungsenergien von $10 < E^* < 45$ MeV angewendet.

4. Ergebnisse

Die Berechnungen wurden mit experimentellen Daten verglichen und zeigen folgende Ergebnisse:

• Massendistribution der Spaltfragmente:
  • Die Berücksichtigung der Neutronenemission verbessert die Übereinstimmung zwischen Theorie und Experiment erheblich.
  • Die Emission von Neutronen verstärkt den Einfluss von Schaleffekten auf die Ausbeuten.
  • Die Dämpfung der Schaleffekte mit steigender Anregungsenergie wird durch das verwendete Modell korrekt wiedergegeben.
• Verteilung der Emissionszeitpunkte:
  • Bei niedrigen Energien ($E^* = 10$ MeV) werden Neutronen fast ausschließlich nach dem Überqueren des Spaltbarriers emittiert. Emission im Grundzustandsbereich führt hier meist zu einer Verdampfungsrückstandsbildung (keine Spaltung), da die Energie nach der Emission unter die Barriere fällt.
  • Bei höheren Energien ($E^* \ge 20$ MeV) wird ein kleinerer Anteil der Neutronen bereits im Grundzustandsbereich emittiert.
  • Der Großteil der prä-scission-Neutronen wird im Bereich großer Deformationen (zwischen Sattelpunkt und Spaltpunkt) emittiert.
• Neutronenspektren:
  • Die berechneten Spektren der prä-scission-Neutronen stimmen gut mit der bekannten Watt-Verteilung überein.
• Prä-scission-Neutronen-Multiplizität ($M_{pre}$):
  • Die berechnete Multiplizität (Anzahl der Neutronen pro Spaltungsereignis vor dem Spaltpunkt) stimmt gut mit den wenigen verfügbaren experimentellen Daten überein.
  • Die Stabilität des Ergebnisses wurde durch Tests mit unterschiedlichen Anzahlen von Trajektorien ($N_{tra}$) verifiziert.

5. Bedeutung und Fazit

Diese Arbeit stellt einen bedeutenden Fortschritt in der theoretischen Beschreibung der Kernspaltung dar. Durch die Integration der Neutronenemission direkt in den mehrdimensionalen Langevin-Prozess konnten die Autoren:

1. Die Diskrepanz zwischen theoretischen Vorhersagen und experimentellen Massendistributionen signifikant verringern.
2. Ein detailliertes Bild davon gewinnen, wann und wo im Spaltprozess Neutronen emittiert werden, was experimentell nicht direkt messbar ist.
3. Zeigen, dass Schaleffekte auch bei der Berücksichtigung von Teilchenemission eine entscheidende Rolle für die Fragmentverteilung spielen.

Das entwickelte Modell bietet eine robuste Grundlage für zukünftige Studien zur Spaltung schwerer Kerne, insbesondere für die Vorhersage von Neutronenmultiplizitäten und Fragmentverteilungen bei verschiedenen Anregungsenergien, was für Anwendungen in der Kernphysik und Reaktorphysik relevant ist.