📚 课程目录第3章噪声的数学表示

3. 噪声的数学表示

在量子纠错编码的语境里，“噪声”不能只停留在“比特翻转”“相位翻转”这样的口头描述上。编码方案能否被证明有效、阈值能否被计算、解码器能否被设计，都依赖于噪声有一套严格的数学表示。本节沿着“状态表示—演化表示—算符分解—物理约束”这条线索，把噪声从直观图像推进到可计算的映射形式，并给出一个可以手算验证的小例子。

一、密度矩阵：把“不确定”写进状态里

纯态用态矢量 |ψ⟩ 描述，但噪声场景下我们经常不知道系统究竟处于哪个纯态，只知道它以概率 p_i 处于 |ψ_i⟩。这时用密度矩阵表示：

ρ = Σ_i p_i |ψ_i⟩⟨ψ_i|,   p_i ≥ 0,  Σ_i p_i = 1

密度矩阵必须满足三条性质：厄米性 ρ = ρ^†；半正定性 ρ ≥ 0；迹为 1，Tr(ρ) = 1。纯态满足 Tr(ρ²) = 1，混合态满足 Tr(ρ²) < 1，因此 Tr(ρ²) 可以作为“纯度”的度量。对单量子比特，ρ 可写成 Bloch 形式 ρ = (I + r_xX + r_yY + r_zZ)/2，其中 Bloch 矢量 r 的长度 |r| ≤ 1，|r| = 1 对应纯态，|r| = 0 对应完全混合态 I/2。这个几何图像在后续理解退相干时非常有用：噪声往往就是把 Bloch 矢量“缩短”或“旋转”的操作。

多比特情形下，密度矩阵的维度随比特数指数增长，n 个量子比特对应 2ⁿ × 2ⁿ 的矩阵。这正是经典模拟量子系统困难的根本原因，也是纠错编码必须依赖结构化假设（如局域噪声）的原因。

二、量子信道：状态的合法演化

封闭系统的演化是酉变换 ρ → UρU^†。但噪声意味着系统与环境发生了不可控的相互作用，系统本身不再封闭，演化变成更一般的映射 ℰ：ρ → ℰ(ρ)。这个映射被称为量子信道（quantum channel）。

并非任意映射都能描述物理过程。一个映射要成为合法的量子信道，必须满足两个条件：完全正性（completely positive，CP）与保迹性（trace-preserving，TP），合称 CPTP 映射。保迹性保证概率守恒：Tr(ℰ(ρ)) = 1。完全正性的含义更微妙：不仅 ℰ 本身要把正定矩阵映为正定矩阵，而且当我们把 ℰ 作用在更大系统的一个子系统上、即 ℰ ⊗ I_k 时，对任意辅助维度 k 都必须保持正性。只要求“正性”而不要求“完全正性”是不够的，存在正但非完全正的映射，它作用在纠缠态上会产生负本征值，从而给出非物理的“负概率”。

记忆要点：酉演化是 CPTP 映射的特例；CPTP 是“物理可实现”的充要刻画，任何真实噪声过程都落在这一集合内。

三、Kraus 算符：CPTP 映射的显式分解

Kraus 定理给出了 CPTP 映射的等价刻画：映射 ℰ 是完全正且保迹的，当且仅当存在一组算符 {K_i}，使得

ℰ(ρ) = Σ_i K_i ρ K_i†,   且   Σ_i K_i† K_i = I

这组 K_i 称为 Kraus 算符，第二个等式称为完备性关系，它正是保迹性的直接体现。Kraus 表示不唯一：若 {K_i} 是一组解，则对任意酉矩阵 u，{K′_j = Σ_i u_jiK_i} 给出同一信道。这种不唯一性对应“环境基矢的选择自由”，在物理上无关紧要，但在数值实现时需要固定一种规范。

Kraus 算符的个数最多为 d²（d 为系统维度），但实际噪声通常只需少数几个算符。下面用三个典型信道说明如何从物理图像写出 Kraus 算符。

信道	物理图像	Kraus 算符	对 Bloch 矢量的作用
比特翻转	以概率 p 施加 X	K₀=√(1−p)I，K₁=√p X	r_x 不变，r_y、r_z 乘 (1−2p)
相位翻转	以概率 p 施加 Z	K₀=√(1−p)I，K₁=√p Z	r_z 不变，r_x、r_y 乘 (1−2p)
退极化	以概率 p 完全随机化	K₀=√(1−3p/4)I，K_1..3=√(p/4){X,Y,Z}	整个 r 乘 (1−p)

验证退极化信道的完备性关系：K₀^†K₀ + Σ_j=1..3 K_j^†K_j = (1−3p/4)I + 3·(p/4)I = I，成立。注意参数化方式：这里 p 是“发生非恒等错误的概率”，当 p = 1 时信道把任意态映为 I/2，即完全混合态，Bloch 矢量长度归零。

四、一个可手算的示例：比特翻转信道作用于 |0⟩

设 p = 0.1，输入 ρ = |0⟩⟨0| = diag(1, 0)。按 Kraus 公式计算：

K_0 ρ K_0† = (1−p) · diag(1, 0) = diag(0.9, 0)
K_1 ρ K_1† = p · X|0⟩⟨0|X = p · diag(0, 1) = diag(0, 0.1)
ℰ(ρ) = diag(0.9, 0.1)

结果是一个对角混合态，纯度 Tr(ρ²) = 0.81 + 0.01 = 0.82 < 1，符合混合态预期。若把同一信道作用在 |+⟩ = (|0⟩+|1⟩)/√2 上，由于 X|+⟩ = |+⟩，信道对它是“透明”的，输出仍为纯态 |+⟩⟨+|。这说明比特翻转噪声在 X 本征基下不产生错误——这正是纠错编码选择测量基、设计稳定子的物理依据：噪声的作用依赖于基，编码就是把信息放进对目标噪声“不敏感”的子空间。

五、从抽象映射到可计算表示

在数值实现中，Kraus 表示可以直接用于模拟：给定 ρ 和算符列表，逐项累加即可。但若要做阈值估计或解码优化，常把信道写成其他等价形式。例如 Choi 矩阵 J(ℰ) = Σ_ij |i⟩⟨j| ⊗ ℰ(|i⟩⟨j|)，它把信道变成一个作用在双倍维度上的密度矩阵，CP 性等价于 J(ℰ) ≥ 0，TP 性等价于对第二个子系统求偏迹得到 I。另一种常用形式是 Pauli 信道（Pauli channel）：任何单比特信道都可展开为 ℰ(ρ) = Σ_{P∈{I,X,Y,Z}} p_P PρP，其中 p_P ≥ 0 且和为 1。Pauli 信道在稳定子码分析中极其方便，因为 Pauli 算符在 Clifford 操作下结构稳定，错误传播可以用“错误模式”的符号演算追踪。

常见误区：把“正性”当成“完全正性”。仅验证 ℰ(ρ) ≥ 0 对单个输入成立，并不能保证 ℰ ⊗ I 在纠缠输入上仍为正。工程上若用自定义算符构造信道，务必检查完备性关系 Σ K_i^†K_i = I，这是最直接的 TP 检验；CP 性则由“存在 Kraus 分解”自动保证。

下图把本节的核心结构串起来：从物理噪声出发，经 CPTP 约束落到 Kraus 分解，再展开为三种可计算表示，并标出各自的检验条件。

把上述结构落到工程实践，通常按以下顺序推进：先根据器件物理确定噪声类型（例如超导比特以退极化与相位退相干为主，离子阱以退相干为主），写出对应的 Kraus 算符或 Pauli 信道概率；再用完备性关系做一次数值校验，确认 Σ K_i^†K_i 与单位矩阵的最大偏差在 10⁻¹² 量级；随后用 Choi 矩阵检查半正定性，若出现负本征值则说明算符构造有误或参数越界；最后把信道接入编码模拟器，观察逻辑错误率随物理错误率 p 的变化曲线，曲线的交叉点即阈值估计的起点。这一流程把“噪声是什么”转化为“噪声如何被检验、被计算、被比较”，是后续稳定子码与阈值分析的必要前提。

需要强调的是，Kraus 表示与 Pauli 信道并非互相替代，而是同一信道的不同“视角”：前者适合直接模拟与严格证明，后者适合稳定子码的符号化错误追踪。理解两者之间的转换关系（例如把 Kraus 算符按 Pauli 基展开并合并同类项），是阅读纠错文献时的一项基本技能。