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 versus coherences?

  • Find the steady state of the Lindblad equation for amplitude damping with operator \(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?

Lecture Notes#

Overview#

For an isolated system, the Schrödinger equation perfectly captures quantum evolution. But real systems interact with their surroundings. The Lindblad master equation describes how open systems evolve, accounting for both unitary dynamics and irreversible environmental effects. It is the most general Markovian (memoryless) evolution compatible with quantum mechanics.

The Lindblad equation decomposes open-system evolution into: (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.

Physical Example: Spontaneously Emitting Atom

Consider a two-level atom in an excited state \(\vert 1\rangle\) coupled to the electromagnetic field. The excited state decays to the ground state \(\vert 0\rangle\) by emitting a photon at rate \(\gamma\). In the Lindblad picture, the upper level’s population decays as \(\dot{\rho}_{11} = -\gamma\rho_{11}\). Remarkably, coherences between \(\vert 0\rangle\) and \(\vert 1\rangle\) decay at the slower rate \(\dot{\rho}_{01} = -(\gamma/2)\rho_{01}\). This asymmetry—coherences outlive populations by a factor of 2—is the simplest Lindblad process.

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)

(259)#\[ \frac{\mathrm{d}\hat{\rho}}{\mathrm{d}t} = -\frac{\mathrm{i}}{\hbar}[\hat{H}, \hat{\rho}] + \sum_k \gamma_k \left( \hat{L}_k \hat{\rho} \hat{L}_k^\dagger - \frac{1}{2}\{\hat{L}_k^\dagger \hat{L}_k, \hat{\rho}\} \right) \]

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:

Properties of Lindblad Evolution#

The Lindblad equation preserves all fundamental properties of quantum states:

1. Trace preservation: \(\frac{\mathrm{d}}{\mathrm{d}t}\operatorname{Tr}(\hat{\rho}) = 0\), so \(\operatorname{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\).

4. CPTP map: The time evolution is a completely positive trace-preserving (quantum) channel: \(\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), while \(-\frac{1}{2}\{\hat{L}_k^\dagger \hat{L}_k,\hat{\rho}\}\) represents loss (probability outflow). Together they preserve trace:

\[ \operatorname{Tr}\left(\hat{L}_k\hat{\rho} \hat{L}_k^\dagger - \frac{1}{2}\{\hat{L}_k^\dagger \hat{L}_k,\hat{\rho}\}\right) = 0 \]

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\) (spontaneous emission) at rate \(\Gamma_1 = 1/T_1\). For the two-level atom, the Lindblad operator \(\hat{L} = \sqrt{\gamma}\,\vert 0\rangle\langle 1\vert\) describes spontaneous emission at rate \(\gamma\).

(260)#\[ \frac{\mathrm{d}\hat{\rho}}{\mathrm{d}t} = \Gamma_1\!\left(\hat{\sigma}_-\hat{\rho}\hat{\sigma}_+ - \frac{1}{2}\{\hat{\sigma}_+\hat{\sigma}_-,\hat{\rho}\}\right) \]

Effect: \(\dot{\hat{\rho}}_{11} = -\Gamma_1\hat{\rho}_{11}\) (decay), \(\dot{\hat{\rho}}_{00} = +\Gamma_1\hat{\rho}_{11}\) (gain), \(\dot{\hat{\rho}}_{01} = -\frac{\Gamma_1}{2}\hat{\rho}_{01}\) (coherences decay at half rate).

2. Phase Damping: \(L = \sqrt{\Gamma_\varphi}\,\hat{Z}\).

Physical process: Random phase fluctuations (dephasing) with no energy loss; the resulting coherence decays at rate \(2\Gamma_\varphi\), so the pure dephasing time is \(T_\varphi = 1/(2\Gamma_\varphi)\).

(261)#\[ \frac{\mathrm{d}\hat{\rho}}{\mathrm{d}t} = \Gamma_\varphi\!\left(\hat{Z}\hat{\rho}\hat{Z} - \hat{\rho}\right) \]

Effect: Only off-diagonal terms decay: \(\hat{\rho}_{01}(t) = \mathrm{e}^{-2\Gamma_\varphi t}\hat{\rho}_{01}(0)\). Diagonal terms unchanged—this is the microscopic realization of the dephasing channel (§6.3.3, §6.4.1).

3. Depolarizing: \(L_\mu = \sqrt{\Gamma_d/3}\,\hat{\sigma}^\mu\) for \(\mu = x,y,z\).

Physical process: Symmetric noise 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\).

Quantum Jump Picture#

The Lindblad equation averages over an ensemble of quantum trajectories. Each trajectory evolves stochastically:

  1. Between jumps: The state evolves under the non-Hermitian effective Hamiltonian

(262)#\[ \hat{H}_\text{eff} = H - \frac{\mathrm{i}\hbar}{2}\sum_k\gamma_k \hat{L}_k^\dagger \hat{L}_k \]

The imaginary part causes the norm to decay, tracking the “no-jump” probability.

  1. 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\,\operatorname{Tr}(\hat{L}_k^\dagger \hat{L}_k\hat{\rho})\,\mathrm{d}t\).

  2. Averaging over all trajectories recovers the Lindblad master equation.

Physical interpretation: Each quantum jump corresponds to a physical event—photon emission, phonon absorption, spin flip. The non-Hermitian Hamiltonian reduces the state’s norm between jumps, tracking the probability that no jump has occurred. When a jump happens, the state is projected and renormalized. Averaging over many trajectories with different jump sequences recovers the Lindblad master equation.

Summary#

  • Lindblad master equation: \(\frac{\mathrm{d}\hat{\rho}}{\mathrm{d}t} = -\mathrm{i}[\hat{H}, \hat{\rho}] + \sum_k \gamma_k (L_k \hat{\rho} L_k^\dagger - \frac{1}{2}\{L_k^\dagger L_k, \hat{\rho}\})\) (GKSL); describes Markovian open-system dynamics.

  • Jump operators \(L_k\): Represent decoherence and dissipation; decay rates \(\gamma_k \geq 0\). Amplitude damping, dephasing, and decay are modeled by specific choices.

  • Steady state: Long-time limit \(\hat{\rho}(t \to \infty) = \hat{\rho}_{ss}\) satisfies \(\frac{\mathrm{d}\hat{\rho}}{\mathrm{d}t} = 0\). Describes equilibrium or thermal state.

  • Positivity and CPTP: Lindblad form guarantees complete positivity and trace preservation for any \(L_k\) and \(\gamma_k \geq 0\).

See Also

  • 6.4.1 Decoherence: Environment-induced superselection and timescales before Markovian master equations.

  • 6.4.3 Error Correction: Fighting channel noise—dynamical decoupling and quantum error correction as open-system control.

  • 6.3.3 Quantum Channels: CPTP structure and complete positivity constraints on Lindblad forms.

Homework#

1. Amplitude Damping Lindblad. Consider a qubit under amplitude damping, \(L = \sqrt{\Gamma}\hat{\sigma}_-\) where \(\hat{\sigma}_- = \vert 0\rangle\langle 1\vert\).

(a) Show that \(\operatorname{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. Decay Rate Dynamics. 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) \neq 0\), show that \(\hat{\rho}_{01}(t) = \hat{\rho}_{01}(0)\,\mathrm{e}^{-t/(2T_1)}\). Coherences decay at half the relaxation rate.

3. Pure Dephasing Lindblad. A qubit undergoes pure dephasing via \(L = \sqrt{\Gamma_\varphi}\hat{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. Thermal Radiation Coupling. 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{Z}/2}\) (Boltzmann distribution).

5. Quantum Problem. In the quantum jump picture, the system evolves under \(\hat{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. Depolarizing Lindblad Equation. 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) = \operatorname{Tr}(\hat{\rho}(t)^{2})\). Using the identity \(\sum_{\mu}\hat{\sigma}^{\mu}\hat{\rho}\hat{\sigma}^{\mu} = 2\hat{I} - \hat{\rho}\) for any qubit state, show that the Bloch vector decays as \(\dot{\boldsymbol{r}} = -(4\Gamma/3)\boldsymbol{r}\), hence \(\dot{\mathcal{P}} = -(8\Gamma/3)(\mathcal{P} - 1/2)\). Solve to get \(\mathcal{P}(t) = \tfrac{1}{2} + \tfrac{1}{2}\mathrm{e}^{-8\Gamma t/3}\).

(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. Multiple Jump Operators. A Lindblad equation has two operators: decay \(L_1 = \sqrt{\Gamma_1}\hat{\sigma}_-\) and dephasing \(L_2 = \sqrt{\Gamma_\varphi}\hat{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?