For the computational modeling of complex chemical reactions, Kohn-Sham density functional theory (KS-DFT) 1 is one of the most widely used electronic structure methods. However, due to its high computational cost, KS-DFT is practically out of reach for larger systems and molecular dynamics on long timescales. In some cases, atomic potentials or force fields offer a low-cost alternative, but the lack of explicit consideration of electronic structure makes them unsuitable for studying chemical reactions or different electronic states.