Hereby, standard nonadiabatic surface-hopping dynamics simulations are used to sample the important regions of the conformational space. Subsequently, an interpolation scheme is used to predict PE surfaces and gradients based on the stored data. Using fitted properties significantly reduces the computational cost, as no electronic structure calculations are performed, which are typically time determining in direct-dynamics simulations. Finally, having global PE surfaces at hand, energy-only calculations are performed, where gradients are calculated directly from the fitted PE surfaces. This opens the possibility to perform nonadiabatic excited-state dynamic simulations, where only energies have to be calculated using electronic structure methods.