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门更丰富。
四、构造案例对比与选择依据
下面用三个典型构造案例说明参数差异。所有数值均为按公式计算得到的示例,用于对比结构特征,不代表特定硬件实测结果。
| 构造方法 | 码参数示例 | 稳定子权重 | 逻辑算符 | 译码特点 |
|---|---|---|---|---|
| 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码构造方案。