In all dynamics calculations a primitive B-spline basis has been employed using B-spline functions of order 10 and a basis density of 0.8. The boundaries of the basis sets used were determined by the reference PESs for each system and used in all computations for consistency in comparisons. These boundaries can be found in the supporting information (SI).