Multiple-time correlation functions for non-markovian interaction: beyond the quantum regression theorem.

We derive the dynamical equation of the reduced propagator, an object that evolves state vectors of the system conditioned to the dynamics of its environment, which is not necessarily in the vacuum state at the initial time. Such a reduced propagator is essential to obtain multiple-time correlation functions (MTCFs). We also study the evolution of MTCFs within the weak-coupling limit and show that the quantum regression theorem is, in general, not satisfied. We illustrate the theory in two different cases: first, solving an exact model, and, second, presenting the results of the numerical integration for a system coupled with a thermal environment through a nondiagonal interaction.