4. 量子比特错误模型建模

在量子纠错循环仿真中，错误模型决定了“错误从哪里来、以多大概率发生、以什么形式作用在量子比特上”。如果错误模型建得过于理想，纠错码的性能会被高估；如果建得过于悲观，又会掩盖编码方案本身的优势。因此，本节的核心任务不是简单调用一个噪声接口，而是把物理噪声转写为可参数化、可复现、可组合的仿真算子，并观察它们如何改变纠错循环的综合征、逻辑错误率与阈值行为。

一、从物理噪声到量子信道

仿真中通常用密度矩阵 ρ 描述量子比特状态，用完全正定保迹映射（CPTP）描述噪声过程。对单比特，最常用的参数化形式是 Kraus 算子和表示：

ρ' = Σ_k E_k ρ E_k†，  约束 Σ_k E_k† E_k = I

其中 E_k 是 Kraus 算子，k 是噪声分支。去极化、振幅阻尼、相位阻尼都可以写成这种形式。仿真实现时，我们并不需要每次显式构造完整密度矩阵，而可以在稳定子仿真或蒙特卡洛轨迹仿真中，按概率抽取一个 Kraus 分支作用到当前态上。这样既保留物理意义，又能控制计算量。

二、五类典型错误模型的参数化

1. 去极化噪声（Depolarizing Noise）。它表示量子比特以概率 p 被随机地施加 X、Y、Z 之一，以概率 1−p 保持不变。单比特去极化信道可写为：

E0 = sqrt(1 - 3p/4) I
E1 = sqrt(p/4) X
E2 = sqrt(p/4) Y
E3 = sqrt(p/4) Z

这里 p 是总错误概率，三个非恒等分支各占 p/4。若仿真中采用“以概率 p 随机选一个 Pauli 错误”的简化模型，则每个 Pauli 的概率为 p/3，这与标准去极化信道在 p 较小时近似一致，但严格形式应使用上述系数。去极化噪声常用于评估纠错码的阈值，因为它对各方向错误对称。

2. 振幅阻尼（Amplitude Damping）。它描述能量从 |1⟩ 衰减到 |0⟩ 的过程，典型来自 T1 弛豫。设衰减概率为 γ，Kraus 算子为：

E0 = [[1, 0], [0, sqrt(1 - γ)]]
E1 = [[0, sqrt(γ)], [0, 0]]

γ 与时间 t 和 T1 的关系为 γ = 1 − exp(−t/T1)。在纠错循环中，每个量子比特在每个周期都会经历一次振幅阻尼，因此 γ 应取“单周期内”的衰减量，而不是整个实验时长。振幅阻尼不是对称错误，它会把 |1⟩ 推向 |0⟩，因此对某些编码会引入偏向性错误。

3. 相位阻尼（Phase Damping）。它描述相干性丢失但能量不衰减的过程，典型来自 T2 退相位。设退相位概率为 λ，Kraus 算子为：

E0 = [[1, 0], [0, sqrt(1 - λ)]]
E1 = [[0, 0], [0, sqrt(λ)]]

λ 与 T2 的关系为 λ = 1 − exp(−t/T2)。相位阻尼只影响密度矩阵的非对角元，对 Z 基测量不直接产生比特翻转，但会降低叠加态的可见度，从而影响需要相干操作的纠错电路。

4. 测量错误（Measurement Error）。它描述测量结果被翻转的概率。若真实测量结果为 m，仿真中输出 m⊕1 的概率为 p_m。测量错误通常建模为经典二进制对称信道，作用在综合征比特上。它不会改变量子态本身，但会让解码器收到错误的综合征，从而施加错误的纠正操作。在多轮纠错中，测量错误会沿时间方向传播，因此需要重复测量和时空解码。

5. 门错误（Gate Error）。它描述量子门操作不完美导致的错误。常见建模方式是在每个门后附加一个去极化信道，错误率为 p_g；也可以按门类型分别设置单比特门错误率 p_1q 和双比特门错误率 p_2q。对于 CNOT 门，还可使用“双比特去极化”模型，即 15 个非恒等 Pauli 组合各以 p/15 的概率出现。门错误是纠错循环中最主要的错误来源之一，因为它直接作用于编码电路和综合征提取电路。

三、在仿真中的实现方式

实现错误模型时，推荐采用“按事件注入”的方式，而不是在每个时间步对整个密度矩阵做一次完整信道作用。具体做法是：先根据当前操作类型和持续时间，计算该操作对应的错误概率；然后在蒙特卡洛轨迹中，以该概率抽取是否发生错误，以及发生哪种错误；最后把抽到的 Pauli 或 Kraus 分支作用到当前稳定子态或态矢量上。

以单周期纠错循环为例，假设一个数据量子比特经历以下事件：

1. 空闲等待 t_idle = 1 μs，T1 = 50 μs，T2 = 70 μs
   γ = 1 - exp(-1/50) ≈ 0.0198
   λ = 1 - exp(-1/70) ≈ 0.0142

2. 单比特门 1 次，p_1q = 0.001
3. CNOT 门 2 次，p_2q = 0.01
4. 测量 1 次，p_m = 0.005

在轨迹仿真中，每个事件独立抽样。例如空闲等待阶段，先以 γ 概率施加振幅阻尼分支，再以 λ 概率施加相位阻尼分支；门操作后以 p_1q 或 p_2q 概率施加去极化错误；测量时以 p_m 概率翻转读出结果。这样得到的错误是“事件级”的，能更真实地反映纠错循环中错误发生的位置和顺序。

如果使用稳定子仿真，Pauli 错误可以直接用二进制辛表示更新；振幅阻尼和相位阻尼则需要转换为等效的 Pauli 近似，或使用“Pauli 通道”近似方法。常见做法是把振幅阻尼近似为：以概率 γ 施加 X 错误，以概率 γ/2 施加 Z 错误，其余保持不变。这种近似在 γ 较小时误差可控，但会丢失振幅阻尼的偏向性，因此在高精度仿真中应谨慎使用。

四、对纠错循环的影响分析

不同错误模型对纠错循环的影响差异很大。去极化噪声是对称的，它同时产生 X、Y、Z 错误，因此对纠错码的阈值要求最高；振幅阻尼是偏向性的，它主要产生 |1⟩→|0⟩ 的衰减，在表面码中会表现为特定方向的错误，可能被解码器部分吸收；相位阻尼不直接翻转比特，但会降低相干性，使需要叠加的综合征提取电路更容易出错；测量错误会直接污染综合征，导致解码器误判，通常需要多轮测量和时空解码来抑制；门错误则是最直接、最普遍的来源，它作用在编码和提取电路上，往往决定整个循环的错误率下限。

在仿真中观察这些影响时，建议固定其他参数，只扫描一个参数。例如固定 p_1q = p_2q = 0.001，只改变 p_m，从 0 到 0.02 扫描，记录逻辑错误率 p_L 随 p_m 的变化。若 p_L 随 p_m 快速上升，说明当前解码方案对测量错误敏感，需要增加测量轮数或改进解码权重。类似地，固定 p_m，只改变 γ，可以观察 T1 对逻辑错误率的贡献。

最后，错误模型的参数化应尽量与实验可测数据对齐。T1、T2、单比特门错误率、双比特门错误率、读出错误率都是量子硬件中常规标定的指标。把这些指标直接映射为仿真参数，可以让纠错循环的性能评估更接近真实设备，也更容易解释“为什么某个码在当前硬件上表现更好”。