计算DFT方法详解：密度泛函理论的计算流程与参数选择策略

密度泛函理论（DFT）将多电子波函数的高维问题映射为电子密度的三维问题，用Kohn-Sham方程替代了无法精确求解的薛定谔方程。计算DFT的核心挑战不在于理解理论本身，而在于在有限的计算资源下做出一系列相互牵制的参数选择。

Kohn-Sham方程将多电子体系的基态性质表达为单电子有效势场中的轨道方程，其中唯一无法精确知道的部分是交换关联能泛函。局域密度近似（LDA）在原子核附近的快速密度变化区域失效严重。广义梯度近似（GGA）通过引入密度梯度修正，在分子键长和晶体晶格常数的预测上取得了显著改进，但仍然系统性地低估半导体和绝缘体的带隙。

研究者在选择交换关联泛函时，往往是在一组系统性的偏差中挑选对当前体系伤害最小的那个。PBE泛函对晶格常数的预测通常比实验值大1%到2%，在比较不同材料的性质趋势时往往可以抵消，但在计算高压相变压力时会导致数个GPa的偏移。SCAN泛函作为元GGA家族的代表，对键长、键角和内聚能的预测精度逼近了杂化泛函的水平，代价是每一自洽步的迭代次数增加约30%到50%。

平面波基组展开是DFT计算中连接实空间与倒易空间的桥梁。截断能决定了平面波基组的数目，直接控制计算精度与成本的权衡。对于包含第三周期元素的体系，400 eV到500 eV的截断能通常能够保证总能量的收敛精度达到1 meV/atom的量级。当体系中包含第一系列过渡金属时，建议将ENCUT设置在500 eV到600 eV之间。这个项目在优化一种铁基钙钛矿的计算参数时，发现ENCUT从500 eV提升到700 eV的过程中，磁矩收敛值从3.2 μB跃升至3.8 μB，这个变化直接改写了整个磁性相图的计算结果。

k点采样密度决定了倒易空间中布里渊区积分的精细程度。对于绝缘体和半导体，设置4×4×4到6×6×6的k点网格通常足够。金属材料因为费米面处态密度的不连续性，需要显著更密的k点采样，8×8×8以上的网格才能将总能量的误差控制在1 meV/atom以内。

DFT计算的心脏是自洽场（SCF）迭代循环。当体系包含部分占据的窄带时，电荷密度的振荡可能导致迭代发散。引入Gaussian smearing通过在费米面附近引入可控的电子态展宽，能够显著改善收敛性，代价是引入了一个需要最终外推至零的物理参数。

电子步收敛标准（EDIFF）的设定需要在数值精度与计算时间之间仔细权衡。设置为1e-6 eV时，每一自洽步需要15到25次子迭代；放宽为1e-4 eV时，子迭代次数降至8到12次。对于单纯的晶格常数优化，1e-4 eV的收敛标准通常足够；对于需要计算力学性质的二阶导数问题，建议将EDIFF收紧至1e-5 eV或更高。离子步收敛标准（EDIFFG）通过原子受力来判断结构优化是否终止，对于后续的声子谱或弹性常数计算，EDIFFG = -0.01 eV/Å是最低要求。

传统的半局域泛函在描述长程色散相互作用时存在根本性缺陷。石墨的层间距在PBE计算中比实验值大约扩大了20%，直接导致层间结合能被低估了一个数量级。为了解决这个问题，过去十五年中发展出了多个范德华修正方案，从最简单的DFT-D两两加和修正，到更为精密的TS方案，再到完全非局域的vdW-DF系列泛函。

TS方案通过在每一个原子对之间引入与极化率相关的C6系数，能够以极小的额外计算成本显著改善分子晶体和有机框架材料的结构预测精度。vdW-DF系列泛函通过在交换关联能中引入显式的非局域关联项，在原理上更为严格，但计算成本通常是PBE的3到5倍，对于包含数百个原子的体系构成实质性的资源壁垒。

密度泛函理论的根本局限性源自交换关联泛函的近似性质。即使是当前最精密的杂化泛函，在描述强关联体系和电荷转移激发时仍然存在系统性偏差。对于这类体系，DFT+U方法通过在特定的局域轨道上引入一个额外的库仑排斥项U，能够以较小的额外成本改善强关联材料的电子结构描述。U值的确定本身就是一个开放的研究问题，不同的确定方案可能导致电子结构的定性改变。

尽管存在这些局限性，DFT计算在材料科学中的位置仍然不可替代。它为分子动力学模拟提供了第一性原理力场参数，为热力学模型提供了形成能和相图数据。在精心选择泛函和充分验证收敛性的前提下，DFT对晶体结构的预测精度能够达到实验分辨率的极限，对弹性性质和热力学稳定性的预测与高压实验的吻合度在近年来持续提升。