Hybrid quantum-mechanics/molecular-mechanics (QM/MM) potentials are the state-of-the-art approach for simulating chemical transformations in condensed phases. By combining a wavefunction or density function theory (DFT) based Hamiltonian for a small portion of the system (the QM region), while describing the rest with a computationally efficient molecular mechanics forcefield, one can consider simulation systems consisting of >10 atoms while still being able to describe electronic rearrangements and bond formation/breaking in the region of interest. However, the performance bottleneck of QM/MM simulations still lies in the QM method of choice. Thus, even for relatively small QM regions (≈10 2 atoms) and using a relatively efficient (non-hybrid) DFT functional, the timescale of QM/MM molecular dynamics (MD) simulations is typically limited to nanoseconds at best.