Closed quantum systems evolve by unitary operators. Open systems exchange information and energy with surroundings, so their reduced states can decohere, relax and approach steady states. The Lindblad equation is the canonical generator of continuous-time Markovian quantum dynamics when complete positivity, trace preservation and semigroup composition are required.
Guide 6 introduced quantum channels through Kraus operators. A Lindblad generator is the infinitesimal version of a continuous channel semigroup. Instead of specifying one map ρ→𝒩(ρ), we specify a superoperator ℒ satisfying
dρ/dt=ℒ(ρ)
so that ρ(t)=e^{tℒ}ρ(0).
This guide develops the GKSL form, amplitude damping, dephasing, vectorisation, Liouvillian spectra, steady-state multiplicity, relaxation gaps, quantum-jump trajectories, dark states, thermal generators and dissipative phase transitions. The Markovian assumption is explicit: not every open quantum process admits a time-independent Lindblad semigroup description.
System + memoryless environment model → CPTP semigroup → Lindblad generator → Liouvillian spectrum → steady states and decay modes.
1. Density matrices are required
A system entangled with an environment is generally mixed after the environment is ignored.
Therefore open-system dynamics acts on density operators ρ rather than only state vectors.
A physical finite-time map must preserve positivity even when the system is entangled with an untouched reference. That stronger requirement is complete positivity.
Trace preservation ensures total probability remains one.
2. Quantum dynamical semigroup
A time-homogeneous Markovian semigroup is a family {𝒯t} for t≥0 satisfying
𝒯_0=I;𝒯_{t+s}=𝒯_t𝒯_s;- each 𝒯t is completely positive and trace preserving;
- the map is suitably continuous in t.
Then
𝒯_t=e^{tℒ}
for a generator ℒ.
The semigroup law expresses absence of explicit memory: evolution over the next interval depends only on the present reduced state, not on how that state was reached.
3. GKSL/Lindblad form
For finite-dimensional systems, the most general generator of a norm-continuous CPTP semigroup can be written
ℒ(ρ)=−i[H,ρ]+Σ_k(L_kρL_k†−(1/2){L_k†L_k,ρ}).
H is Hermitian. The Lk are Lindblad or jump operators. Positive rates can be absorbed into Lk by writing Lk=√γkAk. [1,2]
The anticommutator is {A,B}=AB+BA.
4. Hamiltonian and dissipative parts
The commutator term
−i[H,ρ]
is ordinary unitary Schrödinger evolution written for density matrices.
The dissipator
𝒟[L](ρ)=LρL†−(1/2){L†L,ρ}
contains both a jump term LρL† and a compensating anticommutator that preserves trace and produces the correct no-jump evolution.
5. Verify trace preservation
Use cyclicity of trace:
Tr(LρL†)=Tr(L†Lρ).
The anticommutator contributes
−(1/2)Tr(L†Lρ)−(1/2)Tr(ρL†L)=−Tr(L†Lρ).
Thus each dissipator has zero trace. The commutator also has zero trace.
Therefore
dTrρ/dt=0.
6. Amplitude damping
For a two-level atom that spontaneously decays from |e⟩ to |g⟩ at rate γ, choose
L=√γ σ_- = √γ |g⟩⟨e|.
With H=0, the master equation gives
dρ_ee/dt=−γρ_ee;dρ_gg/dt=γρ_ee;dρ_eg/dt=−(γ/2)ρ_eg.
Population decays at rate γ while coherence decays at rate γ/2 when this is the only noise process.
7. Solve amplitude damping
The equations integrate to
ρ_ee(t)=e^{-γt}ρ_ee(0)
and
ρ_eg(t)=e^{-γt/2}ρ_eg(0).
Trace preservation gives
ρ_gg(t)=1−ρ_ee(t)
for an initially normalised qubit.
The unique steady state is |g⟩⟨g|.
8. Worked damping time
Let γ=1/(20 μs)=0.05 μs−1.
After t=40 μs,
ρ_ee(t)/ρ_ee(0)=e^{-2}≈0.1353.
The coherence factor is
e^{-1}≈0.3679.
Population and coherence therefore have different decay constants even within one simple Lindblad process.
9. Pure dephasing
Choose one jump operator
L=√κ Z.
Since Z²=I, the dissipator becomes
κ(ZρZ−ρ).
Diagonal populations are unchanged. Off-diagonal elements satisfy
dρ_eg/dt=−2κρ_eg.
This guide therefore defines dephasing rate of the coherence itself as 2κ for this exact jump-operator convention. Other texts choose L=√(γφ/2)Z so the displayed coherence-decay rate is γφ.
10. T1 and T2 relation
If amplitude damping at rate γ₁ and independent pure coherence dephasing at rate γφ act together in the convention where off-diagonal pure dephasing is −γφρeg, then
1/T₁=γ₁
and
1/T₂=γ₁/2+γφ.
Therefore T₂≤2T₁ for this standard Markovian qubit model.
A measured violation of a simple relation may indicate additional physics, non-Markovianity, model mismatch or inconsistent rate conventions rather than a failure of quantum mechanics.
11. Finite-time Kraus form of amplitude damping
Let η=e−γt. The amplitude-damping channel can be written with Kraus operators
K_0=|g⟩⟨g|+√η|e⟩⟨e|
and
K_1=√(1−η)|g⟩⟨e|.
They satisfy K_0†K_0+K_1†K_1=I.
Expanding these Kraus operators for small dt recovers the Lindblad generator with jump √γσ−.
12. Lindblad representation is not unique
Different sets of jump operators can generate the same ℒ.
Unitary mixing among jumps leaves the total dissipator invariant. Shifting a jump by a multiple of identity can be compensated by changing H.
Therefore an individual Lk should not automatically be interpreted as a uniquely real physical microscopic event unless a particular unraveling or environment measurement scheme has been specified.
13. Liouvillian as a linear operator on operator space
A d×d density matrix has d² operator components.
Vectorise ρ by stacking its matrix entries into a d²-component vector |ρ⟩⟩.
Then the master equation becomes an ordinary linear differential equation
d|ρ⟩⟩/dt=𝕃|ρ⟩⟩
where 𝕃 is a d²×d² matrix representation of the Liouvillian superoperator.
Unlike a Hamiltonian, 𝕃 is generally non-Hermitian.
14. Vectorisation identity
For column-stacking vectorisation,
vec(AρB)=(B^T⊗A)vec(ρ).
Therefore the Hamiltonian commutator becomes
−i(I⊗H−H^T⊗I)
under this stacking convention.
Each dissipator similarly becomes a Kronecker-product matrix. This makes exact diagonalisation of small Liouvillians straightforward with ordinary linear algebra.
15. Steady state
A steady state satisfies
ℒ(ρ_ss)=0.
In vector form, |ρss⟩⟩ is a right eigenvector of 𝕃 with eigenvalue zero.
Trace preservation guarantees a corresponding left zero mode associated with the identity operator:
⟨⟨I|𝕃=0.
A physical steady-state solution must additionally be Hermitian, positive and trace one.
16. Multiple steady states
The zero eigenvalue can be degenerate.
This occurs with conserved quantities, disconnected sectors, decoherence-free subspaces or symmetry-protected steady-state manifolds.
When steady states are not unique, the long-time state can retain memory of the initial condition through its projection onto the stationary manifold.
One should not speak of “the steady state” until uniqueness has been checked.
17. Liouvillian eigenmodes
Suppose
ℒ(R_j)=λ_jR_j.
Then an eigenmode contributes
e^{λ_jt}R_j
to the time evolution.
For a stable finite-dimensional CPTP semigroup, nonstationary eigenvalues have nonpositive real parts. Negative real parts produce decay; imaginary parts produce damped oscillation.
18. Liouvillian spectrum of amplitude damping
For the H=0 amplitude-damping qubit above, convenient operator eigenmodes have eigenvalues
- 0 for the stationary ground-state mode;
- −γ for population relaxation;
- −γ/2 and −γ/2 for the two coherence quadratures.
The slowest exponential decay rate is γ/2.
The Liouvillian gap under the simple spectral definition is therefore γ/2.
19. Liouvillian spectral gap
For a unique steady state, define a common asymptotic gap
Δ_L=−max_{λ_j≠0}Re λ_j.
A positive ΔL means every nonstationary eigenmode decays exponentially asymptotically.
In simple normal Liouvillians, 1/ΔL is a natural longest relaxation timescale.
For non-normal Liouvillians, eigenvectors can be highly nonorthogonal and large transient amplification or long crossover times can make the actual approach to steady state much slower than the naive inverse-gap estimate. [5]
20. Non-normality
A matrix is normal if it commutes with its adjoint.
Liouvillians need not be normal. Their right eigenmodes can be strongly nonorthogonal, so large coefficients can cancel initially and reveal slow behaviour only after a long transient.
Pseudospectra, singular values or symmetrised generators can then give useful information not visible from eigenvalues alone.
Thus “Liouvillian gap equals relaxation time” is a heuristic requiring structural conditions, not a universal identity.
21. Dark states
A pure state |ψ⟩ is dark to jump Lk if
L_k|ψ⟩=0
for all relevant jumps.
If it is also invariant under H up to phase, then |ψ⟩⟨ψ| is a steady state.
Engineered dissipation can deliberately make a desired entangled or bosonic code state dark, causing environmental coupling to stabilise rather than destroy the target manifold.
22. Decoherence-free subspace
If every state in a subspace experiences the same trivial action from all relevant noise operators, information encoded inside that subspace can be immune to the corresponding decoherence channel.
For collective dephasing L∝ΣZj, states sharing the same total Z eigenvalue acquire the same environmental phase information and can form decoherence-free encodings.
This is passive error avoidance, distinct from active syndrome extraction in Guides 12, 25, 26, 41 and 42.
23. Quantum-jump unraveling
The same Lindblad master equation can be represented as an ensemble average over stochastic pure-state trajectories.
Between jumps, evolve with the non-Hermitian effective Hamiltonian
H_eff=H−(i/2)Σ_kL_k†L_k.
A jump k occurs with infinitesimal probability
dp_k=dt⟨ψ|L_k†L_k|ψ⟩
and updates the state proportionally to Lk|ψ⟩.
Averaging trajectories reproduces ρ(t).
24. Unraveling depends on environment monitoring
Photon counting suggests discrete quantum jumps. Homodyne monitoring of the same emitted field can produce a diffusive stochastic trajectory instead.
The unconditional density matrix may obey the same master equation while individual conditioned records differ.
Trajectories are therefore physical relative to a specified measurement scheme on the environment, not unique hidden paths demanded by the Lindblad equation itself.
25. Thermal Lindblad dynamics
A two-level system coupled weakly to a thermal bath can have downward and upward jumps
L_-=√γ_- σ_-;L_+=√γ_+ σ_+.
Detailed balance at transition frequency ω requires
γ_+/γ_-=e^{-βω}
for a bath satisfying the corresponding equilibrium assumptions.
The stationary excited/ground population ratio is then e−βω, producing the Gibbs state for the qubit.
26. Davies generators
Under weak system–bath coupling, rapid bath-correlation decay and secular approximations, one can derive a Davies-type Lindblad generator whose jump operators connect system energy eigenspaces at Bohr frequencies.
Kubo–Martin–Schwinger/detailed-balance relations in the bath set the upward/downward rates.
This provides a microscopic route to Gibbs thermalisation, but the approximations can fail under strong coupling, dense spectra, structured environments or significant memory.
27. Non-Markovian dynamics
If the environment retains information and later returns it to the system, a time-independent semigroup can be inadequate.
One may encounter time-local master equations with time-dependent rates, memory-kernel equations or enlarged-system Markovian embeddings.
Negative instantaneous decay rates in some time-local representations can signal information backflow, but precise non-Markovianity definitions vary.
“Not Lindblad semigroup” does not automatically mean “unphysical”. It may mean the reduced dynamics has memory.
28. Dissipative state engineering
Choose jump operators so an intended target is the unique dark steady state.
Instead of fighting every environmental coupling, engineer a reservoir that removes entropy from unwanted modes while leaving the target invariant.
Examples include optical pumping, autonomous stabilisation of cat manifolds, dissipative entanglement generation and preparation of stabilizer states.
The design problem becomes: choose ℒ whose stationary manifold has the desired structure and whose nonzero spectrum gives sufficiently rapid convergence.
29. Dissipative phase transitions
In a sequence of many-body open systems of increasing size, a qualitative steady-state change can sharpen into a dissipative phase transition.
The Liouvillian gap often closes in the thermodynamic limit, making relaxation increasingly slow near the transition.
Steady-state observables, correlations and fluctuations can become nonanalytic as control parameters cross critical values.
Finite systems remain smooth except at special degeneracies; scaling with system size is required to establish a thermodynamic dissipative transition.
30. Exceptional points
Because Liouvillians are non-Hermitian, eigenvalues and eigenvectors can coalesce at exceptional points.
At such a point the Liouvillian may become nondiagonalizable and acquire Jordan blocks.
Time evolution can then contain polynomial factors multiplying exponentials, such as t e^{λt}, rather than a sum of independent pure exponentials.
Exceptional points are a specifically non-Hermitian spectral feature and should not be confused with ordinary Hermitian degeneracies.
31. Common misconception: every environment produces Lindblad dynamics
The GKSL form exactly characterises continuous Markovian CPTP semigroups. Structured reservoirs, strong coupling and memory can require more general reduced dynamics.
32. Common misconception: jump operators are unique microscopic events
The Lindblad representation has gauge/unitary freedom. A particular jump interpretation becomes operationally meaningful only after an environment-monitoring unraveling is specified.
33. Common misconception: inverse Liouvillian gap always equals mixing time
It captures asymptotic eigenmode decay in simple settings. Strong non-normality, poor eigenvector conditioning and long transient regimes can make finite-time relaxation much slower than 1/ΔL.
34. Worked synthesis problem
A qubit undergoes amplitude damping with γ=0.04 ns−1 and no Hamiltonian.
Step 1: T₁. T₁=1/γ=25 ns.
Step 2: Coherence lifetime from damping alone. Off-diagonal elements decay at γ/2, so T₂=2/γ=50 ns if no additional dephasing exists.
Step 3: Liouvillian eigenvalues. The decay spectrum is {0,−0.04,−0.02,−0.02} ns^{-1}.
Step 4: Spectral gap. Δ_L=0.02 ns^{-1}, so the asymptotic slow-mode timescale is 50 ns.
Step 5: Steady state. The unique stationary density matrix is |g⟩⟨g|. Adding a thermal excitation jump would change both the steady state and the Liouvillian spectrum.
35. Practice set
- What mathematical properties define a quantum dynamical semigroup here?
- Write the GKSL/Lindblad generator.
- Why does the dissipator preserve trace?
- What jump describes spontaneous decay of a two-level system?
- How fast does coherence decay under pure amplitude damping?
- What does a Z jump do to populations?
- What does vectorisation accomplish?
- Define a steady state.
- What is a Liouvillian eigenmode?
- Define the common Liouvillian spectral gap for a unique steady state.
- Why can non-normality invalidate a naive mixing-time estimate?
- What is a dark state?
Answers
- CPTP maps with identity at t=0, semigroup composition and suitable continuity.
−i[H,ρ]+Σ(LρL†−{L†L,ρ}/2).- The positive jump trace exactly cancels the two anticommutator trace terms; commutators also have zero trace.
√γ|g⟩⟨e|.- At rate γ/2 for off-diagonal density-matrix elements.
- Leaves them unchanged in a pure-dephasing model while suppressing coherences.
- Turns the superoperator equation on d×d matrices into an ordinary d²-dimensional linear differential equation.
- A physical density operator satisfying
ℒ(ρ_ss)=0. - An operator R_j satisfying
ℒ(R_j)=λ_jR_j, contributinge^{λ_jt}. −max_{λ≠0}Reλ.- Nonorthogonal eigenmodes can generate large transient/crossover effects not captured by eigenvalue real parts alone.
- A state annihilated by the relevant jumps and invariant under the coherent dynamics strongly enough to remain stationary.
Sources and further study
[1] Göran Lindblad, On the Generators of Quantum Dynamical Semigroups, Communications in Mathematical Physics 48, 119–130 (1976). The foundational generator theorem.
[2] Vittorio Gorini, Andrzej Kossakowski and E. C. G. Sudarshan, Completely Positive Dynamical Semigroups of N-Level Systems, Journal of Mathematical Physics 17, 821–825 (1976). The finite-dimensional GKSL structure theorem.
[3] Heinz-Peter Breuer and Francesco Petruccione, The Theory of Open Quantum Systems, Oxford University Press (2002). A standard reference on microscopic master-equation derivations, Markovian dynamics and open-system techniques.
[4] Daniel Manzano, A short introduction to the Lindblad master equation, AIP Advances 10, 025106 (2020). An accessible mathematical introduction including vectorised Liouville-space methods.
[5] Takashi Mori and Tatsuhiko Shirai, Symmetrized Liouvillian Gap in Markovian Open Quantum Systems, Physical Review Letters 130, 230404 (2023). A rigorous demonstration that the standard Liouvillian gap need not control finite transient relaxation in nonequilibrium non-normal settings.
Batch 13 series navigation
- Guide 49: Berry Phases, Berry Curvature, Chern Numbers and Quantum Geometry
- Guide 50: Topological Bands, Quantum Hall Mathematics, Edge States and Bulk–Boundary Correspondence
- Guide 51: Quantum Field Theory, Second Quantisation, Fock Space, Wick’s Theorem and Propagators
- Guide 52: Lindblad Master Equations, Liouvillian Spectra, Steady States and Dissipative Quantum Dynamics
- Return to the BTT Mathematics Learning Hub
Educational note: Lindblad equations are exact generators of Markovian CPTP semigroups. Applying one to laboratory data requires a model of coupling, timescales and environmental memory; a convenient fit does not by itself prove microscopic Markovianity.