The conformational space was explored, the data was stored in the database environment and afterwards used to produce fitted PE surfaces. With these fitted surfaces so-called energy-only surface-hopping simulations, where gradients are calculated numerically from the PE surfaces, were performed. For the pyrazine model system as well as the S 1 /S 2 conical intersection of the SO 2 molecule, the energy-only simulations predicted the dynamics correctly. The proposed protocol allows to perform surface-hopping simulations using only adiabatic energies. Specifically for large molecules and electronic structure methods, where