6.4.2 Lindblad Master Equation#
Prompts
Write down the Lindblad equation and explain the physical role of each term — the Hamiltonian part, the gain term \(\hat{L}_k\hat{\rho} \hat{L}_k^\dagger\), and the loss term \(-\frac{1}{2}\{\hat{L}_k^\dagger \hat{L}_k,\hat{\rho}\}\).
What are the three canonical Lindblad operators for amplitude damping, phase damping, and depolarizing noise? What does each one do to populations vs. coherences?
Find the steady state of the Lindblad equation for amplitude damping (\(L = \sqrt{\Gamma}\hat{\sigma}_-\)). Why is it the ground state?
Describe the quantum jump picture. How does averaging over stochastic quantum trajectories reproduce the deterministic Lindblad equation?
How does the Lindblad equation relate to quantum channels from §6.3.3? Show they give equivalent descriptions of open-system evolution.
Lecture Notes#
For an isolated system, the Schrödinger equation perfectly captures quantum evolution. But real systems interact with their surroundings — heat baths, photon fields, uncontrolled degrees of freedom. The Lindblad master equation describes how open systems evolve, accounting for both unitary dynamics and the irreversible effects of environmental coupling. It is the most general form of Markovian (memoryless) evolution compatible with quantum mechanics.
Overview#
The Lindblad equation decomposes open-system evolution into two parts: (1) unitary evolution under the system Hamiltonian \(\hat{H}\), and (2) dissipation described by jump operators \(\{\hat{L}_k\}\) with rates \(\{\gamma_k\}\). The equation preserves trace, Hermiticity, and positivity — necessary conditions for valid quantum dynamics. Every quantum channel (§6.3.3) is obtained from a Lindblad equation integrated over a fixed time interval.
Discussion: Markovian Approximation
The Lindblad equation assumes the system and environment are uncorrelated at each instant — a Markovian (memoryless) approximation. But real environments can have memory: the environment’s response depends on the system’s past state. How do we know when Markovian is justified?
The approximation breaks down if the environment correlation time \(\tau_E\) is comparable to the system timescale. For typical quantum optics setups, \(\tau_E \sim 10^{-15}\) s while atomic evolution is \(\sim 10^{-9}\) s — safely Markovian.
Some systems show “non-Markovian revivals” — coherence temporarily returns after decoherence. Does this invalidate the Lindblad equation, or is it simply outside its regime?
Lindblad Master Equation#
For an open system coupled to an environment, the reduced density matrix \(\hat{\rho}_S\) evolves as:
Lindblad Master Equation (GKSL Form)
where:
\(H\) = system Hamiltonian (Hermitian), governing unitary evolution.
\(\hat{L}_k\) = Lindblad operators (jump operators); generally non-Hermitian.
\(\gamma_k \geq 0\) = decay rates (non-negative real numbers).
\(\{A, B\} = AB + BA\) = anticommutator.
This is the most general Markovian master equation compatible with complete positivity and trace preservation (Lindblad 1976; Kossakowski, Sudarshan — GKSL form).
The Lindblad equation is derived by tracing out the environment from unitary system+environment evolution:
Derivation: Lindblad Equation
Start with total Hamiltonian:
The full system+environment evolves unitarily. Tracing out the environment and assuming the Markovian approximation (short environment memory, \(\tau_E \ll \tau_S\)) gives:
The dissipation superoperator \(\mathcal{D}\) must preserve trace and positivity. The most general form satisfying these constraints is precisely the Lindblad form \(\sum_k\gamma_k(\hat{L}_k\hat{\rho} \hat{L}_k^\dagger - \frac{1}{2}\{\hat{L}_k^\dagger \hat{L}_k,\hat{\rho}\})\).
Properties of Lindblad Evolution#
Properties of the Lindblad Equation
The Lindblad equation preserves all fundamental properties of quantum states:
1. Trace preservation: \(\frac{\mathrm{d}}{\mathrm{d}t}\text{Tr}(\hat{\rho}) = 0\), so \(\text{Tr}(\hat{\rho}(t)) = 1\) for all \(t\).
2. Hermiticity preservation: If \(\hat{\rho}(0) = \hat{\rho}(0)^\dagger\), then \(\hat{\rho}(t) = \hat{\rho}(t)^\dagger\).
3. Positivity: If \(\hat{\rho}(0) \geq 0\), then \(\hat{\rho}(t) \geq 0\) for all \(t\) (valid quantum state).
4. CPTP map: The time evolution \(\hat{\rho}(0) \mapsto \hat{\rho}(t)\) is a completely positive trace-preserving (quantum) channel, equivalent to the Kraus representation \(\hat{\rho}(t) = \sum_\mu K_\mu \hat{\rho}(0) K_\mu^\dagger\) with \(\sum_\mu K_\mu^\dagger K_\mu = I\).
Gain and loss balance: The \(\hat{L}_k\hat{\rho} \hat{L}_k^\dagger\) term represents gain (probability inflow from quantum jumps), while \(-\frac{1}{2}\{\hat{L}_k^\dagger \hat{L}_k,\hat{\rho}\}\) represents loss (probability outflow). Their combined action preserves trace:
Common Lindblad Operators#
1. Amplitude Damping: \(L = \sqrt{\Gamma_1}\,\hat{\sigma}_-\)#
Physical process: Excited state \(\vert 1\rangle\) decays to ground state \(\vert 0\rangle\) (energy dissipation, spontaneous emission) at rate \(\Gamma_1 = 1/T_1\).
Effect on populations: \(\dot{\hat{\rho}}_{11} = -\Gamma_1\hat{\rho}_{11}\) (exponential decay), \(\dot{\hat{\rho}}_{00} = +\Gamma_1\hat{\rho}_{11}\) (gain).
Effect on coherences: \(\dot{\hat{\rho}}_{01} = -\frac{\Gamma_1}{2}\hat{\rho}_{01}\) (decay at half the rate).
2. Phase Damping: \(L = \sqrt{\Gamma_\varphi}\,\hat{\sigma}^z\)#
Physical process: Random phase fluctuations (dephasing) with no energy loss at rate \(\Gamma_\varphi = 1/T_\varphi\).
Effect: Only off-diagonal terms decay: \(\hat{\rho}_{01}(t) = \mathrm{e}^{-2\Gamma_\varphi t}\hat{\rho}_{01}(0)\). Diagonal (population) terms unchanged. This is the microscopic realization of the dephasing channel (§6.3.3 and §6.4.1).
3. Depolarizing: \(L_\mu = \sqrt{\Gamma_d/3}\,\hat{\sigma}^\mu\) for \(\mu = x,y,z\)#
Physical process: Symmetric noise on all three axes; generic random error channel at rate \(\Gamma_d\).
Effect: Drives any state toward the maximally mixed state \(I/2\) on timescale \(1/\Gamma_d\).
Steady States#
A steady state \(\hat{\rho}_\infty\) satisfies \(\mathcal{L}(\hat{\rho}_\infty) = 0\), i.e., \(\frac{\mathrm{d}\hat{\rho}_\infty}{\mathrm{d}t} = 0\).
Example: Amplitude Damping
Problem. Find the steady state for amplitude damping \(L = \sqrt{\Gamma_1}\hat{\sigma}_-\).
Solution. From \(\dot{\hat{\rho}}_{11} = -\Gamma_1\hat{\rho}_{11} = 0\) we get \(\hat{\rho}_{11}^\infty = 0\), hence \(\hat{\rho}_{00}^\infty = 1\).
The steady state is \(\hat{\rho}_\infty = \vert 0\rangle\langle 0\vert\) — all population has decayed to the ground state.
Example: Thermal Steady State
Problem. For an atom coupled to a thermal bath at temperature \(T\), include both decay (\(L_1 = \sqrt{\Gamma\bar{n}_\text{th}+\Gamma}\,\hat{\sigma}_-\)) and absorption (\(L_2 = \sqrt{\Gamma\bar{n}_\text{th}}\,\hat{\sigma}_+\)). Find the steady state.
Solution. Balancing gain and loss: \(\hat{\rho}_{11}^\infty/\hat{\rho}_{00}^\infty = \bar{n}_\text{th}/(\bar{n}_\text{th}+1) = \mathrm{e}^{-\hbar\omega/k_BT}\), which is the Boltzmann factor. The steady state is the thermal density matrix \(\hat{\rho}_\infty = \mathrm{e}^{-\beta \hat{H}}/Z\), consistent with §6.1.2.
Quantum Jump Picture#
The Lindblad equation averages over an ensemble of quantum trajectories. Each individual trajectory evolves stochastically:
Between jumps: The state evolves under the non-Hermitian effective Hamiltonian
(non-Hermitian: the norm decreases, tracking the “no-jump” probability.)
At random times (Poisson process), a quantum jump \(\hat{L}_k\) occurs: \(\hat{\rho} \to \hat{L}_k\hat{\rho} \hat{L}_k^\dagger/p_k\) with probability \(p_k \propto \gamma_k\,\text{Tr}(\hat{L}_k^\dagger \hat{L}_k\hat{\rho})\,\mathrm{d}t\).
Averaging over all trajectories (all possible jump sequences and times) recovers the Lindblad master equation.
Physical interpretation: Each quantum jump corresponds to a physical event — photon emission, phonon absorption, spin flip. Decoherence and dissipation arise from random jumps; the deterministic master equation is the statistical average.
Summary#
The Lindblad equation \(\dot{\hat{\rho}} = -\frac{\mathrm{i}}{\hbar}[\hat{H},\hat{\rho}] + \sum_k\gamma_k(\hat{L}_k\hat{\rho} \hat{L}_k^\dagger - \frac{1}{2}\{\hat{L}_k^\dagger \hat{L}_k,\hat{\rho}\})\) is the most general Markovian master equation for open quantum systems.
Lindblad operators \(\hat{L}_k\) represent quantum jumps; rates \(\gamma_k \geq 0\); the equation preserves trace, Hermiticity, and positivity (CPTP).
Amplitude damping (\(L = \hat{\sigma}_-\)): decay to ground state. Phase damping (\(L = \hat{\sigma}^z\)): coherence decay without energy loss. Depolarizing (\(L_\mu = \hat{\sigma}^\mu\)): drives state to \(I/2\).
Steady states satisfy \(\mathcal{L}(\hat{\rho}_\infty) = 0\); for thermal baths the steady state is the Boltzmann distribution \(\mathrm{e}^{-\beta \hat{H}}/Z\).
The quantum jump picture: Lindblad evolution is the average over stochastic trajectories where jumps occur at random times.
See Also
6.4.1 Decoherence: Microscopic origin of dephasing and decoherence
6.3.3 Quantum Channels: Channels as the integrated form of Lindblad evolution
6.1.2 Entropy: Thermal states as Lindblad steady states
Homework#
1. Consider a qubit under amplitude damping, \(L = \sqrt{\Gamma}\hat{\sigma}_-\) where \(\hat{\sigma}_- = \vert 0\rangle\langle 1\vert\).
(a) Show that \(\text{Tr}(\mathrm{d}\hat{\rho}/\mathrm{d}t) = 0\) (trace is conserved). Use the cyclic property of trace.
(b) Verify \(\hat{\sigma}_+\hat{\sigma}_- = \vert 1\rangle\langle 1\vert\) (projector onto excited state). Use this to write the loss term explicitly.
(c) Show that if \(\hat{\rho}(0)\) is Hermitian, then \(\mathrm{d}\hat{\rho}/\mathrm{d}t\) is also Hermitian (Hermiticity is preserved).
2. A qubit decays via amplitude damping at rate \(\Gamma_1 = 1/T_1\). Start with \(\hat{\rho}(0) = \vert 1\rangle\langle 1\vert\).
(a) Write the Lindblad master equation explicitly for populations \(\hat{\rho}_{00}(t)\) and \(\hat{\rho}_{11}(t)\).
(b) Solve to find \(\hat{\rho}_{11}(t) = \mathrm{e}^{-t/T_1}\) and \(\hat{\rho}_{00}(t) = 1 - \mathrm{e}^{-t/T_1}\).
(c) Show the steady state is \(\hat{\rho}_\infty = \vert 0\rangle\langle 0\vert\) (all population in the ground state).
(d) Starting from \(\hat{\rho}(0)\) with coherence \(\hat{\rho}_{01}(0) eq 0\), show that \(\hat{\rho}_{01}(t) = \hat{\rho}_{01}(0)\,\mathrm{e}^{-t/(2T_1)}\). Coherences decay at half the relaxation rate.
3. A qubit undergoes pure dephasing via \(L = \sqrt{\Gamma_\varphi}\hat{\sigma}^z\). Start with \(\hat{\rho}(0) = \vert+\rangle\langle+\vert\).
(a) Write \(\hat{\rho}(0)\) explicitly as a \(2\times 2\) matrix.
(b) Apply the Lindblad superoperator: show that diagonal elements (populations) are unchanged and off-diagonal elements decay as \(\hat{\rho}_{01}(t) = \mathrm{e}^{-2\Gamma_\varphi t}\hat{\rho}_{01}(0)\).
(c) As \(t \to \infty\), show that \(\hat{\rho}_\infty = I/2\) (maximally mixed state). Explain physically: why does dephasing destroy coherence while preserving energy?
4. An atom couples to a thermal radiation field at temperature \(T\). The relevant processes are spontaneous emission at rate \(A\), absorption at rate \(Bn_\text{th}\) (where \(n_\text{th} = 1/(\mathrm{e}^{\hbar\omega/k_BT}-1)\)), and dephasing at rate \(\gamma_\varphi\).
(a) Write the three Lindblad operators \(L_1\), \(L_2\), \(L_3\) for these processes.
(b) At \(T = 0\): \(n_\text{th} \to 0\), so only spontaneous emission survives. Write the simplified master equation and find the steady state.
(c) At high temperature (\(k_BT \gg \hbar\omega\)): \(n_\text{th} \approx k_BT/\hbar\omega\). Show that the steady state approaches \(\hat{\rho}_\infty \propto \mathrm{e}^{-\beta\hbar\omega\hat{\sigma}^z/2}\) (Boltzmann distribution).
5. In the quantum jump picture, the system evolves under \(H_\text{eff} = -\frac{\mathrm{i}\hbar\Gamma}{2}\vert 1\rangle\langle 1\vert\) (for amplitude damping) between jumps.
(a) For initial state \(\vert\psi(0)\rangle = \vert 1\rangle\), show that the no-jump evolution gives \(\vert\psi(t)\rangle \propto \mathrm{e}^{-\Gamma t/2}\vert 1\rangle\). The squared norm \(\Vert\psi(t)\Vert^2 = \mathrm{e}^{-\Gamma t}\) is the probability that no jump has occurred.
(b) If a jump occurs at time \(t\), the state becomes \(\hat{\sigma}_-\vert\psi(t)\rangle \propto \vert 0\rangle\) (ground state). Interpret this physically: what is the jump event?
(c) Show that averaging over all possible trajectories (no jump with probability \(\mathrm{e}^{-\Gamma t}\), one jump at time \(\tau\) with probability \(\Gamma\mathrm{e}^{-\Gamma\tau}\mathrm{d}\tau\)) reproduces the Lindblad master equation solution \(\hat{\rho}(t) = \mathrm{e}^{-\Gamma t}\vert 1\rangle\langle 1\vert + (1 - \mathrm{e}^{-\Gamma t})\vert 0\rangle\langle 0\vert\).
6. A qubit undergoes depolarizing noise via \(L_\mu = \sqrt{\Gamma/3}\hat{\sigma}^\mu\) for \(\mu = x,y,z\).
(a) Show that \(\hat{\rho}_\infty = I/2\) is a steady state: \(\mathcal{L}(I/2) = 0\). (Hint: Use \(\hat{\sigma}^\mu(I/2)\hat{\sigma}^\mu = I/2\) for each Pauli.)
(b) Starting from \(\hat{\rho}(0) = \vert 0\rangle\langle 0\vert\), define purity \(\mathcal{P}(t) = \text{Tr}(\hat{\rho}(t)^2)\). Show that \(\frac{\mathrm{d}\mathcal{P}}{\mathrm{d}t} = -2\Gamma(\mathcal{P} - 1/2)\). Solve to get \(\mathcal{P}(t) = \frac{1}{2} + \frac{1}{2}\mathrm{e}^{-2\Gamma t}\).
(c) Identify the Kraus operators \(\{K_\mu\}\) for the depolarizing channel (integrated Lindblad evolution for short time \(\mathrm{d}t\)). Show the channel matches the form from §6.3.3.
7. A Lindblad equation has two operators: decay \(L_1 = \sqrt{\Gamma_1}\hat{\sigma}_-\) and dephasing \(L_2 = \sqrt{\Gamma_\varphi}\hat{\sigma}^z\).
(a) Write the full master equation for \(\dot{\hat{\rho}}_{11}\) and \(\dot{\hat{\rho}}_{01}\). Show that populations decay at rate \(\Gamma_1\) and coherences decay at rate \(\Gamma_1/2 + 2\Gamma_\varphi\).
(b) Identify the dephasing time \(T_2\) from the coherence decay rate. Show that \(\frac{1}{T_2} = \frac{1}{2T_1} + \frac{1}{T_\varphi}\) where \(T_1 = 1/\Gamma_1\) and \(T_\varphi = 1/(2\Gamma_\varphi)\).
(c) Explain physically why \(T_2 \leq 2T_1\). Under what condition is this bound saturated?