Once a satisfactory PES has been obtained, based on the usual precision metrics, it important to go beyond these metrics. One way to do this is via molecular dynamics (MD) calculations, which we did using our own software. For the fragmentation method we used F-PES, and we used MB-PES for the many-body approach. In all applications trajectories were initiated at the global minimum of the potential and the total energy was specified and distributed uniformly among the kinetic energies of all atoms. The total angular momentum was adjusted to be zero for each trajectory.