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 个量子比特。
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 的初态制备和测量需要额外的基变换线路,这部分在浅层线路中占比可能不可忽略。
一个可复用的对比检查表
下面给出一个可直接用于项目初筛的对比表,字段按实际估算需要设置。示例数据为 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 数”“估计线路深度”“目标硬件”三列的版本,并在每次更换分子或硬件时重新填写。映射选择不是一次性决定,而是随体系规模和硬件条件动态调整的工程决策。