The Nakajima–Zwanzig generalized quantum master equation provides a general, and formally exact, prescription for simulating the reduced dynamics of a quantum system coupled to a quantum bath. In this equation, the memory kernel accounts for the influence of the bath on the system’s dynamics. The standard approach is based on using a perturbative treatment of the system–bath coupling for calculating this kernel, and is therefore restricted to systems weakly coupled to the bath. In this paper, we propose a new approach for calculating the memory kernel for an arbitrary system–bath coupling. The memory kernel is obtained by solving a set of two coupled integral equations that relate it to a new type of two-time system-dependent bath correlation functions. The feasibility of the method is demonstrated in the case of an asymetrical two-level system linearly coupled to a harmonic bath.
The Nakajima–Zwanzig generalized quantum master equation provides a general, and formally exact, prescription for simulating the reduced dynamics of a quantum system coupled to a quantum bath. In this equation, the memory kernel accounts for the influence of the bath on the system’s dynamics. The standard approach is based on using a perturbative treatment of the system–bath coupling for calculating this kernel, and is therefore restricted to systems weakly coupled to the bath. In this paper, we propose a new approach for calculating the memory kernel for an arbitrary system–bath coupling. The memory kernel is obtained by solving a set of two coupled integral equations that relate it to a new type of two-time system-dependent bath correlation functions. The feasibility of the method is demonstrated in the case of an asymetrical two-level system linearly coupled to a harmonic bath.
The Nakajima-Zwanzig generalized quantum master equation (GQME) provides a general, and formally exact, prescription for simulating the reduced dynamics of a quantum system coupled to a, possibly anharmonic, quantum bath. In this equation, a memory kernel superoperator accounts for the influence of the bath on the dynamics of the system. In a previous paper [Q. Shi and E. Geva, J. Chem. Phys. 119, 12045 (2003)] we proposed a new approach to calculating the memory kernel, in the case of arbitrary system-bath coupling. Within this approach, the memory kernel is obtained by solving a set of two integral equations, which requires a new type of two-time system-dependent bath correlation functions as input. In the present paper, we consider the application of the linearized semiclassical (LSC) approximation for calculating those correlation functions, and subsequently the memory kernel. The new approach is tested on a benchmark spin-boson model. Application of the LSC approximation for calculating the relatively short-lived memory kernel, followed by a numerically exact solution of the GQME, is found to provide an accurate description of the relaxation dynamics. The success of the proposed LSC-GQME methodology is contrasted with the failure of both the direct application of the LSC approximation and the weak coupling treatment to provide an accurate description of the dynamics, for the same model, except at very short times. The feasibility of the new methodology to anharmonic systems is also demonstrated in the case of a two level system coupled to a chain of Lennard-Jones atoms.
The Nakajima-Zwanzig generalized quantum master equation provides a general, and formally exact, prescription for simulating the reduced dynamics of a quantum system coupled to a quantum bath. In this equation, the memory kernel accounts for the influence of the bath on the system's dynamics, and the inhomogeneous term accounts for initial system-bath correlations. In this paper, we propose a new approach for calculating the memory kernel and inhomogeneous term for arbitrary initial state and system-bath coupling. The memory kernel and inhomogeneous term are obtained by numerically solving a single inhomogeneous Volterra equation of the second kind for each. The new approach can accommodate a very wide range of projection operators, and requires projection-free two-time correlation functions as input. An application to the case of a two-state system with diagonal coupling to an arbitrary bath is described in detail. Finally, the utility and self-consistency of the formalism are demonstrated by an explicit calculation on a spin-boson model.
The Nakajima-Zwanzig generalized master equation provides a formally exact framework to simulate quantum dynamics in condensed phases. Yet, the exact memory kernel is hard to obtain and calculations based on perturbative expansions are often employed. By using the spin-boson model as an example, we assess the convergence of high order memory kernels in the Nakajima-Zwanzig generalized master equation. The exact memory kernels are calculated by combining the hierarchical equation of motion approach and the Dyson expansion of the exact memory kernel. High order expansions of the memory kernels are obtained by extending our previous work to calculate perturbative expansions of open system quantum dynamics [M. Xu et al., J. Chem. Phys. 146, 064102 (2017)]. It is found that the high order expansions do not necessarily converge in certain parameter regimes where the exact kernel show a long memory time, especially in cases of slow bath, weak system-bath coupling, and low temperature. Effectiveness of the Padé and Landau-Zener resummation approaches is tested, and the convergence of higher order rate constants beyond Fermi's golden rule is investigated.
We present a modified approach for simulating electronically nonadiabatic dynamics based on the Nakajima-Zwanzig generalized quantum master equation (GQME). The modified approach utilizes the fact that the Nakajima-Zwanzig formalism does not require casting the overall Hamiltonian in system-bath form, which is arguably neither natural nor convenient in the case of the Hamiltonian that governs nonadiabatic dynamics. Within the modified approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density operator is fully captured by a memory kernel super-operator. A methodology for calculating the memory kernel from projection-free inputs is developed. Simulating the electronic dynamics via the modified approach, with a memory kernel obtained using exact or approximate methods, can be more cost effective and/or lead to more accurate results than direct application of those methods. The modified approach is compared to previously proposed GQME-based approaches, and its robustness and accuracy are demonstrated on a benchmark spin-boson model with a memory kernel which is calculated within the Ehrenfest method.
The generalized quantum master equation provides an exact equation of motion for the reduced density operator of a quantum system coupled to one or more baths, with all bath-induced effects encoded in a memory kernel superoperator. In this paper, we investigate whether the memory kernel for a system coupled simultaneously to multiple independent baths can be written as a sum of memory kernels associated with the corresponding single-bath problems. We show that such additivity holds in the weak system-bath coupling limit, where the second-order memory kernel separates into independent bath contributions. Beyond this limit, however, the baths become dynamically correlated through their mutual coupling to the system, and the exact multi-bath memory kernel is no longer equal to the sum of exact single-bath memory kernels. We demonstrate this non-additivity for a two-level system linearly coupled to hot and cold harmonic baths, using the hierarchical equations of motion method to construct the memory kernels. The results show that non-additivity increases with system-bath coupling strength and is process-dependent: coherence damping remains comparatively well described by an additive kernel, whereas population-coherence and coherence-coherence transfer exhibit pronounced non-additive behavior. These findings clarify the limitations of additive multi-bath descriptions and highlight the ambiguity of decomposing heat flow into independent bath-resolved contributions beyond weak system-bath coupling.
Abstract The dynamics of quantum systems coupled to baths are typically studied using the Nakajima-Zwanzig memory kernel ( K K ) or the influence functions (I), particularly when memory effects are present. Despite their significance, formal connections between the two have not been explicitly known. We establish their connections by examining the system propagator for a N-level system linearly coupled to Gaussian baths with various types of system-bath coupling. For a certain class of problems, we devised a non-perturbative, diagrammatic approach to construct K K from I for (driven) systems interacting with Gaussian baths, bypassing conventional projection-free dynamics inputs. Our work provides a way to interpret approximate path integral methods in terms of approximate memory kernels. Moreover, it offers a Hamiltonian learning procedure to extract the bath spectral density from reduced system trajectories, opening new avenues in quantum sensing and engineering. The insights we provide advance our understanding of non-Markovian dynamics and will serve as a stepping stone for future theoretical and experimental developments in this area.