A new approximate method for the calculation of quantum time correlation functions is derived. The new theory combines the linearized path integral approach and the mapping Hamiltonian formulation so that time correlation functions can be calculated for systems in which non-adiabatic electronic transitions are important. (c) 2005 Elsevier B.V. All rights reserved.