4.1 基组选择：平面波、局域原子轨道与实空间网格

4.1 基组选择：平面波、局域原子轨道与实空间网格

本节摘要：KS 方程的波函数要在计算机里表示，必须用一组基函数展开——这就是基组。平面波基组天然适配周期性体系、收敛系统；局域原子轨道基组紧凑、直观、适合分子；实空间网格基组概念简单、适合时间演化。本节讲清三类基组的原理、优缺点与 BSSE 问题，并给出选基组的决策框架。

学习目标

阅读完本节，你应当能够：

1. 解释为什么需要基组：把无限维波函数空间投影到有限维
2. 说出平面波基组的系统收敛性（截断能 E_cut）与无 BSSE 的优点
3. 理解 LCAO 基组的 zeta 层级、极化函数与弥散函数的作用
4. 解释 BSSE 的来源与反作用校正（Counterpoise）思路
5. 知道实空间网格基组的适用场景（TDDFT、复杂边界）
6. 根据体系类型选择基组：周期性→平面波，分子→LCAO

问题与直觉

电子的波函数理论上是个连续函数，在无穷多维的函数空间里。但计算机只能存有限个数字。怎么把一个连续函数「装进」计算机？

答案就是基组：选一组已知的、彼此独立的函数 \{\phi_i\}，用它们的线性组合去逼近真实波函数：

Ψ(r) ≈ Σ_i c_i φ_i(r)
​

系数 c_i 是计算机里真正存的东西。基组选择就像选「调色板」：函数选得好，少量系数就能描述复杂波函数；选得不好，就得堆很多系数。

不同研究领域自然而然地发展出了不同的「调色板」：

• 做固体的人面对的是周期体系，波函数天然周期，平面波就是天作之合
• 做分子的人面对的是局域在原子周围的电子云，原子轨道（局域函数）更合适
• 做时间演化的人需要直接在实空间追踪电子运动，网格更顺手

三条路各有胜负手，没有哪条绝对更好。理解它们的权衡，是配置计算的第一步。

还需要澄清一个常见误解：「基组越大越好」。基组变大的确逼近「完备极限」，但收益递减——从双 Zeta 到三 Zeta 的提升明显，从三 Zeta 到四 Zeta 就小了，代价却一路涨。工程上「够用就行」的基组选择，背后是对「误差预算」的规划：基组误差要和泛函误差、数值误差大致匹配，单方面砸基组纯属浪费算力。真正专业的做法是先做一次收敛扫描，再选「平台区起点」的基组，把省下的算力花在别处。

核心原理

2.1 平面波基组：周期世界的傅里叶之舞

平面波基组把波函数写成平面波叠加：

Ψ_k(r) = Σ_G c_G e^(i(k+G)·r)
​

k 是布里渊区波矢，G 是倒格矢。这个形式天然满足布洛赫定理——周期势场中的电子波函数就该是这个形状。对晶体、表面这类周期性体系，它是「原生的」基组。

系统收敛性是它最大的优点：只需要一个参数——截断能 E_{cut}（满足 |\mathbf{k}+\mathbf{G}|^2/2m \le E_{cut} 的平面波才保留）——就能系统性地提高精度。截断能翻倍，基函数变多，结果单调收敛。这个「一个旋钮控制精度」的特性让结果极具可比性。

无 BSSE是另一个大优点：平面波均匀铺满全空间，分子不会「借用」邻居的基函数，基组叠加误差（BSSE）天然不存在。这对算弱相互作用很重要。

代价：原子核附近波函数剧烈振荡，要很多平面波才能描述。所以平面波几乎总是配赝势或 PAW（4.2 节）——只处理价电子，把核附近「糊掉」。

2.2 局域原子轨道基组（LCAO）

LCAO 的哲学是「以小见大」：分子轨道 ≈ 原子轨道的线性组合：

Ψ(r) = Σ_A Σ_μ c_Aμ φ_Aμ(r)
​

基函数以原子核为中心。历史上先有斯莱特型轨道（STO，指数衰减，像真实原子轨道），后有高斯型轨道（GTO，两个高斯乘积还是高斯，积分高效），现代程序几乎都用 GTO 的组合。

基组有清晰的「升级路径」：

• 最小基组：每个价轨道一个函数，如 STO-3G，定性用
• 双 Zeta（DZ）：每个价轨道两个函数，如 6-31G，允许轨道在不同环境伸缩
• 三 Zeta（TZ）：每个价轨道三个函数，如 6-311G，更高精度
• 极化函数：加高角动量函数（如碳上加 d、氢上加 p），描述键合时电子云变形
• 弥散函数：加「胖」的高斯，描述阴离子、弱束缚电子、激发态

LCAO 的最大优点是紧凑与化学直观：基函数少（比平面波少几个量级），能做全电子计算（不依赖赝势），轨道图像直接对应化学家的直觉。最大缺点有二：BSSE（两个分子靠近时互相「借用」基函数，虚假降低能量）和非系统收敛（没有像截断能那样的单一旋钮，只能一级一级换更大的基组）。

关于「非系统收敛」多说一句：这并不意味着 LCAO 精度不可控，而是说它的收敛路径是「离散台阶」——从 DZ 到 TZ 再到 QZ，每一步都跳变，没有连续旋钮。实践中人们用「外推」技巧：算 DZ、TZ、QZ 三个点，外推到完备基极限。这套「基组外推」是高精度分子计算的标准操作，也是「系统收敛」的一种补偿方案。

2.3 BSSE 与反作用校正

BSSE 的机制：A 分子和 B 分子靠近时，A 的电子可以「借用」B 的基函数来描述自己，能量被虚假拉低——相互作用能因此被算得更负（更稳定）。就像两个学生共享答案，分数虚高。

最常用的修正是反作用校正（Counterpoise）：算 A 和 B 各自的能量时，把对方的「幽灵轨道」（空基函数，无电子）也加上。这样 A 在复合物里能用的基函数集合，和单独算时一致，BSSE 被系统性扣除。

E_bind = E_AB(AB基组) - E_A(AB基组) - E_B(AB基组)
​

BSSE 的严重程度有个规律值得记住：基组越小，BSSE 越严重。因为小基组下每个分子都「饿着」，更急于借用邻居的函数。反之基组越大越接近完备，借用的动机就越小。所以用大基组（如带弥散的三 Zeta）本身就能大幅缓解 BSSE，Counterpoise 是对小基组/弱相互作用的兜底手段。工程上两条路都走：用尽量完备的基组，同时对弱作用体系做 Counterpoise 复核。

2.4 实空间网格基组

实空间方法不走「基函数展开」的路，直接把空间切成网格，在网格点上定义波函数/密度的数值。动能算符用有限差分或 FFT 处理。

动能算符作用于 Ψ ≈ [Ψ(x+Δx) - 2Ψ(x) + Ψ(x-Δx)] / Δx²
​

优点：概念简单、无 BSSE、边界条件灵活（周期性/开放都行）、系统收敛（减小 Δx）。特别适合实时 TDDFT——直接在实空间追踪电子在光场下的演化，是模拟超快过程、强场物理的理想工具。

缺点：网格点随 Δx 缩小立方增长（三维），内存需求巨大；非局域势（如 HF 交换）在实空间处理困难；描述芯电子同样需要赝势。

2.5 一个直观的「网格间距」例子

实空间网格的收敛逻辑值得用数字感受一下：把网格间距 Δx 减半，一维上点数翻倍，三维就是八倍。所以「更精确的网格」不是线性贵，是立方贵。这也解释了为什么实空间方法的内存开销常成为瓶颈——它把「计算精度」直接兑换成了「存储量」。工程上做实空间计算，第一步永远是确定「多大的 Δx 能让能量收敛到我要的精度」，而不是盲目加密。这个「找最粗但足够」的哲学，其实对所有基组都适用。

工程实践要点

1. 周期性体系（晶体、表面、二维材料）→ 平面波 + PAW：VASP、QE 的标准路线。
2. 孤立分子、团簇 → LCAO：Gaussian、ORCA 的标准路线。大分子时 LCAO 紧凑性优势明显。
3. 时间演化、复杂边界 → 实空间：OCTOPUS、GPAW 等。做 TDDFT 动力学、纳米结构时考虑。
4. 收敛测试是铁律：无论哪种基组，都必须做收敛测试——平面波扫截断能、LCAO 逐级升基组、网格减小间距，直到关键量（能量、结构、带隙）稳定。
5. 弱相互作用一定要处理 BSSE：LCAO 算分子间作用时，Counterpoise 校正几乎是必须的。

三类基组对比矩阵

FAQ：几个高频疑问

问：6-31G(d) 和 cc-pVTZ 这种名字怎么读？

答：Pople 系（6-31G 等）和 Dunning 相关一致系（cc-pVDZ/TZ 等）是两套主流命名。以 6-31G(d) 为例：「6-31」表示芯函数用 6 个高斯、价层用 3+1 分裂，「(d)」表示加 d 极化函数。cc-pVTZ 表示相关一致三 Zeta（带极化）。读法不重要，重要的是知道「越大越贵越准」以及「对什么体系用多大合适」——经验值可以查同类文献。

问：为什么平面波不用做 BSSE 校正，LCAO 却要？

答：BSSE 的前提是「基函数局域在原子/分子周围，邻居可以借用」。平面波铺满全空间，不存在「谁借谁的」——任何位置都由同一组平面波描述。LCAO 的基函数锚在原子核上，两个分子靠近时基函数空间重叠，借用发生了。这是两类基组结构性差异的直接后果。

问：可不可以平面波和 LCAO 混用？

答：可以，而且混合基组（原子核附近 LCAO、远处平面波）是活跃的研究方向。工程上更常见的是「分步混用」：固体用 VASP（平面波）优化结构，截取片段用 Gaussian（LCAO）做高精度化学计算。第 5 章催化、生物案例会反复出现这种组合。

⚠️ 常见坑：平面波算孤立分子，结果收敛极慢又浪费资源——周期性基组描述真空里一个分子，等于在倒空间铺了一层用不上的平面波。要么用 LCAO，要么用大真空盒配截断能测试。

💡 关键直觉：基组选择本质是「精度—成本—直观性」三角权衡。平面波胜在系统收敛与无 BSSE，LCAO 胜在紧凑与直观，实空间胜在边界灵活与时间演化。选基组前先想清楚「我算什么体系、需要什么性质」，而不是「哪个基组最流行」。

概念性伪代码：基组收敛测试

def basis_convergence_test(basis_series, system):
    results = []
    for basis in basis_series:      # 如截断能 300→400→500→600
        energy = run_scf(system, basis)
        results.append((basis, energy))
    # 找到能量进入平台区的拐点，用最小的基组
    return plateau_point(results)
​

本节要点回顾

• 基组作用：把无限维波函数投影到有限维，系数是真正的存储对象
• 平面波：周期体系原生、系统收敛（截断能）、无 BSSE，需配赝势
• LCAO：紧凑直观、能做全电子，有 BSSE、收敛非系统
• STO vs GTO：STO 像真轨道但积分难，GTO 积分高效成主流
• 基组升级路径：最小→DZ→TZ→+极化→+弥散
• BSSE：分子互相借用基函数虚低能量，Counterpoise 修正
• 实空间网格：概念简单、无 BSSE、边界灵活，适合 TDDFT，内存大
• 误差预算：基组误差要与泛函、数值误差匹配，别单方面砸基组
• 收敛测试：无论哪种基组，都必须验证关键量稳定
• 选型铁律：周期→平面波，分子→LCAO，时间演化→实空间