and the nuclear coordinates evolve according to classical equations of motion with the above force. More details on the derivation and implementation of the MFE approach can be found in the literature. The MFE approach can also be equivalently written down in terms of the reduced density matrix of the quantum subsystem (Eq. 8) as follows