As proof of principle for surface-hopping simulations using fitted properties as well as the energy-only algorithm, dynamics simulations of a pyrazine model are performed. To analyse the energy-only surface hopping simulations and to validate the accuracy of the fitting algorithm our results are compared with previous simulations. Moreover, the population transfer at the S 1 /S 2 conical intersection of SO 2 was simulated based on ab-initio calculations and compared with the dynamics on interpolated surfaces. The simulations followed a general protocol; First, an independent set of 100 trajectories is propagated to sample the conformational space and the data is stored in a database. Second, the data (i.e. energies and gradients) is used as training data for the fitting algorithm (for SO 2 , additional data from a grid sampling supplements the training data). Subsequently, the quality of the fitted properties has to be checked in the validation step. For this, another set of 100 trajectories, based on an independent Wigner sampling, is propagated and the data is stored. The latter datapoints are compared with the predictions of the interpolator. Finally, the fitted properties are used to perform dynamics simulations with 1000 trajectories (production runs), which are compared with reference calculations. The reference calculations are a set of 1000 trajectories, starting with the same initial conditions as the production calculations, but taking energies and gradients directly from the model or ab-initio calculations without interpolation.