Skip to content

A Lindblad operator identifies which states a dissipative term connects or distinguishes. Its normalization determines the observable decay rate. This reference translates common operators into population and coherence dynamics; a detector interpretation additionally requires a monitoring scheme.

Required background. The Lindblad–GKLS Equation defines the generator, and Lindblad Operators explains its nonunique representations.

Use

D[L]ρ=LρL†−12{L†L,ρ}.\mathcal D[L]\rho=L\rho L^\dagger-\frac12\{L^\dagger L,\rho\}.

Rates may be separate, γD[L]\gamma\mathcal D[L], or absorbed, D[γL]\mathcal D[\sqrt\gamma L]. When LL is dimensionless, γ\gamma has units of inverse time; the absorbed operator has units of inverse square root of time. If LL itself is dimensional, the coefficient must compensate. For example, a coefficient of D[x]\mathcal D[x] has units 1/(length2 time)1/(\text{length}^2\,\text{time}).

The formulas below give only the dissipative contribution. A Hamiltonian can add rotations, change populations in the displayed basis, or move states out of a protected subspace.

Let σ−=∣g⟩⟨e∣\sigma_-=|g\rangle\langle e|, σ+=∣e⟩⟨g∣\sigma_+=|e\rangle\langle g|, and Z=∣g⟩⟨g∣−∣e⟩⟨e∣Z=|g\rangle\langle g|-|e\rangle\langle e|.

Term Population change Coherence decay
Γ↓D[σ−]\Gamma_\downarrow\mathcal D[\sigma_-] p˙e=−Γ↓pe\dot p_e=-\Gamma_\downarrow p_e ρ˙eg=−Γ↓ρeg/2\dot\rho_{eg}=-\Gamma_\downarrow\rho_{eg}/2
Γ↑D[σ+]\Gamma_\uparrow\mathcal D[\sigma_+] p˙e=Γ↑(1−pe)\dot p_e=\Gamma_\uparrow(1-p_e) ρ˙eg=−Γ↑ρeg/2\dot\rho_{eg}=-\Gamma_\uparrow\rho_{eg}/2
(Γϕ/2)D[Z](\Gamma_\phi/2)\mathcal D[Z] Populations fixed ρ˙eg=−Γϕρeg\dot\rho_{eg}=-\Gamma_\phi\rho_{eg}

With both transitions and pure dephasing,

1T1=Γ↓+Γ↑,1T2=12T1+Γϕ.\begin{gathered} \frac1{T_1}=\Gamma_\downarrow+\Gamma_\uparrow,\\ \frac1{T_2}=\frac1{2T_1}+\Gamma_\phi. \end{gathered}

These are the constant nonnegative-rate, two-level Markov formulas. They should not be applied indiscriminately to nonexponential decay, ensemble inhomogeneity, or echo times. An equilibrium bath at inverse temperature β\beta gives Γ↑/Γ↓=e−βℏω0\Gamma_\uparrow/\Gamma_\downarrow=e^{-\beta\hbar\omega_0} for gap ℏω0\hbar\omega_0. Arbitrary pumping rates need not represent such a bath. See Amplitude Damping Master Equation and Thermal Master Equations.

For Pe=∣e⟩⟨e∣P_e=|e\rangle\langle e| there is the exact two-level identity

2ΓϕD[Pe]=Γϕ2D[Z].2\Gamma_\phi\mathcal D[P_e] =\frac{\Gamma_\phi}{2}\mathcal D[Z].

It explains a common factor-of-four difference between projector and Pauli operator coefficients. The coherence-decay rate is the same. See Pure Dephasing Master Equation.

Oscillator loss, gain, and number monitoring

Section titled “Oscillator loss, gain, and number monitoring”

For a bosonic mode, a∣n⟩=n ∣n−1⟩a|n\rangle=\sqrt n\,|n-1\rangle and N=a†aN=a^\dagger a.

Term Mean occupation Mean field
κD[a]\kappa\mathcal D[a] d⟨N⟩/dt=−κ⟨N⟩d\langle N\rangle/dt=-\kappa\langle N\rangle d⟨a⟩/dt=−κ⟨a⟩/2d\langle a\rangle/dt=-\kappa\langle a\rangle/2
gD[a†]g\mathcal D[a^\dagger] d⟨N⟩/dt=g(⟨N⟩+1)d\langle N\rangle/dt=g(\langle N\rangle+1) d⟨a⟩/dt=g⟨a⟩/2d\langle a\rangle/dt=g\langle a\rangle/2

The +1+1 in the gain law is spontaneous excitation. The ideal linear amplifier has no saturation and its occupation grows without bound at long times. A thermal reservoir combines loss and gain with coefficients κ(nˉ+1)\kappa(\bar n+1) and κnˉ\kappa\bar n; then d⟨N⟩/dt=−κ(⟨N⟩−nˉ)d\langle N\rangle/dt=-\kappa(\langle N\rangle-\bar n).

Number dephasing leaves all Fock populations fixed:

(γND[N]ρ)mn=−γN2(m−n)2ρmn.\bigl(\gamma_N\mathcal D[N]\rho\bigr)_{mn} =-\frac{\gamma_N}{2}(m-n)^2\rho_{mn}.

Large number separations lose coherence faster. This term describes phase diffusion or unread number monitoring, not photon loss. The identities assume the required moments exist; a finite Fock cutoff changes oscillator commutators at its upper boundary. Check occupation near that boundary as described in Solving Lindblad Equations.

Projector, Pauli, and collective operators

Section titled “Projector, Pauli, and collective operators”

For any orthogonal projector PP, γPD[P]\gamma_P\mathcal D[P] damps coherences between the PP and I−PI-P sectors at rate γP/2\gamma_P/2, while preserving each diagonal block. Such a generator can represent continuous unread monitoring. One instantaneous nonselective projective measurement is the channel ρ↦PρP+(I−P)ρ(I−P)\rho\mapsto P\rho P+(I-P)\rho(I-P), rather than this finite-rate generator. A POVM effect alone specifies neither process; see Quantum Instruments.

Independent Pauli terms ∑j=x,y,zγjD[σj]\sum_{j=x,y,z}\gamma_j\mathcal D[\sigma_j] give

r˙x=−2(γy+γz)rx,\dot r_x=-2(\gamma_y+\gamma_z)r_x,

with cyclic versions for ry,rzr_y,r_z. Equal coefficients γj=γ\gamma_j=\gamma therefore shrink the Bloch vector at rate 4γ4\gamma, not γ\gamma. See Pauli Channels for the corresponding finite-time map.

Collective emission uses J−=∑jσ−(j)J_-=\sum_j\sigma_-^{(j)}. Its dissipator contains cross terms between emitters and differs from ∑jD[σ−(j)]\sum_j\mathcal D[\sigma_-^{(j)}]. For two identical qubits, the antisymmetric single-excitation state is dark under J−J_- but decays under independent loss. Remaining dark during the full evolution also requires that the Hamiltonian preserve the dark sector; see Decoherence-Free Subspaces.

Equivalent operators and different records

Section titled “Equivalent operators and different records”

For rate-absorbed operators, a unitary mixing Mα=∑kUαkLkM_\alpha=\sum_k U_{\alpha k}L_k leaves ∑αD[Mα]=∑kD[Lk]\sum_\alpha\mathcal D[M_\alpha]=\sum_k\mathcal D[L_k] exactly unchanged. Identity shifts can also be compensated by a Hamiltonian change. Thus an operator list does not uniquely identify microscopic couplings or individual clicks. Unravelings explains how different monitoring choices can produce the same averaged equation.

Read a coefficient. A model contains 3 ms−1D[Z]3\,\mathrm{ms}^{-1}\mathcal D[Z]. What is the pure coherence-decay time?

Solution

The off-diagonal entry of ZρZ−ρZ\rho Z-\rho is −2ρeg-2\rho_{eg}. Its decay rate is 6 ms−16\,\mathrm{ms}^{-1}, so Tϕ=1/6 msT_\phi=1/6\,\mathrm{ms}. The dissipator coefficient is not itself the coherence-decay rate.

Number separation. Compare the decay rates of ρ01\rho_{01} and ρ03\rho_{03} under γND[N]\gamma_N\mathcal D[N].

Solution

They decay at γN/2\gamma_N/2 and 9γN/29\gamma_N/2, respectively. Their populations remain unchanged. A single quoted dephasing time cannot describe all number coherences without specifying the pair.

  • Breuer, H.-P., and F. Petruccione. The Theory of Open Quantum Systems. Oxford University Press, 2002. Qubit and oscillator generators.
  • Gardiner, C. W., and P. Zoller. Quantum Noise. 3rd edition, Springer, 2004. Loss, amplification, and phase diffusion.
  • Lindblad, G. “On the generators of quantum dynamical semigroups.” Communications in Mathematical Physics 48, 119–130 (1976). doi:10.1007/BF01608499.
  • Wiseman, H. M., and G. J. Milburn. Quantum Measurement and Control. Cambridge University Press, 2010. Monitoring and alternative unravelings.