We note here that we choose the time interval of the calculation of the one-mode densities t dens such that new one-mode densities are always calculated after the same time-step as the TD-ADGA iterations are also carried out. This is not a necessity for the algorithm but it makes the most sense to use information from the wave function at the current time in the TD-ADGA iterations.