空位缺陷的第一性原理建模：结构定义、超胞尺寸效应、形成能与电子结构校验

说明：本文主要介绍空位缺陷模型的结构定义、超胞尺寸、形成能参照、电荷态处理和电子结构校验。

空位缺陷模型的结构对象是什么

空位缺陷模型的起点是一个明确的晶格位点。晶体中的某个原子从原本的Wyckoff位置、表面配位位或二维层内位点移除，局部化学计量和近邻配位随之改变。空位不是一个抽象空洞，而是带有元素类型、坐标、近邻原子和局部弛豫的结构模型。写成V_O、V_S或V_Cl时，符号已经包含“移走哪一种原子”的信息；若同一种元素存在不等价位点，需区分桥氧、顶位氧、面内硫、边缘硫或体相氧。

构建缺陷前，完整晶胞通常先经过几何优化。优化后的晶格常数、原子坐标和磁态构成无缺陷参照；随后在同一晶胞或扩展超胞中移除目标原子。参照结构和缺陷结构应来自同一套泛函、赝势、U值、自旋设置和色散修正，否则形成能、态密度和局域磁矩会混入方法差异。对氧化物、卤化物和过渡金属硫属化物，移除原子后常出现近邻阳离子位移、键长重排和局域电子重新分布，缺陷结构要重新弛豫到力收敛。若弛豫后邻近原子重新成键，空位标签仍应对应最初移除的晶格位点，不能把重构后的空腔位置当作新的缺陷定义。

空位类型还受研究对象约束。体相空位强调三维周期晶格中的缺失位点和体相化学势；表面空位还牵涉slab厚度、真空层、表面终止和上下表面对称；二维材料中的空位会受到面外真空、库仑截断和缺陷周期阵列影响。模型名称应写到具体位点和计算场景，例如“MoS₂单层4×4超胞中的单个V_S”比“一个硫空位模型”更利于后续能量和图谱比较。缺陷浓度也随超胞写法改变，同样移除一个S原子，3×3和6×6超胞对应的空位间距与电子态耦合并不相同。

超胞和周期条件怎样影响空位模型

一个V_S放进4×4 MoS₂超胞后，周期DFT实际计算的是沿三个方向重复的硫空位阵列；一个V_O放进小氧化物晶胞后，缺陷浓度可能远高于实验样品。若研究稀缺陷极限，超胞内相邻空位的最短距离应大于局域弛豫区和缺陷态尾部。超胞尺寸决定缺陷之间的距离，也决定缺陷能级、弛豫应变和形成能误差。超胞过小会放大缺陷-缺陷相互作用，带隙内局域态可能展宽成窄带，局部弛豫也会受到相邻镜像限制。

选择超胞时，目标不是把原子数做得无限大，而是让待比较的量进入稳定区间。形成能可测试2×2、3×3、4×4或更大超胞的变化，缺陷态可查看带隙内峰位和带宽，局域磁矩可检查是否集中在空位近邻。二维材料应让真空层容纳面外电荷尾部，表面模型则要检查slab厚度和偶极修正。k点密度也要随超胞放大同步调整，否则总能差中的布里渊区采样误差会伪装成缺陷浓度效应。形成能收敛、带隙内缺陷峰宽和近邻键长变化可作为超胞尺寸筛选的三类结果，三者若同时随尺寸变化变小，模型才接近孤立空位假设。

缺陷有效区域由结构弛豫和缺陷态空间分布共同限定。若缺陷波函数集中在空位周围几层近邻原子，较远区域可以作为近似体相环境；若缺陷态沿某一方向延展，超胞沿该方向的长度就要增加。局域态、长程极化和弹性弛豫对应不同尺度：局域态常由近邻轨道决定，极化和带电缺陷静电项可能跨越整个超胞，弹性位移在低维材料中衰减更慢。固定原子层、真空厚度和偶极修正应在同一组模型中保持一致，表面空位与体相空位的能量参照也不宜交叉比较。

形成能和电荷态怎样限定空位模型

空位模型完成几何弛豫后，形成能把缺失原子和电子数变化写入同一个能量参照。常用形式可写为E_form(q)=E_defect(q)-E_perfect+Σn_iμ_i+q(E_F+E_VBM)+E_corr。式中μ_i记录被移走原子的化学势，q记录缺陷电荷态，E_F记录费米能级位置，E_corr处理周期带电模型的静电误差。同一个空位在O-rich、O-poor、S-rich或S-poor条件下会得到不同形成能。

带电空位模型比中性空位多一个电子数约束。周期程序通常用均匀背景电荷保持总电荷中性，背景电荷会和缺陷镜像相互作用，尤其在slab、二维材料和真空层较大的模型中更明显。裸的带电超胞总能差会携带镜像电荷和势能对齐误差，校正项应和缺陷电荷分布、介电常数、材料厚度及真空层设置相匹配。

形成能图中的不同斜率对应不同电荷态，交点对应热力学转变能级。若研究V_O⁰、V_O⁺和V_O²⁺，各电荷态要分别优化结构，并比较它们在费米能级范围内的低能区间。稳定电荷态来自能量线的相对位置，不来自手动指定的电荷数本身。对含过渡金属的氧化物，自旋多重度、局域磁矩和U值还会改变缺陷电子在邻近金属原子上的分布。

低维材料的静电校正还要区分面内和面外介电响应。MoS₂这类单层材料只有有限厚度，面外真空不能按三维体相处理；表面空位也有类似问题，因为缺陷电荷一侧面对材料，一侧面对真空。二维空位模型的真空层、介电剖面和电荷分布形状会共同进入校正能，形成能比较前应固定模型几何和校正方案。若比较不同空位位点，化学势、费米能级零点和电荷校正必须沿用同一套参照，能量差才反映位点环境，而不是计算设定变化。

构建后的空位模型要用哪些图谱校验

结构弛豫后的空位模型先看局部几何：近邻键长是否异常拉长或压缩，配位数是否符合预期，表面是否发生重构，二维层是否出现明显翘曲。随后检查总能、力和电子收敛，带磁体系并比较不同初始磁矩给出的终态能量。几何结构、形成能和电子态要指向同一个缺陷位点，否则模型可能只是数值上收敛，物理图像并未闭合。局部结构表可记录空位近邻键长、配位数、磁矩和Bader电荷，这些量能直接对应后续DOS、PDOS和自旋密度图。

DOS和PDOS用于识别空位引入的带隙态、费米能级附近态密度和轨道来源。以MoS₂硫空位为例，总态密度中出现带隙附近缺陷能级，PDOS进一步把这些峰归属到空位邻近Mo的d轨道贡献。缺陷态的位置要和空位周围原子的投影轨道相互对应；若峰位随超胞尺寸明显漂移，或缺陷态在空间上不再局域，模型浓度和周期相互作用还要重新评估。

电荷密度、差分电荷密度、自旋密度和Bader电荷可补充DOS/PDOS。电荷密度查看缺陷电子是否局域在空位近邻，差分电荷密度比较缺陷形成前后的电子积累与耗尽，自旋密度识别未成对电子中心，Bader电荷给出原子盆地内电子数变化。这些图谱分别记录空间分布、轨道贡献和分区电荷，它们服务于同一结构模型时，空位缺陷的位点、形成能、电荷态和局域电子结构才具有可比较的计算参照。若缺陷态主要分布在远离空位的原子上，模型中可能存在表面重构、额外吸附物或电荷补偿方式改变，需回查原子坐标和总电荷设置。