|q⟩ Bad Qubits

intermediate · Programming · Noise, Decoherence & Density-Matrix Simulation

Amplitude Damping (T₁)

Every physical qubit is weakly coupled to its surroundings. An excited qubit (1|1\rangle) can spontaneously emit a photon and relax to the ground state (0|0\rangle). This energy dissipation is called amplitude damping, and its characteristic timescale is T1T_1 (the longitudinal or energy-relaxation time).

Kraus operators of the amplitude-damping channel

For a damping strength γ[0,1]\gamma \in [0, 1], the amplitude-damping channel is defined by two Kraus operators

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}.

You can verify trace-preservation directly:

K0K0+K1K1=(1001γ)+(000γ)=I.K_0^\dagger K_0 + K_1^\dagger K_1 = \begin{pmatrix}1 & 0 \\ 0 & 1-\gamma\end{pmatrix} + \begin{pmatrix}0 & 0 \\ 0 & \gamma\end{pmatrix} = I. \checkmark

The channel acts on a density matrix ρ\rho as

E(ρ)=K0ρK0+K1ρK1.\mathcal{E}(\rho) = K_0 \rho K_0^\dagger + K_1 \rho K_1^\dagger.

K1K_1 is a lowering operator: it maps 1|1\rangle to 0|0\rangle and annihilates 0|0\rangle. Its presence means that excited-state population drains into the ground state at rate γ\gamma.

What happens to each element of the density matrix

Let ρ=(ρ00ρ01ρ10ρ11)\rho = \begin{pmatrix}\rho_{00} & \rho_{01} \\ \rho_{10} & \rho_{11}\end{pmatrix}. Carrying out the matrix products gives

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

Two effects are visible:

At γ=1\gamma = 1 the channel is maximally damping: the qubit always ends in 0|0\rangle regardless of the input, and all coherences vanish.

Effect on the excited state

For the pure excited state ρ=11=(0001)\rho = |1\rangle\langle 1| = \begin{pmatrix}0 & 0 \\ 0 & 1\end{pmatrix}, the formula above simplifies to

E(11)=(γ001γ).\mathcal{E}(|1\rangle\langle 1|) = \begin{pmatrix}\gamma & 0 \\ 0 & 1-\gamma\end{pmatrix}.

The qubit has relaxed to a classical mixture: with probability γ\gamma it is now in 0|0\rangle and with probability 1γ1-\gamma it remains in 1|1\rangle. There is no superposition because the off-diagonal entries of 11|1\rangle\langle 1| were zero to begin with.

Effect on a superposition

A superposition +=12(0+1)|+\rangle = \tfrac{1}{\sqrt{2}}(|0\rangle + |1\rangle) has the density matrix

++=12(1111).|{+}\rangle\langle{+}| = \frac{1}{2}\begin{pmatrix}1 & 1 \\ 1 & 1\end{pmatrix}.

After the amplitude-damping channel with strength γ\gamma:

E(++)=12(1+γ1γ1γ1γ).\mathcal{E}(|{+}\rangle\langle{+}|) = \frac{1}{2}\begin{pmatrix}1 + \gamma & \sqrt{1-\gamma} \\ \sqrt{1-\gamma} & 1-\gamma\end{pmatrix}.

The probability of measuring 0|0\rangle is the (0,0)(0,0) entry:

P(0)=1+γ2=ρ00+γρ11.P(0) = \frac{1 + \gamma}{2} = \rho_{00} + \gamma \cdot \rho_{11}.

For γ=0.5\gamma = 0.5 this gives P(0)=0.75P(0) = 0.75: the excited-state component has half a chance of falling to the ground state, raising the ground-state probability above the ideal 0.50.5.

Connection to T1T_1 on real hardware

On a physical device the damping strength γ\gamma grows with idle time tt:

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

A qubit left idle for tT1t \ll T_1 suffers little energy loss (γt/T1\gamma \approx t/T_1). For t=T1t = T_1, γ0.63\gamma \approx 0.63; for tT1t \gg T_1 the qubit is almost certain to be found in 0|0\rangle. Typical superconducting qubits have T1T_1 in the range of tens to hundreds of microseconds (Krantz et al., 2019; cited in the Qiskit Textbook), which constrains the maximum depth of useful circuits.

Try it

The amplitude-damping channel with γ=0.5\gamma = 0.5 is applied to +|{+}\rangle. Using the formula P(0)=ρ00+γρ11P(0) = \rho_{00} + \gamma \cdot \rho_{11}, compute and return the probability of measuring 0|0\rangle after the channel acts.

Run your code to see the quantum state.

Sign in on the full site to ask questions and join the discussion.