Abstract
This article describes the Ehrenfest method and our second-order implementation (with approximate gradient and Hessian) within a CASSCF formalism. We demonstrate that the second-order implementation with the predictor-corrector integration method improves the accuracy of the simulation significantly in terms of energy conservation. Although the method is general and can be used to study any coupled electron-nuclear dynamics, we apply it to investigate charge migration upon ionization of small organic molecules, focusing on benzene cation. Using this approach, we can study the evolution of a non-stationary electronic wavefunction for fixed atomic nuclei, and where the nuclei are allowed to move, to investigate the interplay between them for the first time. Analysis methods for the interpretation of the electronic and nuclear dynamics are suggested: we monitor the electronic dynamics by calculating the spin density of the system as a function of time.