Abstract

The bonding and antibonding character of individual molecular orbitals has been previously shown to be related to their orbital energy derivatives with respect to nuclear coordinates, known as dynamical orbital forces. Albeit usually derived from Koopmans' theorem, in this work we show a more general derivation from conceptual DFT, which justifies application in a broader context. The consistency of the approach is validated numerically for valence orbitals in Kohn–Sham DFT. Then, we illustrate its usefulness by showcasing applications in aromatic and antiaromatic systems and in excited state chemistry. Overall, dynamical orbital forces can be used to interpret the results of routine ab initio calculations, be it wavefunction or density based, in terms of forces and occupations.