Recent work, particularly by Cederbaum and co-workers, has identified the phenomenon of charge migration, whereby charge flow occurs over a static molecular framework after the creation of an electronic wavepacket. In a real molecule, this charge migration competes with charge transfer, whereby the nuclear motion also results in the re-distribution of charge. To study this competition, quantum dynamics simulations need to be performed. To break the exponential scaling of standard grid-based algorithms, approximate methods need to be developed that are efficient yet able to follow the coupled electron-nuclear motion of these systems. Using a simple model Hamiltonian based on the ionisation of the allene molecule, the performance of different methods based on Gaussian Wavepackets is demonstrated.