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.
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.
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.
Fourth-order quantum master equations (FQMEs) are derived in both time nonlocal and local forms for a general system Hamiltonian, with new detailed expressions for the fourth-order kernel, where the bath correlation functions are explicitly decoupled from the system superoperators. Further simplifications can be made for the model of linearly coupled harmonic oscillator bath. Consideration of the high temperature Ohmic bath limit leads to a general Markovian FQME with compact forms of time independent superoperators. Two examples of this equation are then considered. For the system of a quantum particle in a continuous potential field, the equation reduces to a known form of the quantum Fokker–Planck equation, except for a fourth-order potential renormalization term that can be neglected only in the weak system-bath interaction regime. For a two-level system with off-diagonal coupling to the bath, fourth-order corrections do not alter the relaxation characteristics of the second-order equation and introduce additional coherence terms in the equations for the off-diagonal elements.
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.
Generalized master equations provide a concise formalism for studying reduced population dynamics. Usually, these master equations require a perturbative expansion of the memory kernels governing the dynamics; in order to prevent divergences, these expansions must be resummed. Resummation techniques of perturbation series are ubiquitous in physics, but they have not been readily studied for the time-dependent memory kernels used in generalized master equations. In this paper, we present a comparison of different resummation techniques for such memory kernels up to fourth order. We study specifically the spin-boson Hamiltonian as a model system bath Hamiltonian, treating the diabatic coupling between the two states as a perturbation. A novel derivation of the fourth-order memory kernel for the spin-boson problem is presented; then, the second- and fourth-order kernels are evaluated numerically for a variety of spin-boson parameter regimes. We find that resumming the kernels through fourth order using a Padé approximant results in divergent populations in the strong electronic coupling regime due to a singularity introduced by the nature of the resummation, and thus recommend a non-divergent exponential resummation (the "Landau-Zener resummation" of previous work). The inclusion of fourth-order effects in a Landau-Zener-resummed kernel is shown to improve both the dephasing rate and the obedience of detailed balance over simpler prescriptions like the non-interacting blip approximation, showing a relatively quick convergence on the exact answer. The results suggest that including higher-order contributions to the memory kernel of a generalized master equation and performing an appropriate resummation can provide a numerically-exact solution to system-bath dynamics for a general spectral density, opening the way to a new class of methods for treating system-bath dynamics.