Skip to content

Pauli Channels

A Pauli channel is a stochastic mixture of conjugations by Pauli operators. For one qubit,

Φ(ρ)=p0ρ+p1XρX+p2YρY+p3ZρZ,pμ≥0,∑μ=03pμ=1.\Phi(\rho) = p_0\rho + p_1X\rho X + p_2Y\rho Y + p_3Z\rho Z, \qquad p_\mu\ge0, \quad \sum_{\mu=0}^3p_\mu=1.

Equivalently,

Φ(ρ)=∑μ=03pμσμρσμ,σ0=I,σ1=X,σ2=Y,σ3=Z.\Phi(\rho) = \sum_{\mu=0}^3 p_\mu\sigma_\mu\rho\sigma_\mu, \qquad \sigma_0=I,\quad \sigma_1=X,\quad \sigma_2=Y,\quad \sigma_3=Z.

Pauli channels are useful because they are completely positive, trace preserving, unital, diagonal in the Pauli operator basis, and compatible with stabilizer and error-correction calculations. They are also limited: energy relaxation, leakage, coherent over-rotations, and most correlated environment dynamics are not Pauli channels unless an additional approximation has been made.

For the Pauli algebra itself, see Pauli Matrices. For where this model sits among other noise channels, see Common Noise Channels.

The Kraus operators may be chosen as

Kμ=pμ σμ,μ=0,1,2,3.K_\mu = \sqrt{p_\mu}\,\sigma_\mu, \qquad \mu=0,1,2,3.

The trace-preserving condition is immediate:

∑μ=03Kμ†Kμ=∑μ=03pμI=I.\sum_{\mu=0}^3 K_\mu^\dagger K_\mu = \sum_{\mu=0}^3 p_\mu I = I.

The channel is also unital:

Φ(I)=∑μ=03pμσμIσμ=I.\Phi(I) = \sum_{\mu=0}^3 p_\mu\sigma_\mu I\sigma_\mu = I.

Thus a Pauli channel is a random unitary channel. It can be represented by a classical environment record:

V∣ψ⟩=∑μ=03pμ σμ∣ψ⟩⊗∣μ⟩E.V\lvert\psi\rangle = \sum_{\mu=0}^3 \sqrt{p_\mu}\, \sigma_\mu\lvert\psi\rangle \otimes \lvert\mu\rangle_E.

If the record μ\mu is known, the receiver can in principle apply σμ\sigma_\mu again to undo the error. If the record is ignored, the reduced system state undergoes the Pauli channel. This is a simple example of a Stinespring representation.

Write a qubit density operator as

ρ=12(I+rxX+ryY+rzZ).\rho = \frac12 \left( I+r_xX+r_yY+r_zZ \right).

Conjugation by XX preserves the XX component and flips the YY and ZZ components:

XρX=12(I+rxX−ryY−rzZ).X\rho X = \frac12 \left( I+r_xX-r_yY-r_zZ \right).

Similarly, YY preserves ryr_y and flips rx,rzr_x,r_z, while ZZ preserves rzr_z and flips rx,ryr_x,r_y. Therefore a Pauli channel acts diagonally on the Bloch vector:

(rx,ry,rz)⟼(λxrx,λyry,λzrz),(r_x,r_y,r_z) \longmapsto (\lambda_xr_x,\lambda_yr_y,\lambda_zr_z),

where

λx=p0+p1−p2−p3,λy=p0−p1+p2−p3,λz=p0−p1−p2+p3.\begin{aligned} \lambda_x&=p_0+p_1-p_2-p_3,\\ \lambda_y&=p_0-p_1+p_2-p_3,\\ \lambda_z&=p_0-p_1-p_2+p_3. \end{aligned}

The origin of the Bloch ball is fixed because the channel is unital. The three axes may contract by different factors and may also acquire sign flips when an eigenvalue is negative.

This diagonal form is one reason Pauli channels are analytically convenient: they turn many channel calculations into three scalar eigenvalue calculations.

Conversely, a unital qubit channel that is diagonal in the Pauli basis has the form

(rx,ry,rz)↦(λxrx,λyry,λzrz).(r_x,r_y,r_z) \mapsto (\lambda_xr_x,\lambda_yr_y,\lambda_zr_z).

It is a Pauli channel exactly when the probabilities obtained by inversion are nonnegative:

p0=1+λx+λy+λz4,p1=1+λx−λy−λz4,p2=1−λx+λy−λz4,p3=1−λx−λy+λz4.\begin{aligned} p_0&=\frac{1+\lambda_x+\lambda_y+\lambda_z}{4},\\ p_1&=\frac{1+\lambda_x-\lambda_y-\lambda_z}{4},\\ p_2&=\frac{1-\lambda_x+\lambda_y-\lambda_z}{4},\\ p_3&=\frac{1-\lambda_x-\lambda_y+\lambda_z}{4}. \end{aligned}

Complete positivity is therefore equivalent to the four inequalities

1+λx+λy+λz≥0,1+λx−λy−λz≥0,1−λx+λy−λz≥0,1−λx−λy+λz≥0.\begin{array}{rcl} 1+\lambda_x+\lambda_y+\lambda_z&\ge&0,\\ 1+\lambda_x-\lambda_y-\lambda_z&\ge&0,\\ 1-\lambda_x+\lambda_y-\lambda_z&\ge&0,\\ 1-\lambda_x-\lambda_y+\lambda_z&\ge&0. \end{array}

These inequalities form a tetrahedron inside the cube [−1,1]3[-1,1]^3. Positivity of the map on isolated qubit states only gives the cube; complete positivity gives the smaller tetrahedron. This is a useful concrete example of the distinction explained in Completely Positive Maps.

Using the output-input convention of Choi Matrix, a Pauli channel has

JΦ=∑μ=03pμ∣σμ⟩⟩⟨⟨σμ∣,J_\Phi = \sum_{\mu=0}^3 p_\mu \lvert\sigma_\mu\rangle\rangle \langle\langle\sigma_\mu\rvert,

where ∣A⟩⟩\lvert A\rangle\rangle is the vectorized operator in the chosen input basis.

For one qubit, the four normalized vectors

12∣σμ⟩⟩\frac{1}{\sqrt2} \lvert\sigma_\mu\rangle\rangle

are Bell states up to phases. Thus the normalized Choi state is Bell diagonal, with eigenvalues pμp_\mu:

τΦ=12JΦ,spec⁡(τΦ)={p0,p1,p2,p3}.\tau_\Phi = \frac12J_\Phi, \qquad \operatorname{spec}(\tau_\Phi) = \{p_0,p_1,p_2,p_3\}.

The Choi matrix makes the complete-positivity condition transparent: the probabilities are the Choi-state eigenvalues and must be nonnegative.

The ZZ phase-flip channel is

Φp(ρ)=(1−p)ρ+pZρZ.\Phi_p(\rho) = (1-p)\rho+pZ\rho Z.

It is a Pauli channel with

p0=1−p,p3=p,p1=p2=0.p_0=1-p,\qquad p_3=p,\qquad p_1=p_2=0.

Its Bloch eigenvalues are

λx=λy=1−2p,λz=1.\lambda_x=\lambda_y=1-2p, \qquad \lambda_z=1.

It preserves populations in the ZZ basis and damps transverse coherences. This is the Pauli-channel form of Dephasing Channel.

The XX bit-flip channel

Φp(ρ)=(1−p)ρ+pXρX\Phi_p(\rho) = (1-p)\rho+pX\rho X

has

λx=1,λy=λz=1−2p.\lambda_x=1, \qquad \lambda_y=\lambda_z=1-2p.

The YY bit-phase-flip channel is analogous, preserving the YY component and damping the XX and ZZ components.

These names are basis dependent. They are useful in qubit error models, but they should not be mistaken for basis-independent physical mechanisms.

The qubit Pauli-error convention for depolarizing noise is

Pp(ρ)=(1−p)ρ+p3(XρX+YρY+ZρZ).\mathcal P_p(\rho) = (1-p)\rho + \frac{p}{3} \left( X\rho X+Y\rho Y+Z\rho Z \right).

Thus

p0=1−p,p1=p2=p3=p3.p_0=1-p, \qquad p_1=p_2=p_3=\frac{p}{3}.

All Bloch eigenvalues are equal:

λx=λy=λz=1−4p3.\lambda_x=\lambda_y=\lambda_z = 1-\frac{4p}{3}.

The replacement-probability convention instead writes

Dq(ρ)=(1−q)ρ+qI2.\mathcal D_q(\rho) = (1-q)\rho+q\frac{I}{2}.

The conversion is

q=4p3.q=\frac{4p}{3}.

This convention trap is common: in the Pauli-error convention, p=3/4p=3/4 gives the completely mixed output, while p=1p=1 gives the valid but negative-shrink channel with λ=−1/3\lambda=-1/3. See Depolarizing Channel for the full parameter map.

For nn qubits, use Pauli strings

Pa=Pa1⊗⋯⊗Pan,Paj∈{I,X,Y,Z},P_{\mathbf a} = P_{a_1}\otimes\cdots\otimes P_{a_n}, \qquad P_{a_j}\in\{I,X,Y,Z\},

ignoring overall phases. A multiqubit Pauli channel is

Φ(ρ)=∑apaPaρPa,pa≥0,∑apa=1.\Phi(\rho) = \sum_{\mathbf a} p_{\mathbf a} P_{\mathbf a}\rho P_{\mathbf a}, \qquad p_{\mathbf a}\ge0, \quad \sum_{\mathbf a}p_{\mathbf a}=1.

An independent one-qubit Pauli noise model is the special case where the distribution factors:

pa=∏j=1npaj(j).p_{\mathbf a} = \prod_{j=1}^n p_{a_j}^{(j)}.

Correlated Pauli noise allows arbitrary probabilities pap_{\mathbf a}. For example, the event Z1Z2Z_1Z_2 can occur with a probability not determined by the marginal Z1Z_1 and Z2Z_2 rates.

This distinction matters for error correction. A code that corrects independent single-qubit errors can behave very differently under spatially correlated Pauli strings. The full stabilizer-code theory belongs to quantum information; the open-system point is that “Pauli channel” does not automatically mean “independent errors.” For the algebraic stabilizer background, see Stabilizer States Preview and Stabilizer Identities.

Pauli Noise and Depolarizing Channels owns the QI-facing probability-to-transfer, Clifford-propagation, syndrome/logical-pushforward, and approximation workflow; this page retains the channel definition, complete-positivity conditions, Bloch eigenvalues, Choi form, and formal Pauli-twirl construction.

Pauli twirling averages a channel over Pauli conjugations. For nn qubits, schematically,

T(Φ)(ρ)=14n∑QQ†Φ(QρQ†)Q,\mathcal T(\Phi)(\rho) = \frac{1}{4^n} \sum_Q Q^\dagger \Phi(Q\rho Q^\dagger) Q,

where the sum is over phase-free nn-qubit Pauli strings. The twirled channel is diagonal in the Pauli operator basis and, when completely positive and trace preserving, can be interpreted as a Pauli channel.

Twirling is useful because it removes coherent and off-diagonal Pauli-transfer structure. It is common in benchmarking, randomized protocols, and analytic threshold estimates.

But twirling is an approximation or protocol transformation, not a statement that the original microscopic noise was Pauli noise. It can hide:

  • coherent over-rotations,
  • nonunital relaxation,
  • leakage out of the intended Hilbert space,
  • correlations in time,
  • correlations across subsystems,
  • error mechanisms that affect error correction differently despite similar average figures of merit.

Use a Pauli approximation only after stating what has been averaged and which quantities the approximation is intended to preserve.

Pauli channels are unital:

Φ(I)=I.\Phi(I)=I.

Therefore they cannot model relaxation toward a pure ground state. The Amplitude-Damping Channel is the standard counterexample: it is trace preserving and completely positive, but for γ>0\gamma>0 it sends I/2I/2 toward ∣0⟩⟨0∣\lvert0\rangle\langle0\rvert.

Pauli channels also preserve the system Hilbert space. They do not describe flagged erasure unless the output space includes a loss flag, and they do not describe leakage unless the Hilbert space is enlarged.

Finally, a coherent unitary error such as

ρ↦e−iϵZ/2ρeiϵZ/2\rho \mapsto e^{-i\epsilon Z/2} \rho e^{i\epsilon Z/2}

is not generally the same as stochastic ZZ flips. Coherent errors can accumulate systematically, while stochastic Pauli errors compose probabilistically.

  • Calling any qubit noise model “Pauli noise” because Pauli matrices are a convenient basis.
  • Comparing depolarizing parameters without checking whether pp means nonidentity Pauli-error probability or qq means replacement probability.
  • Treating a twirled Pauli channel as microscopically equivalent to the untwirled channel.
  • Forgetting that Pauli channels are unital and therefore cannot represent zero-temperature relaxation.
  • Assuming multiqubit Pauli noise is independent; correlated Pauli-string distributions are still Pauli channels.
  • Reading a Pauli error label as an observed event without specifying a record or syndrome extraction process.
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press (2010).
  • J. Preskill, Lecture Notes on Quantum Computation, California Institute of Technology.
  • D. Gottesman, “Stabilizer codes and quantum error correction,” PhD thesis, California Institute of Technology (1997).
  • J. Watrous, The Theory of Quantum Information, Cambridge University Press (2018).
  • M. M. Wilde, Quantum Information Theory, Cambridge University Press, 2nd ed. (2017).
  1. Bloch eigenvalues. Derive the formulas for λx,λy,λz\lambda_x,\lambda_y,\lambda_z in terms of p0,p1,p2,p3p_0,p_1,p_2,p_3.
Solution

Use

XXX=X,XYX=−Y,XZX=−Z,X X X=X,\qquad X Y X=-Y,\qquad X Z X=-Z,

and the analogous cyclic identities. The XX component of the output receives sign +1+1 from the II and XX branches and sign −1-1 from the YY and ZZ branches:

λx=p0+p1−p2−p3.\lambda_x=p_0+p_1-p_2-p_3.

Similarly,

λy=p0−p1+p2−p3,λz=p0−p1−p2+p3.\lambda_y=p_0-p_1+p_2-p_3, \qquad \lambda_z=p_0-p_1-p_2+p_3.
  1. Invert the probabilities. Starting from the three λi\lambda_i equations and ∑μpμ=1\sum_\mu p_\mu=1, derive the inverse formulas for pμp_\mu.
Solution

Add the normalization equation to all three eigenvalue equations:

1+λx+λy+λz=4p0.1+\lambda_x+\lambda_y+\lambda_z = 4p_0.

Changing signs selects the other probabilities:

1+λx−λy−λz=4p1,1−λx+λy−λz=4p2,1−λx−λy+λz=4p3.\begin{aligned} 1+\lambda_x-\lambda_y-\lambda_z&=4p_1,\\ 1-\lambda_x+\lambda_y-\lambda_z&=4p_2,\\ 1-\lambda_x-\lambda_y+\lambda_z&=4p_3. \end{aligned}

Dividing by 44 gives the stated inverse formulas. Complete positivity is then pμ≥0p_\mu\ge0 for all μ\mu.

  1. Dephasing as a Pauli channel. Show that the phase-flip channel (1−p)ρ+pZρZ(1-p)\rho+pZ\rho Z maps ρ01\rho_{01} to (1−2p)ρ01(1-2p)\rho_{01}.
Solution

For

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

one has

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

Therefore

ρ01′=(1−p)ρ01−pρ01=(1−2p)ρ01.\rho_{01}' = (1-p)\rho_{01} -p\rho_{01} = (1-2p)\rho_{01}.
  1. Why amplitude damping is not Pauli. Use unitality to show that the amplitude-damping channel with γ>0\gamma>0 cannot be a Pauli channel.
Solution

Every Pauli channel is unital:

Φ(I)=I.\Phi(I)=I.

For amplitude damping with 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},

one finds

K0K0†+K1K1†=(1+γ001−γ).K_0K_0^\dagger+K_1K_1^\dagger = \begin{pmatrix} 1+\gamma&0\\ 0&1-\gamma \end{pmatrix}.

This equals II only when γ=0\gamma=0. Hence nontrivial amplitude damping is not a Pauli channel.