A variety of theoretical approaches of different complexity and accuracy are available for the nonadiabatic excited-state wavepacket propagation in complex molecular systems. On the one hand, a fully quantum evaluation of the Hamiltonian dynamics or the dissipative dynamics is certainly possible, but these methods are suitable only for the excitonic or linearvibronic-coupling Hamiltonians. On the other hand, great efforts have been devoted to development of mixed quantum-classical trajectory-based or Gaussian-wavepacket-based propagation schemes . In the commonly used trajectory surface hopping (TSH) method , for example, the nuclear propagation is governed by classical mechanics and the electronic evolution is described by quantum mechanics. These trajectory-based approaches are especially attractive: Firstly, they usually require only the local information on the potential energy surface along the trajectory.