The adiabatic approximation in time-dependent density matrix functional theory: response properties from dynamics of phase-including natural orbitals.