Can we derive Tully's surface-hopping algorithm from the semiclassical quantum Liouville equation? Almost, but only with decoherence We implemented a method for the treatment of field induced transitions in trajectory surface hopping simulations, in the framework of the local diabatization scheme, especially suited for on-the-fly dynamics. The method is applied to a simple one-dimensional model with an avoided crossing and compared with quantum wavepacket dynamics. The results show the

more »

... e of introducing a proper decoherence correction to surface hopping, in order to obtain meaningful results. Also the energy conservation policy of standard surface hopping must be revised: in fact, the quantum wavepacket energetics is well reproduced if energy absorption/emission is allowed for in the hops determined by radiation-molecule coupling. To our knowledge, this is the first time the issues of decoherence and energy conservation have been analyzed in depth to devise a mixed quantum-classical method for dynamics with molecule-field interactions. © 2014 AIP Publishing LLC. [http://dx.