2. 量子化学问题的二次量子化建模

要把一个分子放进量子计算机里求解，第一步不是写量子线路，而是把化学问题翻译成量子比特能承载的形式。这条翻译链有四个关键环节：从分子结构出发构建二次量子化哈密顿量，把费米子算符映射成泡利算符，选择基组决定轨道数与资源规模，再通过活性空间截取把问题压缩到可执行范围。本节按这条链路逐层展开，并给出可操作的参数与判断依据。

一、分子哈密顿量的二次量子化形式

在 Born-Oppenheimer 近似下，电子在固定核骨架中运动，分子电子哈密顿量写成二次量子化形式：

H = Σ_{pq,σ} h_{pq} a†_{pσ} a_{qσ}
  + ½ Σ_{pqrs,στ} h_{pqrs} a†_{pσ} a†_{qτ} a_{sτ} a_{rσ}
  + E_nuc

其中 a† 与 a 是费米子产生、湮灭算符，满足反对易关系 {a_p, a†_q} = δ_pq。单电子积分 h_{pq} 描述电子在核势场中的动能与势能，双电子积分 h_{pqrs} 描述电子间库仑排斥，E_nuc 是核-核排斥常数。注意指标顺序：h_{pqrs} 对应 a†_p a†_q a_s a_r，这个顺序决定了后续映射中泡利串的符号，写错会导致能量偏差。

实际计算中，这些积分由经典量子化学程序（PySCF、Psi4、OpenFermion 等）在给定基组下算出。以 STO-3G 基组下的 H₂ 分子（键长 0.735 Å）为例，只有 2 个空间轨道、4 个自旋轨道，积分矩阵规模很小；而同样的 H₂ 在 6-31G 下空间轨道数翻倍，积分矩阵维度按 N⁴ 增长。这就是基组选择直接影响资源的第一层原因。

二、费米子到量子比特的映射

量子比特是玻色型二能级系统，不天然满足反对易关系，因此必须通过映射把费米子算符“编码”进泡利算符。常用方法有三类：

• Jordan-Wigner（JW）：每个自旋轨道对应一个量子比特，用 Z 串表示占据数奇偶性。实现简单、局域性好，但算符权重随轨道数线性增长，线路深度大。
• Bravyi-Kitaev（BK）：用部分和结构降低算符权重至 O(log N)，适合中等规模分子，但线路结构更复杂。
• Parity / 对称性约化映射：利用粒子数与自旋对称性，把 N 个自旋轨道压缩到 N−2 个量子比特，直接减少比特数。

以 JW 映射为例，H₂ 的 4 个自旋轨道对应 4 个量子比特，哈密顿量展开后约 15 个泡利项；若换成 6-31G 的 8 个自旋轨道，泡利项数会显著增加。映射方式的选择不是“越高级越好”，而是要在比特数、线路深度和硬件连通性之间权衡。

三、基组选择对资源的影响

基组决定空间轨道数 N，而 N 直接决定自旋轨道数 2N、量子比特数和哈密顿量项数。以 H₂O 为例：

注意：比特数只是最直观的指标，泡利项数和线路深度增长更快。STO-3G 适合教学演示和硬件验证，6-31G 及以上适合需要定量精度的场景，但必须配合活性空间截取，否则当前硬件无法承载。

四、活性空间选取策略

活性空间的核心思想是：把化学上“不活跃”的轨道冻结掉，只让真正参与成键、断键和电子相关的轨道进入量子计算。常用策略有三条：

1. 冻结芯轨道：内层电子几乎不参与化学反应，直接冻结。例如 H₂O 的 O 1s 轨道可冻结，节省 2 个自旋轨道。
2. 冻结高能虚轨道：能量远高于费米能级的空轨道贡献极小，可截断。
3. 保留价层与近简并轨道：价层轨道、孤对电子、以及能量接近 HOMO/LUMO 的轨道必须保留，它们主导电子相关效应。

以 H₂O / STO-3G 为例，7 个空间轨道中，O 1s 冻结后剩 6 个；若再冻结 2 个高能虚轨道，活性空间为 4 个空间轨道、8 个自旋轨道，对应 8 个量子比特（JW）。这个规模在含噪声硬件上仍可执行 VQE。判断活性空间是否合理，可以看两个指标：一是冻结前后基态能量差是否小于化学精度（约 1.6 mHa），二是活性空间内的自然轨道占据数是否接近 0 或 2（说明该轨道可冻结）。