2. 量子LDPC码的构造方法

经典LDPC码的核心是稀疏校验矩阵 H:每行代表一个校验约束,每列代表一个码字比特,行重与列重远小于码长。量子LDPC码要继承这种稀疏性,但必须额外满足量子纠错的两个硬约束:稳定子之间必须对易,且错误不能同时被X校验和Z校验“看见”而破坏可测量性。因此,从经典到量子的扩展不是简单替换,而是围绕“如何把稀疏校验矩阵拆成互相兼容的X部分与Z部分”展开。

一、CSS构造:把两个经典码拼成一个量子码

Calderbank-Shor-Steane(CSS)构造是最直接的桥梁。给定两个经典线性码 C1 和 C2,若满足 C2 ⊂ C1 且 C1 的校验矩阵与 C2 的生成矩阵正交,则可定义量子CSS码:X型稳定子来自 C2 的校验矩阵 H2,Z型稳定子来自 C1 的校验矩阵 H1,并要求 H1·H2^T = 0(模2)。这个正交条件是CSS码可构造的判据,也是后续所有量子LDPC构造的出发点。

实际操作中,常见做法是取 C1 = C2^⊥ 的自正交码,或取 C1 为 C2 的扩展码。例如经典 [7,4,3] Hamming码的对偶是 [7,3,4] 单纯码,二者可构成 [[7,1,3]] Steane码。其校验矩阵行重为3或4,列重为3,已经具备LDPC的稀疏雏形。构造时先写出 H1 与 H2,逐行检查内积是否为0,再计算稳定子数量 n-k 与逻辑量子比特数 k = n - rank(H1) - rank(H2)。

二、超图乘积码:从两个经典码生成量子LDPC码

超图乘积码(Hypergraph Product Code, HPC)把CSS构造推广到更一般的稀疏矩阵。给定经典校验矩阵 H1(m1×n1)和 H2(m2×n2),构造量子校验矩阵:

H_X = [ H1 ⊗ I_n2 ,  I_m1 ⊗ H2^T ]
H_Z = [ I_n1 ⊗ H2 ,  H1^T ⊗ I_m2 ]

其中 ⊗ 为Kronecker积,I 为单位矩阵。可以验证 H_X·H_Z^T = H1⊗H2^T + H1⊗H2^T = 0(模2),因此自动满足CSS对易条件。若 H1、H2 是稀疏的,则 H_X、H_Z 也稀疏,码长约为 n1·n2 + m1·m2,逻辑量子比特数约为 k1·k2 + k1^T·k2^T。HPC的优点是参数可解析计算,且当两个经典码都接近容量时,量子码也能获得线性距离增长。

三、表面码与颜色码的LDPC视角

表面码可以看作HPC的特例:取 H1 = H2 = H 为二维网格的关联矩阵,则 H_X 对应顶点稳定子,H_Z 对应面稳定子,每个稳定子作用在4个数据量子比特上,行重列重均为4,是典型的稀疏校验结构。颜色码则把网格三色化,每个面同时承载X和Z稳定子,其校验矩阵行重为6,列重为3,同样满足LDPC的稀疏性。二者区别在于:表面码的逻辑算符是边界上的链,颜色码的逻辑算符是可穿过任意颜色的膜,这导致颜色码的译码更复杂但支持横向Clifford门更丰富。

从经典LDPC到量子CSS码的构造路径 经典LDPC码 稀疏校验矩阵 H 行重/列重 ≪ n Tanner图:变量节点+校验节点 CSS构造 取 C2 ⊂ C1,H1·H2^T = 0 X稳定子 ← H2 Z稳定子 ← H1 k = n − rank(H1) − rank(H2) 超图乘积码 HPC H_X = [H1⊗I, I⊗H2^T] H_Z = [I⊗H2, H1^T⊗I] 自动满足 H_X·H_Z^T = 0 n ≈ n1n2 + m1m2 表面码(Surface Code) HPC特例:H1 = H2 = 网格关联矩阵 顶点稳定子(X型)+ 面稳定子(Z型) 行重=列重=4,逻辑算符为边界链 译码:最小权匹配(MWPM) 颜色码(Color Code) 三色化网格,每面同时承载X/Z稳定子 行重=6,列重=3,稀疏LDPC结构 逻辑算符为可穿色膜,支持横向Clifford 译码:投影/色码匹配,复杂度更高 共同点:稀疏校验矩阵、稳定子对易、可迭代译码 | 差异:稳定子权重、逻辑算符几何、译码算法与门实现方式 构造选择依据:目标码率、距离、译码延迟、硬件连接度

四、构造案例对比与选择依据

下面用三个典型构造案例说明参数差异。所有数值均为按公式计算得到的示例,用于对比结构特征,不代表特定硬件实测结果。

构造方法 码参数示例 稳定子权重 逻辑算符 译码特点
Steane码(CSS) [[7,1,3]] 3或4 全局X/Z算符 可查表译码,距离小
HPC(两个[7,4,3]码) [[25,1,3]]量级 4~6 乘积链 迭代置信传播+OSD
旋转表面码 [[d²,1,d]],d=5时[[25,1,5]] 4 边界链 MWPM,阈值约1%
颜色码 [[7,1,3]](三角面) 6 穿色膜 投影译码,阈值约0.5%

选择构造时,先看目标:若需要高码率且硬件连接度有限,优先HPC或表面码;若需要横向Clifford门简化容错逻辑,颜色码更有优势;若只是教学验证,Steane码最小可测。实际实现中,还要检查校验矩阵的稀疏度是否匹配硬件耦合图——例如表面码要求数据量子比特排成二维网格,每个稳定子只连接最近邻,这直接决定了超导芯片的布线方案。

构造检查清单:
  • 写出 H1 与 H2,确认 H1·H2^T = 0(模2)。
  • 计算 rank(H1)、rank(H2),得到 k = n − rank(H1) − rank(H2)。
  • 检查每行/每列非零元个数,确认稀疏性(通常 ≤ 6)。
  • 确认稳定子作用范围与硬件耦合图一致,避免长程连接。
  • 选择译码器:表面码用MWPM,HPC用BP+OSD,颜色码用投影匹配。

从经典LDPC到量子CSS码,再到HPC、表面码和颜色码,构造方法的核心始终是“稀疏校验矩阵 + 对易条件 + 可译码结构”。理解这三者的权衡,就能根据码率、距离、译码复杂度和硬件约束,选出合适的量子LDPC码构造方案。