2. 量子化学基础与二次量子化
要理解量子化学计算与量子算法之间的关系，必须先建立一套统一的“语言”：分子体系如何被近似、电子状态如何被表示、能量如何被计算。本节围绕四个核心概念展开——Born-Oppenheimer近似、分子轨道理论、电子相关、二次量子化与Fock空间。它们不是彼此独立的工具，而是一条从物理近似到可计算表示的链条。
一、Born-Oppenheimer近似：把核与电子的运动分开
分子中原子核的质量比电子大三个数量级以上，因此电子运动远快于核运动。Born-Oppenheimer近似（BO近似）据此把总波函数写成电子波函数与核波函数的乘积，并先固定核坐标求解电子问题。此时电子哈密顿量（原子单位制）为：
H_el = -Σ_i (1/2)∇_i² - Σ_{i,A} Z_A / |r_i - R_A| + Σ_{i<j} 1 / |r_i - r_j|其中第一项是电子动能，第二项是电子与核的吸引，第三项是电子间排斥。核坐标 R_A 在此作为参数出现，电子能量 E_el(R) 随核构型变化，形成势能面。BO近似是绝大多数量子化学计算的起点，它把“求解分子”转化为“在给定核构型下求解电子问题”。
二、分子轨道理论：从原子轨道到分子轨道
分子轨道理论（MO）用原子轨道的线性组合（LCAO）构造分子轨道：
φ_i = Σ_μ C_{μi} χ_μ其中 χ_μ 是原子轨道基函数，C_{μi} 是组合系数。以H₂为例，两个1s轨道可组合成成键轨道 σ = (χ_A + χ_B)/√2 和反键轨道 σ* = (χ_A − χ_B)/√2。成键轨道能量低于两个孤立原子轨道，反键轨道能量升高。电子按泡利原理和能量最低原则填充，形成电子组态。
在Hartree-Fock（HF）框架下，每个电子在其余电子平均场中运动，分子轨道通过自洽场（SCF）迭代求解。HF给出的是单行列式波函数，是后续相关能计算的参考态。需要强调：HF能量并不是精确能量，它忽略了电子瞬时相关。
三、电子相关：HF之后缺失的那部分能量
电子相关能定义为精确非相对论能量与HF极限能量之差：
E_corr = E_exact − E_HF相关能通常为负值，量级在每对电子约0.02–0.05 Hartree。它来源于电子间的瞬时排斥，使电子彼此“避让”。相关能可分为动态相关和静态相关：动态相关来自电子短程避让，静态相关出现在近简并轨道（如键断裂、过渡金属配合物）中，单行列式描述失效。
处理相关能的方法包括：组态相互作用（CI）、耦合簇（CCSD、CCSD(T)）、多体微扰（MP2）、密度泛函（DFT，通过交换关联泛函近似）。以H₂O为例，HF/STO-3G能量约为−74.96 Hartree，加入MP2相关修正后约降低0.02–0.03 Hartree。这些数值属于示例量级，用于说明相关能不可忽略。
四、二次量子化与Fock空间：把电子状态写成算符作用
一次量子化的波函数在粒子数变化时不方便，二次量子化引入产生算符 a†_p 和湮灭算符 a_p，作用在Fock空间的态上。Fock空间是所有粒子数扇区的直和：
F = ⊕_{N=0}^{∞} H_N其中 H_N 是N电子希尔伯特空间。真空态 |vac⟩ 满足 a_p|vac⟩=0。单电子态写作 a†_p|vac⟩，多电子态写作 a†_p a†_q … |vac⟩。反对易关系为：
{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其中 h_{pq} 是单电子积分，h_{pqrs} 是双电子积分。这一形式直接对应量子算法中的费米子哈密顿量，是VQE、量子相位估计等算法的输入。
五、从概念到计算：一条可操作的路径
把上述概念串起来，一个典型的量子化学计算流程如下：
- 选择核构型与基组：例如H₂O采用STO-3G基组，核坐标取实验几何。
- BO近似下求解HF方程：SCF迭代至能量收敛，得到分子轨道系数与轨道能。
- 构造二次量子化哈密顿量：计算单电子积分 h_pq 和双电子积分 h_pqrs，写成费米子算符形式。
- 选择相关处理方法：小体系可用FCI或CCSD(T)作为基准；量子算法则用VQE或QPE。
- 映射到量子比特：通过Jordan-Wigner或Bravyi-Kitaev变换，将费米子算符转为泡利算符。
- 测量与判断：比较计算能量与化学精度（约1.6 mHa），判断是否达到目标。
六、常见误区与检查点
- BO近似不是永远成立：在锥形交叉、非绝热过程附近，核与电子运动耦合显著，需超越BO近似。
- 基组不是越大越好：基组需与相关方法匹配，过大基组在小体系上可能浪费资源，过小基组则引入基组误差。
- Fock空间不是“无限大”的借口：实际计算中必须截断到有限轨道数，截断方式直接影响精度。
- 二次量子化不改变物理：它只是换了一种表示，使粒子数变化和算符代数更自然，便于量子算法实现。
掌握这条从BO近似到Fock空间的链条，后续学习VQE、量子相位估计、费米子映射时，就能清楚每个算符和参数的物理来源，而不是把它们当作黑箱公式。