The Redfield equation describes the dynamics of a quantum system weakly coupled to one or more reservoirs and is widely used in theory of open quantum system. However, the assumption of weak system-reservoir coupling is often not fully adequate and higher-order corrections to the Redfield equation with respect to the system-bath coupling is required. Here we propose a general method of derivation of higher-order corrections to the Redfield quantum master equation based on the hierarchical equations of motion (HEOM). Also we derive conditions of validity of the Redfield equation as well as the additional secular approximation for it.
The Redfield equation describes the dynamics of a quantum system weakly coupled to one or more reservoirs and is widely used in theory of open quantum system. However, the assumption of weak system-reservoir coupling is often not fully adequate and higher-order corrections to the Redfield equation with respect to the system-bath coupling is required. Here we propose a general method of derivation of higher-order corrections to the Redfield quantum master equation based on the hierarchical equations of motion (HEOM). Also we derive conditions of validity of the Redfield equation as well as the additional secular approximation for it.
Abstract Following the ideas N. N. Bogoliubov used to derive the classical and quantum nonlinear kinetic equations, we derive the Redfield quantum linear master equation, which is widely used in the theory of open quantum systems. This method allows one to calculate the joint state of a system and a reservoir at every instant of time, as well as to obtain quantum master equations as autonomous differential equations in an arbitrary order of perturbation theory. We prove that under certain conditions the expressions for the corrections of all orders are well defined. We also discuss the question of whether the quantum dynamics described by the derived quantum master equations is Markovian.
An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.
The reduced dynamics of an open quantum system obtained from an underlying microscopic Hamiltonian can in general only approximately be described by a time-local master equation. The quality of that approximation depends primarily on the coupling strength and the structure of the environment. Various such master equations have been proposed with different aims. Choosing the most suitable one for a specific system is not straightforward. By focusing on the accuracy of the reduced dynamics we provide a thorough assessment for a selection of methods (Redfield equation, quantum optical master equation, coarse-grained master equation, a related dynamical map approach, and a partial-secular approximation). We consider secondary, here, whether or not an approach guarantees positivity. We use two qubits coupled to a Lorentzian environment in a spin-boson-like fashion modeling a generic situation with various system and bath time scales. We see that, independent of the initial state, the simple Redfield equation with time-dependent coefficients is significantly more accurate than all other methods under consideration. We emphasize that positivity violation in the Redfield equation formalism becomes relevant only in a regime where any of the perturbative master equations considered here are rendered invalid anyway. This implies that the loss of positivity should in fact be welcomed as an important feature: it indicates the breakdown of the weak-coupling assumption. In addition we present the various approaches in a self-contained way and use the behavior of their errors to provide further insight into the range of validity of each method.
Abstract Quantum master equations are widely used to describe the dynamics of open quantum systems. All these different master equations rely on specific approximations that may or may not be justified. Starting from a microscopic model, applying the standard Born and Markov approximations results in the Redfield equation that does not guarantee to preserve positivity. The latter is typically achieved by additionally applying the secular approximation resulting in a quantum master equation in Lindblad form. There are other ways to obtain an equation in Lindblad form, one of which is the recently proposed universal Lindblad equation (ULE). It has been shown that it is in the same equivalence class of approximations as the Redfield master equation although avoiding the heuristic secular approximation (Nathan and Rudner 2020 Phys. Rev. B 102 115109). In this work, we prove that the solutions of the master equation obtained by applying the secular approximation are also obtained by an approximation of the same order as the one performed to obtain the Redfield equation. We hereby provide a mathematical justification for the secular approximation. We show that the result of applying the secular approximation is obtained naturally by applying a self-consistency argument. This shows that the resulting master equation is also in the same equivalence class of approximations as the Redfield master equation and the ULE. We furthermore give numerical considerations about the accuracy of the secular master equation compared to the ULE.
We put forth a new class of quantum master equations that correctly reproduce the asymptotic state of an open quantum system beyond the infinitesimally weak system-bath coupling limit. Our method is based on incorporating the knowledge of the reduced steady state into its dynamics. The correction not only steers the reduced system toward a correct steady state but also improves the accuracy of the dynamics, thereby refining the archetypal Born-Markov weak-coupling second-order master equations. In case of equilibrium, we use the exact mean-force Gibbs state to correct the Redfield quantum master equation. By benchmarking it with the exact solution of the damped harmonic oscillator, we show that our method also helps correct the long-standing issue of positivity violation, albeit without complete positivity. Our method of a canonically consistent quantum master equation opens a new perspective in the theory of open quantum systems leading to a reduced density matrix accurate beyond the commonly used Redfield and Lindblad equations, while retaining the same conceptual and numerical complexity.
It is very common in the literature to write a Markovian quantum master equation in Lindblad form to describe a system with multiple degrees of freedom and weakly connected to multiple thermal baths which can, in general, be at different temperatures and chemical potentials. However, the microscopically derived quantum master equation up to leading order in a system-bath coupling is of the so-called Redfield form, which is known to not preserve complete positivity in most cases. Additional approximations to the Redfield equation are required to obtain a Lindblad form. We lay down some fundamental requirements for any further approximations to the Redfield equation, which, if violated, leads to physical inconsistencies such as inaccuracies in the leading order populations and coherences in the energy eigenbasis, violation of thermalization, and violation of local conservation laws at the nonequilibrium steady state. We argue that one or more of these conditions will generically be violated in all the weak system-bath-coupling Lindblad descriptions existing in the literature to our knowledge. As an example, we study the recently derived universal Lindblad equation and use these conditions to show the violation of local conservation laws due to inaccurate coherences but accurate populations in the energy eigenbasis. Finally, we exemplify our analytical results numerically in an interacting open quantum spin system.