Skip to content

Amplitude-Damping Channel

The amplitude-damping channel is the standard one-qubit channel for irreversible energy relaxation from an excited state ∣1⟩\lvert1\rangle to a lower state ∣0⟩\lvert0\rangle. It is the finite-time channel version of zero-temperature spontaneous emission or qubit T1T_1 decay.

With decay probability γ\gamma, the channel is

Aγ(ρ)=K0ρK0†+K1ρK1†,0≤γ≤1,\mathcal A_\gamma(\rho) = K_0\rho K_0^\dagger + K_1\rho K_1^\dagger, \qquad 0\le\gamma\le1,

where

K0=(1001−γ),K1=(0γ00).K_0 = \begin{pmatrix} 1&0\\ 0&\sqrt{1-\gamma} \end{pmatrix}, \qquad K_1 = \begin{pmatrix} 0&\sqrt{\gamma}\\ 0&0 \end{pmatrix}.

The basis is part of the model. The state ∣1⟩\lvert1\rangle is the state that can decay, and ∣0⟩\lvert0\rangle is the state that receives the population.

Amplitude damping is population relaxation with energy flowing into unobserved degrees of freedom. It is not pure dephasing and not isotropic randomization.

For a density matrix

ρ=(ρ00ρ01ρ10ρ11),\rho = \begin{pmatrix} \rho_{00}&\rho_{01}\\ \rho_{10}&\rho_{11} \end{pmatrix},

the channel gives

Aγ(ρ)=(ρ00+γρ111−γ ρ011−γ ρ10(1−γ)ρ11).\mathcal A_\gamma(\rho) = \begin{pmatrix} \rho_{00}+\gamma\rho_{11} & \sqrt{1-\gamma}\,\rho_{01}\\ \sqrt{1-\gamma}\,\rho_{10} & (1-\gamma)\rho_{11} \end{pmatrix}.

Thus

ρ11⟼(1−γ)ρ11,ρ00⟼ρ00+γρ11,\rho_{11}\longmapsto(1-\gamma)\rho_{11}, \qquad \rho_{00}\longmapsto\rho_{00}+\gamma\rho_{11},

and

ρ01⟼1−γ ρ01.\rho_{01}\longmapsto\sqrt{1-\gamma}\,\rho_{01}.

The population decay and coherence damping are tied together by complete positivity. The same parameter that transfers population from ∣1⟩\lvert1\rangle to ∣0⟩\lvert0\rangle also reduces the coherence involving ∣1⟩\lvert1\rangle.

For the broader model-selection context, see Common Noise Channels.

For repeated damping steps, the survival probabilities multiply; see Channel Composition and Fixed Points for the general composition rule and fixed-state language.

The channel is trace preserving because

K0†K0=(1001−γ),K1†K1=(000γ),K_0^\dagger K_0 = \begin{pmatrix} 1&0\\ 0&1-\gamma \end{pmatrix}, \qquad K_1^\dagger K_1 = \begin{pmatrix} 0&0\\ 0&\gamma \end{pmatrix},

so

K0†K0+K1†K1=I.K_0^\dagger K_0+K_1^\dagger K_1=I.

It is completely positive because it is written in Kraus form. See Kraus Representation for the general structure.

The branch K1K_1 is trace nonincreasing. It represents the decay event in a quantum-jump interpretation. The branch K0K_0 represents the no-jump operation. Their sum is the deterministic channel obtained when the environment record is not kept.

The channel can be realized by an isometry from the system to system plus environment:

∣0⟩∣0⟩E⟼∣0⟩∣0⟩E,∣1⟩∣0⟩E⟼1−γ ∣1⟩∣0⟩E+γ ∣0⟩∣1⟩E.\begin{aligned} \lvert0\rangle\lvert0\rangle_E &\longmapsto \lvert0\rangle\lvert0\rangle_E,\\ \lvert1\rangle\lvert0\rangle_E &\longmapsto \sqrt{1-\gamma}\, \lvert1\rangle\lvert0\rangle_E + \sqrt{\gamma}\, \lvert0\rangle\lvert1\rangle_E. \end{aligned}

The environment state ∣1⟩E\lvert1\rangle_E records that a decay quantum was emitted. Tracing out the environment gives the amplitude-damping channel.

For an input superposition

∣ψ⟩=α∣0⟩+β∣1⟩,\lvert\psi\rangle = \alpha\lvert0\rangle+\beta\lvert1\rangle,

the joint output is

α∣0⟩∣0⟩E+β1−γ ∣1⟩∣0⟩E+βγ ∣0⟩∣1⟩E.\alpha\lvert0\rangle\lvert0\rangle_E + \beta\sqrt{1-\gamma}\, \lvert1\rangle\lvert0\rangle_E + \beta\sqrt{\gamma}\, \lvert0\rangle\lvert1\rangle_E.

The two environment alternatives are orthogonal, so the reduced system state loses coherence with the emitted-decay record. This is a concrete example of reduced dynamics and of a Stinespring dilation.

If the environment is monitored, the two Kraus operators define two possible conditional updates.

The jump probability is

p1=Tr⁡(K1ρK1†)=γρ11.p_1 = \operatorname{Tr}(K_1\rho K_1^\dagger) = \gamma\rho_{11}.

When the jump is observed and p1>0p_1>0, the normalized state is

ρ1=K1ρK1†p1=∣0⟩⟨0∣.\rho_1 = \frac{K_1\rho K_1^\dagger}{p_1} = \lvert0\rangle\langle0\rvert.

The no-jump probability is

p0=Tr⁡(K0ρK0†)=1−γρ11.p_0 = \operatorname{Tr}(K_0\rho K_0^\dagger) = 1-\gamma\rho_{11}.

The corresponding normalized state is

ρ0=K0ρK0†p0.\rho_0 = \frac{K_0\rho K_0^\dagger}{p_0}.

Even the absence of a detected jump changes the state: it reduces the relative amplitude of ∣1⟩\lvert1\rangle. This is a measurement-backaction statement, not just a dissipative one. Without access to the record, the observer uses Aγ(ρ)\mathcal A_\gamma(\rho).

Write

ρ=12(I+xX+yY+zZ),\rho = \frac12 \left( I+xX+yY+zZ \right),

with ∣0⟩\lvert0\rangle at z=+1z=+1 and ∣1⟩\lvert1\rangle at z=−1z=-1. Then amplitude damping acts as

(x,y,z)⟼(1−γ x,1−γ y,(1−γ)z+γ).(x,y,z) \longmapsto \left( \sqrt{1-\gamma}\,x, \sqrt{1-\gamma}\,y, (1-\gamma)z+\gamma \right).

This is an affine map, not a simple contraction about the origin. The Bloch ball is pulled toward the ground-state pole.

In particular,

Aγ ⁣(I2)=12(1+γ001−γ).\mathcal A_\gamma\!\left(\frac{I}{2}\right) = \frac12 \begin{pmatrix} 1+\gamma&0\\ 0&1-\gamma \end{pmatrix}.

Therefore Aγ\mathcal A_\gamma is not unital for γ>0\gamma>0. Since every mixture of unitary conjugations is unital, amplitude damping is not a random-unitary channel unless γ=0\gamma=0.

Using the output-input convention of Choi Matrix, and the ordered basis

∣00⟩, ∣01⟩, ∣10⟩, ∣11⟩,\lvert00\rangle,\, \lvert01\rangle,\, \lvert10\rangle,\, \lvert11\rangle,

where the first slot is the output system and the second slot is the input reference, the Choi matrix is

JAγ=(1001−γ0γ0000001−γ001−γ).J_{\mathcal A_\gamma} = \begin{pmatrix} 1&0&0&\sqrt{1-\gamma}\\ 0&\gamma&0&0\\ 0&0&0&0\\ \sqrt{1-\gamma}&0&0&1-\gamma \end{pmatrix}.

Its nonzero eigenvalues are

2−γ,γ.2-\gamma, \qquad \gamma.

Thus JAγ≥0J_{\mathcal A_\gamma}\ge0 for 0≤γ≤10\le\gamma\le1, and the channel has Kraus rank 22 for 0<γ≤10\lt\gamma\le1. At γ=0\gamma=0 it reduces to the identity channel and has Kraus rank 11.

The trace-preservation condition is

Tr⁡outJAγ=Iin.\operatorname{Tr}_{\mathrm{out}}J_{\mathcal A_\gamma}=I_{\mathrm{in}}.

The unitality condition fails for γ>0\gamma>0:

Tr⁡inJAγ=(1+γ001−γ)≠I.\operatorname{Tr}_{\mathrm{in}}J_{\mathcal A_\gamma} = \begin{pmatrix} 1+\gamma&0\\ 0&1-\gamma \end{pmatrix} \ne I.

For Markovian zero-temperature relaxation with decay rate Γ\Gamma, the finite-time channel has

γ(t)=1−e−Γt.\gamma(t)=1-e^{-\Gamma t}.

The excited-state population obeys

ρ11(t)=e−Γtρ11(0),\rho_{11}(t) = e^{-\Gamma t}\rho_{11}(0),

and the coherence obeys

ρ01(t)=e−Γt/2ρ01(0)\rho_{01}(t) = e^{-\Gamma t/2}\rho_{01}(0)

apart from any Hamiltonian phase rotation or additional pure dephasing.

The corresponding Lindblad equation is

dρdt=Γ(σ−ρσ+−12{σ+σ−,ρ}),\frac{d\rho}{dt} = \Gamma \left( \sigma_-\rho\sigma_+ - \frac12 \{\sigma_+\sigma_-,\rho\} \right),

where

σ−=∣0⟩⟨1∣,σ+=∣1⟩⟨0∣.\sigma_-=\lvert0\rangle\langle1\rvert, \qquad \sigma_+=\lvert1\rangle\langle0\rvert.

In this convention T1=1/ΓT_1=1/\Gamma. Additional pure dephasing changes the transverse coherence time but not the population decay. See Lindblad–GKSL Equation and the Formula Sheet for the general master-equation conventions.

For the generator-level treatment, including jump interpretation, finite-temperature rates, and oscillator damping, see Amplitude Damping Master Equation.

Relation to Dephasing and Depolarizing Noise

Section titled “Relation to Dephasing and Depolarizing Noise”

Amplitude damping differs from dephasing because it changes populations in the energy basis. Pure dephasing leaves ρ00\rho_{00} and ρ11\rho_{11} fixed.

Dephasing and Amplitude Damping owns the QI-facing conversion from protocol-qualified T1–T2 and equilibrium-population records to one physical, composable idle channel with held-out checks; this page retains the zero-temperature channel definition, Kraus and Choi forms, dilation, fixed point, and formal finite-temperature extension.

It differs from depolarizing noise because it is not isotropic and not unital. Depolarizing noise pulls every state toward I/2I/2; amplitude damping pulls states toward ∣0⟩⟨0∣\lvert0\rangle\langle0\rvert.

It differs from a Pauli channel because Pauli channels are unital mixtures of Pauli conjugations. Amplitude damping has an affine shift in the Bloch ball, so it cannot be written as a Pauli channel for γ>0\gamma>0.

Finite Temperature and Generalized Amplitude Damping

Section titled “Finite Temperature and Generalized Amplitude Damping”

The channel above is a zero-temperature model: it allows downward transitions but no upward thermal excitation. If the bath can excite the system, the correct qubit model includes both downward and upward processes.

At the master-equation level, populations then satisfy

ρ˙11=−Γ↓ρ11+Γ↑ρ00,\dot\rho_{11} = -\Gamma_\downarrow\rho_{11} + \Gamma_\uparrow\rho_{00},

with steady excited population

ρ11ss=Γ↑Γ↑+Γ↓.\rho_{11}^{\mathrm{ss}} = \frac{\Gamma_\uparrow} {\Gamma_\uparrow+\Gamma_\downarrow}.

Thermal detailed balance gives

Γ↑Γ↓=e−βℏω\frac{\Gamma_\uparrow}{\Gamma_\downarrow} = e^{-\beta\hbar\omega}

for a two-level spacing ℏω\hbar\omega coupled to a thermal bath. Using zero-temperature amplitude damping when Γ↑\Gamma_\uparrow is appreciable is a physical error, not just a convention choice.

Treating the computational basis as arbitrary

Section titled “Treating the computational basis as arbitrary”

Amplitude damping is tied to an ordered energy basis. Swapping the labels ∣0⟩\lvert0\rangle and ∣1⟩\lvert1\rangle changes the channel.

The population survival factor is 1−γ1-\gamma, while the coherence factor is 1−γ\sqrt{1-\gamma}. Forgetting the square root gives the wrong T1T_1 contribution to coherence decay.

Amplitude damping is nonunital for γ>0\gamma>0. A Pauli channel is unital. They are not the same model.

Using zero-temperature damping at finite temperature

Section titled “Using zero-temperature damping at finite temperature”

If the environment can drive upward transitions, use a thermal master equation or generalized amplitude-damping model rather than the two-Kraus zero-temperature channel.

Treating Kraus operators as unique histories

Section titled “Treating Kraus operators as unique histories”

The two-Kraus representation matches a natural jump/no-jump monitoring picture, but Kraus representations are not unique. The physical trajectory interpretation depends on the measurement record and unraveling.

Amplitude damping keeps the system inside the two-dimensional Hilbert space. Leakage out of the computational subspace requires a different output space or a larger system model; see Erasure and Loss Channels.

Show that the Kraus operators K0K_0 and K1K_1 satisfy the trace-preserving completeness relation.

Solution

Compute

K0†K0=(1001−γ),K_0^\dagger K_0 = \begin{pmatrix} 1&0\\ 0&1-\gamma \end{pmatrix},

and

K1†K1=(000γ).K_1^\dagger K_1 = \begin{pmatrix} 0&0\\ 0&\gamma \end{pmatrix}.

Adding gives

K0†K0+K1†K1=(1001)=I.K_0^\dagger K_0+K_1^\dagger K_1 = \begin{pmatrix} 1&0\\ 0&1 \end{pmatrix} =I.

Starting from the Kraus representation, derive the matrix action of Aγ\mathcal A_\gamma on

ρ=(ρ00ρ01ρ10ρ11).\rho = \begin{pmatrix} \rho_{00}&\rho_{01}\\ \rho_{10}&\rho_{11} \end{pmatrix}.
Solution

The no-jump term is

K0ρK0†=(ρ001−γ ρ011−γ ρ10(1−γ)ρ11).K_0\rho K_0^\dagger = \begin{pmatrix} \rho_{00}&\sqrt{1-\gamma}\,\rho_{01}\\ \sqrt{1-\gamma}\,\rho_{10}&(1-\gamma)\rho_{11} \end{pmatrix}.

The jump term is

K1ρK1†=(γρ11000).K_1\rho K_1^\dagger = \begin{pmatrix} \gamma\rho_{11}&0\\ 0&0 \end{pmatrix}.

Adding them gives

Aγ(ρ)=(ρ00+γρ111−γ ρ011−γ ρ10(1−γ)ρ11).\mathcal A_\gamma(\rho) = \begin{pmatrix} \rho_{00}+\gamma\rho_{11} & \sqrt{1-\gamma}\,\rho_{01}\\ \sqrt{1-\gamma}\,\rho_{10} & (1-\gamma)\rho_{11} \end{pmatrix}.

Compute Aγ(I/2)\mathcal A_\gamma(I/2) and show that the channel is not unital for γ>0\gamma>0.

Solution

Set ρ00=ρ11=1/2\rho_{00}=\rho_{11}=1/2 and ρ01=0\rho_{01}=0 in the matrix formula:

Aγ(I/2)=((1+γ)/200(1−γ)/2).\mathcal A_\gamma(I/2) = \begin{pmatrix} (1+\gamma)/2&0\\ 0&(1-\gamma)/2 \end{pmatrix}.

This equals I/2I/2 only when γ=0\gamma=0. Therefore the channel is not unital for nonzero damping.

For Markovian amplitude damping, use γ(t)=1−e−t/T1\gamma(t)=1-e^{-t/T_1} to show that the coherence factor is e−t/(2T1)e^{-t/(2T_1)}.

Solution

The finite-time channel multiplies the coherence by

1−γ(t).\sqrt{1-\gamma(t)}.

With

γ(t)=1−e−t/T1,\gamma(t)=1-e^{-t/T_1},

we have

1−γ(t)=e−t/T1=e−t/(2T1).\sqrt{1-\gamma(t)} = \sqrt{e^{-t/T_1}} = e^{-t/(2T_1)}.

Thus population relaxation alone contributes a transverse coherence decay rate 1/(2T1)1/(2T_1).

Use the Choi matrix to determine the Kraus rank of the amplitude-damping channel for 0<γ≤10\lt\gamma\le1.

Solution

The Choi matrix has nonzero eigenvalues

2−γ,γ.2-\gamma, \qquad \gamma.

For 0<γ≤10\lt\gamma\le1, both are positive. Therefore the Choi rank is 22, which equals the minimal number of Kraus operators. At γ=0\gamma=0, the eigenvalue γ\gamma vanishes and the channel reduces to the identity, whose Choi rank is 11.

  • K. Kraus, States, Effects, and Operations: Fundamental Notions of Quantum Theory, Springer (1983).
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press (2010).
  • H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press (2002).
  • H. M. Wiseman and G. J. Milburn, Quantum Measurement and Control, Cambridge University Press (2010).
  • C. W. Gardiner and P. Zoller, Quantum Noise, Springer, 3rd ed. (2004).
  • J. Watrous, The Theory of Quantum Information, Cambridge University Press (2018).