These calculations provide useful and non-trivial tests to systematically investigate the numerical performance of various recently developed diabatic quantum dynamics approaches, going beyond the simple diabatic model systems that have historically been the major workhorse in the quantum dynamics field. At the same time, we hope these available benchmark studies will also foster the development of new quantum dynamics approaches. lation dynamics for the AIMS approach with 14 nuclear samples (solid blue curve) and the spin-LSC approach with 1000 nuclear samples (solid orange curve). Both of these results are duplicated from Fig. . The dotted black curve shows the spin-LSC result using the same 14 initial geometries and momenta as the AIMS result. A vast majority of the features in the population dynamics of AIMS are recovered and can be directly attributed to a lack of convergence in the nuclear sampling scheme. Since similar numbers of nuclear samplings were performed for the other three molecular models, the AIMS results are to be taken as a qualitative result for the population dynamics.