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.
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.
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.
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.
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.
Starting from a microscopic system-baths description, we derive the general conditions for a time-local quantum master equation (QME) to satisfy the first and second laws of thermodynamics at the fluctuating level. Using counting statistics, we show that the fluctuating second law can be rephrased as a generalized quantum detailed balance condition (GQDB), i.e., a symmetry of the time-local generators which ensures the validity of the fluctuation theorem. When requiring in addition a strict system-bath energy conservation, the GQDB reduces to the usual notion of detailed balance which characterizes QMEs with Gibbsian steady states. However, if energy conservation is only required on average, QMEs with non-Gibbsian steady states can still maintain a certain level of thermodynamic consistency. Applying our theory to commonly used QMEs, we show that the Redfield equation breaks the GQDB, and that some recently derived approximation schemes based on the Redfield equation (which hold beyond the secular approximation and allow one to derive a QME of Lindblad form) satisfy the GQDB and the average first law. We find that performing the secular approximation is the only way to ensure the first and second laws at the fluctuating level.
The perturbative master equation (Bloch–Redfield) is used extensively to study dissipative quantum mechanics—particularly for qubits—despite the 25-year-old criticism that it violates positivity (generating negative probabilities). We take an arbitrary system coupled to an environment containing many degrees-of-freedom and cast its perturbative master equation (derived from a perturbative treatment of Nakajima–Zwanzig or Schoeller–Schön equations) in the form of a Lindblad master equation. We find that the equation's parameters are time dependent . This time dependence is rarely accounted for and invalidates Lindblad's dynamical semigroup analysis. We analyse one such Bloch–Redfield master equation (for a two-level system coupled to an environment with a short but non-vanishing memory time), which apparently violates positivity. We analytically show that, once the time dependence of the parameters is accounted for, positivity is preserved.
While quantum master equations (QMEs) are the primary workhorse in quantum information science, quantum optics, spectroscopy, and quantum thermodynamics, bringing an arbitrary N-level QME into Lindbladian form and verifying complete positivity of the associated quantum dynamical map remain open challenges for N≥3. We explore and implement two independent methods to accomplish these tasks, which enable one to directly compute the Kossakowski matrix of an arbitrary Markovian QME from its Liouvillian. In the first method, due to Hall, Cresser, Li, and Andersson, the Kossakowski matrix elements are obtained by evaluating the action of the Liouvillian on the orthonormal SU(N) basis matrices and then computing a sum of matrix-product traces. The second method, developed in this work, is based on the real N-level coherence vector and relies on the Moore-Penrose pseudo-inverse of a rectangular matrix composed of the structure constants of SU(N). We show that both methods give identical results, and apply them to establish the complete positivity of the partial secular Bloch-Redfield QME for the Λ and V-systems driven by incoherent light. We find that the eigenvalues of the Kossakowski matrix of these seemingly different three-level systems are identical, implying close similarities of their dissipative dynamics. By facilitating the expression of multilevel Markovian QMEs in Lindblad form, our results enable testing the QMEs for complete positivity without solving them, as well as restoring complete positivity by keeping only non-negative eigenvalues of the Kossakowski matrix.