说明:本文华算科技主要介绍密度泛函理论中电子密度、Kohn-Sham 映射、交换相关泛函和材料计算结果的使用边界。


密度泛函理论为什么从电子密度开始
DFT 的起点不是逐个追踪每个电子的完整波函数,而是用三维空间中的电子密度 n(r) 描述体系基态。多电子波函数需要 3N 个空间坐标,电子数增加时,问题规模会急剧增加;电子密度只依赖三维坐标,仍然保留原子、键、表面和缺陷周围电子分布的核心信息。
Hohenberg-Kohn 定理奠定了关键理论基础:在给定原子核外势的基态问题中,电子密度可以唯一确定体系的基态能量和其他基态性质。DFT 的核心逻辑是由原子位置构建外势,再由外势确定基态密度,最后由密度计算得到总能量 E[n]。这条映射关系让材料计算可以绕开完整多体波函数。
外势主要由原子核位置、元素类型和赝势决定。晶胞、表面 slab、缺陷超胞和吸附构型一旦确定,程序首先构建电子感受到的外部环境,再寻找与这个环境匹配的基态密度。DFT 将结构问题转化为外势和密度问题,再由能量泛函判断哪个结构或电子分布更有利。

电子密度不是简单的电子云图像。在材料计算里,它参与构建总能、静电势、交换相关势、力和应力,也决定后续的能带、态密度、电荷分析等结果。n(r) 是 DFT 的核心变量,也是许多图谱和能量差的共同来源,因此理解 DFT 首先要理解密度的来源。

密度变化通常对应成键、极化和电荷重新分布,但不能直接等同于价态、活性或稳定性。比如残余密度或差分密度能显示空间分布变化,定量电荷转移还需要 Bader 电荷、积分区域或其他布居方法辅助;能量判断还需要参考态和总能量差。


Kohn-Sham 方程到底在算什么
Kohn-Sham 方法引入一个辅助的非相互作用电子体系,让它产生和真实相互作用电子体系相同的基态密度。这样做并非假设电子真的不存在相互作用,而是将难处理的多体效应纳入交换相关泛函 Exc[n],同时保留可求解的单电子方程形式。
一次 DFT 计算通常从原子位置和初始密度初猜开始,构造有效势,求解 Kohn-Sham 方程,得到轨道和占据数,再更新电子密度。新的密度用于重新构建有效势,计算过程反复迭代直至自洽。自洽之后,程序输出总能量 Etot、原子受力、应力、能级、波函数和电荷密度等结果。

自洽的目的不是让某个数值在表面上收敛,而是使密度、有效势和能量达到自洽一致。收敛标准控制着密度、能量和力的收敛精度,因此电子步收敛不充分会影响后续离子步优化,离子步未收敛也会使总能量差、吸附能和声子计算失去可靠的基准。
Kohn-Sham 轨道和能级具有重要的计算参考价值,但必须明确其物理意义。Kohn-Sham 轨道主要是计算密度和能量的工具,不能把每一个本征值都直接等同于实验激发能。带隙、光吸收、激子、强关联激发等问题,通常需要 HSE、GW、BSE、TDDFT 或更合适的模型补充。

原子结构通过外势和哈密顿量影响电子结构计算。晶格常数、缺陷位置、吸附构型、磁态、U 值、SOC 和 vdW 修正都会改变有效势或能量泛函的组成。DFT 计算结果不是脱离模型的材料真相,而是在给定结构、边界条件和近似泛函下得到的基态电子结构与能量。


交换相关泛函为何决定 DFT 的边界
精确交换相关泛函未知,这是 DFT 精度边界的核心来源。Kohn-Sham 框架将电子间复杂的交换、相关和部分动能差异纳入 Exc[n]。LDA 只使用局域密度,GGA 加入密度梯度,meta-GGA 继续引入动能密度等信息,杂化泛函再加入一定比例的精确交换。
泛函不是可有可无的参数。它会改变平衡结构、键长、相对能量、带隙、磁矩、吸附能、缺陷形成能和反应自由能。同一个原子模型换用不同泛函,收敛密度和总能量都可能改变,这也是文献中常见的 PBE、PBEsol、SCAN、HSE06、PBE+U 等方法差异的来源。

泛函阶梯可以帮助理解不同泛函的精度差异,但不能作为简单排序。LDA、GGA、meta-GGA 和杂化泛函各自使用的变量不同,计算代价和适用场景也不同。更高阶泛函的结果不代表对所有材料都更准确,尤其在强关联、弱相互作用、金属表面和带隙问题中,误差抵消和目标性质会同时影响选择。
交换相关势 vxc(r) 不是简单常数修正,它随局域密度、密度变化和材料环境改变。低密度区、原子附近、强离域或强局域电子区域,对泛函形式更敏感。材料体系一旦包含过渡金属 d 电子、稀土 f 电子、层间弱相互作用或电荷转移界面,泛函选择会显著影响计算结果。

所以,DFT 计算结果需要明确说明所用的泛函、赝势、基组、k 点、截断能、色散修正、磁性设置和收敛标准。可比较性来自同一套参考条件,而不是只比较最后的 eV 数值。当研究目标换成带隙、反应路径、吸附强度或缺陷电荷态时,可靠的计算设置也要随目标性质调整。


材料计算里该怎样使用 DFT 结果
DFT 最常用的直接结果是基态总能量、电子密度和原子受力。结构优化通过力的最小化实现,形成能、吸附能、表面能和反应能均来自总能量差,能带和 DOS 来自 Kohn-Sham 本征值与态密度投影,差分电荷密度、Bader 电荷和功函数则继续使用收敛密度或静电势。

在材料筛选中,DFT 常作为统一参考数据源。研究者可以用同一套计算协议计算得到能量、带隙、力学常数、吸附能或反应自由能,再建立材料库和机器学习模型。高通量 DFT 的价值在于能够提供一致的趋势对比,而不在于每个单点数值都与实验值完全一致。
许多常见物理量其实是 DFT 输出的二次处理结果。形成能来自化学势参考态,吸附能来自表面和吸附物参考态,反应自由能还要加入零点能、熵和电位修正。总能量差需要先定义初态、终态和归一化方式,否则同样是 eV,比较的科学对象可能完全不同。

使用 DFT 时首先要明确科学问题属于基态、近基态还是激发态。晶体结构、相对稳定性、吸附构型和许多反应热力学问题适合从 DFT 起步;光学激发、强关联绝缘体、长时间动力学、溶剂界面和有限温熵效应,则需要额外方法或修正。明确这个边界,才能避免 DFT 结果的混用。
DFT 的核心可以概括为一句话:在给定原子结构和近似泛函下,DFT 通过自洽求解电子密度计算得到基态能量、力和电子结构。它提供了材料计算的基础数据框架;后续的吸附能、形成能、能带、态密度、电荷分析和反应自由能,都要回到这个框架里理解。
