3. 量子比特映射与费米子编码

在量子化学的量子算法实现中，电子结构问题天然定义在费米子二次量子化框架下：哈密顿量由产生算符 \(a_p^\dagger\) 和湮灭算符 \(a_q\) 构成，它们满足反对易关系 \(\{a_p, a_q^\dagger\} = \delta_{pq}\)。而量子计算机的硬件语言是量子比特上的泡利算符 \(X, Y, Z\) 及其张量积。因此，从费米子到量子比特的映射（fermion-to-qubit mapping）是连接化学问题与量子线路的第一道关键桥梁。映射方式的选择不仅决定了哈密顿量的项数与形式，还直接影响量子线路的深度、测量方案以及算法整体复杂度。

为什么需要费米子编码

费米子态具有两个核心特征：一是占据数只能为 0 或 1（泡利不相容原理），二是交换两个费米子会产生 -1 的相位（反对称性）。量子比特天然满足第一个条件——\(|0\rangle\) 表示未占据，\(|1\rangle\) 表示占据。但第二个条件需要额外编码：当我们在某个轨道上作用产生或湮灭算符时，必须根据该轨道之前有多少个轨道被占据，附加相应的符号。这正是不同编码方案的核心差异所在。

以 \(N\) 个自旋轨道为例，一个费米子态可写为 \(|n_0 n_1 \cdots n_{N-1}\rangle\)，其中 \(n_i \in \{0,1\}\)。产生算符的作用为：

\[ a_p^\dagger |n_0 \cdots n_{N-1}\rangle = (-1)^{\sum_{i

关键在于那个 \((-1)^{\sum_{i

Jordan-Wigner 编码：直观但非局域

Jordan-Wigner（JW）编码是最直接的方案。它将第 \(p\) 个轨道的占据数直接映射到第 \(p\) 个量子比特，并用一条从比特 0 到比特 \(p-1\) 的 \(Z\) 算符链来产生宇称符号：

\[ a_p = \left(\prod_{i=0}^{p-1} Z_i\right) \frac{X_p + iY_p}{2}, \quad a_p^\dagger = \left(\prod_{i=0}^{p-1} Z_i\right) \frac{X_p - iY_p}{2} \]

JW 编码的优点是物理意义清晰、占据数直接可读。但它的致命弱点是 \(Z\) 链的长度随轨道序号线性增长。对于 \(N\) 个轨道，最长的算符链涉及 \(O(N)\) 个量子比特。当我们将哈密顿量中的费米子项展开为泡利项时，一个双体项 \(a_p^\dagger a_q^\dagger a_r a_s\) 可能产生权重高达 \(O(N)\) 的泡利串。这意味着量子线路中需要 \(O(N)\) 个受控门来实现一个项，整体线路深度可能达到 \(O(N^4)\) 甚至更高。

Bravyi-Kitaev 编码：用树结构压缩局域性

Bravyi-Kitaev（BK）编码的核心思想是：不直接存储占据数，而是存储部分和（宇称信息）与占据信息的混合。它将轨道索引组织成一棵二叉树，每个量子比特存储一个“部分和”或“占据数”。具体来说，BK 编码将第 \(j\) 个量子比特定义为：

\[ |b_j\rangle = \left| \left(\sum_{i \in S_j} n_i\right) \bmod 2 \right\rangle \]

其中 \(S_j\) 是由二叉树结构决定的轨道子集。这种设计使得产生/湮灭算符的更新只需 \(O(\log N)\) 个量子比特操作，而非 JW 的 \(O(N)\)。对于哈密顿量中的双体项，BK 编码产生的泡利串权重通常为 \(O(\log N)\)，显著降低了线路深度。

但 BK 编码也有代价：占据数不再直接可读，需要额外的解码步骤；且对于某些特定项，其泡利展开可能比 JW 更复杂。实际选择时需根据哈密顿量结构和硬件连接性综合评估。

Parity 编码：以宇称为基

Parity 编码可以看作 JW 的“对偶”方案。它不存储单个轨道的占据数，而是存储从第 0 个到第 \(p\) 个轨道的累积宇称：

\[ |q_p\rangle = \left| \left(\sum_{i=0}^{p} n_i\right) \bmod 2 \right\rangle \]

在这种编码下，宇称算符是局域的（单个 \(Z_p\)），而单个轨道的占据数需要长链 \(Z\) 算符来提取。Parity 编码特别适合需要频繁测量宇称或粒子数的场景，例如在量子化学的粒子数守恒模拟中。它与 JW 编码之间存在简单的基变换关系，可以互相转换。

三种编码的复杂度对比

下表从几个关键维度对比三种编码方案（以 \(N\) 个自旋轨道为参考）：

需要强调的是，项数本身在三种编码下通常都是 \(O(N^4)\) 量级（因为双体相互作用有 \(O(N^4)\) 个组合），但每个项展开为泡利串后的权重和线路实现代价差异巨大。JW 编码下，一个双体项可能展开为多个权重为 \(O(N)\) 的泡利串；而 BK 编码下，同样项展开的泡利串权重降至 \(O(\log N)\)。在含噪声中等规模量子（NISQ）设备上，这意味着 BK 编码可以显著减少 CNOT 门数量和线路深度，从而降低噪声累积。

实际选择与工程考量

在实际量子化学模拟中，编码选择需要结合具体问题。对于小分子（如 \(H_2\)、\(LiH\)）的 NISQ 模拟，JW 编码因其简单直观常被优先采用，且 \(N\) 较小（\(N \leq 12\)）时 \(O(N)\) 的链长尚可接受。对于更大体系（如 \(FeMoco\) 活性空间，\(N \approx 50\)），BK 编码可将泡利串权重从 50 降至约 6，线路深度降低一个数量级，这对当前硬件的相干时间至关重要。

此外，还有更高级的编码方案如对称性保持编码（symmetry-preserving encoding）和低秩编码（low-rank encoding），它们利用粒子数守恒和自旋对称性进一步压缩希尔伯特空间。但这些方案通常以增加线路复杂度为代价，需根据具体硬件和算法目标权衡。

总结来说，费米子编码不是单纯的数学变换，而是决定量子算法资源开销的关键设计决策。JW 编码提供了最直观的起点，BK 编码通过树结构实现了对数级局域性，Parity 编码则在宇称测量场景中展现优势。理解三者的结构差异与复杂度特征，是进行量子化学算法工程优化的基础能力。