4. 基组与轨道表示

在量子化学计算中，哈密顿量最初定义在连续的三维空间里，电子坐标可以取任意实数值。要把这样一个算符交给计算机处理，必须先把连续空间“投影”到一个有限维的函数空间上——这就是基组的作用。基组的选择直接决定了哈密顿量矩阵的维度、矩阵元的计算方式，以及最终能量的精度。本节围绕四个核心概念展开：Slater轨道、高斯型轨道、收缩基组，以及基组重叠与正交化。

一、从连续空间到有限基组

假设我们选取一组已知的、固定的单电子函数 {χ_μ(r)}，称为基函数。分子轨道 φ_i 用它们的线性组合表示：

φ_i(r) = Σ_μ C_{μi} χ_μ(r)

把这个展开式代入电子哈密顿量，并利用变分原理，连续空间的微分方程就转化为矩阵本征值问题：

H C = S C ε

其中 H 是哈密顿量矩阵，S 是重叠矩阵，C 是分子轨道系数矩阵，ε 是轨道能。这一步就是“投影”：连续算符被离散化为有限维矩阵。基函数选得越接近真实电子分布，投影误差越小。

二、Slater轨道：物理上正确的形状

Slater型轨道（STO）的形式为：

χ_STO(r) = N · r^(n-1) · exp(-ζ r) · Y_lm(θ, φ)

它的关键特征是 exp(-ζr) 型径向衰减。对于氢原子1s轨道，精确解正是这种指数形式，因此STO在原子核附近有正确的尖点（cusp），在远距离也有正确的指数尾巴。用STO做基函数，通常只需要很少几个函数就能得到不错的原子能量。

但STO有一个致命缺点：当两个STO分别位于不同原子中心时，它们之间的双电子积分没有解析公式，必须做数值积分。对于含有几十个原子的分子，积分数量随基函数数目四次方增长，数值积分会变得极其昂贵。这就是为什么现代量子化学程序几乎不用纯STO做分子计算。

三、高斯型轨道：用解析性换效率

高斯型轨道（GTO）把指数部分换成 exp(-α r²)：

χ_GTO(r) = N · x^a y^b z^c · exp(-α r²)

GTO最大的优势是“高斯乘积定理”：两个不同中心的高斯函数相乘，可以解析地合并成一个位于新中心的高斯函数。这使得所有双电子积分都能写成解析表达式，计算速度比STO快几个数量级。

代价是形状不对。GTO在核附近是平的（导数为零），没有尖点；远距离衰减太快，像 exp(-αr²) 而不是 exp(-ζr)。单个GTO对真实原子轨道的描述很差。解决办法是用多个GTO线性组合去拟合一个STO的形状，这就是收缩基组的思想。

四、收缩基组：精度与效率的折中

收缩基组的基本做法是：选定一组原高斯函数（primitive GTO），用固定的线性组合系数把它们“收缩”成一个基函数。例如STO-3G基组，就是用3个高斯函数拟合一个Slater轨道。收缩系数和指数在基组定义文件中预先给定，不参与SCF迭代优化。

以氢原子的1s为例，STO-3G的收缩形式为：

χ_1s = 0.154329 · g(3.42525091) 
     + 0.535328 · g(0.62391373) 
     + 0.444635 · g(0.16885540)

三个高斯函数的指数从大到小排列，大指数负责核附近区域，小指数负责尾部区域，组合起来逼近 exp(-ζr) 的形状。收缩基组把“形状拟合”和“分子轨道优化”分成两步：前者用固定系数完成，后者只优化C矩阵。这样既保留了GTO的解析积分优势，又改善了单电子函数的物理形状。

常见的收缩策略包括：

• 最小基组：每个原子轨道只用一个收缩函数，如STO-3G。
• 分裂价基组：价层轨道用两个或更多收缩函数，如3-21G、6-31G，给价电子更多灵活性。
• 极化基组：在分裂价基础上加入更高角动量的函数，如6-31G(d)给重原子加d函数。
• 弥散基组：加入指数很小的函数，描述远离核的电子密度，如6-31+G(d)。

五、基组重叠与正交化

当基函数位于不同原子中心时，它们之间一般不正交。重叠矩阵S的定义为：

S_{μν} = ∫ χ_μ*(r) χ_ν(r) dr

对角元S_μμ=1（归一化），非对角元S_μν衡量两个基函数在空间中的重叠程度。如果两个基函数相距很远，S_μν接近0；如果它们在同一原子附近，S_μν可能达到0.5甚至更高。

重叠矩阵的存在使广义本征值问题 H C = S C ε 不能直接对角化。标准做法是先把基组正交化。常用的正交化方案有两种：

• 对称正交化：令 X = S^-1/2，变换后的基组保持最大相似性。S^-1/2 通过S的本征分解得到：S = U s U^T，则 S^-1/2 = U s^-1/2 U^T。
• 正则正交化：令 X = U s^-1/2，变换后基函数按重叠矩阵本征值排序，数值上更稳定。

正交化之后，广义本征值问题化为标准形式：

H' C' = C' ε,  其中 H' = X^T H X,  C = X C'

这一步在SCF迭代之前完成一次即可。需要注意的是，如果基组中存在近线性相关的函数（例如两个指数非常接近的高斯函数），S矩阵会出现很小的本征值，S^-1/2会放大数值误差。实际程序中通常会设定一个阈值（如10^-6），丢弃对应的小本征值方向，保证数值稳定性。

六、从基组到哈密顿量矩阵的完整链路

把以上环节串起来，哈密顿量从连续空间到可计算矩阵的流程是：

1. 选择基组类型（STO-3G、6-31G(d)等），确定每个原子的收缩函数集合。
2. 计算所有基函数之间的单电子积分：重叠积分S_μν、动能积分T_μν、核吸引积分V_μν。
3. 计算双电子积分 (μν|λσ)，这是计算量最大的部分。
4. 组装哈密顿量核心矩阵 H^core = T + V。
5. 构造重叠矩阵S，做正交化得到变换矩阵X。
6. 在SCF迭代中，用密度矩阵P构建Fock矩阵F = H^core + G(P)，再变换为F' = X^T F X，对角化得到轨道能和新系数。

基组的选择贯穿整个流程：它决定了积分计算的难度、矩阵的维度、S矩阵的条件数，以及最终能量的系统误差。一个经验是：对于平衡几何和常规反应能，6-31G(d)通常能给出定性正确的结果；如果需要定量精度（误差小于1 kcal/mol），则需要加到6-311++G(2df,2pd)或更大。基组收敛性测试——用同一方法、不同基组计算同一性质，观察结果随基组增大的变化——是判断基组是否足够的最直接手段。