During the time-evolution of the system, the time-dependent modals of the Hartree product are changed, carrying the time-dependency of the TDH wave function. The one-mode density matrix in Eq. 20 reduces to unity, meaning that the TDH one-mode densities for each mode are equivalent to the norm square of the time-dependent modal function