3. 费米子到量子比特的映射工程

在分子电子结构量子计算中,我们真正要处理的对象是电子哈密顿量。它天然写在二次量子化的费米子代数里:产生算符 \(a_p^\dagger\) 与湮灭算符 \(a_q\) 满足反对易关系,而量子计算机的硬件语言却是量子比特上的泡利算符 \(X,Y,Z\) 与受控门。映射工程要解决的,就是如何把费米子算符“翻译”成量子比特算符,同时尽量少花量子比特、少用门、并让结果可测量。本节围绕 Jordan-Wigner(JW)、Bravyi-Kitaev(BK)、Parity 三类主流映射,讲清原理、资源开销、局域性差异,并给出工程选型准则。

一、为什么不能直接“一个轨道一个比特”

最朴素的想法是:每个自旋轨道对应一个量子比特,比特为 1 表示该轨道被占据,为 0 表示空。但费米子波函数在交换两个电子时会整体变号,这个“费米子统计”必须被编码进去。JW 的做法是:用单比特的 \(Z\) 算符记录该轨道之前的占据数奇偶性,形成一条“宇称链”。例如在 JW 下,

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_0\cdots Z_{p-1}\) 就是那条链。它保证了反对易关系,但代价是:作用在第 \(p\) 个轨道上的算符,会牵动前面所有比特。这就是 JW 的局域性问题。

二、三种映射的核心思想

Jordan-Wigner:最直观,占据数直接对应比特值,宇称信息用一条 \(Z\) 链串起来。它的优点是物理意义清晰、测量方便;缺点是长链导致高权重泡利项,门深度大。

Bravyi-Kitaev:不直接存“每个轨道的占据数”,而是把占据信息拆成“部分和”与“部分奇偶性”,用一棵平衡二叉树组织。这样,单个费米子算符的 \(Z\) 链长度从 \(O(N)\) 降到 \(O(\log N)\),但代价是每个比特不再直接对应一个轨道的占据数,测量时需要额外的基变换。

Parity:把比特编码为“从某个参考点到当前轨道的累积宇称”。它的最大优势是:分子哈密顿量中常见的粒子数守恒项、以及某些双激发项,在 Parity 下会变成非常局域的算符,甚至直接消失。但代价是,单粒子项反而变得非局域,且需要处理“宇称约束”带来的冗余比特。

三、量子比特资源开销对比

设分子有 \(N\) 个自旋轨道。JW 和 BK 通常使用 \(N\) 个量子比特;Parity 可以利用粒子数与自旋对称性,把比特数降到 \(N-2\) 甚至更少。但“少比特”不等于“少门”。下表给出典型量级对比(以 \(N\) 个轨道、粒子数守恒的分子哈密顿量为例):

映射量子比特数单算符 Z 链长度典型泡利项权重测量复杂度
Jordan-WignerNO(N)高低(直接测量)
Bravyi-KitaevNO(log N)中中(需基变换)
ParityN-2 或更少O(N)(单粒子项)低(双激发项)高(需解码)

需要强调:这里的“权重”指泡利字符串中非恒等算符的个数。权重越高,分解为受控门时所需的 CNOT 越多,电路越深。JW 的权重随轨道序号线性增长,BK 把它压到对数级,Parity 则把双激发项压到常数级,但把单粒子项推高。

四、局域性差异:谁更“好做”

局域性直接决定量子线路深度和测量方案。JW 的算符作用范围大,但结构规整,适合小分子和教学演示;BK 的算符作用范围小,适合中等规模、需要控制线路深度的场景;Parity 在双激发主导的关联能计算中非常高效,但单粒子项和测量解码会带来额外开销。实际工程中,常把三种映射混合使用:对双激发项用 Parity,对单粒子项用 JW 或 BK,以取得折中。

三种费米子映射的局域性与资源开销对比 Jordan-Wigner 比特:N 个 Z 链:O(N) 权重:高 测量:直接 适用:小分子、教学 占据数 ↔ 比特值 Bravyi-Kitaev 比特:N 个 Z 链:O(log N) 权重:中 测量:需基变换 适用:中等规模 二叉树部分和编码 Parity 比特:N-2 或更少 Z 链:单粒子 O(N) 权重:双激发低 测量:需解码 适用:关联能计算 累积宇称编码 局域性示意:作用在第 p 个轨道上的算符牵动多少比特 JW p Z 链覆盖 0 到 p-1,长度 O(N) BK p 二叉树路径长度 O(log N) Parity p 累积宇称,双激发项局域化

五、工程选型准则

选型不是“哪个最好”,而是“哪个最适合当前任务”。可以从四个维度判断:

  • 量子比特预算:若硬件比特数紧张,优先 Parity,利用对称性减少比特;若比特充足但门深度受限,优先 BK。
  • 线路深度:JW 的 Z 链长,深线路下噪声累积严重;BK 的 \(O(\log N)\) 链更适合含噪中等规模设备。
  • 测量方案:JW 可直接测量,适合需要频繁读取占据数的场景;Parity 需要解码,适合一次性测量关联能。
  • 哈密顿量结构:双激发主导时 Parity 优势明显;单粒子项主导时 JW 或 BK 更直接。
工程建议:在真实项目中,先用 JW 快速验证哈密顿量构造与测量逻辑,再根据线路深度和比特数决定是否切换到 BK 或 Parity。切换时务必重新检查费米子算符的符号约定与比特顺序,避免因映射不一致导致能量偏差。
常见陷阱:不同映射下,同一个费米子算符的矩阵表示不同,但物理结果必须一致。若发现能量随映射变化,通常不是映射本身的问题,而是比特顺序、宇称约束或测量基变换写错了。建议用 \(H_2\) 或 \(LiH\) 的小基组做交叉验证。

总结来说,JW 是“直译”,BK 是“压缩链长”,Parity 是“压缩比特并局域化双激发”。工程选型时,先看比特预算,再看线路深度,最后看测量与哈密顿量结构。三者并非互斥,混合映射往往是实际项目中的最优解。