Non-adiabatic dynamics around a conical intersection with surface-hopping coupled coherent states.

Abstract

A surface-hopping extension of the coupled coherent states-method [D. Shalashilin and M. Child, Chem. Phys. 304, 103-120 (2004)] for simulating non-adiabatic dynamics with quantum effects of the nuclei is put forward. The time-dependent Schrödinger equation for the motion of the nuclei is solved in a moving basis set. The basis set is guided by classical trajectories, which can hop stochastically between different electronic potential energy surfaces. The non-adiabatic transitions are modelled by a modified version of Tully's fewest switches algorithm. The trajectories consist of Gaussians in the phase space of the nuclei (coherent states) combined with amplitudes for an electronic wave function. The time-dependent matrix elements between different coherent states determine the amplitude of each trajectory in the total multistate wave function; the diagonal matrix elements determine the hopping probabilities and gradients. In this way, both interference effects and non-adiabatic transitions can be described in a very compact fashion, leading to the exact solution if convergence with respect to the number of trajectories is achieved and the potential energy surfaces are known globally. The method is tested on a 2D model for a conical intersection [A. Ferretti, J. Chem. Phys. 104, 5517 (1996)], where a nuclear wavepacket encircles the point of degeneracy between two potential energy surfaces and interferes with itself. These interference effects are absent in classical trajectory-based molecular dynamics but can be fully incorpo rated if trajectories are replaced by surface hopping coupled coherent states.