1. DFT基础概念解析
密度泛函理论(Density Functional Theory, DFT)是量子力学计算建模方法，广泛应用于物理、化学和材料科学领域。它通过电子密度而非波函数来描述多粒子系统，大大简化了计算复杂度。DFT的核心思想是：系统的所有基态性质都可以由电子密度n(r)唯一确定。
1.1 理论框架与基本方程
DFT的数学基础建立在Hohenberg-Kohn定理上：
- 定理1：外势(因而总能量)是电子密度的唯一泛函
- 定理2：给出基态能量的泛函在输入密度为真实基态密度时取得最小值
Kohn-Sham方程是DFT的核心计算工具：
code复制[-ħ²/2m ∇² + V_s(r)]φ_i(r) = ε_i φ_i(r)
其中有效单粒子势V_s(r)包含三部分：
- 外势V(r)
- Hartree项(描述电子-电子库仑排斥)
- 交换关联势V_XC[n(r)]
1.2 交换关联泛函
精确的交换关联泛函未知，常用近似包括：
- 局域密度近似(LDA)：仅依赖坐标处的密度
- 广义梯度近似(GGA)：考虑密度梯度
- 杂化泛函：结合精确交换和DFT关联
2. DFT计算流程详解
2.1 自洽场迭代过程
- 初始猜测电子密度n(r)
- 构建Kohn-Sham势V_s(r)
- 求解Kohn-Sham方程得轨道φ_i
- 计算新密度n'(r)=Σ|φ_i(r)|²
- 混合新旧密度，检查收敛
- 未收敛则返回步骤2
实际计算中常采用Pulay混合等加速收敛技术，混合参数通常取0.1-0.3
2.2 基组选择策略
| 基组类型 | 特点 | 适用场景 |
|---|---|---|
| 平面波 | 收敛性好，需赝势 | 固体计算 |
| 高斯型 | 计算效率高 | 分子体系 |
| 数值原子轨道 | 局域性好 | 大体系计算 |
3. 实用计算技巧与参数设置
3.1 关键参数优化
-
截断能(平面波)：
- 一般材料：400-600 eV
- 含d电子体系：≥800 eV
- 收敛测试应使总能变化<1 meV/atom
-
k点网格：
- 金属：至少12×12×12
- 半导体