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

Cette étude résout les équations de Langevin dans un modèle multidimensionnel pour simuler l'évolution temporelle de la fission nucléaire, en intégrant l'émission de neutrons à chaque étape afin de calculer la multiplicité pré-scission et d'autres observables comparables aux données expérimentales.

L'essentiel

Le Voyage d'un Atome qui se Fend : Une Histoire de Langevin et de Neutrons

Imaginez que vous observez un atome lourd (comme l'uranium) qui est sur le point de se briser en deux. C'est ce qu'on appelle la fission nucléaire. Mais ce n'est pas une simple cassure instantanée comme un verre qui tombe. C'est un voyage lent et complexe, un peu comme une goutte d'eau qui s'étire, tremble, et finit par se séparer en deux.

Les auteurs de ce papier, F.A. Ivanyuk et ses collègues, ont créé un simulateur très précis pour raconter cette histoire. Voici comment ils ont fait, expliqué simplement :

1. Le Modèle : Une Balle dans un Paysage de Collines

Pour comprendre comment l'atome se déforme, les scientifiques utilisent une équation appelée l'équation de Langevin.

• L'analogie : Imaginez une bille lourde roulant sur un paysage montagneux.
  • Les collines représentent l'énergie qui maintient l'atome ensemble.
  • Les vallées sont les formes que l'atome peut prendre.
  • La bille, c'est l'atome lui-même.
  • Le vent et les secousses aléatoires (représentés par le "bruit" dans l'équation) poussent la bille, la faisant osciller et avancer vers le bas de la montagne jusqu'au point de rupture (le "point de scission").

Jusqu'à présent, les modèles étaient comme des cartes statiques. Ils disaient : "Voici le chemin, l'atome va tomber ici." Mais ils ne disaient pas ce qui se passait pendant le voyage.

2. Le Problème : Les "Éclats" de la Balle (Les Neutrons)

Pendant que l'atome roule vers sa destruction, il est très chaud (excité). Comme une casserole d'eau bouillante, il perd de la chaleur en émettant de la vapeur. Dans le monde nucléaire, cette "vapeur", ce sont des neutrons.

Le problème, c'est que les anciens modèles ignoraient souvent ces pertes de neutrons pendant le voyage. Ils pensaient que l'atome gardait toute son énergie jusqu'à la fin. Or, chaque fois qu'un neutron s'échappe, l'atome devient plus léger et moins chaud. Cela change sa trajectoire sur la carte montagneuse !

3. La Nouvelle Méthode : Un Guide de Voyage en Temps Réel

Dans ce papier, les auteurs ont ajouté une règle intelligente à leur simulation : "À chaque instant, vérifiez si un neutron s'échappe."

Voici comment leur algorithme fonctionne, étape par étape :

1. Le pas de temps : À chaque micro-seconde du voyage de la bille (l'atome), le programme se demande : "Est-ce qu'un neutron va sauter hors de la bille maintenant ?"
2. Le calcul de probabilité : Ils utilisent une formule mathématique (basée sur la physique des gaz et la mécanique quantique) pour estimer cette chance. C'est comme regarder la température de la casserole pour prédire si une bulle va éclater.
3. L'action : Si la réponse est "OUI" :
   • Ils retirent un neutron de la bille.
   • Ils refroidissent la bille (on retire l'énergie nécessaire pour arracher le neutron).
   • Ils changent la carte : la bille roule maintenant sur un terrain légèrement différent (car elle a un neutron de moins).
4. La suite : Le voyage continue avec la nouvelle configuration jusqu'à ce que l'atome se brise.

4. Les Découvertes Surprenantes

En faisant cela, ils ont découvert des choses intéressantes sur le timing de l'émission des neutrons :

• À basse température (10 MeV) : L'atome est si froid au début qu'il ne perd aucun neutron tant qu'il n'a pas franchi la crête de la montagne (la barrière de fission). Les neutrons ne partent que lorsque l'atome est déjà bien engagé dans sa chute, juste avant de se briser.
• À haute température (40 MeV) : L'atome est si chaud qu'il perd des neutrons dès le début, même au fond de la vallée (près de sa forme normale).
• Le résultat final : En tenant compte de ces pertes, la forme des deux morceaux qui se séparent (les fragments de fission) correspond beaucoup mieux à ce que l'on observe dans les expériences réelles. C'est comme si, en tenant compte du poids perdu par la bille pendant sa course, on pouvait prédire exactement où elle atterrirait.

5. Pourquoi c'est Important ?

Imaginez que vous essayez de prédire la météo. Si vous ignorez l'évaporation de l'océan, votre prévision sera fausse. De même, pour comprendre la fission nucléaire (cruciale pour les réacteurs ou la sécurité des déchets), il faut savoir exactement combien de neutrons sont émis et à quel moment.

Ce papier montre que si l'on simule le voyage de l'atome en temps réel, en comptant chaque neutron perdu comme un passager qui descend du bus, on obtient une image beaucoup plus fidèle de la réalité.

En résumé : Les auteurs ont transformé une simulation statique en un film dynamique où l'atome perd du poids et de la chaleur en cours de route. Grâce à cette approche, leurs prédictions sur la façon dont l'atome se brise et sur les neutrons qu'il libère correspondent parfaitement à la réalité observée en laboratoire.

Résumé technique

1. Problématique et Contexte

La fission nucléaire est un processus complexe impliquant l'évolution dynamique de la forme du noyau depuis l'état composé jusqu'au point de scission. Bien que l'approche basée sur les équations de Langevin ait été largement utilisée pour décrire les distributions de masse et les énergies cinétiques des fragments, la plupart des études antérieures se sont concentrées sur la « fission à première chance » (first chance fission), négligeant souvent l'évaporation de particules légères (neutrons, protons, alpha) et de rayons gamma durant l'évolution vers la scission.

Le défi principal réside dans la modélisation simultanée de l'évolution de la forme du noyau et de l'émission neutronique pré-scission. Les modèles existants fournissent souvent des probabilités d'émission pour l'ensemble du processus de fission, ce qui est incompatible avec la nature temporelle et détaillée de la simulation de Langevin, où la trajectoire peut ne jamais atteindre la scission (retour à l'état fondamental) ou être abandonnée. Il est donc nécessaire de développer une méthode permettant de décider, à chaque pas d'intégration temporelle, si un neutron est émis, et d'ajuster dynamiquement l'énergie d'excitation et la surface d'énergie potentielle en conséquence.

2. Méthodologie

Les auteurs ont développé une extension du modèle de Langevin multi-dimensionnel (réduit ici à 4 dimensions pour des raisons de temps de calcul) pour inclure l'émission neutronique en temps réel.

• Équations de Langevin : Le système est décrit par l'évolution temporelle de variables collectives $q_\mu$ (paramètres de forme du modèle à deux centres, TCSM) et de leurs moments conjugués $p_\mu$. Les équations incluent les forces de frottement, l'inertie, et une force aléatoire thermique.
• Surface d'énergie potentielle : L'énergie libre $F(q, T)$ est calculée via un modèle macroscopique-microscopique (correction de coquille de Strutinsky + modèle de Yukawa plié). Les effets de coquille sont atténués avec l'augmentation de l'énergie d'excitation $E^*$.
• Modélisation de l'émission neutronique ($dN/dt$) :
  • Une approximation analytique du taux d'émission instantané $dN/dt$ a été dérivée en combinant l'équation de continuité et le modèle du gaz de Fermi.
  • Le taux dépend de la surface intégrale de la forme du noyau, de la température locale, et du potentiel chimique $\mu$ (énergie du dernier état occupé).
  • La formule finale simplifiée (Eq. 29) permet un calcul rapide sans pénalité excessive sur le temps de calcul.
• Algorithme d'intégration couplée :
  1. À chaque pas de temps $\Delta t$, le taux d'émission est calculé.
  2. Un nombre aléatoire détermine si un neutron est émis.
  3. Si oui :
     • L'énergie d'excitation locale est réduite de l'énergie de séparation neutronique ($S_n = -\mu$) plus l'énergie cinétique moyenne du neutron émis.
     • Le système bascule vers la couche de la surface d'énergie potentielle correspondant à un nombre de neutrons inférieur ($A-1$).
     • L'intégration des équations de Langevin continue sur cette nouvelle surface.
  4. Le processus s'arrête lorsque le rayon du col (neck) atteint une valeur critique ($r_{neck} < 2$ fm), définissant le point de scission.

3. Contributions Clés

• Couplage dynamique temps-réel : Contrairement aux approches statiques, cette méthode permet de suivre l'émission neutronique à chaque étape de la trajectoire de fission, offrant une résolution temporelle fine du processus.
• Approximation analytique du taux d'émission : Développement d'une formule pour $dN/dt$ basée sur le gaz de Fermi, adaptée aux conditions locales (déformation, température) sans dépendre de barrières de fission globales ou de sections efficaces inverses.
• Étude de la dépendance en énergie d'excitation : Application du modèle à la fission de l'uranium-236 ($^{236}\text{U}$) pour une gamme d'énergies d'excitation de 10 à 45 MeV.
• Analyse de la phase d'émission : Identification précise de la déformation (étirement $R_{12}$) à laquelle les neutrons sont émis, distinguant les émissions depuis le puits de potentiel fondamental, la barrière de fission, et la région post-barrière.

4. Résultats Principaux

Les calculs ont été effectués pour $^{236}\text{U}$ et comparés aux données expérimentales disponibles.

• Distributions de masse des fragments :
  • L'inclusion de l'émission neutronique améliore significativement l'accord entre les prédictions théoriques et les données expérimentales.
  • L'émission de neutrons renforce l'impact des effets de coquille sur les rendements de fission, en particulier à des énergies d'excitation modérées.
  • L'amortissement des effets de coquille avec l'augmentation de $E^*$ est correctement reproduit.
• Émission neutronique selon la phase de fission :
  • À basse énergie ($E^* = 10$ MeV), aucune émission n'est observée dans la région du fond du puits de potentiel ; tous les neutrons sont émis après le franchissement de la barrière de fission. Si un neutron était émis trop tôt à cette énergie, l'énergie résiduelle serait insuffisante pour franchir la barrière, menant à un résidu d'évaporation plutôt qu'à une fission.
  • À des énergies plus élevées ($E^* \ge 20$ MeV), une fraction de neutrons est émise depuis la région du fond du puits.
  • La majorité des neutrons pré-scission sont émis dans la région de grande déformation, entre le point selle et le point de scission.
• Spectres et multiplicité neutronique :
  • Les spectres d'énergie des neutrons pré-scission calculés correspondent raisonnablement bien à une distribution de Watt (standard pour les neutrons de fission).
  • La multiplicité neutronique pré-scission ($M_{pre}$) calculée est en bon accord avec les rares données expérimentales disponibles (issues de réactions $\alpha + ^{232}\text{Th}$ et $p + ^{236,238}\text{U}$).

5. Signification et Conclusion

Ce travail démontre que l'intégration explicite de l'émission neutronique dans le cadre dynamique de Langevin est non seulement possible mais essentielle pour une description précise de la fission nucléaire.

• Précision accrue : Le modèle corrige les écarts entre théorie et expérience observés dans les modèles statiques, notamment pour les distributions de masse.
• Compréhension physique : Il fournit des informations inédites sur quand et où les neutrons sont émis durant le processus de fission, confirmant que l'émission pré-scission est dominée par la région de grande déformation, sauf à très haute énergie d'excitation.
• Validité du modèle : La capacité du modèle à reproduire simultanément les distributions de masse, les spectres neutroniques et les multiplicités valide l'approche analytique utilisée pour le taux d'émission et la prise en compte des effets de coquille dépendants de la température.

En résumé, cette étude établit un cadre robuste pour simuler la fission en tenant compte de la désexcitation compétitive par émission de neutrons, ouvrant la voie à des prédictions plus fiables pour des noyaux exotiques ou des conditions extrêmes.