  The presence of pairing correlations within the time-dependent density
functional theory (TDDFT) extension to superfluid systems, is tantamount to the
presence of a quantum collision integral in the evolution equations, which
leads to an obviously non-Markovian behavior of the single-particle occupation
probabilities, unexpected in a traditional quantum extension of kinetic
equations. The quantum generalization of the Boltzmann equation, based on a
collision integral in terms of phase-space occupation probabilities, is the
most used approach to describe nuclear dynamics and which by construction has a
Markovian character. By contrast, the extension of TDDFT to superfluid systems
has similarities with the Baym and Kadanoff kinetic formalism, which however is
formulated with much more complicated evolution equations with long-time memory
terms and non-local interactions.
  The irreversibility of quantum dynamics is properly characterized using the
canonical wave functions/natural orbitals, and the associated canonical
occupation probabilities, which provide the smallest possible representation of
any fermionic many-body wave function. In this basis, one can evaluate the
orbital entanglement entropy, which is an excellent measure of the
non-equilibrium dynamics of an isolated system.
  To explore the phenomena of memory effects and irreversibility, we
investigate the use of canonical wave functions/natural orbitals in nuclear
many-body calculations, assessing their utility for static calculations,
dynamics, and symmetry restoration. As the number of single-particle states is
generally quite large, it is highly desirable to work in the canonical basis
whenever possible, preferably with a cutoff. We show that truncating the number
of canonical wave functions can be a valid approach in the case of static
calculations, but that such a truncation is not valid for time-dependent
calculations, as it leads to the violation of continuity equation, energy
conservation, and other observables, as well as an inaccurate representation of
the dynamics. Indeed, in order to describe the dynamics of a fissioning system
within a TDDFT framework all the canonical states must be included in the
calculation. Finally, we demonstrate that the canonical representation provides
a very efficient basis for performing symmetry projections.
