4. 分子哈密顿量的量子线路构造

要把一个分子的电子结构问题交给量子计算机求解，核心任务是把连续空间的薛定谔方程“翻译”成量子比特能读懂的离散形式。这条翻译链大致是：选定基组 → 确定活性空间 → 计算单/双电子积分 → 二次量子化写出哈密顿量 → 映射为 Pauli 算符串 → 构造量子线路。本节以氢分子 H₂ 与锂氢 LiH 为例，走通这条链路。

一、STO-3G 基组：最小但可用的起点

STO-3G 是 Slater 型轨道的 3 个高斯函数拟合，属于最小基组：每个原子只提供与占据轨道数相同的基函数。对 H 是 1 个 1s，对 Li 是 1s、2s、2p_x、2p_y、2p_z 共 5 个。因此 H₂ 有 2 个空间轨道，LiH 有 6 个空间轨道。空间轨道数 N 对应自旋轨道数 2N，也就对应量子比特数（在 Jordan-Wigner 映射下）。

STO-3G 精度有限，但它的价值在于：哈密顿量规模小、积分可解析计算、便于验证算法正确性。工程上通常先用 STO-3G 跑通流程，再升级到 6-31G、cc-pVDZ 等更大基组。

二、活性空间选择：冻结核心，只算关键轨道

全空间计算代价随轨道数指数增长。化学直觉告诉我们：内层核心轨道几乎不参与成键，可以“冻结”。活性空间记作 CAS(n, m)，n 为活性电子数，m 为活性空间轨道数。冻结轨道对哈密顿量的贡献被折算成一个常数能量项（核心能量 E_core），不进入量子线路。

以 LiH 为例：Li 的 1s 是核心轨道，冻结后剩 2 个价电子进入活性空间。若取 CAS(2, 3)，即 2 个电子分布在 3 个活性轨道（Li 的 2s、2p 与 H 的 1s 组合）中，对应 6 个自旋轨道、6 个量子比特。若取 CAS(2, 5) 则保留全部价轨道，量子比特数升到 10。选择原则是：把占据数接近 0 或 2 的轨道冻结，保留占据数在 0~2 之间显著波动的轨道。

三、积分计算：从原子轨道到分子轨道

量子化学程序（如 PySCF、Psi4）先计算原子轨道基下的单电子积分 h_pq 与双电子积分 (pq|rs)，再通过分子轨道系数 C 做四指标变换，得到分子轨道基下的 h_pq 和 ⟨pq|rs⟩。二次量子化形式下，电子哈密顿量为：

H = E_core + Σ_pq h_pq a†_p a_q
    + ½ Σ_pqrs ⟨pq|rs⟩ a†_p a†_q a_s a_r

其中 p,q,r,s 为自旋轨道指标，a† 与 a 为产生、湮灭算符。E_core 是冻结核心的贡献。这一步的关键是积分索引的对称性与自旋守恒：非零积分必须满足自旋匹配，这能大幅削减项数。

四、Pauli 项分解：Jordan-Wigner 映射

费米子算符不能直接作用于量子比特，需要映射。Jordan-Wigner 变换把每个自旋轨道对应一个量子比特，并引入 Z 串处理反对易关系：

a_j = (X_j + iY_j)/2 ⊗ Z_0 ⊗ ... ⊗ Z_{j-1}
a†_j = (X_j - iY_j)/2 ⊗ Z_0 ⊗ ... ⊗ Z_{j-1}

代入哈密顿量后，H 被展开成 Pauli 串的线性组合 H = Σ_k c_k P_k，P_k ∈ {I, X, Y, Z}^⊗n。例如 H₂ 在 STO-3G、CAS(2,2) 下，经对称性约化后通常得到约 15 个 Pauli 项；LiH 的 CAS(2,3) 则可能产生上百项。项数直接决定后续变分量子本征求解器（VQE）的测量次数与线路深度。

五、从 Pauli 串到量子线路

每个 Pauli 项 P_k 对应一个可测量的观测量。在 VQE 中，用参数化线路 |ψ(θ)⟩ 制备试探态，再对每个 P_k 测量期望值 ⟨P_k⟩，加权求和得到能量 E(θ) = Σ c_k⟨P_k⟩。含 X、Y 的项需在测量前加基变换门（如 H 门、S†H 门）把 X/Y 旋转到 Z 基。线路构造的实质，就是把每个 Pauli 串翻译成“基变换 + CNOT 链 + 旋转 + 逆 CNOT 链”的测量子线路。

六、工程要点与常见陷阱

• 积分阈值截断：小于 10⁻¹⁰ 的积分可丢弃，能显著减少 Pauli 项数而不影响化学精度（约 1.6 mHartree）。
• 对称性约化：利用粒子数守恒、自旋守恒、点群对称性，可把 Pauli 项数压缩数倍。H₂ 的 15 项中，很多项系数为零。
• 量子比特数估算：n_qubit = 2 × m（活性轨道数）。CAS(2,3) → 6 qubit，CAS(2,5) → 10 qubit，CAS(6,6) → 12 qubit。这是当前 NISQ 设备的可行区间。
• 测量分组：对易的 Pauli 项可同时测量，减少线路执行次数。例如所有只含 Z 的项可一次测量完成。

走完这条链路，你就得到了一个可执行的量子线路输入：一组 Pauli 串及其系数、一个参数化试探态线路、一套测量方案。下一步就是选择 VQE 或 QPE 等算法，在优化器驱动下逼近基态能量。STO-3G + 小活性空间是验证整条链路正确性的最佳试验台。