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.
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.
Lindbladian quantum master equations (LEs) are the most popular descriptions for quantum systems weakly coupled to baths. But, recent works have established that in many situations such Markovian descriptions are fundamentally limited: they cannot simultaneously capture populations and coherences even to the leading-order in system-bath couplings. This can cause violation of fundamental properties like thermalization and continuity equations associated with local conservation laws, even when such properties are expected in the actual setting. This begs the question: given a physical situation, how do we know if there exists an LE that describes it to a desired accuracy? Here we show that, for both equilibrium and non-equilibrium steady states (NESS), this question can be succinctly formulated as a semidefinite program (SDP), a convex optimization technique. If a solution to the SDP can be found to a desired accuracy, then an LE description is possible for the chosen setting. If not, no LE description is fundamentally attainable, showing that a consistent Markovian treatment is impossible even at weak system-bath coupling for that particular setting. Considering few qubit isotropic XXZ-type models coupled to multiple baths, we find that in most parameter regimes, LE description giving accurate populations and coherences to leading-order is unattainable, leading to rigorous no-go results. However, in some cases, LE description having correct populations but inaccurate coherences, and satisfying local conservation laws, is possible over some of the parameter regimes. Our work highlights the power of semidefinite programming in the analysis of physically consistent LEs, thereby, in understanding the limits of Markovian descriptions at weak system-bath couplings.
We investigate the possibility of a Markovian quantum master equation (QME) that consistently describes a finite-dimensional system, a part of which is weakly coupled to a thermal bath. To preserve complete positivity and trace, such a QME must be of Lindblad form. For physical consistency, it should additionally preserve local conservation laws and be able to show thermalization. We search for Lindblad equations satisfying these additional criteria. First, we show that the microscopically derived Bloch-Redfield equation violates complete positivity unless in extremely special cases. We then prove that imposing complete positivity and demanding preservation of local conservation laws enforces the Lindblad operators and the lamb-shift Hamiltonian to be local, i.e., to be supported only on the part of the system directly coupled to the bath. We then cast the problem of finding local Lindblad QMEs which can show thermalization into a semidefinite program. We call this the thermalization optimization problem (TOP). For given system parameters and temperature, the solution of the TOP conclusively shows whether the desired type of QME is possible up to a given precision. Whenever possible, it also outputs a form for such a QME. For an XXZ chain of a few qubits, fixing a reasonably high precision, we find that such a QME is impossible over a considerably wide parameter regime when only the first qubit is coupled to the bath. Remarkably, we find that when the first two qubits are attached to the bath, such a QME becomes possible over much of the same parameter regime, including a wide range of temperatures.
The Lindblad form guarantees complete positivity of a Markovian quantum master equation (QME). However, its microscopic derivation for a quantum system weakly interacting with a thermal bath requires several approximations, which may result in inaccuracies in the QME. Recently, various Lindbladian QMEs were derived without resorting to the secular approximation from the Redfield equation which does not guarantee the complete positivity. Here we explicitly calculate, in a perturbative manner, the equilibrium steady states of these Lindbladian QMEs. We compare the results with the steady state of the Redfield equation obtained from an analytic continuation method, which coincides with the so-called mean force Gibbs (MFG) state. The MFG state is obtained by integrating out the bath degrees of freedom for the Gibbs state of the total Hamiltonian. We explicitly show that the steady states of the Lindbladian QMEs are different from the MFG state. Our results indicate that manipulations of the Redfield equation needed to enforce complete positivity of a QME drives its steady state away from the MFG state. We also find that, in the high-temperature regime, both the steady states of the Lindbladian QMEs and MFG state reduce to the same Gibbs state of a system Hamiltonian under certain conditions.
We develop a Markovian master equation in the Lindblad form that enables the efficient study of a wide range of open quantum many-body systems that would be inaccessible with existing methods. The validity of the master equation is based entirely on properties of the bath and the system-bath coupling, without any requirements on the level structure within the system itself. The master equation is derived using a Markov approximation that is distinct from that used in earlier approaches. We provide a rigorous bound for the error induced by this Markov approximation; the error is controlled by a dimensionless combination of intrinsic correlation and relaxation timescales of the bath. Our master equation is accurate on the same level of approximation as the Bloch-Redfield equation. In contrast to the Bloch-Redfield approach, our approach ensures preservation of the positivity of the density matrix. As a result, our method is robust, and can be solved efficiently using stochastic evolution of pure states (rather than density matrices). We discuss how our method can be applied to static or driven quantum many-body systems, and illustrate its power through numerical simulation of a spin chain that would be challenging to treat by existing methods.
We present the Born-Markov approximated Redfield quantum master equation (RQME) description for an open system of noninteracting particles (bosons or fermions) on an arbitrary lattice of N sites in any dimension and weakly connected to multiple reservoirs at different temperatures and chemical potentials. The RQME can be reduced to the Lindblad equation, of various forms, by making further approximations. By studying the N=2 case, we show that RQME gives results which agree with exact analytical results for steady-state properties and with exact numerics for time-dependent properties over a wide range of parameters. In comparison, the Lindblad equations have a limited domain of validity in nonequilibrium. We conclude that it is indeed justified to use microscopically derived full RQME to go beyond the limitations of Lindblad equations in out-of-equilibrium systems. We also derive closed-form analytical results for out-of-equilibrium time dynamics of two-point correlation functions. These results explicitly show the approach to steady state and thermalization. These results are experimentally relevant for cold atoms, cavity QED, and far-from-equilibrium quantum dot experiments.
Away from equilibrium, the properties of open quantum systems depend on the details of their environment. A microscopic derivation of a master equation (ME) is therefore crucial. Of particular interest are Lindblad-type equations, not only because they provide the most general class of Markovian MEs, but also since they are the starting point for efficient quantum trajectory simulations. Lindblad-type MEs are commonly derived from the Born-Markov-Redfield equation via a rotating-wave approximation (RWA). However the RWA is valid only for ultraweak system-bath coupling and often fails to accurately describe nonequilibrium processes. Here we derive an alternative Lindbladian approximation to the Redfield equation, which does not rely on ultraweak system-bath coupling. Applying it to an extended Hubbard model coupled to Ohmic baths, we show that, especially away from equilibrium, it provides a good approximation in large parameter regimes where the RWA fails.