We study the coupled electronic-nuclear dynamics in a model system to compare numerically exact calculations of electronic and nuclear flux densities with those obtained from the Born-Oppenheimer (BO) approximation. Within the adiabatic expansion of the total wave function, we identify the terms which contribute to the flux densities. It is found that only off-diagonal elements that involve the interaction between different electronic states contribute to the electronic flux whereas in the nuclear case the major contribution belongs to the BO electronic state. New flux densities are introduced where in both, the electronic and the nuclear case, the main contribution is contained in the component corresponding to the BO state. As a consequence, they can be determined within the BO approximation, and an excellent agreement with the exact results is found. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Keywords: Born–Oppenheimer approximation; Ehrenfest theorem; flux densities.