where \(\rho (t)\)  is the time-dependent one-body density matrix and \(h[\rho (t)]\) 
is the time-dependent mean field, which is a functional of the time-dependent
density matrix. Although the TDHF equation is much simpler to solve numerically
and was, historically, the first version of TDDFT applied to fission studies
[1]}, [2]}, [3]}, it has strong limitations: unless the
calculation is initialized rather close to the scission point, or with an
artificial boost in energy, the time-evolution of fissile nuclei such as
\(^{240}\) Pu does not lead to fission [4]}, [5]}. In order
to achieve actual scission, pairing correlations play an essential role
[6]}, [7]}. For this reason, it is preferable to use either
the TDHF+BCS theory [8]}, [9]} or the fully-fledged TDHFB
formalism. The TDHFB equation is formally similar to (REF ) as it
reads
\(i\hbar \frac{\partial \mathcal {R}}{\partial t} = \big [ \mathcal {H}(t), \mathcal {R}(t) \big ],\) 
