Not yet reviewed by Pith; the record is open.

This work combines for the first time ab initio molecular dynamics (AIMD) within the Born-Oppenheimer approximation, with a global natural orbital functional (GNOF), an approximate functional of the one-particle reduced density matrix. The most prominent feature of GNOF-AIMD is the ability to display the real-time evolution of natural orbitals, providing detailed information on the time-dependent electronic structure of complex systems and processes, including reactive collisions. The quartet ground-state reaction N($^4$S) + H$_2$($^1\Sigma$) $\rightarrow$ NH($^3\Sigma$) + H($^2$S) is taken as validation test. Collision energy influences on integral cross sections for different initial ro-vibrational states of H$_2$ and rotational-state distributions of NH product are discussed, showing a good agreement with previous high-quality theoretical results.