“…Alternatively, the potential energy surfaces can be calculated on-the-fly and coupled with a fully quantum treatment of the vibrations in methods such as variational multi-configurational Gaussian, , ab initio multiple spawning, multi-configurational Ehrenfest, , etc. More approximate approaches include semiclassical methods that propagate frozen or thawed Gaussian wavepackets on-the-fly. − Instead of obtaining fitted potential energy surfaces or calculating them on-the-fly, the potential energy surfaces can be represented by local quadratic approximations, provided there are no complications such as surface crossings, double wells, or large anharmonicities. Arnold, Vendrell, and Santra have found this approach satisfactory for modeling decoherence in water, paraxylene, and phenylalanine. , Yang and co-workers − examined this approach in detail for the nuclear dynamics of a coherent superposition of the X and A states of iodoacetylene cation and successfully modeled the decoherence of charge migration.…”