Double Chooz publie la première mesure quantitative du flux résiduel d'antineutrinos

Un réacteur nucléaire arrêté n'est pas silencieux. Des produits de fission accumulés dans son cœur et dans ses piscines continuent de décroître pendant des mois ou des années, et leurs chaînes bêta émettent des antineutrinos. La collaboration Double Chooz publie la première mesure quantitative de ce flux résiduel, obtenue à la faveur d'arrêts simultanés des deux réacteurs de la centrale de Chooz.

Source : journals.aps.org

En clair

Quand on éteint un feu de bois, les braises rougeoient encore longtemps. Un réacteur nucléaire fait la même chose, mais avec une lueur qu'aucun œil ne peut voir : un flux de particules fantômes, les antineutrinos, produit par les chaînes de désintégration bêta des produits de fission présents dans le combustible irradié restant dans les cœurs et dans les assemblages usés entreposés en piscine.

Pendant 17,2 jours de temps vif où les deux réacteurs de Chooz, dans les Ardennes, étaient à l'arrêt, le détecteur proche — situé à environ 400 mètres des cœurs et protégé par 115 mètres d'équivalent-eau — a isolé 106 ± 18 candidats résiduels après soustraction des fonds. C'est environ six par jour. Le modèle détaillé prévoit que 44 % du signal vient des piscines où le combustible usé refroidit ; cette répartition n'est pas séparée directement par le détecteur.

L'intérêt pratique envisagé est le contrôle des matières nucléaires : les antineutrinos traversent la matière avec une très faible probabilité d'interaction, ce qui permet d'observer un réacteur sans accès à son cœur. Attention toutefois : ce résultat repose sur plusieurs dizaines de mètres cubes de scintillateur et une protection souterraine de 115 mètres d'équivalent-eau. Il valide un calcul de flux résiduel ; il ne démontre pas qu'un appareil transportable fonctionnerait à la surface d'une centrale.

Explication technique

1. Comment on attrape un antineutrino, et pourquoi il en faut beaucoup. La réaction utilisée est la désintégration bêta inverse (IBD) : $\bar{\nu}_e + p \rightarrow e^+ + n$νˉe​+p→e++n. Un antineutrino électronique interagit avec un proton libre — un noyau d'hydrogène du scintillateur — et produit un neutron et un positon. La signature est double, et c'est elle qui rend la mesure possible malgré un taux de quelques événements par jour. Le positon dépose immédiatement son énergie cinétique puis s'annihile : c'est le signal prompt. Le neutron diffuse avant d'être capturé ; la méthode « total neutron capture » de Double Chooz exploite les captures sur le gadolinium, l'hydrogène et le carbone dans les volumes internes pour former le signal retardé. Exiger les deux signaux dans une fenêtre de 0,5 à 800 µs et à moins de 1,2 m l'un de l'autre élimine l'essentiel de la radioactivité ambiante. L'énergie de l'antineutrino se reconstruit à partir du prompt par $E_{\bar{\nu}_e} \simeq E_{\text{prompt}} + 0{,}78\ \text{MeV}$Eνˉe​​≃Eprompt​+0,78 MeV, et la réaction possède un seuil cinématique d'environ 1,8 MeV : en dessous, ce canal ne détecte pas l'antineutrino.

2. D'où vient la lueur résiduelle : trois chaînes, trois horloges. Un réacteur en marche tire l'essentiel de son flux d'antineutrinos des produits de fission à vie courte, dont la contribution s'effondre dans les heures suivant l'arrêt. Le reliquat détectable vient surtout de trois couples père-fils : pendant l'irradiation puis le refroidissement, un père à vie longue alimente continuellement un fils à vie beaucoup plus courte. Le fils disparaissant rapidement après chaque production, son activité suit principalement celle du père ; l'arrêt retire surtout le flux dominant qui masquait ces chaînes persistantes.

   Couple	Demi-vie du père	Fils : demi-vie et énergie maximale $Q_\beta$Qβ​	Rôle
   $^{144}\mathrm{Ce} \rightarrow {}^{144}\mathrm{Pr}$144Ce→144Pr	285 jours	17,3 min · 3,00 MeV	≈ 54 % du signal 1–3 MeV attendu par le modèle
   $^{106}\mathrm{Ru} \rightarrow {}^{106}\mathrm{Rh}$106Ru→106Rh	372 jours	30,1 s · 3,54 MeV	≈ 38 % du signal attendu par le modèle
   $^{90}\mathrm{Sr} \rightarrow {}^{90}\mathrm{Y}$90Sr→90Y	28,9 ans	3,19 h · 2,28 MeV	devient dominant dans la prédiction du modèle à long terme (> 90 % après une décennie)

   Ces énergies commandent tout le reste : avec des $Q_\beta$Qβ​ de 2 à 3,5 MeV et un seuil IBD à 1,8 MeV, seule la queue haute de chaque spectre bêta est détectable. Les auteurs attendent en conséquence 98,7 % du flux résiduel détectable en dessous de 3 MeV — c'est ce qui justifie de restreindre l'analyse à la fenêtre 1–3 MeV en énergie prompte, et c'est aussi ce qui la rend difficile : c'est là que les bruits de fond sont les plus denses.

3. Ce que mesure exactement la fenêtre 1–3 MeV, et le contrôle qui la valide. Restreindre l'analyse à une fenêtre où l'on attend le signal expose à une objection classique : on trouve ce qu'on cherche là où on le cherche. Les auteurs y répondent par un contrôle négatif. Au-dessus de 3 MeV, les rendements après soustraction du bruit de fond sont compatibles avec zéro. Le modèle prévoit bien une petite queue résiduelle, mais sous 0,1 événement par jour, donc sous la sensibilité actuelle. L'absence d'excès mesurable dans cette zone soutient le modèle de fond ; elle ne prouve pas que le flux physique y soit strictement nul.

4. Les bruits de fond, et pourquoi la protection souterraine compte. À six événements par jour, tout compte. Trois familles sont traitées séparément. Les accidentels — deux dépôts d'énergie sans lien, l'un pris pour un prompt, l'autre pour un retardé — sont supprimés par un réseau de neurones dont le pouvoir de réjection dépasse 400, puis ce qui reste est mesuré in situ à partir du taux de dépôts isolés. Les corrélés viennent des muons cosmiques : soit des neutrons rapides produits par spallation dans la roche, qui miment un prompt puis sont capturés, soit des isotopes β-n comme le $^{9}\mathrm{Li}$9Li, engendré par spallation muonique sur le $^{12}\mathrm{C}$12C du scintillateur, dont la désintégration fournit à elle seule le couple électron-neutron. Ces derniers sont étiquetés par corrélation temporelle et spatiale avec le muon parent. Les 115 m d'équivalent-eau de protection du détecteur proche réduisent le flux muonique, donc la production de ces isotopes, qu'aucun veto n'élimine parfaitement.

5. La prédiction : un calcul d'inventaire, assemblage par assemblage. Le nombre attendu n'est pas un ordre de grandeur mais une simulation détaillée. Le spectre IBD prédit s'écrit
   $\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,
   où $N_p$Np​ est le nombre de protons cibles, $L_i$Li​ la distance à la source $i$i, $\varphi_i$φi​ le spectre d'émission dépendant du temps, $\sigma_{\text{IBD}}$σIBD​ la section efficace, $P_{ee}$Pee​ la probabilité de survie par oscillation et $\mathcal{D}$D la réponse du détecteur. Les activités des produits de fission sortent des codes APOLLO-2.8.4 et DARWIN-3 avec la bibliothèque de données nucléaires JEFF-3.1.1 ; les spectres bêta viennent de la bibliothèque BESTIOLE, qui applique la théorie de Fermi avec corrections électromagnétiques, taille finie du noyau, écrantage atomique et facteurs de forme pour les transitions permises et interdites. Chaque assemblage de combustible a vu son histoire d'irradiation simulée individuellement, avec le suivi de sa puissance thermique sur toute sa vie, jusqu'au début de chaque période d'arrêt.

6. Pourquoi le modèle attribue 44 % aux piscines. Les deux piscines d'entreposage sont proches des cœurs — donc, vues d'un détecteur à environ 400 m, leur facteur géométrique $1/L^2$1/L2 reste comparable. Chaque assemblage y arrive après plusieurs cycles d'irradiation et y refroidit pendant des années. Un cœur à l'arrêt contient un inventaire récent qui décroît vite ; une piscine cumule des assemblages d'âges échelonnés, dont certains isotopes persistent longtemps. Sur un arrêt prolongé, la somme prédite des contributions des piscines devient comparable à celle des cœurs. Le calcul intégré donne 56 % pour les cœurs et 44 % pour les piscines ; la géométrie actuelle ne permet pas de mesurer séparément ces deux fractions.

7. Ce qui limite la prédiction, et ce que cela interdit de conclure. L'incertitude de normalisation du nombre prédit atteint 7,4 %, décomposée en 6,0 % de modélisation spectrale, 2,9 % de géométrie des lignes de base, 2,1 % d'inventaire des produits de fission, 2,0 % d'inventaire des assemblages en piscine et moins de 0,8 % pour le détecteur. Le terme dominant tient à la forme de la transition bêta du $^{144}\mathrm{Pr}$144Pr, dont l'incertitude est prise comme l'écart entre un calcul de structure nucléaire détaillé et l'approximation dite $\xi$ξ. Ces 7,4 % ne mesurent ni la précision de l'accord données-modèle ni une performance d'inférence d'inventaire : cette étude compare un comptage à une prédiction, sans estimer en retour une masse ou une composition de combustible.

Notation du schéma : $N_{\rm obs}$Nobs​ est le comptage résiduel mesuré, $N_{\rm pred}$Npred​ la prédiction du modèle et $H_0(B)$H0​(B) l'hypothèse du fond seul. Les $5{,}9\,\sigma$5,9σ comparent $N_{\rm obs}$Nobs​ à $H_0(B)$H0​(B), pas à $N_{\rm pred}$Npred​.

Pourquoi ça a marché

La mesure tient à une circonstance et à une géométrie. La circonstance : en 2017, les deux réacteurs de Chooz se sont trouvés simultanément à l'arrêt quatre fois, dont une période de 20,8 jours — condition indispensable, puisque le détecteur proche compte environ 900 événements IBD par jour à pleine puissance, contre 6,2 ± 1,1 événements résiduels par jour ici. La géométrie : un détecteur déjà construit, calibré et compris, installé à 400 m, dont les systématiques avaient été établies par des années de mesure d'oscillation.

La comparaison chiffrée est directe. Une analyse antérieure, sur les données 2011–2012, n'avait réuni qu'une vingtaine de candidats — trop peu pour extraire un spectre résiduel. Ici, après soustraction des fonds, 106 ± 18 candidats résiduels sont comparés à 88 ± 7 attendus : l'écart vaut $\Delta = 18 \pm 19$Δ=18±19 événements, c'est-à-dire un accord. Au détecteur lointain, 27 ± 13 pour 14 ± 1 attendus, avec une barre d'erreur qui laisse la conclusion bien plus lâche.

Les limites sont explicites. L'analyse de signal porte sur la fenêtre 1–3 MeV. Au-dessus de 3 MeV, les rendements soustraits sont compatibles avec zéro, tandis que le modèle conserve une petite queue résiduelle jusque vers 4,5 MeV, inférieure à 0,1 événement par jour et non résolue par cet échantillon. L'incertitude de normalisation de 7,4 % caractérise la prédiction ; l'expérience ne démontre pas la précision d'une inférence d'inventaire.

L'écart entre ce qui est démontré et ce qui est envisagé mérite d'être posé. Ce qui est démontré : un flux résiduel est distingué du fond à $5{,}9\,\sigma$5,9σ, et le comptage $106\pm18$106±18 est compatible avec la prédiction $88\pm7$88±7 puisque $\Delta=18\pm19$Δ=18±19. Ce qui ne l'est pas : une surveillance opérationnelle du combustible. La démonstration mobilise plus de 30 m³ de scintillateur sous 115 m d'équivalent-eau, pour environ six événements par jour à 400 m. Un instrument compact en surface devrait encore maîtriser le bruit cosmique et démontrer sa sensibilité à des changements d'inventaire.

Chaîne causale

Fission dans le cœur → accumulation, pendant l'irradiation, de parents à vie longue ($^{144}\mathrm{Ce}$144Ce, $^{106}\mathrm{Ru}$106Ru, $^{90}\mathrm{Sr}$90Sr) alimentant des filles bêta à vie courte → arrêt des réacteurs et disparition rapide du flux dominant → persistance des chaînes déjà établies dans les cœurs et les assemblages transférés en piscine → quatre arrêts simultanés totalisant 24,4 jours en 2017 → 17,2 jours de temps vif au détecteur proche → $106\pm18$106±18 candidats résiduels entre 1 et 3 MeV, incompatibles avec le fond seul à $5{,}9\,\sigma$5,9σ → petite queue prédite au-dessus de 3 MeV, inférieure à 0,1 événement par jour et non résolue → compatibilité avec le modèle assemblage par assemblage ($88\pm7$88±7, $\Delta=18\pm19$Δ=18±19) → première référence quantitative pour tester les calculs de flux résiduel → études futures de séparation cœur-piscine et de détecteurs dédiés.

Anecdote

Double Chooz avait déjà tenté cette mesure avec ses données 2011–2012. La collaboration n'avait alors réuni qu'une vingtaine de candidats résiduels : le détecteur proche n'était pas encore disponible et le temps d'arrêt simultané était trop court pour extraire un spectre. En 2017, le détecteur proche fonctionnait et quatre arrêts simultanés, dont un de 20,8 jours, ont fourni 17,2 jours de temps vif ; ces deux changements ont rendu possible la mesure quantitative.

Héritage et données actuelles

La mesure fournit le premier point de comparaison quantitatif publié entre un calcul de sommation du flux résiduel et un comptage réel. La prédiction possède une incertitude de normalisation de 7,4 % ; la comparaison observée est, elle, $\Delta=18\pm19$Δ=18±19 événements. Ces deux quantités ne doivent pas être confondues.

Sur le versant applicatif, la répartition prédite est l'information neuve la plus directement exploitable : le modèle attribue 44 % du signal du site à l'arrêt aux piscines, pas aux cœurs. Une méthode de surveillance qui ne modéliserait que le combustible encore dans les cœurs sous-estimerait donc fortement le signal attendu. Ce que la publication ne permet pas de chiffrer, en revanche, c'est la sensibilité qu'atteindrait un détecteur compact en surface : aucune donnée de ce type ne figure dans l'article, et l'extrapolation depuis les volumes actifs du dispositif souterrain n'est pas faite par les auteurs.

Le regard du chercheur — questions ouvertes

Ces pistes sont des prolongements possibles, non des résultats de l'étude.

• Le test spectral décisif. Avec plus de statistique, suivre plusieurs mois l'évolution relative du $^{144}\mathrm{Ce}$144Ce et du $^{106}\mathrm{Ru}$106Ru testerait la forme temporelle du modèle. Observer la domination du couple $^{90}\mathrm{Sr}$90Sr–$^{90}\mathrm{Y}$90Y demanderait des années : la prédiction ne dépasse 90 % qu'après plus de dix ans.
• Séparer les piscines des cœurs par la géométrie. Les 56 %/44 % sont ici une prédiction du modèle, pas une séparation mesurée : à 400 m, 38 m d'écart ne se résolvent pas. Un détecteur placé nettement plus près, ou deux détecteurs à des angles différents, pourrait dissocier expérimentalement les deux contributions au lieu de les répartir par calcul.
• Le verrou de la structure nucléaire. L'incertitude dominante ne vient pas du détecteur mais de la forme de la transition bêta du $^{144}\mathrm{Pr}$144Pr. Une mesure dédiée de ce spectre bêta pourrait réduire le terme de 6,0 %, mais l'étude ne quantifie pas la précision finale qu'atteindrait alors le modèle.

Sources

Références vérifiées lors de l'audit de fact-checking du 16 août 2026.

• Double Chooz Collaboration, « First Measurement of Neutrino Emissions from Spent Nuclear Fuel by the Double Chooz Experiment », Physical Review Letters 137, 061803, 4 août 2026 — article évalué par les pairs. DOI : 10.1103/dr26-j19g
• Version préprint du même travail : arXiv:2510.04869 — texte intégral consulté (résumé, méthodes, tableaux de systématiques).

Références de fond

1. P. Vogel et J. F. Beacom — Angular distribution of neutron inverse beta decay — cinématique et section efficace de référence pour la désintégration bêta inverse.

Déclaration de confiance

Confiance élevée dans les comptages, significativités, fenêtres d'énergie, géométrie et composantes d'incertitude reproduits depuis l'article évalué par les pairs. Confiance moyenne dans la portée applicative : la publication évoque les garanties nucléaires, mais ne teste ni détecteur compact en surface ni inférence d'inventaire. La part de 44 % attribuée aux piscines est une sortie du modèle, pas une séparation mesurée. Cette première analyse quantitative n'a pas encore de réplication indépendante. Contrôle d'unicité : aucun DOI, identifiant arXiv ni résultat central équivalent n'apparaît dans le registre A001–A018 ; verdict UNIQUE. DOI, dates, formules, liens et absence de doublon ont été vérifiés.