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

Cette étude utilise la théorie de la fonctionnelle de la densité dépendante du temps pour démontrer que les neutrons de scission constituent une composante significative à haute énergie du spectre de neutrons de fission prompt dans la fission de l' $^{235}\mathrm{U}$, de l' $^{239}\mathrm{Pu}$ et du $^{252}\mathrm{Cf}$, fournissant ainsi la preuve directe que leur inclusion résout la sous-estimation systématique des rendements à haute énergie observée dans les modèles traditionnels basés uniquement sur l'évaporation.

L'essentiel

Imaginez un ballon géant et instable (un noyau atomique) qui se déchire soudainement en deux. C'est la fission nucléaire. Depuis des décennies, les scientifiques savent que lorsqu'un tel événement se produit, les deux morceaux résultants (appelés fragments) s'écartent à des vitesses incroyables et recrachent de minuscules particules appelées neutrons.

Pendant longtemps, les physiciens ont pensé que tous ces neutrons étaient « évaporés » plus tard, comme la vapeur s'élevant d'une tasse de café chaud après que l'eau a bouilli. Ils supposaient que les fragments étaient pleinement formés et se déplaçaient de manière constante avant de commencer à lâcher ces neutrons.

Cependant, ce nouvel article suggère que certains neutrons sont en réalité « expulsés » au moment même où le ballon se déchire. Ce sont des neutrons de scission. Ils naissent dans le chaos de la rupture, dans l'instant précis de la séparation, et non d'un fragment calme et en train de refroidir plus tard.

Voici comment les chercheurs ont trouvé la preuve de ces neutrons du « moment de la rupture », expliquée simplement :

1. La simulation par supercalculateur

Pour observer ce qui se passe lors de la division, les scientifiques n'ont pas utilisé de microscope ; ils ont utilisé un supercalculateur pour lancer un film de l'événement en utilisant une théorie appelée Théorie de la Fonctionnelle de la Densité Dépendante du Temps (TDDFT).

Considérez cela comme un jeu vidéo 3D à haute vitesse où ils simulent la danse des atomes et leur décomposition. Dans les versions précédentes de ce « jeu », le monde virtuel était trop petit. Les neutrons heurtaient le bord de l'écran avant que les scientifiques ne puissent déterminer exactement à quelle vitesse ils allaient ou dans quelle direction ils volaient.

Dans cette étude, ils ont construit un monde virtuel beaucoup plus grand (environ 3 fois plus grand qu'auparavant). Cela a donné aux neutrons assez de place pour s'envoler et se stabiliser afin que les scientifiques puissent les mesurer avec précision sans que les « murs » de la simulation ne faussent les données.

2. La découverte de la « limite de vitesse »

Une fois qu'ils ont eu une vue claire, ils ont observé les neutrons s'échappant selon des angles spécifiques (principalement sur les côtés et légèrement vers l'arrière par rapport à la division). Ils ont découvert quelque chose de surprenant :

• La zone d'interdiction : Il n'y a aucun neutron de scission à basse énergie (en dessous d'environ 1,5 à 2 millions d'électron-volts). C'est comme s'il y avait une limite de vitesse ; rien de lent n'est autorisé à être un neutron de « moment de la rupture ».
• La foule à haute vitesse : Au lieu de cela, ces neutrons sont tous rapides. Ils se regroupent autour d'une vitesse spécifique élevée (3–3,5 MeV) puis s'étirent en une longue traîne de particules encore plus rapides.

C'est comme une foule de gens sautant d'un plongeoir. Les neutrons « évaporés » sont comme des gens descendant calmement du bord de la piscine plus tard. Les neutrons de « scission » sont comme des gens violemment projetés du plongeoir à la seconde exacte où celui-ci se brise. Ceux qui sont projetés du plongeoir sont toujours en mouvement rapide ; on ne voit jamais un mouvement lent provenant de cet événement spécifique.

3. Résoudre le mystère de l'énergie « manquante »

Les scientifiques essaient depuis des années de faire correspondre leurs modèles informatiques avec les expériences réelles. Ils avaient un problème :

• L'ancien modèle : Si l'on ne compte que la « vapeur » (les neutrons évaporés), votre modèle informatique prédit trop peu de neutrons à haute énergie. C'est comme essayer de remplir un seau avec une petite tasse, mais le seceau semble toujours avoir besoin de plus d'eau que ce que la tasse peut fournir.
• Le nouveau modèle : Lorsque les chercheurs ont ajouté les neutrons du « moment de la rupture » (ceux qu'ils ont trouvés dans leur grande simulation) aux neutrons de la « vapeur », le calcul a enfin fonctionné. Le modèle combiné correspond parfaitement aux données de haute énergie mesurées lors d'expériences réelles sur l'Uranium et le Californium.

4. Pourquoi c'est important

C'est un événement majeur car c'est la première fois qu'une théorie purement microscopique (qui ne se contente pas de deviner ou de supposer l'existence de choses) a prédit ces neutrons du « moment de la rupture » et a prouvé qu'ils sont réels.

• Avant : Les scientifiques devaient supposer que ces neutrons existaient parce que les calculs ne concordaient pas.
• Maintenant : La simulation informatique a naturellement produit ces neutrons sans qu'on lui dise d'exister. C'est comme prédire une tempête en observant le mouvement des nuages, plutôt que de simplement supposer qu'une tempête arrive parce que le bulletin météo le dit.

L'essentiel

L'article conclut que lorsqu'un atome se divise, une petite mais importante partie des neutrons (environ 6 % à 10 % du total) naît dans le moment violent de la rupture. Ces neutrons sont distincts car ils sont toujours rapides et jamais lents selon certains angles.

En trouvant cette « empreinte digitale » dans les données, les chercheurs ont enfin réussi à séparer les neutrons du « moment de la rupture » des neutrons de la « vapeur », offrant ainsi une image plus claire et plus précise de la façon dont la fission nucléaire fonctionne réellement. Cela aide à affiner notre compréhension des forces fondamentales qui maintiennent la matière ensemble et la déchirent.

Résumé technique

Résumé technique : Propriétés angulaires et cinétiques des neutrons de scission au sein de la théorie de la densité fonctionnelle dépendante du temps

Énoncé du problème
Malgré la découverte de la fission nucléaire il y a près d'un siècle, une description microscopique complète du processus demeure élusive, particulièrement en ce qui concerne la transition du point de selle extérieur vers la scission et la formation subséquente des fragments de fission (FF) primaires. Les codes d'évaporation statistique actuels supposent que tous les neutrons prompts sont émis isotropement à partir de fragments pleinement accélérés. Cependant, ce cadre ne peut rendre compte des neutrons émis pendant le processus dynamique hautement hors équilibre de la rupture du col, connu sous le nom de neutrons de scission (SN). Bien que les SN aient été conjecturés depuis 1939 et étudiés via divers modèles, ils restent non résolus expérimentalement car les observables de neutrons prompts seuls ne permettent pas de distinguer les neutrons émis lors de la scission de ceux évaporés ultérieurement. Les tentatives précédentes pour prédire les SN en utilisant la théorie de la densité fonctionnelle dépendante du temps (TDDFT) étaient limitées par de petits domaines de simulation, empêchant la détermination simultanée des propriétés angulaires et cinétiques avant que les neutrons n'atteignent les limites de calcul.

Méthodologie
Cette étude emploie un cadre entièrement microscopique basé sur la TDDFT étendue aux systèmes superfluides pour étudier l'émission de SN dans $^{235}\text{U}(n_{th}, f)$, $^{239}\text{Pu}(n_{th}, f)$ et $^{252}\text{Cf}(sf)$.

• Configuration de simulation : Les auteurs ont utilisé le solveur d'oscillateur harmonique à symétrie axiale AxialHOHFB avec la fonctionnelle de densité d'énergie SkM*. Une avancée critique a été l'utilisation d'un domaine de simulation nettement plus grand (un ellipsoïde avec des demi-axes $Z_{max} = 70$ fm et $R_{max} = 55$ fm), environ 3,2 fois plus grand que dans les études précédentes.
• Intégration en deux étapes : Pour gérer les contraintes de calcul tout en maintenant la précision, une procédure d'intégration temporelle en deux étapes a été employée. Le système est évolué de la solution statique de Hartree-Fock-Bogoliubov (HFB) contrainte vers le moment de scission dans une base plus petite ($B_1$), puis transformé vers la base plus grande ($B_2$) pour poursuivre l'évolution jusqu'à ce que les neutrons approchent de la limite.
• Extraction des propriétés des SN : Les SN ont été identifiés en définissant une région $V^{(SN)}_\eta$ où le rapport de la densité d'énergie de flux collectif sur la densité d'énergie cinétique intrinsèque dépasse un seuil $\eta = 3$. Ce critère isole les neutrons dominés par un flux collectif sortant.
• Analyse angulaire et cinétique : À l'intérieur de cette région, les champs de vitesse locale et l'énergie de flux collectif par neutron ont été calculés. L'étude s'est concentrée sur l'intervalle angulaire $110^\circ \leq \theta \leq 150^\circ$ (mesuré par rapport à la direction du fragment léger), où le rendement des SN sature avant que les réflexions de bord ne surviennent.
• Comparaison avec l'expérience : Le spectre d'énergie cinétique des SN calculé a été combiné avec un modèle d'évaporation maxwellien pour la composante de neutron évaporé (EN). Les paramètres du modèle d'évaporation (température du fragment $T$ et rapport de multiplicité $\bar{\nu}_L/\bar{\nu}_H$) ont été contraints exclusivement par les données expérimentales de basse énergie ($E \leq E_{thr}$). Ce modèle combiné a ensuite été comparé aux données de spectre de neutrons de fission prompt (PFNS) à haute énergie pour $^{239}\text{Pu}$ et $^{252}\text{Cf}$.

Résultats clés

• Caractéristiques du spectre des SN : Les spectres de SN extraits pour les trois systèmes présentent un comportement de seuil distinct. Les SN sont absents en dessous d'un seuil d'énergie spécifique ($E_{thr}$), qui est d'environ 1,5–2 MeV (spécifiquement $\approx 2$ MeV pour la fission induite et $\approx 1,5$ MeV pour la fission spontanée). Les spectres culminent près de 3–3,5 MeV et présentent une queue de haute énergie exponentielle s'étendant jusqu'à environ 18 MeV.
• Estimations de rendement : Pour $^{239}\text{Pu}(n_{th}, f)$, le rendement total de SN est calculé à $0,19 \pm 0,01$ (environ 6,6 % des neutrons prompts totaux). Pour $^{252}\text{Cf}(sf)$, le rendement est de $0,39 \pm 0,02$ (environ 10,4 %).
• Validation du modèle : Lorsque la composante SN calculée est ajoutée au modèle d'évaporation (contraint uniquement par les données de basse énergie), le PFNS total résultant reproduit avec précision le rendement mesuré à haute énergie pour $^{239}\text{Pu}$ et $^{252}\text{Cf}$. En revanche, le modèle d'évaporation seule sous-estime systématiquement le rendement à haute énergie, même dans l'intervalle de confiance de 95 %.
• Stabilité angulaire : L'étude confirme que des propriétés asymptotiques fiables ne peuvent être extraites que dans la plage angulaire $110^\circ \leq \theta \leq 150^\circ$. En dehors de cette plage, les effets d'interférence transitoires et les réflexions de bord empêchent le rendement de saturer dans le domaine de simulation actuel.

Signification et affirmations
L'article affirme fournir la première identification d'une signature de neutron de scission au sein des spectres de neutrons de fission prompt mesurés existants en utilisant un cadre entièrement microscopique. En démontrant que la TDDFT produit naturellement des SN sans hypothèses a priori, et que ces SN résolvent la divergence entre les modèles d'évaporation seule et les données expérimentales à haute énergie, l'étude apporte une preuve directe pour une composante de neutrons de scission non négligeable dans la fission prompt.

Les auteurs postulent que ce travail établit un fondement microscopique pour l'interprétation des spectres de neutrons de fission au-delà du schéma statistique d'évaporation standard. De plus, ils identifient une signature expérimentale spécifique : la séparation des contributions SN et EN en dessous du seuil d'énergie ($E_{thr}$). Cela implique que les mesures de basse énergie dans la plage angulaire spécifiée fournissent les contraintes les plus robustes pour la composante d'évaporation, tandis que la queue de haute énergie sert d'indicateur primaire de la contribution SN. L'étude conclut que, bien que les simulations actuelles soient limitées à des plages angulaires spécifiques, la méthodologie isole avec succès une signature mesurable des neutrons de scission qui a été théoriquement conjecturée depuis près d'un siècle.