4. 逻辑错误率与码距关系
在容错量子计算中,码距 d 是决定逻辑错误率的核心参数。物理错误率 p 一旦低于某个阈值 p_th,逻辑错误率 P_L 会随码距增大而呈指数下降,这一现象是量子纠错能够“用更多物理比特换取更可靠逻辑比特”的根本依据。本节围绕 P_L(p,d) 的定量关系展开,重点讲清指数抑制的标度形式、阈值行为的判据、有限尺寸效应带来的偏差,以及工程上常用的经验拟合与解析近似方法。
一、指数抑制的基本标度
对于距离为 d 的量子纠错码,当物理错误率 p 远低于阈值时,逻辑错误率通常可写成如下标度形式:
P_L(p, d) ≈ A · (p / p_th)^((d+1)/2)
其中 A 是与码型、解码器和噪声模型相关的常数,指数 (d+1)/2 来自“最少需要 (d+1)/2 个物理错误才能造成一个不可纠正的逻辑错误”这一组合事实。以表面码为例,d=3 时指数为 2,d=5 时指数为 3,d=7 时指数为 4。若 p/p_th = 0.1,则 d=3 的逻辑错误率约为 A×10^-2,d=5 约为 A×10^-3,d=7 约为 A×10^-4。每增加 2 个码距,逻辑错误率大约下降一个数量级。这就是“指数抑制”的直观含义:码距每增加 2,错误率乘以一个小于 1 的因子 p/p_th。
需要强调的是,该标度只在 p < p_th 时成立。当 p > p_th 时,p/p_th > 1,增大码距反而使逻辑错误率上升,纠错码失去保护作用。因此阈值行为是判断“能否通过增大码距获益”的分水岭。
二、阈值行为与有限尺寸效应
阈值 p_th 的定义是:当 p = p_th 时,不同码距的逻辑错误率曲线相交于同一点,即 P_L 与 d 无关。在 p < p_th 区域,曲线随 d 增大而下降;在 p > p_th 区域,曲线随 d 增大而上升。实际数值实验中,由于解码器不完美、边界效应和有限采样,曲线不会严格交于一点,而是形成一个交叉带,交叉带中心可视为有效阈值。
有限尺寸效应主要体现在三个方面。第一,当 d 较小时,(d+1)/2 指数较小,指数抑制尚未充分显现,P_L 对 d 的敏感度低于渐近预期。第二,当 p 接近 p_th 时,高阶修正项不可忽略,单纯幂律拟合会高估或低估 P_L。第三,解码器在有限码距下可能无法达到最优解码性能,导致实测 P_L 高于理论下界。因此,在评估资源开销时,不能只用渐近公式外推,必须结合目标码距范围内的实际拟合。
三、经验拟合方法
在工程评估中,常用经验拟合从少量码距数据外推更大码距的 P_L。最常用的形式是:
log P_L = log A + ((d+1)/2) · log(p / p_th)
对 log P_L 与 d 做线性回归,斜率即为 (1/2)·log(p/p_th),截距为 log A。实际操作中,先固定 p,取 d=3,5,7,9 的 P_L 数据,做线性拟合得到斜率 k,再由 k = (1/2)·log(p/p_th) 反解 p_th。若拟合残差较大,说明有限尺寸效应显著,应改用带修正项的拟合式:
P_L(p, d) ≈ A · (p / p_th)^((d+1)/2) + B · (p / p_th)^((d+3)/2)
其中第二项代表次主导的“双逻辑错误”路径。对于表面码,B 通常比 A 小一个量级左右,但在 p 接近 p_th 时贡献不可忽略。
四、解析近似与工程判据
解析近似方面,表面码在独立去极化噪声下,阈值约为 p_th ≈ 0.5%~1%(取决于解码器与噪声模型)。在 p 远低于阈值时,常用近似为:
P_L ≈ 0.03 · (p / 0.0057)^((d+1)/2)
该式中的系数和阈值来自特定解码器与噪声模型的拟合,仅作为示例参考,实际项目应基于自身噪声模型重新拟合。工程上判断“是否值得增大码距”的判据是:若当前 p 低于有效阈值,且目标逻辑错误率 P_L^target 满足 P_L(d) ≤ P_L^target,则可通过增大 d 实现;若 p 已接近或超过阈值,增大 d 不会带来收益,此时应优先降低物理错误率或改进解码器。
p < p_th,再用 P_L ≈ A·(p/p_th)^((d+1)/2) 估算达到目标 P_L 所需的最小码距 d_min。若 d_min 超过硬件可承受的物理比特数,则说明当前物理错误率下该目标不可达,需要先降低 p。
五、示例:从目标逻辑错误率反推码距
假设某表面码方案在 p = 0.1% 下测得 p_th = 0.6%,A = 0.03。目标逻辑错误率 P_L^target = 10^-6。代入标度式:
10^-6 = 0.03 · (0.001 / 0.006)^((d+1)/2)
10^-6 / 0.03 = (1/6)^((d+1)/2)
3.33×10^-5 = (1/6)^((d+1)/2)
ln(3.33×10^-5) = ((d+1)/2) · ln(1/6)
-10.31 = ((d+1)/2) · (-1.792)
(d+1)/2 ≈ 5.75
d ≈ 10.5 → 取 d = 11
即需要码距约 11 的表面码。若物理错误率降至 p = 0.05%,则 p/p_th = 1/12,同样目标下 (d+1)/2 ≈ 4.0,d ≈ 7。可见物理错误率减半,所需码距从 11 降到 7,物理比特数从约 2d²=242 降到约 2d²=98,资源开销下降超过一半。这正是“降低物理错误率比增大码距更划算”的定量依据。
六、有限尺寸效应的修正策略
当 d 较小(如 d ≤ 5)或 p 接近阈值时,建议采用以下修正策略:第一,使用带次主导项的拟合式,避免单幂律外推偏差;第二,对每个 d 单独拟合 A 和 p_th,观察其随 d 的变化趋势,若 p_th 随 d 增大而收敛,说明渐近区已接近;第三,在资源评估中预留 1~2 个码距的余量,以覆盖有限尺寸效应和解码器非理想性带来的额外开销。
综上,逻辑错误率与码距的关系可概括为:p < p_th 时 P_L 随 d 指数下降,下降速率由 p/p_th 决定;p > p_th 时关系反转;有限尺寸效应使小码距和近阈值区域偏离理想标度,需用经验拟合与修正项处理。掌握这些关系后,才能在资源开销评估中正确回答“需要多大码距、多少物理比特、物理错误率要降到多少”这三个核心问题。