Multi-configurational Ehrenfest simulations of ultrafast nonadiabatic dynamics in a charge-transfer complex

Multi-configurational Ehrenfest (MCE) approaches, which are intended to remedy the lack of correlations in the standard mean-field Ehrenfest method, have been proposed as coherent-state based ansatze for quantum propagation [D. V. Shalashilin, J. Chem. Phys. 130, 244101 (2009)] and also as the classical limit of the variational Gaussian-based multiconfiguration time dependent Hartree (G-MCTDH) method [S. Romer and I. Burghardt, Mol. Phys. 111, 3618 (2013)]. In the present paper, we establish the formal connection between these schemes and assess the performance of MCE for a coherent-state representation of the classical-limit subsystem. As a representative model system, we address the ultrafast, coherent charge transfer dynamics in an oligothiophene-fullerene donor acceptor complex described by a two-state linear vibronic coupling model. MCE calculations are compared with reference calculations performed with the MC IDH method, for 10-40 vibrational modes. Beyond a dimensionality of 10 modes, it is shown that the correct representation of electronic coherence depends crucially on the sampling of initially unoccupied Gaussians.