Stochastic Schrödinger equation (SSE) methods provide a powerful approach for simulating open quantum system dynamics by unraveling the influence functional into an ensemble average of pure-state trajectories. However, long-time simulations using SSE methods are hindered by the inherent growth of ensemble variance. An alternative strategy is to extract essential memory kernels or transfer tensors from short-time dynamics and extrapolate to arbitrarily long times, as in the transfer tensor method (TTM). Conventional TTM implementations, however, require dynamical maps constructed from an informationally complete set of N2 initial states─where N is the system dimension─limiting their applicability to small-scale problems. Here, we show that the SSE naturally provides an overcomplete set of Kraus operators that captures the essential information needed to construct dynamical maps for TTM. This connection enables the extension of TTM to high-dimensional systems using results obtained from a single converged SSE simulation with a sufficient number of stochastic realizations. We validate our approach with two representative SSE methods applied to a symmetric spin-boson model and a 24-site Fenna–Matthews–Olson (FMO) complex. The results demonstrate that SSE-TTM offers an accurate and computationally efficient strategy for simulating long-time quantum dynamics in relatively large-scale systems.
Gao et al. (Mon,) studied this question.