Discretization of Errors

Classical error correction is intrinsically digital. A bit is flipped or it is not; on nn bits there are only 2n2^n error patterns, and a code that handles each one handles them all. Quantum information looks hopelessly worse off. A qubit lives on the Bloch sphere, and an unintended rotation exp⁡(−iθ n^⋅σ⃗/2)\exp(-i\theta\,\hat n\cdot\vec\sigma/2) can nudge its state by an arbitrarily small, continuous amount; real decoherence is a whole continuum of completely positive maps, none of them a clean bit flip. The naive fear is fatal: correcting a continuum of errors would seem to demand infinite precision and infinite resources, and quantum computing would be impossible before it began. The discretization of errors is the single observation that dissolves this fear, the hinge on which fault-tolerant quantum computing turns. It rests on two facts that look almost too cheap to matter. First, the four Pauli matrices {I,X,Y,Z}\{I,X,Y,Z\} are a basis for all 2×22\times2 operators, so every single-qubit error — unitary or not — is a linear combination of these four. Second, error correction is linear: the Knill–Laflamme conditions that decide correctability are bilinear in the error operators, so a code correcting each element of a set automatically corrects every operator in that set's complex linear span. Together, a code built to defeat only XX, YY, and ZZ on each qubit defeats every continuous single-qubit error for free. The engine that performs the digitization is measurement: syndrome extraction projects the entangled superposition of error branches onto one definite Pauli branch that recovery undoes with no residual angle — the analog error has become a digital event with a definite probability. This article makes that precise, from the Pauli basis through span closure to the collapse of a full decohering channel.

Key takeaways
  • A single-qubit error is a point in a continuum — an arbitrary 2×22\times2 operator or a whole CPTP channel — yet a code corrects only a finite, discrete set of errors and still wins.
  • The escape is algebraic: {I,X,Y,Z}\{I,X,Y,Z\} spans every single-qubit operator and correctability is linear, so correcting these four basis errors corrects their entire complex span.
  • Syndrome measurement physically digitizes the noise: it projects the entangled error superposition onto one definite Pauli branch that recovery undoes exactly. The continuous parameter survives only as an outcome probability.
  • A coherent over-rotation by θ=0.1\theta=0.1 rad on the three-qubit code collapses to a clean bit flip with probability sin⁡2(θ/2)=2.50×10−3\sin^2(\theta/2)=2.50\times10^{-3}, and correction leaves zero residual rotation.
  • Digitization is the precondition for the threshold theorem: analog noise becomes discrete failure events with definite probabilities, so error rates compose like classical ones.
  • Key idea: because {I,X,Y,Z}\{I,X,Y,Z\} spans all single-qubit operators and error correction is linear, correcting these four discrete errors corrects every continuous error; syndrome measurement physically collapses the continuous error amplitude onto one definite Pauli branch.

Where this lives today

Discretization is running, right now, inside every quantum error-correction experiment. When Google's 2023–2024 surface-code processors showed a distance-5 code suppressing logical errors below distance 3, the decoder was fed a stream of discrete syndrome bits, extracted by measuring stabilizers on superconducting qubits whose real noise is a continuous mix of relaxation, dephasing, and control error. That analog noise never reaches the decoder: syndrome extraction digitizes it into Pauli events, and a classical decoder — minimum-weight matching, or a neural network — does the rest. Quantum error correction is, thanks to discretization, a classical inference problem in a quantum coat. The same principle governs the frontier. Randomized benchmarking reports a coherent over-rotation as a fidelity number precisely because, inside a code, that error becomes an incoherent Pauli flip — the N2→NN^2\to N conversion of Worked Example 1 is what stops control errors accumulating across a deep circuit. Where discretization stops is equally live: leakage out of the computational subspace is not a Pauli and demands dedicated leakage-reduction units, and bosonic cat and GKP codes engineer the hardware so the dominant error — photon loss — maps onto a small, correctable, discrete set. The design question is always Shor's: find a finite correctable set whose span swallows the physical noise, and let measurement do the digitizing.

The physical picture

Picture the space of single-qubit errors as a four-dimensional complex vector space with axes II, XX, YY, ZZ. Any error you can name — a feather-light coherent rotation, a violent dephasing kick, a Kraus operator of some exotic channel — is just a point (or sum of points) in this space, and a code that "covers" the four basis directions automatically covers their entire span. Continuity lives in the coefficients, and the coefficients are exactly what linearity handles. The decisive move is the collapse. Before measurement the corrupted codeword is a coherent superposition over error branches, a0 ∣ψL⟩+aX X∣ψL⟩+⋯a_0\,|\psi_L\rangle + a_X\,X|\psi_L\rangle + \cdots — the continuous error has entangled the data with a continuum of error amplitudes. Syndrome measurement asks a sharp, discrete question — which error eigenspace are you in? — and forces the state onto one branch. The tiny amplitude aXa_X never appears in the surviving state; it sets only the probability ∣aX∣2|a_X|^2 of landing in the XX-branch. The first figure follows this collapse for a real over-rotation on the three-qubit bit-flip code.

Figure
One continuous over-rotation on the three-qubit code becomes one of two discrete branches: outcome (+1,+1) leaves the exact codeword with probability cos²(θ/2), and outcome (−1,+1) leaves a clean bit flip X₁ with probability sin²(θ/2), which recovery undoes
A continuous error entangles the codeword; measurement collapses it to one discrete Pauli branch encoded state |ψ_L⟩ = α|000⟩+β|111⟩ apply R_X₁(θ) continuous angle θ coherent superposition cos(θ/2) · |ψ_L⟩ + −i sin(θ/2) · X₁|ψ_L⟩ measure Z₁Z₂, Z₂Z₃ outcome (+1,+1) state = |ψ_L⟩ exactly prob = cos²(θ/2) outcome (−1,+1) state = X₁|ψ_L⟩ prob = sin²(θ/2) recover X₁ → |ψ_L⟩ cos²(θ/2) + sin²(θ/2) = 1: θ vanishes from both final states, surviving only in the probabilities.

Discretization lets quantum error correction inherit the power of classical coding theory: each round of syndrome extraction stamps the noise into one of finitely many discrete Pauli outcomes whose probabilities compose like classical failure rates, so the threshold theorem can suppress them toward zero — the continuous accumulation that doomed analog computers never occurs. And the principle is not limited to unitary rotations: a general decohering channel has every Kraus operator expanded in Paulis, and syndrome measurement leaves an effective Pauli channel with definite per-error probabilities. The second figure shows this for amplitude damping — the paradigm of messy analog decoherence — reduced to four numbers.

Figure
Amplitude damping with γ = 0.1, after encoding and syndrome measurement, becomes an effective Pauli channel: the digitized Pauli weights are p_I = 0.949, p_X = 0.025, p_Y = 0.025, p_Z = 6.6×10⁻⁴, and they sum to 1 (probability shown on a logarithmic axis)
The punchline: analog amplitude damping (γ = 0.1) digitized into four Pauli probabilities weight p_a = Σ_k of the squared coefficient of P_a in each Kraus operator 10⁰ 10⁻¹ 10⁻² 10⁻³ probability p_a 0.949 0.025 0.025 6.6×10⁻⁴ I X Y Z p_I + p_X + p_Y + p_Z = 1 — a continuous channel reduced to a discrete Pauli channel

The core equations

The whole principle is powered by one fact: the Pauli matrices are a basis. Any single-qubit error operator — Hermitian, unitary, or neither — is a unique combination of the four, its coefficients read off by an inner product.

E=c0 I+c1 X+c2 Y+c3 Z,ck=12 Tr(σk†E)(1)E = c_0\, I + c_1\, X + c_2\, Y + c_3\, Z,\qquad c_k=\tfrac12\,\mathrm{Tr}(\sigma_k^\dagger E) \qquad (1)

where: EE — an arbitrary 2×22\times2 single-qubit error operator (need not be unitary or Hermitian); I,X,Y,ZI,X,Y,Z — the four single-qubit Pauli matrices; ckc_k — complex expansion coefficients. Reading: the four Paulis carry all the information in any single-qubit map, so handling them means handling everything. That the coefficients are unique — that {I,X,Y,Z}\{I,X,Y,Z\} is a genuine basis — follows from orthonormality under the Hilbert–Schmidt inner product ⟨A,B⟩=12Tr(A†B)\langle A,B\rangle=\tfrac12\mathrm{Tr}(A^\dagger B).

12 Tr(σj†σk)=δjk(2)\tfrac12\,\mathrm{Tr}(\sigma_j^\dagger\sigma_k)=\delta_{jk} \qquad (2)

where: σj∈{I,X,Y,Z}\sigma_j\in\{I,X,Y,Z\} — the Pauli operators; δjk\delta_{jk} — the Kronecker delta. Reading: four orthonormal vectors in the four-dimensional matrix space form a complete basis, so exactly four errors suffice. Real noise is a channel: the operator-sum (Kraus) form writes any CPTP map as a sum over Kraus operators, each expanded in Paulis.

E(ρ)=∑kEk ρ Ek†,∑kEk†Ek=I,Ek=∑acka Pa(3)\mathcal{E}(\rho)=\sum_k E_k\,\rho\,E_k^\dagger,\quad \sum_k E_k^\dagger E_k=I,\qquad E_k=\sum_a c_{ka}\,P_a \qquad (3)

where: E\mathcal{E} — the CPTP noise map; EkE_k — its Kraus operators; ρ\rho — the qubit's density operator; Pa∈{I,X,Y,Z}P_a\in\{I,X,Y,Z\}; ckac_{ka} — the Pauli coefficient of PaP_a in EkE_k. Reading: every EkE_k is a Pauli combination, so the whole channel lives in the Pauli span — no continuous part hides elsewhere. Correctability is decided by the Knill–Laflamme conditions: a code with orthonormal codewords ∣ψLi⟩|\psi_L^i\rangle corrects an error set {Ea}\{E_a\} exactly when these hold.

⟨ψLi∣Ea†Eb∣ψLj⟩=Cab δij(4)\langle \psi_L^i | E_a^\dagger E_b | \psi_L^j\rangle = C_{ab}\,\delta_{ij} \qquad (4)

where: ∣ψLi⟩|\psi_L^i\rangle — orthonormal logical codewords; Ea,EbE_a,E_b — errors in the correctable set (e.g. weight-≤t\le t Paulis); CabC_{ab} — a Hermitian matrix independent of i,ji,j. Reading: the left side is bilinear in the errors, and that bilinearity — not anything special about Paulis — makes the span argument work. The generic continuous error is the coherent rotation, and (1) makes its Pauli content explicit.

Rn^(θ)=e−iθ(n^⋅σ⃗)/2=cos⁡θ2 I−isin⁡θ2 (n^⋅σ⃗)(5)R_{\hat n}(\theta)=e^{-i\theta(\hat n\cdot\vec\sigma)/2}=\cos\tfrac{\theta}{2}\,I - i\sin\tfrac{\theta}{2}\,(\hat n\cdot\vec\sigma) \qquad (5)

where: n^\hat n — unit rotation axis; θ\theta — continuous rotation angle; σ⃗=(X,Y,Z)\vec\sigma=(X,Y,Z). Reading: however small θ\theta is, Rn^(θ)R_{\hat n}(\theta) sits in span{I,X,Y,Z}\mathrm{span}\{I,X,Y,Z\}, so a code for the four Paulis corrects it. Measurement is the step that digitizes: projecting onto the syndrome-ss subspace collapses the superposition, the Born rule setting the outcome probability.

p(s)=∥Ms∣Ψ⟩∥2=∣as∣2,∣ψpost⟩=Ms∣Ψ⟩p(s)=Pa(s)∣ψL⟩(6)p(s)=\big\|M_s|\Psi\rangle\big\|^2=|a_s|^2,\qquad |\psi_{\text{post}}\rangle=\frac{M_s|\Psi\rangle}{\sqrt{p(s)}}=P_{a(s)}|\psi_L\rangle \qquad (6)

where: MsM_s — projector onto the syndrome-ss subspace; ∣Ψ⟩=∑aasPa∣ψL⟩|\Psi\rangle=\sum_a a_s P_a|\psi_L\rangle — the corrupted superposition; asa_s — amplitude of the branch carrying syndrome ss; Pa(s)P_{a(s)} — the unique Pauli tagged by syndrome ss. Reading: the continuous amplitude asa_s appears only inside p(s)p(s); the survivor is a definite Pauli error. Applying the matched recovery to every syndrome branch returns the exact logical state.

ρrec=∑sRs Ms E(ρL) Ms Rs†=ρL,Rs=Pa(s)†(7)\rho_{\text{rec}}=\sum_s R_s\,M_s\,\mathcal{E}(\rho_L)\,M_s\,R_s^\dagger=\rho_L,\qquad R_s=P_{a(s)}^\dagger \qquad (7)

where: RsR_s — recovery for syndrome ss, the inverse of the tagged Pauli; ρL\rho_L — the recovered logical state. Reading: the formal statement that "correct the four Paulis" implies "correct the whole channel" (Derivation C).

Derivation ACorrecting a spanning set corrects its entire linear span

The goal is to prove the algebraic heart of discretization: if a code satisfies the Knill–Laflamme conditions on a finite set {Ea}\{E_a\}, it automatically satisfies them on every operator in that set's complex linear span — so a code built against the four Paulis is a code against every continuous single-qubit error, since equation (1) puts every such error inside the span. We feed a general span element through the Knill–Laflamme test and watch the codeword-index dependence reduce to one constant.

Step 1. Let F=∑aαaEaF=\sum_a \alpha_a E_a be an arbitrary span element with complex weights αa\alpha_a — exactly the shape of a continuous error like Rn^(θ)R_{\hat n}(\theta). Why: writing the continuous error as a combination of correctable basis errors is the only input needed.

Step 2. Form the Knill–Laflamme element and expand bilinearly, ⟨ψLi∣F†F∣ψLj⟩=∑a,bαa∗αb ⟨ψLi∣Ea†Eb∣ψLj⟩\langle\psi_L^i|F^\dagger F|\psi_L^j\rangle=\sum_{a,b}\alpha_a^*\alpha_b\,\langle\psi_L^i|E_a^\dagger E_b|\psi_L^j\rangle. Why: sesquilinearity pulls the coefficients outside the codeword bracket.

Step 3. Insert the basis conditions (4), giving ∑a,bαa∗αb Cab δij=(∑a,bαa∗αbCab)δij\sum_{a,b}\alpha_a^*\alpha_b\,C_{ab}\,\delta_{ij}=\big(\sum_{a,b}\alpha_a^*\alpha_b C_{ab}\big)\delta_{ij}. Why: all codeword-index dependence funnels through the single δij\delta_{ij}.

Step 4. So FF satisfies Knill–Laflamme with a new constant C′=∑a,bαa∗αbCabC'=\sum_{a,b}\alpha_a^*\alpha_b C_{ab} independent of i,ji,j — it is correctable. Since the four Paulis span all single-qubit operators, every single-qubit error is such an FF. Why: a finite correctable set implies a continuous correctable span — the whole content of "discretization."

Go deeper · Derivation B · 1 minCollapse of a small over-rotation on the three-qubit bit-flip code
Derivation BCollapse of a small over-rotation on the three-qubit bit-flip code

The goal is to watch the span argument in a concrete circuit: we apply a continuous over-rotation RX1(θ)R_{X_1}(\theta) — a systematic control error leaking amplitude into a bit-flip branch — to qubit 1 of an encoded state and show that measuring the two stabilizers projects onto either the exact codeword or a single clean bit flip, so recovery removes the error with no residual angle. We use ∣ψL⟩=α∣000⟩+β∣111⟩|\psi_L\rangle=\alpha|000\rangle+\beta|111\rangle with stabilizers S1=Z1Z2S_1=Z_1Z_2 and S2=Z2Z3S_2=Z_2Z_3.

Step 1. Expand the rotation with (5): RX1(θ)=cos⁡θ2 I−isin⁡θ2 X1R_{X_1}(\theta)=\cos\tfrac{\theta}{2}\,I-i\sin\tfrac{\theta}{2}\,X_1, an arbitrarily small error parametrized by θ\theta. Why: the "analog" object we must digitize.

Step 2. Act on the codeword: RX1(θ)∣ψL⟩=cos⁡θ2 (α∣000⟩+β∣111⟩)−isin⁡θ2 (α∣100⟩+β∣011⟩)R_{X_1}(\theta)|\psi_L\rangle=\cos\tfrac{\theta}{2}\,(\alpha|000\rangle+\beta|111\rangle)-i\sin\tfrac{\theta}{2}\,(\alpha|100\rangle+\beta|011\rangle), a superposition of the no-error and X1X_1 branches. Why: the continuous amplitude is now physically in the state.

Step 3. The no-error branch is a +1+1 eigenstate of both stabilizers, syndrome (+1,+1)(+1,+1); the X1X_1 branch flips Z1Z2Z_1Z_2 but not Z2Z3Z_2Z_3, syndrome (−1,+1)(-1,+1) — orthogonal eigenspaces. Why: distinct syndromes are what let measurement tell the branches apart.

Step 4. Outcome (+1,+1)(+1,+1) occurs with probability cos⁡2θ2\cos^2\tfrac{\theta}{2}, collapsing to α∣000⟩+β∣111⟩=∣ψL⟩\alpha|000\rangle+\beta|111\rangle=|\psi_L\rangle exactly — no correction, no leftover angle. Why: the state has forgotten θ\theta.

Step 5. Outcome (−1,+1)(-1,+1) occurs with probability sin⁡2θ2\sin^2\tfrac{\theta}{2}, collapsing to α∣100⟩+β∣011⟩=X1∣ψL⟩\alpha|100\rangle+\beta|011\rangle=X_1|\psi_L\rangle; recovery applies X1X_1, and X1(X1∣ψL⟩)=∣ψL⟩X_1(X_1|\psi_L\rangle)=|\psi_L\rangle. Why: the error is now one discrete Pauli the code undoes.

Step 6. Check: cos⁡2θ2+sin⁡2θ2=1\cos^2\tfrac\theta2+\sin^2\tfrac\theta2=1, both branches recover ∣ψL⟩|\psi_L\rangle, and the residual rotation is zero for any θ\theta. Why: the continuous parameter survives only as the branching probability.

Go deeper · Derivation C · 1 minA general decohering channel is digitized into a Pauli channel
Derivation CA general decohering channel is digitized into a Pauli channel

The goal is to lift the picture from a single rotation to an arbitrary CPTP channel — amplitude damping, dephasing, anything — and show that measurement plus recovery makes the code experience an effective Pauli channel with definite per-error probabilities, the strongest form of discretization. The mechanism: the projectors annihilate the coherences between distinct Pauli branches, leaving diagonal Pauli-error terms each undone by its matched recovery. We assume a nondegenerate code, so distinct Paulis carry distinct syndromes.

Step 1. Expand each Kraus operator, Ek=∑ackaPaE_k=\sum_a c_{ka}P_a, so E(ρL)=∑k∑a,bckackb∗ PaρLPb†\mathcal{E}(\rho_L)=\sum_k\sum_{a,b}c_{ka}c_{kb}^*\,P_a\rho_L P_b^\dagger. Why: the cross-terms a≠ba\ne b are coherences between error branches — the analog residue.

Step 2. The syndrome action is MsPa∣ψL⟩=δs,s(a) Pa∣ψL⟩M_s P_a|\psi_L\rangle=\delta_{s,s(a)}\,P_a|\psi_L\rangle: each Pauli sends the code space into one distinct syndrome subspace, and the projectors are orthogonal. Why: distinct Paulis, distinct orthogonal subspaces.

Step 3. Sandwiching between MsM_s, the cross-term MsPaρLPb†MsM_s P_a\rho_L P_b^\dagger M_s vanishes unless s(a)=s(b)s(a)=s(b), which for a nondegenerate code forces a=ba=b. Why: measurement erases the interference between branches — the act of digitization.

Step 4. Apply recovery Rs(a)=Pa†R_{s(a)}=P_a^\dagger and sum; since Rs(a)Pa=IR_{s(a)}P_a=I, the survivor is ∑a(∑k∣cka∣2)ρL=∑apa ρL\sum_a\big(\sum_k|c_{ka}|^2\big)\rho_L=\sum_a p_a\,\rho_L with pa=∑k∣cka∣2p_a=\sum_k|c_{ka}|^2. Why: a discrete Pauli channel applying PaP_a with probability pap_a — the second figure.

Step 5. Trace preservation ∑kEk†Ek=I\sum_k E_k^\dagger E_k=I gives ∑apa=1\sum_a p_a=1, so ρrec=ρL\rho_{\text{rec}}=\rho_L exactly. Why: the probabilities sum to one — perfect correction of an arbitrary continuous channel.

Limiting cases

The principle is powerful but not unconditional: it holds when the error lies in the span of a correctable set, and outside that it weakens or fails.

CaseResultMeaning
Single-qubit, correctable span{Ea}={I,X,Y,Z}i\{E_a\}=\{I,X,Y,Z\}_iIdeal case: the four Paulis span all single-qubit operators and are mutually correctable, so every continuous single-qubit error is corrected.
Coherent over-rotationR=cos⁡θ2I−isin⁡θ2ZR=\cos\tfrac{\theta}{2}I-i\sin\tfrac{\theta}{2}ZLives in span{I,Z}\mathrm{span}\{I,Z\}; both correctable, so collapse fixes it exactly however small θ\theta is.
General CPTP channelE(ρ)=∑kEkρEk†\mathcal{E}(\rho)=\sum_k E_k\rho E_k^\daggerEach EkE_k expands in Paulis; measurement yields an effective Pauli channel — discretization survives full decoherence.
Leakage out of qubit space∣0⟩→∣2⟩|0\rangle\to|2\rangleFails: a third-level jump is not in span{I,X,Y,Z}\mathrm{span}\{I,X,Y,Z\} and needs a dedicated leakage-reduction unit.
Correlated / high-weight errorE=ZiZjE=Z_iZ_j beyond distanceNeeds the two-qubit Pauli space (16 operators); a distance-3 code correcting weight-1 errors need not span it.
Non-correctable spanning set⟨ψLi∣Ea†Eb∣ψLj⟩≠Cabδij\langle\psi_L^i|E_a^\dagger E_b|\psi_L^j\rangle\ne C_{ab}\delta_{ij}If Knill–Laflamme fails on the basis, no measurement helps: the discrete errors are themselves uncorrectable.

Worked examples

1 · A coherent over-rotation, corrected to zero residual (confidence builder)

A control miscalibration applies an unwanted RX1(θ)R_{X_1}(\theta) with θ=0.1\theta=0.1 rad to qubit 1 of a state encoded in the three-qubit bit-flip code. What is the probability that syndrome extraction flags an error, and how large is the residual coherent rotation after recovery?

Step 1. By Derivation B the state splits into a no-error branch with amplitude cos⁡(θ/2)\cos(\theta/2) and an X1X_1 branch with amplitude −isin⁡(θ/2)-i\sin(\theta/2); the half-angle is θ/2=0.05\theta/2=0.05 rad.

Step 2. The error is flagged (syndrome (−1,+1)(-1,+1)) with probability sin⁡2(θ/2)=sin⁡2(0.05)=(0.049979)2=2.50×10−3\sin^2(\theta/2)=\sin^2(0.05)=(0.049979)^2=\mathbf{2.50\times10^{-3}}, and left uncorrected (syndrome (+1,+1)(+1,+1)) with probability cos⁡2(0.05)=0.99750\cos^2(0.05)=\mathbf{0.99750}. These sum to 11.

Step 3. In the flagged branch recovery applies X1X_1, giving X1(X1∣ψL⟩)=∣ψL⟩X_1(X_1|\psi_L\rangle)=|\psi_L\rangle; in the unflagged branch the state already is ∣ψL⟩|\psi_L\rangle. Either way the output is the exact codeword, and the residual coherent rotation is exactly zero — θ\theta appears nowhere in the final state, only inside the 2.50×10−32.50\times10^{-3} probability.

Step 4. The contrast with an uncorrected qubit is the lesson. Left alone, NN over-rotations compose coherently to angle NθN\theta, with infidelity ∼N2θ2/4\sim N^2\theta^2/4; digitized by per-round measurement, each contributes an independent flip probability θ2/4\theta^2/4, so failure grows only as Nθ2/4N\theta^2/4. Discretization converts a dangerous N2N^2 coherent build-up into a benign linear NN accumulation — exactly what the threshold theorem needs.

2 · Amplitude damping digitized into a Pauli channel (fully worked)

Amplitude damping with rate γ=0.1\gamma=0.1 — a qubit relaxing toward ∣0⟩|0\rangle, the archetype of analog decoherence — has Kraus operators E0=(1001−γ)E_0=\begin{pmatrix}1&0\\0&\sqrt{1-\gamma}\end{pmatrix} and E1=(0γ00)E_1=\begin{pmatrix}0&\sqrt\gamma\\0&0\end{pmatrix}. Expand it in the Pauli basis and find the effective Pauli-error probabilities a single-error-correcting code sees after syndrome measurement.

Step 1. Expand E0E_0 with (1). It is diagonal, so only II and ZZ appear: E0=1+1−γ2 I+1−1−γ2 ZE_0=\tfrac{1+\sqrt{1-\gamma}}{2}\,I+\tfrac{1-\sqrt{1-\gamma}}{2}\,Z. With 0.9=0.948683\sqrt{0.9}=0.948683, the coefficients are c0,I=0.974342c_{0,I}=0.974342 and c0,Z=0.025658c_{0,Z}=0.025658.

Step 2. Expand E1E_1. Since X+iY=(0200)X+iY=\begin{pmatrix}0&2\\0&0\end{pmatrix}, we have E1=γ2(X+iY)E_1=\tfrac{\sqrt\gamma}{2}(X+iY), giving c1,X=γ2=0.158114c_{1,X}=\tfrac{\sqrt\gamma}{2}=0.158114 and c1,Y=0.158114 ic_{1,Y}=0.158114\,i. So amplitude damping spreads over all four Paulis — a point in span{I,X,Y,Z}\mathrm{span}\{I,X,Y,Z\}.

Step 3. By Derivation C, after syndrome measurement the code sees an effective Pauli channel with weights pa=∑k∣cka∣2p_a=\sum_k|c_{ka}|^2: pI=(0.974342)2=0.94934p_I=(0.974342)^2=\mathbf{0.94934}, pZ=(0.025658)2=6.58×10−4p_Z=(0.025658)^2=\mathbf{6.58\times10^{-4}}, and pX=pY=(0.158114)2=0.02500p_X=p_Y=(0.158114)^2=\mathbf{0.02500}.

Step 4. Normalization checks: pI+pX+pY+pZ=0.94934+0.02500+0.02500+0.00066=1.00000p_I+p_X+p_Y+p_Z=0.94934+0.02500+0.02500+0.00066=\mathbf{1.00000}. The continuous, coherence-carrying amplitude-damping channel has become a discrete Pauli channel. A single-error-correcting code corrects each Pauli, hence corrects amplitude damping to leading order in γ\gamma — the four numbers plotted in the second figure.

Physics of the portrait

Discretization made quantum error correction thinkable, and it entered the world inside Peter Shor's 1995 nine-qubit code — a code that corrects an arbitrary single-qubit error by defeating only the discrete Paulis.

Origins · 1 minPeter Shor and how this was discovered
PS
Peter Shor
born 1959 · quantum error correction and the first quantum code

Peter Williston Shor is the American mathematician whose 1994 factoring algorithm proved a quantum computer could break public-key cryptography — and who, a year later, showed such a machine could survive its own noise. The 1994 result drew an immediate objection: qubits decohere continuously, classical redundancy seemed barred by the no-cloning theorem, and the errors seemed too rich to enumerate. Shor's 1995 paper answered both at once. His nine-qubit code — a bit-flip code nested inside a phase-flip code — encoded one logical qubit so that measuring a set of commuting check operators revealed which discrete Pauli error had occurred without disturbing the logical information. The breakthrough underneath was discretization: because any single-qubit error is a combination of I,X,Y,ZI,X,Y,Z, a code correcting just those four corrects every continuous rotation and decoherence process, since measurement collapses the continuous error onto a discrete Pauli branch. Within a year Steane found a smaller seven-qubit code and the CSS family appeared; Knill and Laflamme then distilled the general correctability conditions that make the span argument a theorem, and Gottesman's stabilizer formalism recast the whole apparatus as Pauli measurement. Discretization became the quiet axiom on which the threshold theorem now rests.

History
1994
— Shor's factoring algorithm makes protecting quantum information urgent
1995
— Shor's nine-qubit code corrects an arbitrary single-qubit error
1996
— Steane's seven-qubit code and the Calderbank–Shor–Steane construction
1996
— Bennett, DiVincenzo, Smolin, Wootters give general correction conditions
1997
— Knill and Laflamme prove the conditions for quantum error correction
1997
— Gottesman's stabilizer formalism casts syndromes as Pauli measurements

Going further

Texts. Nielsen & Chuang, Quantum Computation and Quantum Information, §10.2–10.3 (discretization and the general conditions) · Preskill, Lecture Notes on Quantum Computation, Ch. 7 · Gottesman, Stabilizer Codes and Quantum Error Correction (PhD thesis, Caltech, 1997) · Knill & Laflamme, Phys. Rev. A 55, 900 (1997) · Wilde, Quantum Information Theory, Ch. 22 · Lidar & Brun (eds.), Quantum Error Correction (2013). The Pauli expansion (1) and orthonormality (2) are standard results, verifiable by direct matrix multiplication; the rotation identity (5) follows from (n^⋅σ⃗)2=I(\hat n\cdot\vec\sigma)^2=I. The three-qubit syndrome of X1X_1 is (−1,+1)(-1,+1), with flip probability sin⁡2(0.05)=2.50×10−3\sin^2(0.05)=2.50\times10^{-3}; the amplitude-damping weights are pI=0.94934p_I=0.94934, pX=pY=0.02500p_X=p_Y=0.02500, pZ=6.58×10−4p_Z=6.58\times10^{-4}, summing to 11.

Prerequisites. Quantum Error Model and Decoherence Channels · Pauli Group and Clifford Group · Three-Qubit Bit-Flip Code · Knill–Laflamme Conditions. Leads to. Stabilizer Formalism · Syndrome Measurement and Decoding · Shor Nine-Qubit Code · CSS Codes.

Ask the tutor: Ask me to walk through the syndrome collapse for a specific channel — depolarizing, dephasing, or your own Kraus operators — and watch the continuous error digitize into Pauli probabilities.

Does discretization mean the physical error was secretly discrete all along?

No — this is the most common misreading. Before measurement a small rotation genuinely places the codeword in a continuous superposition of error branches; the amplitude −isin⁡(θ/2)-i\sin(\theta/2) is a real, continuous quantity physically present in the state, and the qubit really has drifted. What is true is that measurement makes it digital: projecting onto a syndrome subspace collapses the superposition onto one Pauli branch, and only then is the error a discrete event with a definite probability. Discretization describes what measurement does to a continuous error, not a claim that it was secretly discrete beforehand.

If a rotation by θ = 0.001 is tiny, why bother correcting it at all?

Because tiny coherent errors are insidious when left uncorrected. A single over-rotation by small θ\theta costs almost nothing, but NN of them add coherently, composing to angle NθN\theta and an infidelity growing like N2θ2/4N^2\theta^2/4 — the quadratic build-up that dooms a long computation. Encoding and measuring each round digitizes every rotation into an independent Pauli flip of probability ≈θ2/4\approx\theta^2/4, so the damage grows only linearly, like Nθ2/4N\theta^2/4. Correcting the tiny error is not about the one rotation; it breaks the coherent accumulation across the circuit.

Why do the coherences between different Pauli branches vanish?

Because distinct correctable Paulis carry distinct syndromes and live in orthogonal stabilizer eigenspaces. A channel produces cross-terms PaρLPb†P_a\rho_L P_b^\dagger with a≠ba\ne b — coherences between error branches — but the projector MsM_s keeps only the branch with syndrome ss: for a nondegenerate code just the a=ba=b term survives, and off-diagonal terms are annihilated because PaP_a and PbP_b sit in different syndrome subspaces. Measurement destroys the interference, leaving a diagonal (incoherent) Pauli channel — the microscopic event that turns analog decoherence into a discrete Pauli mixture.

Does discretization work for full decoherence, or only for coherent rotations?

It works for any completely positive trace-preserving channel — the general model of decoherence. Each Kraus operator expands in the Pauli basis by (1), so the whole channel lives in the Pauli span; syndrome measurement kills the branch coherences and leaves an effective Pauli channel with probabilities pa=∑k∣cka∣2p_a=\sum_k|c_{ka}|^2, as Worked Example 2 shows for amplitude damping. So a code correcting the four Paulis corrects amplitude damping, dephasing, depolarizing, and every other single-qubit channel — coherent rotations are just the simplest case.

What kind of error does discretization fail to handle?

Two kinds. First, anything outside the correctable span: leakage out of the computational subspace — a transmon excited to ∣2⟩|2\rangle — is not a 2×22\times2 Pauli combination, slips past every syndrome, and needs a dedicated leakage-reduction unit. Second, correlated or high-weight errors beyond the code's distance: a distance-3 code spanning only weight-1 Paulis does not span ZiZjZ_iZ_j. Discretization is a theorem about a correctable span; errors that leave it are genuinely uncorrected.

Is measurement really necessary, or could recovery digitize the error by itself?

Measurement — or something physically equivalent — is essential, because digitization is the collapse of a superposition, which unitary recovery alone cannot do: a linear recovery fed a continuous superposition of error branches simply outputs a continuous superposition, so the analog character survives. What breaks the linearity is projective syndrome measurement (or coupling the syndrome to an ancilla that is then measured). Autonomous error correction swaps explicit measurement for engineered dissipation, but that dissipation plays the same role: it irreversibly funnels the continuous error onto discrete pointer states. Either way, an entropy-removing, non-unitary step is what turns the analog error digital.

Before this

Where this leads