2. 分子电子结构问题建模

在把分子基态能量问题交给 VQE 之前，必须先把它从化学语言翻译成量子计算可执行的数学形式。这一步通常称为分子电子结构建模，它决定了后续量子线路要处理多少个量子比特、哈密顿量有多少项、以及最终能量精度能到哪一档。本节沿着“分离核与电子运动 → 写出二次量子化哈密顿量 → 选择基组 → 划分活性空间”的顺序展开，每一步都直接影响 VQE 的资源开销与可解性。

一、Born-Oppenheimer 近似：先把核的运动冻结

分子中原子核质量比电子大三个数量级以上，电子运动的时间尺度远快于核振动。Born-Oppenheimer 近似据此把总波函数写成电子波函数与核波函数的乘积，并在求解电子问题时把核坐标视为固定参数。于是分子电子哈密顿量（原子单位制）为：

H_elec = -Σ_i (1/2)∇_i²  - Σ_{i,A} Z_A / |r_i - R_A|
         + Σ_{i<j} 1/|r_i - r_j|  + Σ_{A<B} Z_A Z_B / |R_A - R_B|

前三项分别是电子动能、电子与核的库仑吸引、电子间库仑排斥；最后一项是核-核排斥能，在固定核坐标下只是一个常数。对 VQE 而言，这一常数不参与量子线路，只需在测量得到电子能量后加上即可。因此“分子基态能量”实际求解的是电子哈密顿量的最低本征值，再加上核排斥常数。

二、二次量子化形式：从波函数到费米子算符

直接处理多电子波函数的反对称性很麻烦。二次量子化把电子态表示为占据数态，用产生算符 a†_p 和湮灭算符 a_p 描述电子在自旋轨道 p 上的产生与消失，并自动满足反对易关系 {a_p, a†_q} = δ_pq。分子电子哈密顿量写成：

H = Σ_{pq} h_pq a†_p a_q
  + (1/2) Σ_{pqrs} h_pqrs a†_p a†_q a_s a_r + E_nuc

其中 h_pq 是单电子积分（动能加核吸引），h_pqrs 是双电子积分（电子间排斥）。单电子项对应电子在轨道间的跃迁，双电子项对应两电子相互作用。这个形式的好处是：只要选定一组自旋轨道，所有积分都可以由经典量子化学程序算出，哈密顿量就变成一组费米子算符的线性组合。

量子计算机不能直接执行费米子算符，需要映射为泡利算符。常用 Jordan-Wigner 变换把每个自旋轨道映射为一个量子比特，占据为 |1⟩、空为 |0⟩，并引入 Z 串保证反对易性。映射后哈密顿量变成泡利串之和：

H = Σ_k c_k P_k,  P_k ∈ {I, X, Y, Z}^⊗n

项数随轨道数增长，例如 4 个自旋轨道（2 个空间轨道）在 Jordan-Wigner 下通常产生几十到上百个泡利项。VQE 的测量次数与泡利项数直接相关，因此控制轨道数是工程上的关键。

三、基组选择：精度与量子比特数的权衡

基组决定用多少空间轨道展开分子轨道。最小基组 STO-3G 每个原子只用价层轨道，例如 H₂ 用 2 个空间轨道、4 个自旋轨道，对应 4 个量子比特；6-31G 会翻倍甚至更多。基组越大，电子相关能描述越准，但量子比特数和哈密顿量项数迅速上升。

实际选择时通常遵循：先用 STO-3G 验证线路与流程，再逐步升级到 6-31G 或 cc-pVDZ。对 VQE 教学与原型验证，STO-3G 已能给出定性正确的键长与能量趋势；若追求化学精度（约 1 kcal/mol），则需要更大基组并配合活性空间方法。

四、活性空间划分：把“重要电子”留给量子处理器

即使基组不大，全空间量子比特数也可能超出当前硬件。活性空间方法把轨道分成三部分：

• 冻结核心轨道：内层电子几乎不参与成键，能量贡献近似常数，可经典处理。
• 活性轨道：价层中参与成键、反键和电子相关的轨道，交给量子线路处理。
• 冻结虚轨道：高能空轨道对基态贡献小，可忽略或经典微扰处理。

常用记号 CAS(n, m) 表示 n 个电子在 m 个空间轨道中完全活跃。例如 H₂ 在 STO-3G 下是 CAS(2, 2)，即 2 个电子、2 个空间轨道，对应 4 个自旋轨道、4 个量子比特；若把 1s 核心冻结，则活性空间进一步缩小。活性空间划分直接决定量子比特数、哈密顿量项数和 VQE 可达到的精度。

五、从建模到 VQE 的衔接

完成上述四步后，VQE 的输入就明确了：一个由泡利串表示的哈密顿量 H = Σ c_k P_k，以及一个作用在 n 个量子比特上的参数化试探态 |ψ(θ)⟩。VQE 通过测量各泡利串的期望值求和得到能量 E(θ) = ⟨ψ(θ)|H|ψ(θ)⟩，再用经典优化器调整 θ 使 E 最小。此时能量下限由活性空间和基组共同决定，线路深度由试探态 ansatz 决定，测量次数由泡利项数决定。

一个常见的工程判断是：如果活性空间选得太小，即使 VQE 优化到收敛，能量也会偏离真实值；如果选得太大，量子比特数和测量次数会超出硬件能力。因此建模阶段需要在“化学精度”和“量子资源”之间做显式权衡，而不是把全部轨道都交给量子处理器。

总结本节的关键映射关系：Born-Oppenheimer 近似把核运动冻结，得到电子哈密顿量加常数 E_nuc；二次量子化把电子问题写成费米子算符；基组决定空间轨道数；活性空间划分决定实际参与量子计算的轨道和电子数；Jordan-Wigner 等映射把自旋轨道变成量子比特，把哈密顿量变成泡利串之和。这四步共同决定了 VQE 的量子比特数、哈密顿量项数和可达到的精度，是后续线路设计与测量方案的基础。