手机版
           

DFT密度泛函理论从基础到实践——一个计算材料研究者的十年认知

发布时间:2026-06-13   来源:科研学术网    
字号:

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实战的经验分享和参数对照案例。

图说天下

×
cp2k计算
dft计算
Gaussian计算
MS计算
VASP计算