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-Wigner | N | O(N) | 高 | 低(直接测量) |
| Bravyi-Kitaev | N | O(log N) | 中 | 中(需基变换) |
| Parity | N-2 或更少 | O(N)(单粒子项) | 低(双激发项) | 高(需解码) |
需要强调:这里的“权重”指泡利字符串中非恒等算符的个数。权重越高,分解为受控门时所需的 CNOT 越多,电路越深。JW 的权重随轨道序号线性增长,BK 把它压到对数级,Parity 则把双激发项压到常数级,但把单粒子项推高。
四、局域性差异:谁更“好做”
局域性直接决定量子线路深度和测量方案。JW 的算符作用范围大,但结构规整,适合小分子和教学演示;BK 的算符作用范围小,适合中等规模、需要控制线路深度的场景;Parity 在双激发主导的关联能计算中非常高效,但单粒子项和测量解码会带来额外开销。实际工程中,常把三种映射混合使用:对双激发项用 Parity,对单粒子项用 JW 或 BK,以取得折中。
五、工程选型准则
选型不是“哪个最好”,而是“哪个最适合当前任务”。可以从四个维度判断:
- 量子比特预算:若硬件比特数紧张,优先 Parity,利用对称性减少比特;若比特充足但门深度受限,优先 BK。
- 线路深度:JW 的 Z 链长,深线路下噪声累积严重;BK 的 \(O(\log N)\) 链更适合含噪中等规模设备。
- 测量方案:JW 可直接测量,适合需要频繁读取占据数的场景;Parity 需要解码,适合一次性测量关联能。
- 哈密顿量结构:双激发主导时 Parity 优势明显;单粒子项主导时 JW 或 BK 更直接。
总结来说,JW 是“直译”,BK 是“压缩链长”,Parity 是“压缩比特并局域化双激发”。工程选型时,先看比特预算,再看线路深度,最后看测量与哈密顿量结构。三者并非互斥,混合映射往往是实际项目中的最优解。