The small matrix path integral (SMatPI) methodology allows fully quantum mechanical calculations on system–bath Hamiltonians under a variety of conditions with minimal array storage and is thus applicable to multistate systems. This work focuses on the SMatPI iteration algorithm, which can be computationally prohibitive in systems with many states. By transforming the SMatPI propagation matrices to the desired initial state, the cost of the iterative multiplication is reduced by a factor of n2, where n is the number of system states. This leads to dramatic savings when n≫10 and allows the simulation of large systems over very long propagation times. Applications are presented on models with up to 50 system states, including exciton transfer dynamics in one-dimensional aggregates and a square lattice.
Nancy Makri (Tue,) studied this question.