Schmidt J R, Tully John C
Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA.
J Chem Phys. 2007 Sep 7;127(9):094103. doi: 10.1063/1.2757170.
Within the adiabatic approximation, it is trivial to generalize existing imaginary time path-integral techniques to the case of multiple electronic surfaces. However, there are many times where nonadiabatic effects can play an important role. To this end, we reformulate the well-known path-integral expressions to incorporate multiple potential surfaces, without necessitating the adiabatic approximation. We show that the resulting expression, like its adiabatic counterpart, can be interpreted in terms of a simple classical isomorphic system and thus is amenable to simulation through Monte Carlo techniques. We derive simple expressions to compute expectation values of a general operator in both the nuclear coordinate and electronic state, and demonstrate the existence of a simple internal diagnostic that can be used to evaluate the magnitude of equilibrium nonadiabatic effects.
