Non-Markovian character and irreversibility of real-time quantum many-body dynamics

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.