Le code FELIX (Finite Element soLver for fissIon dynamiX) résout une équation de Schrödinger locale et dépendant du temps dans des espaces possédant un petit nombre de dimensions. $$ i \hbar \frac{\partial}{\partial t} g(\mathbf{q},t)= \left[ -\frac{\hbar^2}{2} \sum_{kl} \frac{1}{\sqrt{\gamma(\mathbf{q})}} \frac{\partial}{\partial q_{k}} \sqrt{\gamma(\mathbf{q})} B_{kl}(\mathbf{q}) \frac{\partial}{\partial q_{l}} + V(\mathbf{q}) \right] \cdot g(\mathbf{q},t) $$

• $ g(\mathbf{q,t}) $ est l’inconnue et correspond à une fonction complexe des variables $\mathbf{q}$ et du temps $t$.
• Les quantités $\gamma(\mathbf{q}), B(\mathbf{q})$ et $V(\mathbf{q})$ s’identifient à une métrique, un tenseur d’inertie et un potentiel. Elles définissent le problème à résoudre.

Une telle équation se présente dans le contexte de la simulation de la dynamique de la fission nucléaire. Les variables $\mathbf{q}$ correspondent alors à des déformations arbitraires du noyau et l’on cherche l’évolution temporelle d’un paquet d’onde dans l’espace des déformations possibles. Ce code se base sur une discretisation en éléments finis spectraux de l’espace des $\mathbf{q}$.

FELIX-2.0 est un code open source, et vous pouvez le télécharger depuis la bibliothèque de codes du journal Computer Physics Communication.