4. 分子哈密顿量的二次量子化表示
在上一节中，我们把分子电子结构问题化归为一个在给定基组下的矩阵本征值问题。但要在量子计算机上求解它，还需要把“电子在分子轨道上的分布”翻译成量子比特能理解的语言。这一节要完成的核心任务，就是建立从分子轨道到量子比特的完整映射链条：分子轨道 → 占据数表象 → 费米子算符 → 量子比特算符。理解这条链条，是后续构造变分量子线路、测量能量期望值的前提。
一、从分子轨道到二次量子化
分子轨道（MO）是单电子波函数，通常由原子轨道线性组合（LCAO）得到。设分子有 N 个空间轨道，考虑自旋后得到 2N 个自旋轨道。在二次量子化框架下，我们不再追踪每个电子的坐标，而是用“每个自旋轨道上有几个电子”来描述状态。对于费米子，每个自旋轨道的占据数只能是 0 或 1（泡利不相容原理）。
引入产生算符 ap† 和湮灭算符 ap，它们分别表示在自旋轨道 p 上增加或移除一个电子。这些算符满足反对易关系：
{a_p, a_q†} = δ_pq
{a_p, a_q} = 0
{a_p†, a_q†} = 0反对易关系是费米子统计的直接体现，也是后续所有变换必须保持的核心结构。分子的电子哈密顿量在二次量子化形式下写成：
H = Σ_pq h_pq a_p† a_q + 1/2 Σ_pqrs h_pqrs a_p† a_q† a_r a_s + E_nuc其中 hpq 是单电子积分（动能加核吸引），hpqrs 是双电子积分（电子间库仑排斥），Enuc 是核间排斥能。对于 Norb 个空间轨道，双电子积分项的数量约为 O(Norb4)，这是分子哈密顿量复杂度的主要来源。
二、占据数表象与费米子算符的作用
占据数表象用一组占据数 |n0n1…n2N-1⟩ 表示多电子态，其中 np∈{0,1}。例如，|1100⟩ 表示前两个自旋轨道被占据，后两个为空。这种表示天然适合量子比特编码，因为每个占据数恰好对应一个量子比特的 |0⟩ 或 |1⟩。
但费米子算符作用在占据数态上时，会带有一个符号因子。具体地：
a_p† |n_0 … n_p …⟩ = (-1)^{Σ_{q<p} n_q} (1-n_p) |n_0 … 1 …⟩
a_p |n_0 … n_p …⟩ = (-1)^{Σ_{q<p} n_q} n_p |n_0 … 0 …⟩这个 (-1) 的符号来自费米子波函数的反对称性，称为宇称因子或 Jordan-Wigner 弦。它意味着：当我们在量子比特上直接实现 ap 或 ap† 时，必须把这种非局域的符号信息编码进去。这正是 Jordan-Wigner 和 Bravyi-Kitaev 变换要解决的问题。
三、Jordan-Wigner 变换
Jordan-Wigner（JW）变换是最直接的费米子-量子比特映射。它把每个自旋轨道的占据数直接对应到一个量子比特，并规定：
a_p = (Π_{q<p} Z_q) ⊗ (X_p + iY_p)/2
a_p† = (Π_{q<p} Z_q) ⊗ (X_p - iY_p)/2其中 Zq 是作用在第 q 个量子比特上的泡利 Z 算符，Xp、Yp 是泡利 X、Y 算符。直观理解：第 p 个量子比特上的 X 和 Y 负责翻转占据数，而前面所有量子比特上的 Z 算符串（即 JW 弦）负责产生正确的宇称符号。
JW 变换的优点是形式简单、易于实现；缺点是 JW 弦的长度随轨道序号线性增长，导致高序号轨道上的算符变成多量子比特的长程相互作用。在量子线路中，这意味着需要更多的 CNOT 门来模拟，线路深度较大。
四、Bravyi-Kitaev 变换
Bravyi-Kitaev（BK）变换采用不同的编码策略：它不把占据数直接放在单个量子比特上，而是把占据数信息分散到一组量子比特中，使得宇称信息可以部分局域化。BK 变换的核心思想是利用二叉树结构，将每个轨道的占据数编码为“部分和”形式。
在 BK 编码下，第 p 个轨道的宇称信息只需要 O(log p) 个量子比特的 Z 算符串，而不是 JW 中的 O(p) 个。因此，对于大规模分子，BK 变换通常能显著降低量子线路的深度和门数。代价是编码和解码逻辑更复杂，测量和初始化时需要额外的处理。
下表对比两种变换的关键特性：
| 特性 | Jordan-Wigner | Bravyi-Kitaev |
|---|---|---|
| 量子比特数 | 2N（自旋轨道数） | 2N |
| 宇称弦长度 | O(p) | O(log p) |
| 线路深度 | 较大 | 较小 |
| 实现复杂度 | 低 | 中高 |
| 适用场景 | 小分子、教学演示 | 较大分子、追求线路效率 |
五、映射链条的整体视图
下面这张 SVG 图展示了从分子轨道到量子比特的完整映射链条，以及 JW 与 BK 两种变换在“宇称弦长度”上的差异。图中左侧是物理对象，中间是二次量子化表示，右侧是量子比特编码，底部用条形长度直观对比两种变换的弦长增长趋势。
六、实际实现中的关键细节
在具体实现时，有几个细节需要特别注意。第一，自旋轨道排序：通常约定先排 α 自旋轨道，再排 β 自旋轨道，或者按空间轨道交替排列。不同的排序会影响 JW 弦的具体形式，但不会改变物理结果。第二，粒子数守恒：分子哈密顿量保持电子数不变，因此 JW 或 BK 变换后的算符应保持粒子数守恒。在变分量子线路中，可以通过只使用粒子数守恒的激发算符来约束搜索空间。第三，积分对称性：hpqrs 具有交换对称性，可以利用这些对称性减少独立积分数量，从而降低哈密顿量项数。
以一个简化示例说明：假设只有 2 个自旋轨道（对应 1 个空间轨道），占据数态为 |10⟩ 和 |01⟩。在 JW 变换下，a0†a1 对应 (X0+iY0)/2 · Z0 · (X1-iY1)/2，展开后得到 X0Y1、Y0X1 等项。这些项可以直接用标准量子门实现。对于更大的分子，项数会迅速增加，因此需要借助 OpenFermion 等工具自动生成变换后的哈密顿量。
jordan_wigner() 和 bravyi_kitaev() 函数直接对费米子算符进行变换。建议先用小分子（如 H₂，4 个自旋轨道）验证变换的正确性，再扩展到更大体系。验证方法：比较变换前后哈密顿量在某个参考态上的期望值是否一致。
总结来说，本节建立了从分子轨道到量子比特的完整映射。分子轨道提供单电子基函数，二次量子化用产生/湮灭算符描述电子关联，占据数表象为量子比特编码提供直接对应，而 JW 和 BK 变换则解决了费米子反对易关系带来的符号问题。理解这条链条，是后续构造变分量子线路、测量分子能量的基础。