In Eq. ( ) we have introduced the electronic wavefunction, ψ L , for an arbitrary adiabatic state, L, which depends on the electronic coordinates r. The parametric dependence of the electronic wavefunction on the nuclear geometry, R, is also given. Generally, a single BOMD trajectory will thus propagate the classical nuclear coordinates on a single PES corresponding to a specific adiabatic electronic state.