Multi-trajectory Ehrenfest (MTE) is based on an uncorrelated Ansatz for the molecular wavefunction Ψ(r, R, t) χ(R, t)Φ(r, t). Therefore, TDSE-like evolution equations for each wavefunction, say the nuclear χ(R, t) (electronic Φ(r, t)), can be determined by averaging over the degrees of freedom of the other wavefunction, say the electronic (nuclear), such that the nuclei evolve in a mean-field created by the electrons and vice versa. To introduce a trajectory-based NAMD scheme, the nuclear dynamics is approximated by a swarm of independent trajectories that are generated by Newton's equation with mean-field forces determined from the instantaneous expectation value of the operator -∇ ν Ĥel (r, R(t)) on the time-dependent electronic state. Similarly to TSH, the electronic wavefunction can be expressed using a BH-like expansion, such that evolution equations for the expansion coefficients, ĊIT l (t), are determined and the mean-field force F MF (t) is expressed in terms of adiabatic PESs and nonadiabatic coupling vectors. Note that in MTE ĊIT l (t) is the same as in TSH. Being a mean-field method, MTE may yield incorrect nuclear dynamics, as for instance it is not able to capture, by construction, reaction pathways occurring along diverging pathways in nuclear configuration space. This is because the dynamics is driven by a mean potential rather than by the adiabatic PESs. Also, MTE suffers from overcoherence as TSH .