Path-integral simulations beyond the adiabatic approximation

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.