Double Chooz veröffentlicht die erste quantitative Messung des Antineutrino-Restflusses

Ein abgeschalteter Kernreaktor ist nicht still. Spaltprodukte, die sich in seinem Kern und seinen Becken angesammelt haben, zerfallen noch Monate oder Jahre weiter, und ihre Betazerfallsketten emittieren Antineutrinos. Die Double-Chooz-Kollaboration veröffentlicht die erste quantitative Messung dieses Restflusses, gewonnen während gleichzeitiger Stillstände beider Reaktoren des Kernkraftwerks Chooz.

Quelle: journals.aps.org

Einfach erklärt

Wird ein Holzfeuer gelöscht, glüht die Glut noch lange weiter. Ein Kernreaktor tut dasselbe, jedoch mit einem Leuchten, das kein Auge sehen kann: einem Strom geisterhafter Teilchen, den Antineutrinos, erzeugt durch die Betazerfallsketten der Spaltprodukte im bestrahlten Brennstoff, der in den Kernen verbleibt, und in den abgebrannten Brennelementen, die in Becken gelagert werden.

Während 17,2 Tagen Messzeit, in denen beide Reaktoren von Chooz in den französischen Ardennen abgeschaltet waren, isolierte der Nahdetektor — rund 400 Meter von den Kernen entfernt und durch 115 Meter Wasseräquivalent abgeschirmt — nach Abzug des Untergrunds 106 ± 18 verbleibende Kandidaten. Das sind etwa sechs pro Tag. Das detaillierte Modell sagt voraus, dass 44 % des Signals aus den Becken stammen, in denen der abgebrannte Brennstoff abkühlt; der Detektor trennt diesen Anteil nicht direkt ab.

Als praktische Anwendung wird die Überwachung von Kernmaterial erwogen: Antineutrinos durchqueren Materie mit sehr geringer Wechselwirkungswahrscheinlichkeit, sodass sich ein Reaktor beobachten lässt, ohne auf seinen Kern zuzugreifen. Allerdings beruht dieses Ergebnis auf mehreren Dutzend Kubikmetern Szintillator und einer unterirdischen Abschirmung von 115 Metern Wasseräquivalent. Es validiert eine Berechnung des Restflusses; es weist nicht nach, dass ein transportables Gerät an der Oberfläche eines Kernkraftwerks funktionieren würde.

Technische Erklärung

1. Wie man ein Antineutrino einfängt und warum so viele nötig sind. Genutzt wird der inverse Betazerfall (IBD): $\bar{\nu}_e + p \rightarrow e^+ + n$νˉe​+p→e++n. Ein Elektron-Antineutrino wechselwirkt mit einem freien Proton — einem Wasserstoffkern im Szintillator — und erzeugt ein Neutron sowie ein Positron. Die Signatur besteht aus zwei Teilen, und genau das macht die Messung trotz einer Rate von nur wenigen Ereignissen pro Tag möglich. Das Positron gibt sofort seine kinetische Energie ab und annihiliert anschließend: Das ist das prompte Signal. Das Neutron streut, bevor es eingefangen wird; die „Total Neutron Capture“-Methode von Double Chooz nutzt Einfänge an Gadolinium, Wasserstoff und Kohlenstoff in den inneren Volumina, um das verzögerte Signal zu bilden. Die Forderung, beide Signale in einem Fenster von 0,5 bis 800 µs und mit weniger als 1,2 m Abstand voneinander zu finden, entfernt den Großteil der Umgebungsradioaktivität. Die Antineutrinoenergie wird aus dem prompten Signal über $E_{\bar{\nu}_e} \simeq E_{\text{prompt}} + 0{,}78\ \text{MeV}$Eνˉe​​≃Eprompt​+0,78 MeV rekonstruiert; die Reaktion besitzt eine kinematische Schwelle von etwa 1,8 MeV: Unterhalb davon kann dieser Kanal das Antineutrino nicht nachweisen.

2. Woher das Restleuchten stammt: drei Ketten, drei Uhren. Im Betrieb bezieht ein Reaktor den Großteil seines Antineutrinoflusses aus kurzlebigen Spaltprodukten, deren Beitrag in den Stunden nach der Abschaltung einbricht. Der nachweisbare Rest stammt vor allem von drei Mutter-Tochter-Paaren: Während der Bestrahlung und der anschließenden Abkühlung speist ein langlebiges Mutternuklid fortlaufend ein wesentlich kurzlebigeres Tochternuklid. Weil die Tochter nach jeder Entstehung rasch wieder verschwindet, folgt ihre Aktivität hauptsächlich derjenigen der Mutter; die Abschaltung beseitigt vor allem den dominanten Fluss, der diese beständigen Ketten verdeckt hatte.

   Paar	Halbwertszeit der Mutter	Tochter: Halbwertszeit und maximale Energie $Q_\beta$Qβ​	Rolle
   $^{144}\mathrm{Ce} \rightarrow {}^{144}\mathrm{Pr}$144Ce→144Pr	285 Tage	17,3 min · 3,00 MeV	≈ 54 % des vom Modell erwarteten Signals von 1–3 MeV
   $^{106}\mathrm{Ru} \rightarrow {}^{106}\mathrm{Rh}$106Ru→106Rh	372 Tage	30,1 s · 3,54 MeV	≈ 38 % des vom Modell erwarteten Signals
   $^{90}\mathrm{Sr} \rightarrow {}^{90}\mathrm{Y}$90Sr→90Y	28,9 Jahre	3,19 h · 2,28 MeV	wird in der langfristigen Modellvorhersage dominant (> 90 % nach einem Jahrzehnt)

   Diese Energien bestimmen alles Weitere: Bei $Q_\beta$Qβ​-Werten von 2 bis 3,5 MeV und einer IBD-Schwelle von 1,8 MeV ist nur der hochenergetische Ausläufer jedes Betaspektrums nachweisbar. Entsprechend erwarten die Autoren 98,7 % des nachweisbaren Restflusses unterhalb von 3 MeV — das rechtfertigt die Beschränkung der Analyse auf das Fenster von 1–3 MeV prompter Energie und macht sie zugleich schwierig: Dort sind die Untergründe am dichtesten.

3. Was das Fenster von 1–3 MeV genau misst und welche Kontrolle es validiert. Eine Analyse auf ein Fenster zu beschränken, in dem das Signal erwartet wird, provoziert einen klassischen Einwand: Man findet das Gesuchte dort, wo man danach sucht. Die Autoren begegnen ihm mit einer Negativkontrolle. Oberhalb von 3 MeV sind die Raten nach Untergrundabzug mit null vereinbar. Das Modell sagt zwar einen kleinen Restausläufer voraus, aber mit weniger als 0,1 Ereignis pro Tag und damit unterhalb der heutigen Empfindlichkeit. Dass in diesem Bereich kein messbarer Überschuss vorliegt, stützt das Untergrundmodell; es beweist nicht, dass der physikalische Fluss dort strikt null ist.

4. Die Untergründe und weshalb die unterirdische Abschirmung zählt. Bei sechs Ereignissen pro Tag ist jedes Detail wichtig. Drei Gruppen werden getrennt behandelt. Zufällige Koinzidenzen — zwei voneinander unabhängige Energiedepositionen, von denen eine für ein promptes und die andere für ein verzögertes Signal gehalten wird — werden durch ein neuronales Netz unterdrückt, dessen Rückweisungsvermögen 400 übersteigt; der verbleibende Anteil wird anschließend in situ aus der Rate isolierter Depositionen gemessen. Korrelierte Untergründe stammen von kosmischen Myonen: entweder von schnellen Neutronen, die durch Spallation im Gestein entstehen, ein promptes Signal nachahmen und danach eingefangen werden, oder von β-n-Isotopen wie $^{9}\mathrm{Li}$9Li, das durch Myonenspallation am $^{12}\mathrm{C}$12C des Szintillators entsteht und bei seinem Zerfall allein schon das Elektron-Neutron-Paar liefert. Letztere werden über ihre zeitliche und räumliche Korrelation mit dem erzeugenden Myon markiert. Die Abschirmung des Nahdetektors durch 115 m Wasseräquivalent reduziert den Myonenfluss und damit die Entstehung dieser Isotope, die kein Veto vollständig beseitigen kann.

5. Die Vorhersage: eine Inventarberechnung, Brennelement für Brennelement. Die erwartete Zahl ist keine Größenordnungsschätzung, sondern eine detaillierte Simulation. Das vorhergesagte IBD-Spektrum lautet
   $\frac{dN_{\bar{\nu}_e}}{dE_{\text{vis}}} = \sum_i \frac{N_p}{4\pi L_i^2}\iint \varphi_i(E,t)\,\sigma_{\text{IBD}}(E)\,P_{ee}(E,L_i)\,\mathcal{D}(E_{\text{vis}}|E)\,dE\,dt,$dEvis​dNνˉe​​​=i∑​4πLi2​Np​​∬φi​(E,t)σIBD​(E)Pee​(E,Li​)D(Evis​∣E)dEdt,
   wobei $N_p$Np​ die Zahl der Targetprotonen, $L_i$Li​ der Abstand zur Quelle $i$i, $\varphi_i$φi​ das zeitabhängige Emissionsspektrum, $\sigma_{\text{IBD}}$σIBD​ der Wirkungsquerschnitt, $P_{ee}$Pee​ die Überlebenswahrscheinlichkeit durch Oszillation und $\mathcal{D}$D die Detektorantwort ist. Die Aktivitäten der Spaltprodukte stammen aus den Codes APOLLO-2.8.4 und DARWIN-3 mit der Kerndatenbibliothek JEFF-3.1.1; die Betaspektren kommen aus der BESTIOLE-Bibliothek, die die Fermi-Theorie mit Korrekturen für elektromagnetische Effekte, endliche Kerngröße, atomare Abschirmung sowie Formfaktoren für erlaubte und verbotene Übergänge anwendet. Die Bestrahlungsgeschichte jedes Brennelements wurde einzeln simuliert und seine thermische Leistung über seine gesamte Lebensdauer bis zum Beginn jedes Stillstandszeitraums verfolgt.

6. Warum das Modell 44 % den Becken zuordnet. Die beiden Lagerbecken liegen nahe bei den Kernen — aus Sicht eines rund 400 m entfernten Detektors bleibt ihr geometrischer Faktor $1/L^2$1/L2 somit vergleichbar. Jedes Brennelement gelangt nach mehreren Bestrahlungszyklen dorthin und kühlt dort über Jahre ab. Ein abgeschalteter Kern enthält ein junges Inventar, das rasch zerfällt; ein Becken sammelt Brennelemente unterschiedlichen Alters, deren Isotope zum Teil lange fortbestehen. Bei einem längeren Stillstand wird die Summe der vorhergesagten Beiträge aus den Becken mit derjenigen aus den Kernen vergleichbar. Die integrierte Berechnung ergibt 56 % für die Kerne und 44 % für die Becken; mit der gegenwärtigen Geometrie lassen sich diese beiden Anteile nicht getrennt messen.

7. Was die Vorhersage begrenzt und welche Schlussfolgerungen dadurch ausgeschlossen sind. Die Normierungsunsicherheit der vorhergesagten Zahl erreicht 7,4 %; sie zerfällt in 6,0 % für die Spektralmodellierung, 2,9 % für die Geometrie der Basislinien, 2,1 % für das Spaltproduktinventar, 2,0 % für das Inventar der Brennelemente in den Becken und weniger als 0,8 % für den Detektor. Der dominante Term beruht auf der Form des Betazerfalls von $^{144}\mathrm{Pr}$144Pr; dessen Unsicherheit wird als Differenz zwischen einer detaillierten Kernstrukturrechnung und der sogenannten $\xi$ξ-Näherung angesetzt. Diese 7,4 % messen weder die Genauigkeit der Übereinstimmung zwischen Daten und Modell noch die Leistungsfähigkeit einer Inventarinferenz: Die Studie vergleicht eine Zählung mit einer Vorhersage, ohne umgekehrt eine Brennstoffmasse oder -zusammensetzung zu schätzen.

Notation der Grafik: $N_{\rm obs}$Nobs​ ist die gemessene Restzählung, $N_{\rm pred}$Npred​ die Modellvorhersage und $H_0(B)$H0​(B) die reine Untergrundhypothese. Die $5{,}9\,\sigma$5,9σ vergleichen $N_{\rm obs}$Nobs​ mit $H_0(B)$H0​(B), nicht mit $N_{\rm pred}$Npred​.

Warum es funktioniert hat

Die Messung verdankt sich einem Umstand und einer Geometrie. Der Umstand: Im Jahr 2017 waren die beiden Reaktoren von Chooz viermal gleichzeitig abgeschaltet, darunter ein Zeitraum von 20,8 Tagen — eine unverzichtbare Voraussetzung, denn der Nahdetektor zählt bei voller Leistung rund 900 IBD-Ereignisse pro Tag, gegenüber 6,2 ± 1,1 verbleibenden Ereignissen pro Tag in dieser Messung. Die Geometrie: ein bereits gebauter, kalibrierter und verstandener Detektor in 400 m Entfernung, dessen systematische Effekte durch jahrelange Oszillationsmessungen bestimmt worden waren.

Der Zahlenvergleich ist direkt. Eine frühere Analyse anhand der Daten von 2011–2012 hatte nur etwa zwanzig Kandidaten erfasst — zu wenige, um ein Restspektrum zu extrahieren. Hier werden nach Untergrundabzug 106 ± 18 verbleibende Kandidaten mit 88 ± 7 erwarteten verglichen: Die Differenz beträgt $\Delta = 18 \pm 19$Δ=18±19 Ereignisse, also eine Übereinstimmung. Am Ferndetektor stehen 27 ± 13 den erwarteten 14 ± 1 gegenüber, mit einem Fehlerbalken, der nur eine wesentlich lockerere Schlussfolgerung zulässt.

Die Grenzen sind ausdrücklich benannt. Die Signalanalyse bezieht sich auf das Fenster von 1–3 MeV. Oberhalb von 3 MeV sind die untergrundbereinigten Raten mit null vereinbar, während das Modell bis etwa 4,5 MeV einen kleinen Restausläufer beibehält, der unter 0,1 Ereignis pro Tag liegt und in dieser Stichprobe nicht aufgelöst wird. Die Normierungsunsicherheit von 7,4 % charakterisiert die Vorhersage; das Experiment belegt nicht die Genauigkeit einer Inventarinferenz.

Der Abstand zwischen dem Nachgewiesenen und dem Erwogenen sollte klar benannt werden. Nachgewiesen ist: Ein Restfluss wird bei $5{,}9\,\sigma$5,9σ vom Untergrund unterschieden, und die Zählung $106\pm18$106±18 ist mit der Vorhersage $88\pm7$88±7 vereinbar, weil $\Delta=18\pm19$Δ=18±19. Nicht nachgewiesen ist: eine operative Brennstoffüberwachung. Der Nachweis nutzt mehr als 30 m³ Szintillator unter 115 m Wasseräquivalent, für etwa sechs Ereignisse pro Tag in 400 m Entfernung. Ein kompaktes Instrument an der Oberfläche müsste den kosmischen Untergrund weiterhin beherrschen und seine Empfindlichkeit für Inventaränderungen nachweisen.

Kausalkette

Spaltung im Kern → Ansammlung langlebiger Mutternuklide während der Bestrahlung ($^{144}\mathrm{Ce}$144Ce, $^{106}\mathrm{Ru}$106Ru, $^{90}\mathrm{Sr}$90Sr), die kurzlebige Betatöchter speisen → Abschaltung der Reaktoren und rasches Verschwinden des dominanten Flusses → Fortbestehen der bereits etablierten Ketten in den Kernen und in den in die Becken überführten Brennelementen → vier gleichzeitige Stillstände mit insgesamt 24,4 Tagen im Jahr 2017 → 17,2 Tage Messzeit am Nahdetektor → $106\pm18$106±18 verbleibende Kandidaten zwischen 1 und 3 MeV, mit reinem Untergrund bei $5{,}9\,\sigma$5,9σ unvereinbar → kleiner vorhergesagter Ausläufer oberhalb von 3 MeV, unter 0,1 Ereignis pro Tag und nicht aufgelöst → Vereinbarkeit mit dem Brennelement-für-Brennelement-Modell ($88\pm7$88±7, $\Delta=18\pm19$Δ=18±19) → erste quantitative Referenz zum Test von Restflussberechnungen → künftige Studien zur Trennung von Kern- und Beckenbeiträgen sowie mit eigens dafür vorgesehenen Detektoren.

Eine Anekdote

Double Chooz hatte diese Messung bereits mit den Daten von 2011–2012 versucht. Die Kollaboration hatte damals nur etwa zwanzig verbleibende Kandidaten erfasst: Der Nahdetektor war noch nicht verfügbar, und die gleichzeitige Stillstandszeit war zu kurz, um ein Spektrum zu extrahieren. Im Jahr 2017 war der Nahdetektor in Betrieb, und vier gleichzeitige Stillstände, darunter einer von 20,8 Tagen, lieferten 17,2 Tage Messzeit; diese beiden Veränderungen machten die quantitative Messung möglich.

Vermächtnis und aktueller Datenstand

Die Messung liefert den ersten veröffentlichten quantitativen Vergleich zwischen einer Summationsrechnung des Restflusses und einer realen Zählung. Die Vorhersage besitzt eine Normierungsunsicherheit von 7,4 %; beim beobachteten Vergleich beträgt die Differenz $\Delta=18\pm19$Δ=18±19 Ereignisse. Diese beiden Größen dürfen nicht verwechselt werden.

Auf der Anwendungsseite ist die vorhergesagte Aufteilung die unmittelbar nutzbarste neue Information: Das Modell ordnet 44 % des Signals des abgeschalteten Standorts den Becken zu, nicht den Kernen. Eine Überwachungsmethode, die nur den noch in den Kernen befindlichen Brennstoff modelliert, würde das erwartete Signal daher deutlich unterschätzen. Was die Veröffentlichung dagegen nicht beziffern lässt, ist die Empfindlichkeit eines kompakten Oberflächendetektors: Der Artikel enthält keine derartigen Daten, und die Autoren extrapolieren nicht von den aktiven Volumina der unterirdischen Anlage.

Aus Sicht der Forschung — offene Fragen

Diese Ansätze sind mögliche Fortsetzungen, keine Ergebnisse der Studie.

• Der entscheidende Spektraltest. Mit mehr Statistik könnte die relative Entwicklung von $^{144}\mathrm{Ce}$144Ce und $^{106}\mathrm{Ru}$106Ru über mehrere Monate verfolgt und so der zeitliche Verlauf des Modells geprüft werden. Die Dominanz des Paares $^{90}\mathrm{Sr}$90Sr–$^{90}\mathrm{Y}$90Y zu beobachten, würde Jahre erfordern: Die Vorhersage übersteigt 90 % erst nach mehr als zehn Jahren.
• Becken und Kerne durch Geometrie trennen. Die 56 %/44 % sind hier eine Modellvorhersage, keine gemessene Trennung: Bei 400 m lässt sich ein Abstandsunterschied von 38 m nicht auflösen. Ein deutlich näher platzierter Detektor oder zwei Detektoren unter verschiedenen Winkeln könnten die beiden Beiträge experimentell auseinanderhalten, statt sie rechnerisch aufzuteilen.
• Der Engpass der Kernstruktur. Die dominante Unsicherheit stammt nicht vom Detektor, sondern von der Form des Betazerfalls von $^{144}\mathrm{Pr}$144Pr. Eine gezielte Messung dieses Betaspektrums könnte den Term von 6,0 % verringern; die Studie beziffert jedoch nicht, welche Endgenauigkeit das Modell dann erreichen würde.

Quellen

Beim Faktencheck-Audit vom 16. August 2026 verifizierte Referenzen.

• Double Chooz Collaboration, „First Measurement of Neutrino Emissions from Spent Nuclear Fuel by the Double Chooz Experiment“, Physical Review Letters 137, 061803, 4. August 2026 — begutachteter Artikel. DOI: 10.1103/dr26-j19g
• Preprint-Version derselben Arbeit: arXiv:2510.04869 — konsultierter Volltext (Zusammenfassung, Methoden, Tabellen der systematischen Unsicherheiten).

Hintergrundliteratur

1. P. Vogel und J. F. Beacom — Angular distribution of neutron inverse beta decay — Referenz für Kinematik und Wirkungsquerschnitt des inversen Betazerfalls.

Vertrauensaussage

Hohes Vertrauen in die Zählungen, Signifikanzen, Energiefenster, Geometrie und Unsicherheitskomponenten, die aus dem begutachteten Artikel übernommen wurden. Mittleres Vertrauen in die praktische Tragweite: Die Veröffentlichung behandelt nukleare Sicherungsmaßnahmen, testet jedoch weder einen kompakten Oberflächendetektor noch eine Inventarinferenz. Der den Becken zugeordnete Anteil von 44 % ist eine Modellausgabe, keine gemessene Trennung. Für diese erste quantitative Analyse gibt es noch keine unabhängige Replikation. Eindeutigkeitsprüfung: Im Register A001–A018 erscheint weder ein DOI noch eine arXiv-Kennung oder ein gleichwertiges zentrales Ergebnis; Urteil UNIQUE. DOI, Daten, Formeln, Links und das Fehlen eines Duplikats wurden geprüft.