3. 二次量子化与哈密顿量构造的资源影响

在量子计算化学中，我们最终要处理的对象并不是薛定谔方程本身，而是经过二次量子化后得到的哈密顿量。二次量子化把 N 电子体系的多体波函数转化为费米子算符在福克空间中的表达式，其一般形式为：

H = Σ_{pq} h_{pq} a†_p a_q
  + ½ Σ_{pqrs} h_{pqrs} a†_p a†_q a_r a_s + E_nuc

其中 h_{pq} 是单电子积分，h_{pqrs} 是双电子积分，a† 和 a 分别是产生与湮灭算符。这个表达式看起来简洁，但它的规模完全由所选轨道集合决定。轨道集合越大，非零积分越多，费米子算符项数越多，后续 Jordan-Wigner 或 Bravyi-Kitaev 映射产生的量子比特数和门数也随之膨胀。因此，分子轨道选择、活性空间和基组大小，是资源估算中最上游、也最具杠杆效应的三个参数。

轨道选择：从全空间到有效模型

一个含 M 个空间轨道的体系，在自旋轨道层面有 2M 个模式。若不做任何裁剪，双电子积分张量 h_{pqrs} 的名义规模是 (2M)^4，即使利用对称性，实际非零项仍随 M 快速增长。以水分子为例，STO-3G 基组下氧的 1s、2s、2p 和两个氢的 1s 共 7 个空间轨道，自旋轨道 14 个；而 cc-pVDZ 下空间轨道增至 24 个，自旋轨道 48 个。仅这一变化，双电子积分的名义数量就从约 3.8×10^4 量级跳到约 5.3×10^6 量级。

但并非所有轨道都同等重要。内层轨道能量极低，电子几乎不参与成键与反应；高能虚轨道在基态中占据数极小。真正决定化学性质的是费米能级附近的价层轨道。因此，轨道选择的核心操作是：先做 Hartree-Fock 计算得到分子轨道及其能量，再按能量排序和占据数分析，把轨道划分为冻结核心、活性占据、活性虚轨道和冻结虚轨道四类。

活性空间：CAS(n, m) 的参数含义

活性空间通常记作 CAS(n, m)，表示 n 个电子在 m 个空间轨道中完全相关。冻结核心轨道不参与相关计算，只贡献一个常数能量和有效核势修正；冻结虚轨道完全不进入哈密顿量。这样，费米子算符只作用在 2m 个自旋轨道上，双电子积分规模降为 (2m)^4 量级。

以水分子为例，若选择 CAS(4, 4)，即 4 个电子在 4 个空间轨道中，自旋轨道数为 8。相比全空间 14 个自旋轨道，双电子积分名义项数从约 3.8×10^4 降到约 4.1×10^3，减少近一个数量级。若进一步选择 CAS(2, 2)，自旋轨道数为 4，积分项数降到约 256，但代价是丢失了大部分动态相关能，总能量误差可能达到数十毫哈特里。

基组大小：从 STO-3G 到 cc-pVTZ 的连锁反应

基组决定了空间轨道数 M。STO-3G 是最小基组，每个原子只用最少函数描述；cc-pVDZ、cc-pVTZ 则逐级增加极化函数和弥散函数。对同一分子，M 随基组增大而增大，而费米子算符项数大致按 M^4 增长，映射后的量子比特数按 2M 线性增长，但门数受算符项数影响更接近高阶增长。

下表给出一个示意性的对比，帮助理解三者的连锁影响。表中数值为按常见对称性约化后的量级估计，用于教学说明，不是某次具体计算的实测结果。

可以看到，基组从 STO-3G 换到 cc-pVDZ，全空间量子比特数从 14 增到 48，积分项数增加两个数量级以上。但如果同时施加 CAS(4,4) 活性空间，量子比特数又回到 8，积分项数回到 4.1×10³。这说明活性空间是控制资源的关键阀门，而基组主要影响活性空间之外的环境描述精度。

从费米子算符到量子线路的映射开销

费米子算符不能直接作用在量子比特上，必须经过映射。Jordan-Wigner 映射把每个自旋轨道对应一个量子比特，但产生湮灭算符会带上 Z 字符串，导致一个双电子项展开成多个泡利字符串。Bravyi-Kitaev 映射通过二叉树结构把 Z 字符串长度从 O(N) 降到 O(log N)，但增加了辅助比特和线路深度。无论哪种映射，映射后的泡利项数都正比于费米子算符项数，而费米子算符项数又正比于活性空间自旋轨道数的四次方。

因此，资源估算的链条可以概括为：基组决定 M，活性空间决定 m，m 决定费米子算符项数，映射方式决定泡利项数和量子比特数，最终决定量子线路的门数和深度。下面这张图展示了这一传导关系。

图中最值得关注的是中间那一层：基组和轨道选择先汇聚成“空间轨道数 M”，再经过活性空间裁剪变成“自旋轨道数 2m”，最后决定积分张量规模和费米子算符项数。很多初学者会误以为基组越大越好，却忽略了基组增大后如果不加活性空间约束，费米子算符项数会以四次方速度膨胀，映射后的量子线路深度可能超出当前硬件的相干时间。

工程上的取舍策略

在实际资源估算中，建议按以下顺序操作。第一步，用 STO-3G 或 def2-SVP 做 Hartree-Fock，得到轨道能量和占据数，确定哪些是核心、哪些是价层。第二步，选定一个初始活性空间，例如 CAS(4,4) 或 CAS(6,6)，计算费米子算符项数和映射后的泡利项数。第三步，逐步扩大活性空间或基组，观察能量误差和资源开销的变化曲线，找到误差下降趋缓而资源仍可接受的拐点。

一个实用的经验是：对于轻元素分子的基态性质，CAS(4,4) 到 CAS(8,8) 通常能捕获大部分静态相关能；若需要更高精度，再考虑增大基组而不是盲目扩大活性空间。因为基组增大主要改善动态相关，而动态相关在量子线路中往往需要更多项才能体现，资源代价更高。

总结来说，二次量子化后的哈密顿量规模由轨道集合直接决定。分子轨道选择决定了哪些轨道进入计算，活性空间决定了费米子算符的实际作用范围，基组大小决定了轨道总数和积分张量的上限。三者共同影响费米子算符项数，进而通过映射影响量子比特数、泡利项数和线路深度。资源估算的第一步，就是把这几个参数明确下来，并理解它们之间的四次方传导关系。