Angular and Kinetic Properties of Scission Neutrons within Time-dependent Density Functional Theory

Diese Studie nutzt die zeitabhängige Dichtefunktionaltheorie, um zu zeigen, dass Spaltneutronen eine signifikante, hochenergetische Komponente des prompten Spaltneutronenspektrums bei der Spaltung von $^{235}\mathrm{U}$, $^{239}\mathrm{Pu}$ und $^{252}\mathrm{Cf}$ darstellen, was einen direkten Beleg dafür liefert, dass deren Einbeziehung die systematische Unterschätzung der Hochenergieausbeuten behebt, die in traditionellen Modellen, die nur auf Verdampfung basieren, beobachtet wurde.

Das Wesentliche

Stellen Sie sich einen riesigen, instabilen Ballon (einen Atomkern) vor, der plötzlich in zwei Teile zerplatzt. Das ist die Kernspaltung. Seit Jahrzehnten wissen Wissenschaftler, dass die beiden resultierenden Teile (sogenannte Fragmente) mit unglaublicher Geschwindigkeit auseinanderfliegen und winzige Teilchen namens Neutronen ausspucken.

Lange Zeit glaubten Physiker, dass all diese Neutronen erst später „verdampfen“, wie Dampf, der von einer heißen Tasse Kaffee aufsteigt, nachdem das Wasser gekocht hat. Sie nahmen an, dass die Fragmente bereits vollständig geformt waren und sich stetig bewegten, bevor sie anfingen, diese Neutronen loszulassen.

Dieses neue Paper legt jedoch nahe, dass einige Neutronen genau in dem Moment, in dem der Ballon zerplatzt, „herausgeschleudert“ werden. Dies sind sogenannte Scissions-Neutronen (Spaltmoment-Neutronen). Sie werden in dem chaotischen, sekundenbruchteillangen Chaos des Bruchs geboren, nicht erst von einem ruhigen, abkühlenden Fragment zu einem späteren Zeitpunkt.

Hier ist die Erklärung, wie die Forscher den Beweis für diese „Snap-Zeit“-Neutronen fanden, vereinfacht dargestellt:

1. Die Supercomputer-Simulation

Um zu sehen, was während der Spaltung geschieht, nutzten die Wissenschaftler keinen Mikroskop; stattdessen verwendeten sie einen Supercomputer, um einen Film des Ereignisses mithilfe einer Theorie namens zeitabhängiger Dichtefunktionaltheorie (TDDFT) laufen zu lassen.

Man kann sich das wie ein hochgeschwindigkeitsfähiges, 3D-Videospiel vorstellen, in dem sie die Atome tanzen und auseinanderbrechen lassen. In früheren Versionen dieses „Spiels“ war die virtuelle Welt zu klein. Die Neutronen stießen gegen den Rand des Bildschirms, bevor die Wissenschaftler genau bestimmen konnten, wie schnell sie waren oder in welche Richtung sie flogen.

In dieser Studie bauten sie eine viel größere virtuelle Welt (etwa 3 Mal größer als zuvor). Dies gab den Neutronen genug Raum, um herauszufliegen und sich zu setzen, damit die Wissenschaftler sie genau messen konnten, ohne dass die „Wände“ der Simulation die Daten verfälschen würden.

2. Die Entdeckung der „Geschwindigkeitsbegrenzung“

Sob-elle hatten eine klare Sicht, beobachteten sie die Neutronen, die in bestimmten Winkeln herausfliegen (hauptsächlich seitlich und leicht rückwärts bezüglich der Spaltung). Dabei fanden sie etwas Überraschendes:

• Die „No-Go“-Zone: Es gibt keine Scissions-Neutronen mit niedriger Energie (unter etwa 1,5 bis 2 Millionen Elektronenvolt). Es ist, als gäbe es eine Geschwindigkeitsbegrenzung; nichts Langsames darf ein „Snap-Zeit“-Neutron sein.
• Die Hochgeschwindigkeits-Gruppe: Stattdessen sind diese Neutronen alle schnell. Sie häufen sich um eine spezifische hohe Geschwindigkeit (3–3,5 MeV) und laufen dann in einem langen Schweif noch schnellerer Teilchen aus.

Es ist wie eine Menge Menschen, die von einem Sprungbrett springt. Die „verdampfenden“ Neutronen sind wie Menschen, die später lässig vom Pooldeck treten. Die „Scissions“-Neutronen sind wie Menschen, die genau in der Sekunde, in der das Brett bricht, heftig vom Brett weggeschleudert werden. Diejenigen, die vom Brett weggeschleudert werden, bewegen sich immer schnell; man sieht niemals ein langsames Neutron von diesem speziellen Ereignis.

3. Die Lösung des Rätsels der „fehlenden“ Energie

Wissenschaftler versuchen schon seit Jahren, ihre Computermodelle mit realen Experimenten in Einklang zu bringen. Sie hatten ein Problem:

• Das alte Modell: Wenn man nur die „Dampf“-Neutronen zählt, sagt dein Computermodell zu wenige hochenergetische Neutronen voraus. Es ist, als versuche man, einen Eimer mit einer kleinen Tasse zu füllen, aber der Eimer benötigt ständig mehr Wasser, als die Tasse liefern kann.
• Das neue Modell: Als die Forscher die „Snap-Zeit“-Neutronen (diejenigen, die sie in ihrer großen Simulation fanden) zu den „Dampf“-Neutronen hinzufügten, stimmte die Mathematik endlich. Das kombinierte Modell entsprach perfekt den hochenergetischen Daten, die in realen Experimenten für Uran und Californium gemessen wurden.

4. Warum das wichtig ist

Dies ist eine große Sache, denn es ist das erste Mal, dass eine rein mikroskopische Theorie (eine, die nicht einfach nur raten oder annehmen, dass Dinge existieren) diese „Snap-Zeit“-Neutronen vorhergesagt und bewiesen hat, dass sie real sind.

• Vorher: Mussten Wissenschaftler raten, dass diese Neutronen existieren, weil die Mathematik nicht aufging.
• Jetzt: Die Computersimulation produzierte sie ganz natürlich, ohne dass man ihr gesagt hatte, dass sie existieren. Es ist, als würde man einen Sturm vorhersagen, indem man die Wolken wandern beobachtet, anstatt nur anzunehmen, dass ein Sturm kommt, weil der Wetterbericht es sagt.

Das Fazeresultat

Das Paper kommt zu dem Schluss, dass beim Zerfall eines Atoms ein kleiner, aber wichtiger Teil der Neutronen (etwa 6 % bis 10 % der Gesamtzahl) im gewaltsamen Moment des Bruchs geboren wird. Diese Neutronen sind dadurch unterscheidbar, dass sie immer schnell und niemals langsam innerhalb bestimmter Winkel sind.

Indem sie diesen „Fingerabdruck“ in den Daten fanden, ist es den Forschern endlich gelungen, die „Snap-Zeit“-Neutronen von den „Dampf“-Neutronen zu trennen, was uns ein klareres und genaueres Bild davon vermittelt, wie Kernspaltung tatsächlich funktioniert. Dies hilft dabei, unser Verständnis der fundamentalen Kräfte, die Materie zusammenhalten oder zerreißen, zu verfeinern.

Technische Zusammenfassung

Technisches Resümee: Winkel- und kinetische Eigenschaften von Scissions-Neutronen innerhalb der zeitabhängigen Dichtefunktionaltheorie

Problemstellung
Trotz der Entdeckung der Kernspaltung vor fast einem Jahrhundert bleibt eine vollständige mikroskopische Beschreibung des Prozesses aus, insbesondere hinsichtlich des Übergangs vom äußeren Sattelpunkt zur Scission (Spaltung) und der anschließenden Bildung primärer Spaltfragmente (FFs). Aktuelle statistische Verdampfungscodes gehen davon aus, dass alle prompten Neutronen isotrop von voll beschleunigten Fragmenten emittiert werden. Dieses Framework kann jedoch nicht die Neutronen erklären, die während des hochgradig nicht-äquilibrierten dynamischen Prozesses des Halsbruchs (Neck Rupture) emittiert werden, bekannt als Scissions-Neutronen (SNs). Während SNs seit 1939 postuliert und mittels verschiedener Modelle untersucht wurden, bleiben sie experimentell ungelöst, da prompt-neutronen-bezogene Observablen allein nicht zwischen bei der Scission emittierten Neutronen und später verdampften Neutronen unterscheiden können. Frühere Versuche, SNs mittels zeitabhängiger Dichtefunktionaltheorie (TDDFT) vorherzusagen, waren durch kleine Simulationsdomänen begrenzt, was eine gleichzeitige Bestimmung der Winkel- und kinetischen Eigenschaften verhinderte, bevor die Neutronen die Rechengrenzen erreichten.

Methodik
Diese Studie verwendet einen vollständig mikroskopischen Rahmen basierend auf TDDFT, erweitert auf superflüssige Systeme, um die SN-Emission in $^{235}\text{U}(n_{th}, f)$, $^{239}\text{Pu}(n_{th}, f)$ und $^{252}\text{Cf}(sf)$ zu untersuchen.

• Simulationsaufbau: Die Autoren nutzten den axial-symmetrischen Harmonischen-Oszillator-Solver AxialHOHFB mit dem SkM*-Energiedichtefunktional (EDF). Ein entscheidender Fortschritt war die Verwendung eines wesentlich größeren Simulationsbereichs (ein Ellipsoid mit den Halbachsen $Z_{max} = 70$ fm und $R_{max} = 55$ fm), der etwa 3,2-mal größer ist als in früheren Studien.
• Zweistufige Integration: Um Rechenbeschränkungen zu bewältigen und gleichzeitig die Genauigkeit zu wahren, wurde ein zweistufiges Zeitintegrationsverfahren angewandt. Das System wurde von der statischen eingeschränkten Hartree-Fock-Bogoliubov-Lösung (HFB) zum Scissions-Moment in einer kleineren Basis ($B_1$) entwickelt, dann in eine größere Basis ($B_2$) transformiert, um die Entwicklung fortzusetzen, bis die Neutronen sich der Grenze näherten.
• Extraktion der SN-Eigenschaften: SNs wurden durch Definition einer Region $V^{(SN)}_\eta$ identifiziert, in der das Verhältnis der kollektiven Fluss-Energiedichte zur intrinsischen kinetischen Energiedichte einen Schwellenwert von $\eta = 3$ überschreitet. Dieses Kriterium isoliert Neutronen, die von einem nach außen gerichteten kollektiven Fluss dominiert werden.
• Winkel- und kinetische Analyse: Innerhalb dieser Region wurden lokale Geschwindigkeitsfelder und die kollektive Fluss-Energie pro Neutron berechnet. Die Studie konzentrierte sich auf das Winkelintervall $110^\circ \leq \theta \leq 150^\circ$ (gemessen relativ zur Richtung des leichten Spaltfragments), in dem die SN-Ausbeute sättigt, bevor Grenzflächenreflexionen auftreten.
• Vergleich mit Experimenten: Das berechnete SN-Kinetikspektrum wurde mit einem Maxwellschen Verdampfungsmodell für die verdampfte Neutronen-Komponente (EN) kombiniert. Die Parameter des Verdampfungsmodells (Fragmenttemperatur $T$ und Multiplizitätsverhältnis $\bar{\nu}_L/\bar{\nu}_H$) wurden ausschließlich durch niederenergetische experimentelle Daten ($E \leq E_{thr}$) eingeschränkt. Dieses kombinierte Modell wurde dann mit hochenergetischen Prompt-Fission-Neutronenspektren (PFNS) für $^{239}\text{Pu}$ und $^{252}\text{Cf}$ verglichen.

Kernergebnisse

• Charakteristika der SN-Spektren: Die extrahierten SN-Spektren für alle drei Systeme zeigen ein ausgeprägtes Schwellenwertverhalten. SNs fehlen unterhalb eines spezifischen Energieschwellenwerts ($E_{thr}$), der etwa 1,5–2 MeV beträgt (speziell $\approx 2$ MeV für induzierte Spaltung und $\approx 1,5$ MeV für spontane Spaltung). Die Spektren peaken bei 3–3,5 MeV und weisen einen exponentiellen hochenergetischen Ausläufer auf, der bis zu etwa 18 MeV reicht.
• Schätzungen der Ausbeute: Für $^{239}\text{Pu}(n_{th}, f)$ wird die gesamte SN-Ausbeute mit $0,19 \pm 0,01$ berechnet (ca. 6,6 % der gesamten prompten Neutronen). Für $^{252}\text{Cf}(sf)$ beträgt die Ausbeute $0,39 \pm 0,02$ (ca. 10,4 %).
• Modellvalidierung: Wenn die berechnete SN-Komponente zum Verdampfungsmodell hinzugefügt wird (das nur durch niederenergetische Daten eingeschränkt ist), reproduziert das resultierende gesamte PFNS die gemessene hochenergetische Ausbeute für sowohl $^{239}\text{Pu}$ als auch $^{252}\text{Cf}$ präzise. Im Gegensatz dazu unterschätzt das reine Verdampfungsmodell die hochenergetische Ausbeute systematisch, selbst innerhalb der 95 %-Konfidenzintervalle.
• Winkelstabilität: Die Studie bestätigt, dass zuverlässige asymptotische Eigenschaften nur innerhalb des Winkelbereichs $110^\circ \leq \theta \leq 150^\circ$ extrahiert werden können. Außerhalb dieses Bereichs verhindern transiente Interferenzeffekte und Grenzflächenreflexionen, dass die Ausbeute innerhalb des aktuellen Simulationsbereichs sättigt.

Bedeutung und Ansprüche
Die Arbeit beansprucht, die erste Identifizierung einer Scissions-Neutronen-Signatur innerhalb bestehender gemessener Prompt-Fission-Neutronenspektren mittels eines vollständig mikroskopischen Rahmens geliefert zu haben. Durch den Nachweis, dass die TDDFT-Methode SNs ohne a priori Annahmen natürlich erzeugt und dass diese SNs die Diskrepanz zwischen reinen Verdampfungsmodellen und hochenergetischen experimentellen Daten auflösen, liefert die Studie einen direkten Beweis für eine nicht vernachlässigbare Scissions-Neutronen-Komponente in der prompten Spaltung.

Die Autoren postulieren, dass diese Arbeit ein mikroskopisches Fundament für die Interpretation von Spaltneutronenspektren jenseits des Standard-Statistik-Verdampfungsbildes schafft. Darüber hinaus identifizieren sie eine spezifische experimentelle Signatur: die Trennung der SN- und EN-Beiträge unterhalb des Energieschwellenwerts ($E_{thr}$). Dies impliziert, dass niederenergetische Messungen im spezifizierten Winkelbereich die robustesten Randbedingungen für die Verdampfungskomponente liefern, während der hochenergetische Ausläufer als primärer Indikator für den SN-Beitrag dient. Die Studie kommt zu dem Schluss, dass die Simulationen zwar auf spezifische Winkelbereiche beschränkt sind, die Methodik jedoch erfolgreich eine messbare Signatur von Scissions-Neutronen isoliert, die theoretisch seit fast einem Jahrhundert postuliert wird.