The Lindblad equation is an evolution equation for the density matrix in quantum theory. It is the general linear, Markovian, form which ensures that the density matrix is Hermitian, trace 1, positive and completely positive. Some elementary examples of the Lindblad equation are given. The derivation of the Lindblad equation presented here is ‘simple’ in that all it uses is the expression of a Hermitian matrix in terms of its orthonormal eigenvectors and real eigenvalues. Thus, it is appropriate for students who have learned the algebra of quantum theory. Where helpful, arguments are first given in a two-dimensional Hilbert space.
The Lindblad equation is an evolution equation for the density matrix in quantum theory. It is the general linear, Markovian, form which ensures that the density matrix is Hermitian, trace 1, positive and completely positive. Some elementary examples of the Lindblad equation are given. The derivation of the Lindblad equation presented here is ‘simple’ in that all it uses is the expression of a Hermitian matrix in terms of its orthonormal eigenvectors and real eigenvalues. Thus, it is appropriate for students who have learned the algebra of quantum theory. Where helpful, arguments are first given in a two-dimensional Hilbert space.
We give a simple argument that, for a large class of jump operators, the Lindblad evolution can be written as a gradient flow in the space of density operators acting on a Hilbert space of dimension D. We give explicit expressions for the (matrix-valued) eigenvectors and eigenvalues of the Lindblad evolution using this formalism. We argue that in many cases the interpretation of the evolution is simplified by passing from the complex D2-dimensional space of density operators to the real ${D}^{2}\ensuremath{-}1$-dimensional space of Bloch vectors. When jump operators are non-Hermitian the evolution is not in general gradient flow, but we show that it nevertheless resembles gradient flow in two particular ways. Importantly, the steady states of Lindbladian evolution are still determined by the potential in all cases.
The exact solution of the Lindblad equation with a quadratic Hamiltonian and linear coupling operators was derived within the chord representation, that is, for the Fourier transform of the Wigner function, also known as the characteristic function. It is here generalized for several degrees of freedom, so as to provide an explicit expression for the reduced density operator of any subsystem, as well as moments expressed as derivatives of this evolving chord function. The Wigner function is then the convolution of its straightforward classical evolution with a widening multidimensional Gaussian window, eventually ensuring its positivity. Futher on, positivity also holds for the Glauber-Sundarshan P function, which guarantees separability of the components. In the context of several degrees of freedom, a full dissipation matrix is defined, whose trace is equal to twice the previously derived dissipation coefficient. This governs the rate at which the phase space volume of the argument of the Wigner function contracts, while that of the chord function expands. Examples of Markovian evolution of a triatomic molecule and of an array of harmonic oscillators are discussed.
The exact solution of the Lindblad equation with a quadratic Hamiltonian and linear coupling operators was derived within the chord representation, that is, for the Fourier transform of the Wigner function, also known as the characteristic function. It is here generalized for several degrees of freedom, so as to provide an explicit expression for the reduced density operator of any subsystem, as well as moments expressed as derivatives of this evolving chord function. The Wigner function is then the convolution of its straightforward classical evolution with a widening multidimensional Gaussian window, eventually ensuring its positivity. Futher on, positivity also holds for the Glauber-Sundarshan P function, which guarantees separability of the components. In the context of several degrees of freedom, a full dissipation matrix is defined, whose trace is equal to twice the previously derived dissipation coefficient. This governs the rate at which the phase space volume of the argument of the Wigner function contracts, while that of the chord function expands. Examples of Markovian evolution of a triatomic molecule and of an array of harmonic oscillators are discussed.
The structure of the initial system-environment state is fundamental to determining the nature and characteristics of the evolution of such an open quantum system. The usual assumption is to consider that the initial system-environment state is separable. Here, we go beyond this simple case and derive the evolution equations, up to second order in a weak-coupling expansion, that describe the evolution of the reduced density matrix of the system for any arbitrary system-environment initial state. The structure of these equations allows us to determine the initial conditions for which the evolution of the reduced density matrix can be written in terms of a set of Lindblad-like equations, once considering the Markov and secular approximations. Moreover, we show that for initial states belonging to a subset of separable states such Lindblad-like equations become actual Lindblad equations and that the evolution of ${\ensuremath{\rho}}_{s}$ thereby obtained is also trace and positive preserving.
In the presence of Lindblad decoherence, i.e., dissipative effects in an open quantum system due to interaction with an environment, we examine the transition probabilities between the eigenstates in the two-level quantum system described by non-Hermitian Hamiltonians with the Lindblad equation, for which the parity-time-reversal (PT) symmetry is conserved. First, the density matrix formalism for PT-symmetric non-Hermitian Hamiltonian systems is developed. It is shown that the Lindblad operators Lj are pseudo-Hermitian, namely, $\ensuremath{\eta}{L}_{j}^{}{\ensuremath{\eta}}^{\ensuremath{-}1}={L}_{j}^{\ifmmode\dagger\else\textdagger\fi{}}$ with $\ensuremath{\eta}$ being a linear and positive-definite metric, and respect the PT symmetry as well. We demonstrate that the generalized density matrix ${\ensuremath{\rho}}_{\mathrm{G}}^{}(t)\ensuremath{\equiv}\ensuremath{\rho}(t)\ensuremath{\eta}$, instead of the normalized density matrix ${\ensuremath{\rho}}_{\mathrm{N}}^{}(t)\ensuremath{\equiv}\ensuremath{\rho}(t)/\mathrm{tr}\left[\ensuremath{\rho}(t)\right]$, should be implemented for the calculation of the transition probabilities in accordance with the linearity requirement. Second, the density matrix formalism is used to derive the transition probabilities in general cases of PT-symmetric non-Hermitian Hamiltonians. In some concrete examples, we calculate compact analytical formulas for the transition probabilities and explore their main features with numerical illustrations. We also make a comparison between our present results and our previous ones using state vectors in the absence of Lindblad decoherence.
A semiclassical approximation for an evolving density operator, driven by a 'closed' Hamiltonian operator and 'open' Markovian Lindblad operators, is obtained. The theory is based on the chord function, i.e. the Fourier transform of the Wigner function. It reduces to an exact solution of the Lindblad master equation if the Hamiltonian operator is a quadratic function and the Lindblad operators are linear functions of positions and momenta. Initially, the semiclassical formulae for the case of Hermitian Lindblad operators are reinterpreted in terms of a (real) double phase space, generated by an appropriate classical double Hamiltonian. An extra 'open' term is added to the double Hamiltonian by the non-Hermitian part of the Lindblad operators in the general case of dissipative Markovian evolution. The particular case of generic Hamiltonian operators, but linear dissipative Lindblad operators, is studied in more detail. A Liouville-type equivariance still holds for the corresponding classical evolution in double phase space, but the centre subspace, which supports the Wigner function, is compressed, along with expansion of its conjugate subspace, which supports the chord function. Decoherence narrows the relevant region of double phase space to the neighbourhood of a caustic for both the Wigner function and the chord function. This difficulty is avoided by a propagator in a mixed representation, so that a further 'small-chord' approximation leads to a simple generalization of the quadratic theory for evolving Wigner functions.