Calculation of reactive flux correlation functions for systems in a condensed phase environment: a multilayer multiconfiguration time-dependent Hartree approach.

A numerically exact quantum mechanical approach is proposed to evaluate thermal rate constants for systems in a model condensed phase environment. Employing the reactive flux correlation function formalism, the approach efficiently combines the multilayer multiconfiguration time-dependent Hartree theory with an importance sampling scheme for thermal distribution of the initial states. The performance of the method is illustrated by applications to two models of condensed phase dynamics: the donor-acceptor electron transfer model also known as the spin-boson model and a model for proton transfer reactions in the condensed phase.