手机版

密度泛函理论DFT计算入门：从Hohenberg-Kohn定理到VASP实操的完整路径

发布时间：2026-06-15 来源：科研学术网

小中大

字号：

密度泛函理论DFT计算的入门门槛不在数学推导，而在理解”每一步近似引入了什么误差”。Hohenberg-Kohn定理告诉你电子密度唯一确定了体系的基态性质，Kohn-Sham方程把多体问题映射到单体问题——但这两个优雅的理论框架在走向实际计算的过程中，要做一系列物理近似，每一个近似都有适用边界。

从多体问题到Kohn-Sham方程：近似的链条

薛定谔方程对N个电子体系的精确求解，计算量随电子数呈指数增长——这就是所谓的”多体问题灾难”。DFT的核心洞察是：基态总能量是电子密度的泛函——不需要知道多电子波函数的全部信息，只需要电子密度就够了。

但Hohenberg-Kohn定理只是存在性证明，它告诉你”存在这样一个泛函”，但没有给出具体的泛函形式。Kohn-Sham引入了一套无相互作用的辅助粒子系统，通过求解单粒子Kohn-Sham方程来得到电子密度。代价是引入了交换关联泛函Exc[n]——这个泛函的精确形式未知，必须做近似。

从精确到近似的链条是这样的：精确的薛定谔方程 → Hohenberg-Kohn定理（精确） → Kohn-Sham方程（引入近似泛函） → LDA/GGA（局域/半局域近似） → 杂化泛函（部分精确交换）。每往下一层，计算量增加，精度提升——但没有哪一层是”完美的”。

泛函层级：LDA、GGA、杂化泛函的取舍

LDA（局域密度近似）假设每一点的交换关联能等于均匀电子气在该密度的值。这对电荷密度变化缓慢的体系（简单金属）效果不错，但对电荷密度变化剧烈的体系（分子、表面、带隙材料）表现较差。

GGA（广义梯度近似）在LDA基础上加入了密度的梯度修正。PBE是目前最常用的GGA泛函，对大部分材料的结构优化和能量计算给出了合理的结果。但PBE有三个已知缺陷：低估带隙（0.3-1.0 eV的系统性低估）、对强关联体系（d/f电子）描述不足、对弱相互作用（vdW力）不敏感。

杂化泛函（HSE06、PBE0）通过混合Hartree-Fock的精确交换来弥补GGA的不足。HSE06用25%的精确交换，并且对长程交换做了屏蔽处理，计算成本约为PBE的10-15倍。它对带隙的修正效果显著——TiO₂的带隙从PBE的2.2 eV修正到HSE06的3.4 eV，实验值3.2 eV，偏差只有6%。

选哪个泛函，取决于你关注什么物理量。如果只做结构优化，PBE就够了；如果要看带隙和光学性质，HSE06几乎是必须的；如果体系有强关联电子（NiO、MnO），还需要加DFT+U。

VASP实操：INCAR参数的物理意义

理解了泛函层级之后，VASP的INCAR参数设置就不只是一份参数清单，而是对物理理解的体现。

ENCUT（截断能）：决定了平面波基组的完备性。截断能越高，波函数展开越精确，但计算成本约∝E^3。每个新体系的第一步测试就是ENCUT收敛性——从350 eV开始，每次加50 eV，看总能量变化<2 meV/atom时的ENCUT就是收敛值。一般来说，PAW势对ENCUT的要求比USPP高，ENCUT通常在400-500 eV。

k点密度：决定了Brillouin zone的采样密度。金属体系的k点密度需要比绝缘体更高（因为费米面的积分需要更密的采样）。一般用Monkhorst-Pack网格，经验规则是每个方向至少5个点（小cell）到1个点（大supercell）。

ISMEAR和SIGMA：控制费米面展宽方法。绝缘体和半导体必须用ISMEAR=0（Gaussian展宽）或-5（四面体法，仅静态计算），金属可以用ISMEAR=1（Methfessel-Paxton）。如果对绝缘体用了ISMEAR=1，能带在带边附近会出现人为的虚假展宽，直接导致带隙被低估。

EDIFF和EDIFFG：电子步收敛阈值和离子步收敛阈值。EDIFF=1E-5是大部分体系的默认值，但如果要做高精度能量差计算（如形成能、缺陷能），建议收紧到1E-6。EDIFFG控制离子步的结构优化判据，一般设-0.01到-0.02 eV/Å。