Our first example is a set of 2-dimensional 2-state linear vibronic coupling models which represent molecular systems with a conical intersection (CI): bis(methylene) adamantyl (BMA) radical cation , butatriene cation and pyrazine. The dynamics through the CIs in these systems was studied with ab initio multiple spawning (AIMS) and SH in Ref. [ 54], where both AIMS and SH performed well in describing the dynamics through the CI in these systems. It is has been argued 53 that the lack of geometric phase effects in methods such as AIMS (SPA0) and SH is compensated by their neglect of the diagonal BO corrections.