3. DFT总能量计算基础

缺陷形成能本质上是一个能量差：把缺陷引入完美晶体后，体系总能量相对参考态的变化。因此，任何影响总能量数值精度的设置，都会直接传递到形成能、电荷态跃迁能级和离化能的计算结果中。本节围绕平面波基组、赝势、交换关联泛函、截断能和k点采样这五个核心要素，说明它们各自在算什么、如何影响缺陷能量，以及实际计算中应如何选择和验证。

平面波基组：把波函数展开成一组“频率”

在周期性晶体中，Kohn-Sham 轨道满足 Bloch 定理，可写成平面波与周期函数的乘积。平面波基组就是把波函数用一组倒格矢 G 对应的平面波展开：

ψ_{n,k}(r) = Σ_G c_{n,k}(G) · exp[i(k+G)·r]

这组基函数天然满足周期性边界条件，不需要人为构造原子轨道，也不存在基组重叠积分问题。代价是：要描述原子核附近剧烈振荡的波函数，需要大量高频率平面波，计算量随截断能迅速上升。这正是赝势存在的理由。

赝势：把芯电子“藏起来”，只算价电子

内层芯电子几乎不参与成键，但其波函数在核附近剧烈振荡，直接展开代价极高。赝势（或投影缀加波 PAW 方法）用一个平滑的等效势替代芯电子与核的联合作用，使价电子波函数在核附近变得平缓，从而可以用较小的截断能描述。常见选择包括超软赝势（USPP）和 PAW，后者在缺陷计算中更常用，因为它对磁性和电荷转移的描述更稳定。

赝势的关键影响在于：不同赝势的“硬度”不同，所需截断能不同；同一元素的不同赝势可能给出略有差异的形成能。因此，缺陷计算中应固定一套赝势，并在方法部分明确说明版本和价电子配置。

交换关联泛函：决定能量精度的“物理近似”

交换关联泛函是 DFT 中唯一需要近似的部分。缺陷计算中最常用的是 GGA-PBE，它计算成本低、对结构弛豫较可靠，但会低估带隙，导致缺陷能级位置和电荷态跃迁能级出现系统偏差。对于强关联体系（如含过渡金属氧化物），可考虑 DFT+U 或杂化泛函 HSE06。杂化泛函能显著改善带隙和局域缺陷态描述，但计算量大幅增加，通常用于对关键缺陷做最终修正，而非全部筛选。

截断能与k点收敛：两个必须验证的数值参数

截断能 E_cut 决定平面波基组的完备程度。E_cut 太低，总能量未收敛，形成能随 E_cut 漂移；E_cut 太高，计算量浪费。k 点采样决定布里渊区积分的精度。对于缺陷计算，超胞越大，k 点可以越稀疏；但缺陷态往往在带边附近局域，k 点太少会导致占据数、磁矩和总能量不稳定。

实际做法是：先对完美超胞做截断能和 k 点收敛测试，观察总能量随参数的变化。通常以总能量变化小于 1 meV/atom 作为收敛判据。缺陷形成能对这两个参数更敏感，因为它是两个大能量之差，误差可能放大。建议在最终形成能计算中，将截断能比收敛值再提高 10%–20%，k 点网格至少覆盖缺陷超胞的倒空间。

参数选择与验证的实际流程

下面给出一套可直接参考的收敛测试与缺陷计算参数设置流程，适用于常见的半导体或氧化物缺陷研究。

1. 确定赝势与泛函：选择 PAW-PBE 赝势，记录版本号与价电子配置。若体系含 d 电子，评估是否需要 DFT+U。
2. 构建完美超胞：例如 2×2×2 的 64 原子超胞，先做结构弛豫，得到平衡晶格常数。
3. 截断能收敛测试：固定 k 点为 4×4×4，令 E_cut 从 300 eV 逐步增至 600 eV，记录总能量。当相邻步能量差小于 1 meV/atom 时，取该值为收敛截断能，再上浮 10%–20%。
4. k 点收敛测试：固定收敛后的 E_cut，令 k 网格从 2×2×2 增至 6×6×6，观察总能量与带隙变化。缺陷超胞通常取 2×2×2 或 3×3×3 即可。
5. 缺陷计算：在收敛参数下，分别计算带电缺陷和中性缺陷的总能量，代入形成能公式。注意修正背景电荷和势能对齐。
6. 结果复核：检查磁矩、缺陷态在带隙中的位置、原子位移是否合理；若形成能随 k 点或 E_cut 仍有明显漂移，需回到第 3 步重新测试。

常见参数设置参考

总结来说，平面波基组和赝势决定了“能不能算得动”，交换关联泛函决定了“算得准不准”，截断能和 k 点决定了“数值是否收敛”。在缺陷形成能计算中，这五个要素不是独立的：赝势越硬，所需截断能越高；泛函越精确，所需 k 点往往也越密。实际工作中，先用 PBE 加中等参数做快速筛选，再对关键缺陷用 HSE06 和更密 k 点做最终确认，是兼顾效率与精度的常见策略。