If we just consider the TDSE for independent trajectories, as in TSH, we can easily obtain the timeevolution of a set of diabatic coefficients we shall call D tsh . By inserting the electronic wavefunction expansion 26 into the TDSE (Eq. 2) we obtain Ḋtsh = -i ℏ H(t)D tsh .