📚 课程目录第4章超胞模型构建与缺陷浓度

4. 超胞模型构建与缺陷浓度

在缺陷形成能计算中，超胞（supercell）是承载缺陷的周期性重复单元。它决定了缺陷浓度、缺陷间距以及周期性镜像相互作用的强弱。超胞尺寸选择不当，会直接导致形成能、电荷态跃迁能级和缺陷浓度预测出现系统性偏差。本节围绕“多大超胞才够用”这一核心问题，梳理尺寸选择、间距控制、镜像相互作用与有限尺寸修正的完整逻辑。

一、超胞尺寸与缺陷浓度的定量关系

设原胞体积为 Ω₀，沿三个方向分别扩胞 n₁、n₂、n₃，则超胞体积 Ω = n₁n₂n₃Ω₀，超胞内原子数为 N = n₁n₂n₃N₀。若每个超胞中放置一个缺陷，则缺陷浓度（原子分数）为：

c_def = 1 / N = 1 / (n₁·n₂·n₃·N₀)

示例：Si 原胞 N₀ = 2
  2×2×2 超胞：N = 16，c_def ≈ 6.25 at.%
  3×3×3 超胞：N = 54，c_def ≈ 1.85 at.%
  4×4×4 超胞：N = 128，c_def ≈ 0.78 at.%
  5×5×5 超胞：N = 250，c_def ≈ 0.40 at.%

可见，超胞越大，缺陷浓度越低，越接近“孤立缺陷”极限。但计算量随 N 近似线性增长，因此需要在精度与成本之间取平衡。经验上，对于共价半导体中的中性缺陷，3×3×3 或 4×4×4 超胞（约 64～128 原子）常能给出收敛到 0.1 eV 以内的形成能；对于带电缺陷或金属体系，往往需要更大超胞。

二、缺陷间距与周期性镜像相互作用

超胞的周期性意味着缺陷在空间中按晶格平移重复排列。相邻缺陷之间的最小距离 d 由超胞尺寸和缺陷位置共同决定。以立方超胞边长 L 为例，若缺陷位于超胞中心，则最近镜像距离约为 L；若缺陷偏离中心，最近镜像距离可能缩短到 L/2 甚至更小。

镜像相互作用主要来自三个方面：

弹性相互作用：缺陷引起的晶格畸变在周期边界上叠加，导致应力场不能完全弛豫。
静电相互作用：带电缺陷之间存在库仑排斥或吸引，在周期性边界下形成无限点电荷晶格，能量发散，需要背景电荷补偿。
能带与波函数重叠：缺陷态波函数在相邻镜像间交叠，使缺陷能级展宽、跃迁能级发生移动。

这些相互作用会使形成能随超胞尺寸变化。通常，带电缺陷的形成能对超胞尺寸更敏感，因为库仑项按 1/(εL) 衰减，收敛较慢；中性缺陷的弹性项衰减更快，一般按 1/L³ 量级。

三、有限尺寸修正策略

当超胞无法无限增大时，需要引入有限尺寸修正（finite-size correction），把超胞计算得到的形成能外推到孤立缺陷极限。常用策略包括：

外推法：对一系列超胞尺寸计算形成能，按 1/L 或 1/L³ 拟合，外推到 L→∞。对带电缺陷，常用 1/L 线性外推；对中性缺陷，常用 1/L³ 外推。
Makov-Payne 修正：针对带电缺陷的库仑相互作用，扣除周期性点电荷晶格的马德隆能。修正项为 q²α/(2εL)，其中 α 为马德隆常数，ε 为介电常数。该修正可显著改善带电缺陷形成能的收敛性。
势对齐修正：带电缺陷会改变超胞的平均静电势，需将缺陷超胞与完美超胞的静电势对齐，否则形成能会出现额外偏移。常用方法是对比远离缺陷区域的静电势平台。
背景电荷补偿：在周期性边界下，带电超胞需要均匀背景电荷中和，否则能量发散。计算中应确保总电荷为零，并在后处理中扣除背景电荷贡献。

实用建议：对于带 +1 或 +2 电荷的缺陷，建议至少使用 3×3×3 超胞，并配合 Makov-Payne 修正；若条件允许，用 4×4×4 与 3×3×3 的结果对比，检查修正后形成能差异是否小于 0.05 eV。若差异仍较大，说明超胞不足或修正参数（介电常数、马德隆常数）需要重新标定。

四、超胞构建的实际操作要点

以 VASP 为例，构建超胞并计算缺陷形成能时，可按以下流程操作：

# 1. 扩胞（以 Si 为例，原胞 2 原子）
# 使用 ASE 或 pymatgen 生成 3×3×3 超胞
python -c "
from ase.build import bulk
si = bulk('Si', 'diamond', a=5.43)
sc = si.repeat((3,3,3))
sc.write('POSCAR_3x3x3.vasp')
"

# 2. 确定缺陷位置：优先选择高对称位点
# 空位：直接删除一个原子
# 替位：替换目标原子
# 间隙：在八面体/四面体间隙插入原子

# 3. 设置 INCAR 关键参数
# ISIF = 2        # 固定体积，仅弛豫离子
# ISMEAR = 0      # 半导体用 Gaussian 展宽
# SIGMA = 0.05
# EDIFF = 1E-6
# NELM = 200

# 4. 带电缺陷：设置 NELECT 调整电子数
# 例如 +1 电荷：NELECT = 原电子数 - 1
# 同时设置 IDIPOL = 3 和 LDIPOL = .TRUE. 进行偶极修正

构建超胞时，缺陷位置应尽量靠近超胞中心，使最近镜像距离最大化。对于各向异性材料，应沿介电常数较小的方向扩胞更多，以更快衰减静电相互作用。此外，超胞形状应尽量接近立方，避免某一方向过短导致镜像作用异常增强。

五、收敛性检查与结果判读

判断超胞是否足够，不能只看单次计算，而应做系统的收敛性检查。具体做法是：固定缺陷类型和电荷态，依次计算 2×2×2、3×3×3、4×4×4 超胞的形成能，观察其随 1/L 的变化。若修正后形成能在最后两个尺寸间变化小于 0.05 eV，可认为收敛；若仍明显漂移，则需继续增大超胞或检查修正参数。

同时应检查缺陷态能级是否随超胞尺寸移动。若缺陷能级在带隙中的位置随超胞增大而显著变化，说明波函数交叠仍不可忽略，需要更大超胞。对于带电缺陷，还应检查静电势对齐是否合理：若对齐后形成能出现非物理的跳变，往往说明势对齐区域选择不当或背景电荷补偿有误。

常见误区：直接用 2×2×2 超胞计算带电缺陷形成能，且不做任何修正，得到的数值可能偏离孤立缺陷极限 0.5 eV 以上。另一个误区是只增大超胞但不做势对齐，导致形成能随超胞尺寸非单调变化，误判为“已收敛”。

综上，超胞模型构建的核心是：用足够大的超胞降低缺陷浓度和镜像相互作用，用有限尺寸修正把结果外推到孤立缺陷极限，并通过多尺寸收敛性检查确认精度。只有把尺寸选择、间距控制、镜像修正和结果判读作为一个整体来考虑，才能得到可靠、可比较的缺陷形成能数据。