第4章 · 计算实现与方法
文档摘要
第 4 章 · 计算实现与方法 章节摘要：理论再好，落不到代码上就是空中楼阁。本章回答「DFT 到底怎么算」：波函数用什么基组展开（平面波/局域原子轨道/实空间网格）、芯电子怎么处理（赝势与 PAW）、自洽场 SCF 如何迭代到收敛、矩阵对角化与几何优化有哪些数值算法，以及 VASP、Quantum ESPRESSO、Gaussian、ORCA 这些软件各自的长处。读完本章，你应当能独立设计并跑通一轮基本的 DFT 计算。
第 4 章 · 计算实现与方法
章节摘要：理论再好，落不到代码上就是空中楼阁。本章回答「DFT 到底怎么算」：波函数用什么基组展开（平面波/局域原子轨道/实空间网格）、芯电子怎么处理（赝势与 PAW）、自洽场 SCF 如何迭代到收敛、矩阵对角化与几何优化有哪些数值算法，以及 VASP、Quantum ESPRESSO、Gaussian、ORCA 这些软件各自的长处。读完本章，你应当能独立设计并跑通一轮基本的 DFT 计算。
学习目标
阅读完本章，你应当能够：
- 对比平面波、局域原子轨道、实空间网格三类基组的收敛性、BSSE 与适用体系
- 说出赝势的「冻结芯电子」思想、范数守恒与超软赝势的区别，以及 PAW 如何兼顾全电子精度与平面波效率
- 复述 SCF 迭代六步：初猜、构建哈密顿量、对角化、更新密度、混合、收敛判断
- 区分简单混合与 Pulay/DIIS 混合，知道金属体系收敛困难的应对（smearing 等）
- 说明满矩阵对角化与迭代求解器（CG、Davidson、LOBPCG）的复杂度差异
- 根据「分子还是周期性体系」在 VASP/QE/ABINIT/Gaussian/ORCA 之间选型
核心概念速览
金句：DFT 计算是一次「参数联调」——基组、赝势、k 点、收敛阈值互相牵制，哪一环没收敛测试，结果就不可信。
计算实现五环节架构
子章节导航
4.1 基组选择
波函数的离散化思想、平面波基组的截断能与系统性收敛、LCAO 的 zeta 层级与极化/弥散函数、实空间网格与 TDDFT 的配合，以及 BSSE 与反作用校正。
4.2 赝势与 PAW 方法
全电子计算的瓶颈、赝势的构造原理与分类（范数守恒/超软）、PAW 的线性变换与原子球重构思想、VASP 等软件中的势库选择。
4.3 自洽场（SCF）迭代过程
SCF 的「鸡生蛋」困境、六步迭代流程、密度混合策略、收敛判据与常见收敛失败的对策。
4.4 数值求解与优化算法
KS 方程的矩阵化、满矩阵对角化与迭代求解器、几何优化中的 BFGS/L-BFGS、过渡态搜索的 NEB、从头算分子动力学与并行计算。
4.5 主流计算软件包：VASP 等
VASP 的 PAW 平面波与商业授权、QE 的开源模块化、ABINIT 的教学友好、Gaussian 的分子化学全栈、ORCA 的大分子效率，以及选型决策逻辑。
子章节之间的逻辑关系
4.1 基组 ──► 4.2 赝势/PAW ──► 4.3 SCF ──► 4.4 数值算法 ──► 4.5 软件 │ │ │ │ │ 用什么展开 芯电子怎么处理 怎么迭代到自洽 怎么算得快 怎么组织参数
五节基本是「从下往上搭积木」：基组和赝势决定方程的离散形式，SCF 决定怎么解方程，数值算法决定解多快，软件则把这些选择打包成输入文件。理解了前四节，第 5 节的软件选型自然顺理成章。
前置知识与后续延伸
- 前置：第 2 章的 KS 方程与 SCF 雏形；第 3 章的泛函选择（泛函决定部分计算参数）
- 为后续铺垫：第 5 章所有应用都建立在「会跑计算」的基础上；第 6 章 TDDFT、DFT+U 都要在特定软件里打开对应开关，本章的软件地图是前提
每节对应的「动手任务」
| 章节 | 最小动手任务 | 验证标准 |
|---|---|---|
| 4.1 基组 | 对一个简单晶体逐步提高截断能 | 总能量收敛到目标精度 |
| 4.2 赝势/PAW | 换用不同势库重复同一次计算 | 结果差异在可接受范围 |
| 4.3 SCF | 观察迭代日志中的能量序列 | 能量单调逼近并达到收敛 |
| 4.4 数值算法 | 用几何优化找最小构型 | 受力收敛到阈值以下 |
| 4.5 软件 | 用目标软件跑通一个最小算例 | 输出可复现、格式正确 |
常见误区
- 「截断能随便设大一点就行」——设太大会浪费计算资源，设太小结果不可信。正确做法是递增测试，找到能量随截断能变化曲线进入平台区的拐点。
- 「赝势选官方库就万无一失」——官方库质量高，但首次用于特定体系时仍应验证：算该元素的晶格常数、体模量，与实验或高精度结果对比。
- 「SCF 不收敛就加大迭代次数」——根源多半是混合策略不当或初猜太差。先换 Pulay/DIIS、降混合因子、尝试不同初猜，比盲目加次数更有效。
一句话记忆
本章的五个环节可以浓缩成一句实操口诀：「基组选型看体系，赝势验证看精度，SCF 收敛靠混合，数值加速靠算法，软件选型看授权与生态」。把这五个短语记住，跑计算时每个环节该检查什么，心里就有数了。
计算成本的一个现实视角
普通 DFT 计算的复杂度通常是体系原子数的三次方量级（O(N^3)），杂化泛函因为引入精确交换会升到四次方。这意味着：体系翻一倍，普通计算慢八倍，杂化泛函慢十六倍。这也是为什么大体系只能用 GGA、小体系才能上杂化或 RPA 的根本原因。本章 4.4 节的迭代求解器、4.5 节的 GPU 并行，都是在跟这个标度搏斗。读第 5 章应用时会反复碰到「成本决定方法」这条线，现在先有个量级感——它会在你日后每次设计计算方案时跳出来提醒你。若条件允许，建议本章读完后立刻在目标软件里跑一遍「硅晶体基态计算」这类最小算例，把四个输入文件的关键参数对照正文逐项核对，实操一遍胜过读十遍。读完本章，你应当能独立设计并跑通一轮基本的 DFT 计算。