Skip to content

The Lindblad theorem characterizes the generators of Markovian completely positive trace-preserving quantum dynamics. In finite dimensions, it says that a norm-continuous quantum dynamical semigroup is CPTP if and only if its generator has Gorini–Kossakowski–Sudarshan–Lindblad form.

The theorem is structural. It does not derive the generator of a particular atom, circuit, molecule, detector, or reservoir. It tells you what form the generator must have once the semigroup, complete-positivity, and trace-preservation assumptions are imposed.

Helpful background. Completely positive maps and the Choi matrix explain positivity with a reference. The Lindblad–GKSL equation gives elementary examples of the generator.

Let H\mathcal H be finite-dimensional and let T(H)\mathcal T(\mathcal H) be the vector space of trace-class operators, which is just the full matrix algebra in finite dimension.

A quantum dynamical semigroup is a family of linear maps

{Φt:t≥0}\{\Phi_t:t\ge0\}

such that

Φ0=id,Φt+s=ΦtΦs,t,s≥0.\Phi_0=\mathrm{id}, \qquad \Phi_{t+s}=\Phi_t\Phi_s, \qquad t,s\ge0.

The family has time-homogeneous reduced composition: the same map Φt\Phi_t advances any modeled state from ss to s+ts+t. This concerns one-time reduced maps; it does not determine arbitrary multi-time experiments with interventions.

The physical requirements are:

  • Φt\Phi_t is trace preserving for every t≥0t\ge0;
  • Φt\Phi_t is completely positive for every t≥0t\ge0;
  • Φt\Phi_t is continuous at t=0t=0 in the operator norm, so a bounded generator exists.

The generator is

L=lim⁡t↓0Φt−idt,\mathcal L = \lim_{t\downarrow0} \frac{\Phi_t-\mathrm{id}}{t},

and the semigroup is

Φt=etL.\Phi_t=e^{t\mathcal L}.

In finite dimensions, a linear map L\mathcal L generates a norm-continuous CPTP semigroup Φt=etL\Phi_t=e^{t\mathcal L} if and only if there exist a Hermitian operator H=H†H=H^\dagger and operators LjL_j such that

L(ρ)=−iℏ[H,ρ]+∑j(LjρLj†−12{Lj†Lj,ρ}).\mathcal L(\rho) = - \frac{i}{\hbar}[H,\rho] + \sum_j \left( L_j\rho L_j^\dagger - \frac12 \{L_j^\dagger L_j,\rho\} \right).

Equivalently, one may write Lj=γjVjL_j=\sqrt{\gamma_j}V_j with γj≥0\gamma_j\ge0:

L(ρ)=−iℏ[H,ρ]+∑jγj(VjρVj†−12{Vj†Vj,ρ}).\mathcal L(\rho) = - \frac{i}{\hbar}[H,\rho] + \sum_j \gamma_j \left( V_j\rho V_j^\dagger - \frac12 \{V_j^\dagger V_j,\rho\} \right).

The first direction says: if the generator has this form, then etLe^{t\mathcal L} is a CPTP semigroup.

The converse says: if etLe^{t\mathcal L} is a CPTP semigroup under the assumptions above, then L\mathcal L can be written in this form.

The adjoint generator acts on observables. With the trace pairing

Tr⁡[A L(ρ)]=Tr⁡[L†(A)ρ],\operatorname{Tr}[A\,\mathcal L(\rho)] = \operatorname{Tr}[\mathcal L^\dagger(A)\rho],

the Heisenberg-picture form is

L†(A)=iℏ[H,A]+∑j(Lj†ALj−12{Lj†Lj,A}).\mathcal L^\dagger(A) = \frac{i}{\hbar}[H,A] + \sum_j \left( L_j^\dagger A L_j - \frac12 \{L_j^\dagger L_j,A\} \right).

Trace preservation of L\mathcal L is equivalent to

L†(I)=0.\mathcal L^\dagger(I)=0.

Thus the observable identity stays fixed:

Φt†(I)=I.\Phi_t^\dagger(I)=I.

The semigroup condition

Φt+s=ΦtΦs\Phi_{t+s}=\Phi_t\Phi_s

imposes time-homogeneous composition. A generic time-dependent drive or finite-bath memory model will not satisfy it. An alternative memory-kernel representation of a semigroup can still be written, so the presence of an integral in some representation is not by itself decisive.

Complete positivity means

(Φt⊗idR)(X)≥0wheneverX≥0(\Phi_t\otimes\mathrm{id}_R)(X)\ge0 \quad \text{whenever} \quad X\ge0

for every finite reference system RR. It is stronger than positivity on the system alone and is essential for consistency when the system may be entangled with a spectator. See Completely Positive Maps.

Trace preservation keeps total probability normalized:

Tr⁡Φt(ρ)=Tr⁡ρ.\operatorname{Tr}\Phi_t(\rho)=\operatorname{Tr}\rho.

At the generator level this becomes

Tr⁡L(ρ)=0.\operatorname{Tr}\mathcal L(\rho)=0.

The anticommutator term in the Lindblad form is what makes each dissipator trace preserving.

Continuity at t=0t=0 ensures that an infinitesimal generator exists. In finite dimensions this is usually technically mild. In infinite dimensions it becomes one of the places where unbounded generators and domains enter.

The full proof belongs in a rigorous treatment of quantum dynamical semigroups. The finite-dimensional strategy can be summarized as follows.

First, choose an orthonormal operator basis consisting of the identity and traceless operators:

F0=Id,Tr⁡(Fa†Fb)=δab,a,b≥1.F_0=\frac{I}{\sqrt d}, \qquad \operatorname{Tr}(F_a^\dagger F_b)=\delta_{ab}, \qquad a,b\ge1.

Trace preservation and Hermiticity preservation constrain the generator. Terms involving the identity basis element combine into a commutator with a Hermitian operator. The remaining dissipative part can be written using a coefficient matrix cabc_{ab}:

Ldiss(ρ)=∑a,b≥1cab(FaρFb†−12{Fb†Fa,ρ}).\mathcal L_{\mathrm{diss}}(\rho) = \sum_{a,b\ge1} c_{ab} \left( F_a\rho F_b^\dagger - \frac12 \{F_b^\dagger F_a,\rho\} \right).

To see where the positivity constraint enters, use the output-first Choi convention

J(L)=∑mnL(∣m⟩⟨n∣)⊗∣m⟩⟨n∣.J(\mathcal L)=\sum_{mn}\mathcal L(|m\rangle\langle n|)\otimes|m\rangle\langle n|.

The identity map has J(id)=∣Ω⟩⟨Ω∣J(\mathrm{id})=|\Omega\rangle\langle\Omega|, with ∣Ω⟩=∑m∣m,m⟩|\Omega\rangle=\sum_m|m,m\rangle. On the orthogonal complement of ∣Ω⟩|\Omega\rangle, this matrix is zero. Since J(Φt)=J(id)+tJ(L)+O(t2)J(\Phi_t)=J(\mathrm{id})+tJ(\mathcal L)+O(t^2) must remain positive, QJ(L)Q≥0QJ(\mathcal L)Q\ge0, where Q=I−∣Ω⟩⟨Ω∣/dQ=I-|\Omega\rangle\langle\Omega|/d. The traceless basis vectors span this complement; their coefficient matrix is cc. Thus c≥0c\ge0. This is conditional positivity of the generator, not the false claim that L\mathcal L itself is a CP map.

Because c≥0c\ge0, it can be diagonalized:

cab=∑jλjujaujb∗,λj≥0.c_{ab} = \sum_j \lambda_j u_{ja}u_{jb}^*, \qquad \lambda_j\ge0.

Defining

Lj=λj∑aujaFaL_j = \sqrt{\lambda_j} \sum_a u_{ja}F_a

puts the dissipative part in Lindblad form. Here ujau_{ja} is component aa of eigenvector jj; the index order of the displayed coefficient formula fixes the conjugation convention.

For the converse, let R=∑jLj†LjR=\sum_jL_j^\dagger L_j and choose h>0h>0 with h∥R∥≤1h\|R\|\le1. The operators

K0=e−iHh/ℏI−hR,Kj=hLjK_0=e^{-iHh/\hbar}\sqrt{I-hR},\qquad K_j=\sqrt h L_j

satisfy ∑jKj†Kj+K0†K0=I\sum_jK_j^\dagger K_j+K_0^\dagger K_0=I exactly. Their channel is Ψh=id+hL+O(h2)\Psi_h=\mathrm{id}+h\mathcal L+O(h^2). In finite dimension, (Ψt/n)n→etL(\Psi_{t/n})^n\to e^{t\mathcal L}; compositions and limits preserve complete positivity and trace. This proves sufficiency without mistaking an Euler step for an exactly physical channel.

This sketch hides the technical work in showing exactly how complete positivity at infinitesimal times constrains cabc_{ab}, but it captures the theorem’s structure: positivity of a coefficient matrix becomes a sum of dissipators.

A complex coefficient matrix: worked example

Section titled “A complex coefficient matrix: worked example”

For a qubit choose F1=X/2F_1=X/\sqrt2, F2=Y/2F_2=Y/\sqrt2 and

c=κ(1i−i1),κ>0.c=\kappa\begin{pmatrix}1&i\\-i&1\end{pmatrix},\qquad\kappa>0.

Its eigenvalues are 2κ2\kappa and zero. An eigenvector for 2κ2\kappa is (1,−i)/2(1,-i)/\sqrt2, so one jump suffices:

L=κ(F1−iF2)=2κ ∣1⟩⟨0∣.L=\sqrt\kappa(F_1-iF_2)=\sqrt{2\kappa}\,|1\rangle\langle0|.

Thus the imaginary off-diagonal entries encode a directed 0→10\to1 transition. Replacing cc by its entrywise real part would instead produce separate X and Y noise and a different stationary state. This example also explains why checking only the diagonal entries of cc is insufficient; positivity concerns the whole Hermitian matrix.

Mathematical validity does not establish microscopic accuracy, a bath temperature, or a measured jump record. Those require the physical derivation or preparation and monitoring model.

Every Φt\Phi_t in a Lindblad semigroup is a quantum channel and therefore has a Kraus representation. The converse is not true: an arbitrary channel Φ\Phi need not sit inside a continuous semigroup

Φ=Φt∗=et∗L\Phi=\Phi_{t_*}=e^{t_*\mathcal L}

with a time-independent Lindblad generator.

The theorem is about continuous-time semigroups, not about isolated finite-time noise maps. Embedding a channel into a Markovian semigroup is a separate problem.

Many useful master equations have

dρdt=Lt(ρ),\frac{d\rho}{dt} = \mathcal L_t(\rho),

with explicitly time-dependent rates and operators. If Lt\mathcal L_t has Lindblad form with nonnegative rates for every time, then the evolution is CP-divisible under suitable regularity assumptions. But it is not a time-homogeneous semigroup unless Lt\mathcal L_t is constant.

A negative canonical rate in a traceless orthonormal basis prevents the instantaneous generator from having GKSL form. A negative coefficient in a redundant operator list need not: 2γD[Z]−γD[Z]=γD[Z]2\gamma\mathcal D[Z]-\gamma\mathcal D[Z]=\gamma\mathcal D[Z]. Even when a canonical rate is negative, maps from a fixed initial time can remain physical; check them separately. This is common in time-convolutionless master equations.

Many physically important systems have infinite-dimensional Hilbert spaces: oscillators, fields, particles on a line, Brownian motion, and scattering continua. In these settings, formal Lindblad-looking equations can involve unbounded HH or LjL_j.

Then several finite-dimensional shortcuts fail:

  • the generator may be unbounded and defined only on a dense domain;
  • the exponential etLe^{t\mathcal L} may require semigroup theory rather than a matrix exponential;
  • trace preservation can fail if domains and boundary terms are mishandled;
  • complete positivity of a formal differential expression may not imply a well-defined CPTP semigroup;
  • different operator topologies lead to different continuity assumptions.

The bounded-generator Lindblad theorem is still the right reference point, but infinite-dimensional master equations require domain and closability checks that are not visible from the compact formula alone.

Trace preservation of the Lindblad generator

Section titled “Trace preservation of the Lindblad generator”

Show that

Tr⁡L(ρ)=0\operatorname{Tr}\mathcal L(\rho)=0

for

L(ρ)=−iℏ[H,ρ]+∑j(LjρLj†−12{Lj†Lj,ρ}).\mathcal L(\rho) = - \frac{i}{\hbar}[H,\rho] + \sum_j \left( L_j\rho L_j^\dagger - \frac12\{L_j^\dagger L_j,\rho\} \right).
Solution

The trace of the commutator vanishes by cyclicity:

Tr⁡[H,ρ]=0.\operatorname{Tr}[H,\rho]=0.

For one dissipator,

Tr⁡(LρL†)=Tr⁡(L†Lρ)=Tr⁡(ρL†L).\operatorname{Tr}(L\rho L^\dagger) = \operatorname{Tr}(L^\dagger L\rho) = \operatorname{Tr}(\rho L^\dagger L).

Therefore

Tr⁡[LρL†−12{L†L,ρ}]=0.\operatorname{Tr} \left[ L\rho L^\dagger - \frac12\{L^\dagger L,\rho\} \right] = 0.

Summing over jj gives Tr⁡L(ρ)=0\operatorname{Tr}\mathcal L(\rho)=0.

Suppose Φt=etL\Phi_t=e^{t\mathcal L} in finite dimensions. Show that

L=ddtΦt∣t=0.\mathcal L = \left. \frac{d}{dt}\Phi_t \right|_{t=0}.
Solution

Using the matrix exponential,

etL=I+tL+O(t2).e^{t\mathcal L} = \mathcal I + t\mathcal L + O(t^2).

Hence

Φt−It=L+O(t),\frac{\Phi_t-\mathcal I}{t} = \mathcal L+O(t),

so the limit as t↓0t\downarrow0 is L\mathcal L.

Consider

ρ˙=γ(t)D[L]ρ,γ(t)≥0.\dot\rho = \gamma(t)\mathcal D[L]\rho, \qquad \gamma(t)\ge0.

Is this automatically a time-homogeneous semigroup?

Solution

No. Nonnegative γ(t)\gamma(t) supports CP-divisible time-local evolution under suitable regularity assumptions, but the maps generally depend on both start and end times:

Φ(t,s)≠Φ(t−s,0)\Phi(t,s) \ne \Phi(t-s,0)

For a fixed nonzero dissipator, time homogeneity requires the integrated coefficient to depend only on interval length, hence γ(t)\gamma(t) must be constant almost everywhere. If the dissipator vanishes, its coefficient has no physical effect. A time-homogeneous semigroup requires Φt+s=ΦtΦs\Phi_{t+s}=\Phi_t\Phi_s with the same map family indexed only by elapsed time.

Let cabc_{ab} be positive semidefinite and let

cab=∑jλjujaujb∗,λj≥0.c_{ab} = \sum_j\lambda_j u_{ja}u_{jb}^*, \qquad \lambda_j\ge0.

Show how the operators

Lj=λj∑aujaFaL_j = \sqrt{\lambda_j}\sum_a u_{ja}F_a

turn the coefficient-matrix dissipator into a sum of Lindblad dissipators.

Solution

Substitute the definition:

∑jLjρLj†=∑j,a,bλjujaujb∗FaρFb†=∑a,bcabFaρFb†.\sum_j L_j\rho L_j^\dagger = \sum_{j,a,b} \lambda_j u_{ja}u_{jb}^* F_a\rho F_b^\dagger = \sum_{a,b} c_{ab}F_a\rho F_b^\dagger.

Similarly,

∑jLj†Lj=∑a,bcabFb†Fa.\sum_j L_j^\dagger L_j = \sum_{a,b} c_{ab}F_b^\dagger F_a.

The same replacement inside the anticommutator gives exactly the coefficient-matrix dissipator.

  • R. Alicki and K. Lendi, Quantum Dynamical Semigroups and Applications, Springer (1987).
  • H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press (2002). Publisher.
  • E. B. Davies, Quantum Theory of Open Systems, Academic Press (1976).
  • V. Gorini, A. Kossakowski, and E. C. G. Sudarshan, “Completely positive dynamical semigroups of N-level systems,” Journal of Mathematical Physics 17, 821–825 (1976). DOI.
  • G. Lindblad, “On the generators of quantum dynamical semigroups,” Communications in Mathematical Physics 48, 119–130 (1976). DOI.
  • Á. Rivas and S. F. Huelga, Open Quantum Systems: An Introduction, Springer (2012). Author manuscript.