Accurate Two-Time Correlation from a Pseudomode-Based Variational Approach with the Multiple Davydov D1 Density...

The pseudomode framework provides an efficient description of open quantum systems in non-Markovian and nonperturbative coupling regimes by utilizing simple Markovian Lindblad master equations that incorporate a set of damped bosonic pseudomodes. Developing an efficient numerical scheme to solve these equations accurately is highly desirable. In this work, we propose a modified multiple Davydov D1 density operator ansatz and combine it with the Dirac-Frenkel time-dependent variational principle to simulate the dynamics governed by these master equations. We apply this variational method to the spin-boson model, both with and without a periodic external driving field, and compute both single- and two-time expectation values in non-Markovian, strong-coupling regimes. The accuracy of our variational approach is validated by comparing its predictions with results obtained from the small matrix path integral method. We demonstrate that the variational and path-integral approaches are in excellent agreement for both single- and two-time expectations across nondriven and driven cases, spanning both low- and high-temperature regimes. Our results establish the reliability and broad applicability of the pseudomode-based variational approach for calculating two-time correlation functions in open quantum systems.