4. 费米子到量子比特的映射方法对比
要把小分子哈密顿量送上量子计算机，第一步不是选量子门，而是把费米子算符“翻译”成量子比特算符。费米子满足反对易关系，而量子比特是张量积结构，二者之间必须通过一套映射规则建立同构。Jordan-Wigner（JW）、Bravyi-Kitaev（BK）与 Parity 是三种最常用的映射，它们的差别不在于“能不能映射”，而在于占用多少量子比特、产生多少泡利项、局域性如何。本节以氢分子 H₂ 在 STO-3G 基组下的 4 自旋轨道模型为实战对象，把三种映射完整走一遍，并给出可对照的资源开销。
一、先明确费米子哈密顿量的输入形式
H₂ 在 STO-3G 下经对称性约化后，活跃空间为 2 个空间轨道、4 个自旋轨道（记作 0,1,2,3，其中 0、1 为成键轨道 α/β，2、3 为反键轨道 α/β）。二次量子化形式为：
H = h00 a0†a0 + h11 a1†a1 + h22 a2†a2 + h33 a3†a3
+ h0110 a0†a1†a1a0 + h2332 a2†a3†a3a2
+ h0330 a0†a3†a3a0 + h1221 a1†a2†a2a1
+ (h0220 a0†a2†a2a0 + h1331 a1†a3†a3a1) 等双激发项教学示例中取一组合理假设的单电子积分与双电子积分（单位 Hartree），例如 h00≈-1.25、h11≈-1.25、h22≈-0.48、h33≈-0.48，双电子项量级在 0.1～0.7 之间。这些数值仅用于演示映射流程与项数统计，不代表真实量化计算结果。映射的任务，就是把每个 a†、a 替换成泡利串，最终得到 H = Σ cᵢ Pᵢ 的形式，Pᵢ ∈ {I,X,Y,Z}^⊗n。
二、三种映射的构造逻辑
Jordan-Wigner 映射最直接：第 j 个费米子模式对应第 j 个量子比特，占据数 nⱼ=(1-Zⱼ)/2。湮灭算符写成 aⱼ = (Xⱼ + iYⱼ)/2 ⊗ Z₀⊗…⊗Zⱼ₋₁，即用一条“Z 串”记录前面所有模式的占据宇称，从而保证反对易。优点是构造简单、占据数直接可读；缺点是 Z 串长度随模式序号线性增长，非局域性最强。
Bravyi-Kitaev 映射把信息分成“占据数”和“宇称”两部分，存放在一棵二叉树的节点上。第 j 个比特不再单纯表示第 j 个模式的占据，而是表示一段区间的宇称。更新一个模式只需 O(log n) 个比特翻转，因此单激发算符的泡利串长度降为 O(log n)。代价是比特含义不再直观，读取某个轨道占据数需要组合多个比特。
Parity 映射则反过来：第 j 个比特直接存储“前 j 个模式占据数的奇偶性”。这样宇称信息是局域的，但单个模式的占据数变成非局域量。它适合以宇称为主要判据的场景，例如某些对称性约化和纠错编码接口。
三、资源开销的定量对照
对 H₂ 这个 4 比特例子，三种映射的比特数相同（都等于自旋轨道数 4），泡利项总数也都在 15 项左右，差别主要体现在串长分布和测量分组上。JW 下最长串含 3 个 Z；BK 下最长串约 2 个因子；Parity 下宇称测量最省，但占据数测量最贵。若把体系放大到 8～12 个自旋轨道（如 LiH 的活跃空间），JW 的 Z 串长度线性增长，BK 的对数优势才明显拉开。
| 映射 | 比特数 | 最长泡利串 | 单模式更新代价 | 占据数可读性 | 典型适用场景 |
|---|---|---|---|---|---|
| Jordan-Wigner | n | O(n) | O(n) | 直接可读 | 小分子、教学、需频繁测占据 |
| Bravyi-Kitaev | n | O(log n) | O(log n) | 需解码 | 中等规模、门深度敏感 |
| Parity | n | O(n)（占据）/ O(1)（宇称） | O(n) | 需组合 | 宇称判据、纠错编码接口 |
四、实战：把 H₂ 哈密顿量映射成泡利串
以 JW 为例，取单电子项 h00 a0†a0。由 a0†a0 = (I − Z₀)/2，该项直接变成 (h00/2)I − (h00/2)Z₀。再看双电子项 h0110 a0†a1†a1a0，按 JW 展开后 a0†a1†a1a0 = (1/4)(I − Z₀)(I − Z₁) 的变体，最终得到只含 I、Z₀、Z₁、Z₀Z₁ 的泡利串。所有项相加后，用字典合并同类泡利串的系数，得到形如：
H_JW = c_I·IIII + c_Z0·ZIII + c_Z1·IZII
+ c_Z0Z1·ZZII + c_X0X1·XXII + c_Y0Y1·YYII + ...实操中不必手推。用 OpenFermion 的 jordan_wigner()、bravyi_kitaev()、parity() 三个函数即可一键转换，再用 InteractionOperator 读入积分。关键检查点是：转换后泡利串的系数应为实数（厄米性），且 I 项系数应接近体系总能量偏移。若出现明显非厄米或系数异常，通常是积分索引顺序（物理学家约定 vs 化学家约定）弄反了。
五、结果如何判读
映射完成后，应核对三件事：一是泡利项总数与理论估计是否一致（H₂ 约 15 项）；二是最长串长度是否符合该映射的渐近规律；三是把映射后的哈密顿量作用在 Hartree-Fock 参考态 |1100⟩ 上，得到的期望值应接近该构型能量。这三步能在进入变分求解前，快速确认映射链路没有接错。至此，费米子哈密顿量已经变成量子比特可执行的泡利测量序列，可以交给后续的拟设构造与测量分组环节。