[142] we focused on the problem of the initialization of the electronic dynamics in trajectory-based simulations. The initial conditions for the classical-like nuclear dynamics {R α ν , P α ν } are often sampled from the (harmonic) Wigner distribution of the ground state nuclear wavepacket assuming vertical excitation or from a classical Boltzmann distribution at some given temperature 177 . However, several options can be proposed for the initialization of the electronic dynamics via the electronic coefficients {C α k (0)} when more than one electronic state needs to be populated to initiate the nonadiabatic simulations. In Ref.