DFT密度泛函理论在计算材料领域几乎无处不在。十年前刚接触时,把Kohn-Sham方程的推导从头推了一遍,觉得”把多体问题约化成单电子在有效势中运动”这套逻辑很漂亮。但后面真正的认知转折都不是在推导里发生的,而是在跑实际计算时被数据和预期之间的差距逼出来的。

Hohenberg-Kohn定理告诉我们什么——以及没告诉我们什么
DFT的数学基础是1964年的两个定理:第一定理说基态电子密度唯一决定了外势(从而决定了哈密顿量),第二定理说能量作为密度的泛函在基态密度处取极小。这两个定理从原理上证明了”用电子密度代替多体波函数做基本变量”是严格可行的。
但定理本身不提供泛函的具体形式。动能泛函T[ρ]不知道怎么写,交换关联能E_xc[ρ]也不知道怎么写。Kohn和Sham在1965年的聪明做法是:把动能中”无相互作用参考系”的部分T_s[ρ]单独拿出来精确计算,剩下的所有不知道怎么写的东西全部塞进E_xc[ρ]——包括真实动能和T_s[ρ]的差、以及所有多体交换关联效应。
这个拆解的精妙之处在于:它把DFT的核心困难全部集中到了一个量(E_xc)上,然后允许用各种近似去逼近它。它没解决困难本身,但把困难包装成了一个可以逐级逼近的模块。接下来几十年DFT发展的主线,说到底就是不断给E_xc找更好的近似。
交换关联泛函的四层阶梯
第一层是LDA(局域密度近似)。假设空间中每一点的交换关联能只依赖于该点的电子密度,等于均匀电子气的值。从头算的固体物理里,LDA对简单金属和半导体的晶格常数预测准确到百分之几,对很多体系的结合能也不错。但LDA低估带隙、高估结合能的倾向是系统性的,对含d/f电子的体系问题更突出。
第二层是GGA(广义梯度近似)。在LDA基础上引入密度梯度修正,常见的有PW91和PBE两种形式。PBE目前是计算材料领域的默认选择。相比LDA,PBE改善了结合能(从高估10-20%优化到高估5-10%)、对氢键和分子晶体的描述也合理了一些。但带隙预测依然系统偏低——这是DFT本身的问题,不是GGA特有的。
第三层是meta-GGA。在密度梯度之上再引入动能密度或密度的二阶导数。SCAN泛函(2015年)是这层的代表作,对很多体系的精度接近杂化泛函,但计算成本只比GGA略高。SCAN对水的结构预测、材料的形成能、以及部分强关联体系的带隙都有明显改善。
第四层是杂化泛函。把一部分精确的Hartree-Fock交换能和GGA交换能混合,HSE06是最常见的版本(混合25%的精确交换,并在短程做屏蔽)。HSE06显著改善了带隙的预测精度,但计算成本比GGA高一个数量级,对超过100个原子的体系用起来吃力。
四层阶梯对应着精度和算力之间的不同折中。实践中没有”最好”的泛函,只有”在给定精度要求和算力约束下最合适”的泛函。
实践中经常被忽略的DFT局限
DFT能做的事情很多,但有三类问题即使最好的泛函也解决不了:
第一,强关联体系。过渡金属氧化物、稀土化合物中d/f电子的局域化效应,LDA和GGA通常给出金属态而实验上是绝缘体。DFT+U可以部分补救——人为加一个Hubbard U项惩罚d/f电子的非物理离域——但U值本身的选取有任意性,不同U值给出的结论可能相反。
第二,范德华体系。标准LDA和GGA不包含非局域的长程色散相互作用。石墨、MoS₂这类层状材料的层间距和层间结合能,不加范德华修正会算错一个数量级。DFT-D3和TS方法是两种流行的经验修正方案,但都引入了新的可调参数。
第三,激发态和带隙问题。DFT本质上是基态理论,Kohn-Sham本征值在严格意义上不等于准粒子能量,所以”带隙”也是没有严格理论保证的量。GW近似和多体微扰理论可以修正这个问题,但计算成本和复杂度跳了一个台阶。
了解这些局限不是为了否定DFT——恰恰相反,知道DFT在哪里会出错,才能避开把错误的结果当真。
从使用者视角看DFT的工具链
不管用VASP、Quantum ESPRESSO还是CP2K,DFT计算的决策链大致相同:选泛函→设ENCUT和k点→跑结构优化→跑静态计算→根据需要做后处理(态密度、能带、Bader电荷、弹性常数等)。
这个链条上每个环节的决策都会影响最终结果。选PBE还是HSE06?ENCUT设400还是520 eV?k点密度对应每原子多少个k点?有没有加自旋极化?有没有做偶极修正?——每个问题背后都有量级在0.01到0.5 eV之间的误差,累加起来就是文章和文章之间”结论不同”的来源。
十年下来最大的体会是:DFT密度泛函理论不是一个”输入结构、输出正确结果”的黑箱,而是一套需要逐个参数验证、逐项误差评估的计算框架。项目的质量不取决于用了多高端的版本或多大的超算,而取决于参数选择背后的物理依据是否经得起推敲。这是在任何环节都无法跳过的。在科研学术网首页上能看到更多关于DFT实战的经验分享和参数对照案例。
材料模拟计算方法全景——从量子力学到有限元的多尺度策略
DFT密度泛函理论从基础到实践——一个计算材料研究者的十年认知
GITT法计算锂离子扩散系数的实验设计与数据处理要点
金属原子间键能计算的方法选择与精度控制——从DFT到EAM势
密度泛函理论DFT计算在催化材料筛选中的实践经验
第一性原理计算层错能:超胞建多大、层错怎么引入、弛豫要不要开
DFT计算表面能:切面、弛豫、化学势——三条线交叉决定表面能算出来是0.8还是1.2 J/m²
DFT计算材料相容性:界面的电子结构不匹配,宏观相容性无从谈起