第一性原理基础：薛定谔方程、密度泛函理论（DFT）、交换关联泛函简介
各位同学，今天咱们来聊聊第一性原理计算的根基。说白了，就是怎么从最基本的物理定律出发，算出材料的性质。我刚开始接触这个领域时，也觉得这些理论离实际应用太远。但做了十几年项目，我越来越觉得，理解这些基础，就像盖楼打地基——地基不牢，楼盖得再高也悬。
1. 薛定谔方程：一切计算的起点
先问大家一个问题：为什么我们需要第一性原理？因为实验太贵、太慢。你想想看，一个新材料从设计到合成，可能要花几个月甚至几年。而计算，只要一台服务器，几天就能给出初步判断。
薛定谔方程，就是描述微观粒子行为的“宪法”。它的形式很简单：
ĤΨ = EΨ其中 Ĥ 是哈密顿算符，Ψ 是波函数，E 是能量。但别被这个简洁的形式骗了。对于多电子体系，这个方程复杂到几乎无法精确求解。
我记得刚读博时，导师让我算一个包含50个原子的体系。我天真地以为直接解薛定谔方程就行。结果算了三天三夜，服务器直接崩溃了。后来我才明白，对于实际材料，我们根本不可能精确求解薛定谔方程——这就是所谓的“指数墙”问题。
核心难点： 电子数量每增加一个，波函数的变量就增加3个（三维空间坐标）。对于N个电子，波函数是3N维的。N=10时，计算量就已经大到无法承受。
2. 密度泛函理论（DFT）：化繁为简的智慧
既然直接求解波函数行不通，那怎么办？密度泛函理论给出了一个巧妙的思路：用电子密度 ρ(r) 代替波函数 Ψ。
为什么可以这样？因为电子密度只依赖于三维空间坐标，不管你有多少个电子。这就把3N维的问题，降到了3维。你想想看，这节省了多少计算量！
Hohenberg-Kohn 定理告诉我们：基态能量是电子密度的唯一泛函。也就是说，只要知道电子密度，就能确定体系的所有性质。这个定理在1964年提出，可以说是DFT的“出生证明”。
但这里有个坑：Hohenberg-Kohn 定理只证明了存在性，没告诉我们怎么构造这个泛函。这就引出了Kohn-Sham方程。
[-½∇² + V_ext(r) + V_H(r) + V_xc(r)] φ_i(r) = ε_i φ_i(r)这个方程把多电子问题，转化成了单电子在有效势场中的运动。我个人习惯把这个过程叫做“拆解”——把复杂的多体问题，拆成一个个独立的单电子问题，然后自洽求解。
我的经验： 做DFT计算时，一定要检查自洽收敛。我曾经因为收敛阈值设得太松，算出来的形成能差了0.3 eV。这个误差在缺陷计算中，足以让你得出完全相反的结论。
3. 交换关联泛函：LDA、GGA、杂化泛函
Kohn-Sham方程里有个神秘项——V_xc(r)，也就是交换关联势。它包含了所有电子-电子相互作用的复杂部分。问题是，我们不知道它的精确形式。所以，只能近似。
这就好比我们知道目的地，但地图不完整。不同的近似，就是不同精度的地图。
3.1 LDA：简单但有效
局域密度近似（LDA）是最早的近似。它假设：在空间某一点的交换关联能，只取决于该点的电子密度。就像你只看脚下的一小块地，不管周围是什么。
LDA的优点：计算快，对某些体系（如金属）结果不错。缺点：对带隙的估计严重偏低，对氢键等弱相互作用描述很差。
我记得有一次做钙钛矿材料的缺陷计算，用LDA算出来的带隙只有0.8 eV，而实验值是1.5 eV。这个误差直接导致缺陷能级的位置完全不对。嗯，从那以后，我对LDA的带隙结果就多留了个心眼。
3.2 GGA：更精细的地图
广义梯度近似（GGA）在LDA基础上，加入了电子密度的梯度信息。也就是说，它不光看脚下，还看周围的地势变化。常用的有PBE、PW91等泛函。
GGA比LDA更准确，尤其是对分子和半导体。但它仍然低估带隙，只是程度稍轻。
| 泛函类型 | 精度 | 计算成本 | 适用场景 |
|---|---|---|---|
| LDA | 一般 | 低 | 金属、简单体系 |
| GGA-PBE | 较好 | 中等 | 半导体、表面、缺陷 |
| 杂化泛函 | 高 | 高 | 带隙、缺陷能级 |
3.3 杂化泛函：精度与成本的权衡
杂化泛函的思路很直接：把DFT的交换能和Hartree-Fock的精确交换能混合起来。常用的有HSE06、PBE0等。
为什么这样做？因为DFT低估交换能，而Hartree-Fock高估它。两者混合，正好取长补短。
我建议做缺陷计算时，如果条件允许，尽量用HSE06。虽然计算量是GPA的5-10倍，但带隙和缺陷能级的精度明显提升。我曾经对比过PBE和HSE06算出的氧空位形成能，差了将近1 eV。这个差距，足以决定一个缺陷是否稳定。
避坑指南： 我曾经用PBE算一个宽禁带半导体的缺陷，结果发现形成能是负的——这意味着缺陷会自发形成。但换了HSE06后，形成能变成了正的。所以，对于带隙大于2 eV的材料，千万别迷信GGA的结果。
4. 知识体系总览
下面这张图，是我自己总结的DFT计算知识框架。你可以把它当作一张“地图”，随时回来对照。
5. 总结与建议
好了，这一章的内容就到这里。总结一下：
- 薛定谔方程是理论基础，但实际计算中必须做近似
- DFT通过电子密度简化了问题，是当前第一性原理计算的主流方法
- 交换关联泛函的选择直接影响计算结果，需要根据体系特点权衡精度和成本
我个人建议：做缺陷形成能计算时，先用GGA-PBE做结构优化（因为快），然后用HSE06算单点能（因为准）。这样既保证了效率，又兼顾了精度。
下一章，我们会具体讲怎么搭建计算模型、设置参数。到时候我会分享一些实际项目中的“血泪教训”，保证让你少走弯路。
记住： 没有完美的泛函，只有合适的泛函。理解每个近似的假设和局限，比盲目追求“最精确”的方法更重要。
公众号:蓝海资料掘金营，微信deep3321