4. 费米子到量子比特映射的资源对比

在量子计算化学中,分子哈密顿量天然写在费米子二次量子化框架里,而量子硬件只认识量子比特上的 Pauli 算符。把费米子算符翻译成量子比特算符的过程,就是费米子—量子比特映射。这一步不是无关紧要的“格式转换”:它直接决定需要多少个量子比特、哈密顿量展开成多少条 Pauli 项、每条项作用在哪些比特上,以及最终编译出的线路深度。资源估算工程中,映射选择往往比后续的线路优化更早、更根本地决定了整个算法的可行性。

三种主流映射的基本构造

设自旋轨道数为 N,则三种映射都使用 N 个量子比特,但编码方式不同。

Jordan-Wigner(JW)映射最直观:第 p 个自旋轨道的占据数直接对应第 p 个量子比特的 Z 本征值,即 np = (1 − Zp)/2。产生/湮灭算符写成

a_p† = (X_p − iY_p)/2 · Z_0 Z_1 ... Z_{p−1}
a_p  = (X_p + iY_p)/2 · Z_0 Z_1 ... Z_{p−1}

其中 Z 串称为 Jordan-Wigner 弦,用来保证不同轨道上的费米子算符满足反对易关系。JW 的优点是局域占据数含义清晰、变分 ansatz 容易构造;缺点是弦的长度随轨道序号线性增长,高序号轨道的算符会拖出很长的 Pauli 串。

Bravyi-Kitaev(BK)映射用“占据数 + 宇称”的混合编码替代纯占据数编码。它借助二叉树结构,把每个轨道的更新信息分散到 O(log N) 个比特上,因此 JW 弦的长度从 O(N) 降到 O(log N)。代价是单个量子比特不再直接对应某个轨道的占据数,测量和初态制备需要额外的基变换。

Parity 映射则反过来:第 p 个量子比特编码前 p 个轨道的累积宇称。它的弦结构同样可以做到 O(log N) 量级,并且在只关心粒子数守恒扇区时,可以利用宇称约束减少一个量子比特(2-qubit reduction),把 N 个自旋轨道压缩到 N−1 个量子比特。

三种映射的编码结构与资源特征对比(N = 8 自旋轨道示例) Jordan-Wigner 编码:n_p ↔ Z_p(占据数) q0 q1 q2 q3 q4 q5 q6 a₆† 的 JW 弦: Z0 Z1 Z2 Z3 Z4 Z5 X6 弦长 = p,最坏 O(N) 资源特征 • 量子比特数:N • Pauli 项数:O(N⁴) • 平均权重:O(N) • 线路深度:较深 • 初态制备:简单 • 测量:单比特 Z 适合小体系、教学与 ansatz 设计 Bravyi-Kitaev 编码:占据数 + 宇称(二叉树) q0 q1 q2 q3 q4 q5 q6 更新集大小 = ⌈log₂(N+1)⌉ = 3 资源特征 • 量子比特数:N • Pauli 项数:O(N⁴) • 平均权重:O(log N) • 线路深度:中等 • 初态制备:需基变换 • 测量:需解码 适合中等规模、追求低权重 Parity 映射 编码:累积宇称 P_p = n₀⊕…⊕n_p P0 P1 P2 P3 P4 P5 宇称约束:P_{N−1} = 粒子数奇偶性 → 可去掉 1 个量子比特 资源特征 • 量子比特数:N−1(2-qubit reduction) • Pauli 项数:O(N⁴) • 平均权重:O(log N) • 线路深度:中等 • 初态制备:需宇称基 • 测量:需解码 适合粒子数守恒体系 省 1 比特,但解码开销增加

Pauli 项数与权重:真正的资源瓶颈

三种映射下,哈密顿量展开后的 Pauli 项总数都随自旋轨道数 N 以 O(N⁴) 增长,因为分子哈密顿量包含双电子积分项 hpqrs,每个双激发算符在映射后会产生若干 Pauli 串。真正拉开差距的是每条 Pauli 项的平均权重(即非恒等 Pauli 因子的个数)。

以 N = 8 为例,JW 映射下高序号轨道的算符会携带长度接近 8 的 Z 串,平均权重约为 O(N);而 BK 和 Parity 映射通过二叉树或宇称结构,把平均权重压到 O(log N),在 N = 8 时约为 3~4。权重直接决定每条 Pauli 项编译成线路时需要多少个 CNOT:权重为 w 的项通常需要 2(w−1) 个 CNOT 来对角化。因此,当 N 增大时,JW 的线路深度会显著膨胀,而 BK/Parity 的增长要平缓得多。

需要强调的是,Pauli 项数本身在三种映射下数量级相同,但项数与权重的乘积才近似反映总 CNOT 消耗。JW 的项数可能略少(因为部分项被合并),但权重高;BK/Parity 项数可能略多,但权重低。实际选择时应对具体分子做一次统计,而不是凭经验拍板。

量子比特数的差异

JW 和 BK 都使用 N 个量子比特。Parity 映射在粒子数守恒的化学体系中,可以利用总宇称守恒把最后一个量子比特消去,实现 2-qubit reduction,即用 N−1 个量子比特。对于 N = 8 的小分子,这只省 1 个比特;但在 N = 50 以上的体系中,省下的比特数虽然比例不大,却可能正好让问题落进某个硬件的连通图容量之内。需要注意的是,减少比特会引入额外的解码线路和测量后处理,这部分开销必须计入总资源。

线路深度与硬件连通性

映射选择还会影响线路与硬件拓扑的匹配程度。JW 的长 Z 串在编译时会被拆成一系列 CNOT,如果硬件是线性连通或重六角连通,长距离 CNOT 需要插入 SWAP,进一步加深线路。BK 和 Parity 的 O(log N) 权重意味着更短的 Pauli 串,编译后 CNOT 层数更少,对连通性较差的硬件更友好。但 BK/Parity 的初态制备和测量需要额外的基变换线路,这部分在浅层线路中占比可能不可忽略。

工程建议:在资源估算流程中,先对目标分子分别用 JW、BK、Parity 生成哈密顿量,统计三个量:量子比特数、Pauli 项总数、加权平均权重(Σ wi / 项数)。再结合目标硬件的连通图和原生门集,估算编译后的 CNOT 数与线路深度。不要只看项数,权重和连通性往往才是决定因素。

一个可复用的对比检查表

下面给出一个可直接用于项目初筛的对比表,字段按实际估算需要设置。示例数据为 N = 8 的假设性统计,仅用于说明填写方式,不代表任何真实分子的实测结果。

映射 量子比特数 Pauli 项数(量级) 平均权重 初态制备 测量解码 适用场景
Jordan-Wigner N O(N⁴) O(N) 简单 单比特 Z 小体系、教学、ansatz 设计
Bravyi-Kitaev N O(N⁴) O(log N) 需基变换 需解码 中等规模、追求低权重
Parity N−1 O(N⁴) O(log N) 需宇称基 需解码 粒子数守恒体系

在实际项目中,建议把这张表扩展为包含“编译后 CNOT 数”“估计线路深度”“目标硬件”三列的版本,并在每次更换分子或硬件时重新填写。映射选择不是一次性决定,而是随体系规模和硬件条件动态调整的工程决策。