半导体缺陷迁移计算以明确的路径对象为单位。空位、间隙位和反位缺陷在不同电荷态下会有不同局域弛豫；电荷态、端点身份和原子映射必须同时固定。后续端点优化、图像插值和能量读取都围绕同一缺陷跳跃通道展开。
误差主要来自三处耦合：带电态通过 `NELECT` 改变电子数，端点优化改变局域构型，NEB 图像把多个非平衡结构放在同一次路径优化中。任一位置混入不同 q 态、不同缺陷中心或不同原子编号，曲线仍可能收敛，但峰值已经偏离目标迁移通道。
q 态由形成能区间决定
带电缺陷迁移计算通常由形成能图确定候选电荷态。图1左侧的形成能折线给出不同 q 态随费米能级变化的稳定区间；形成能筛选 q 态，NEB 计算鞍点。右侧的 NEB/CI-NEB 示意说明迁移势垒来自同一势能面上的鞍点，二者对应不同层级的问题。
目标费米能级附近稳定的是 VI+ 时，迁移计算围绕 VI+ 建端点；稳定区转到 VI0 时，中性路径单独建立。中性端点的 NEB 曲线若被用于带电缺陷迁移归因，结构弛豫、电子占据和势垒来源会同时混入结果。
图1 半导体缺陷形成能与 NEB/CI-NEB 能量路径示意。左图说明电荷态稳定区随费米能级变化，右图说明迁移势垒来自反应坐标上的鞍点。DOI: 10.1021/acsmaterialsau.4c00095。
式(1)给出形成能表达式。迁移势垒来自 NEB 鞍点，但式中的 q、EF、Ecorr 决定了哪一个带电缺陷进入迁移计算。外源掺杂或光照下的非平衡缺陷需要把所选 q 态写入模型前提，后续势垒讨论才能对应同一电子数条件。
Ef(Dq) = Etot(Dq) – Etot(bulk) – Σniμi + q(EF + EVBM + ΔV) + Ecorr 式(1)
NELECT 怎样在整条路径中保持一致
VASP 中带电态通常通过 `NELECT` 控制总电子数。以同一个缺陷超胞为参照，q=+1 表示少一个电子，q=-1 表示多一个电子；同一条 NEB 路径的 00、01、02 到终态目录必须使用相同 NELECT。实际输入应以中性缺陷超胞的总电子数为基准核对，否则图像之间已经不在同一电子数条件下比较能量。
半导体体系常用较小展宽，例如 `ISMEAR = 0`、`SIGMA = 0.05` 作为起点；含局域磁矩的缺陷要显式设置 `ISPIN` 和合理的 `MAGMOM`。若缺陷态在带隙中局域，端点优化和 NEB 静态检查阶段还应比较总磁矩、局域磁矩、带隙内占据和缺陷附近电荷分布，确认多余电子或空穴没有转移到无关表面态或远端缺陷。
/// INCAR 参考示例：带电缺陷端点优化
SYSTEM = charged defect endpoint relax
ISPIN = 2
ISMEAR = 0
SIGMA = 0.05
ENCUT = same_as_bulk_test
EDIFF = 1E-5
EDIFFG = -0.02
IBRION = 2
NSW = 200
ISIF = 2
ISYM = 0
NELECT = neutral_defect_electrons – q
LASPH = .TRUE.
初态和末态要先分别优化
端点构型来自目标 q 态下的两个局域极小值。初态 A 和末态 B 都在同一带电态下单独弛豫；NEB 的端点应是同一电荷态下的两个局域极小值。同一超胞、同一元素顺序、同一赝势、同一 k 点和同一约束都要保留，人工拖拽结构只作为末态初猜。
图2中的 CsPbI3 碘间隙和碘空位路径说明，Ii–、Ii0、Ii+ 以及 VI–、VI0、VI+ 的局域位移并不相同。电荷态改变后，缺陷附近键长、配位和迁移通道都会变化，每个 q 态都需要独立端点。
图2 CsPbI3 中碘间隙和碘空位在不同电荷态下的迁移构型。不同 q 态对应不同局域位移和缺陷态占据，不能混用同一条 NEB 路径。DOI: 10.1021/acs.jpclett.5c01139。
端点 A、端点 B、NEB 路径、鞍点静态检查和频率计算分别保留 OUTCAR、CONTCAR、OSZICAR 和关键后处理结果。每个能量值同时绑定 q 态、端点编号、图像编号和后处理口径，后续回看时可以直接定位到对应计算目录。
原子编号和空位表示怎样核对
NEB 插值依赖初末态中原子顺序完全一致。空位迁移的常见偏差来自空位没有原子编号：真正移动的是邻近原子跳入空位，因此初态和末态的 POSCAR 仍然少同一个元素总数，但参与跳跃的那个原子编号要保持可追踪。ASE、pymatgen 或 VTST 工具插值前，元素顺序和原子顺序需要锁定。
|
|
|
|
|---|---|---|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
图像数和优化分两步处理
半导体缺陷迁移的图像数由通道形状决定，不能只由初末态距离决定。间隙位绕过配位多面体、空位跨越瓶颈、反位伴随局域换位时，反应坐标可能有弯曲。3 到 5 个图像可用于初始通道测试，峰附近再按曲率和力收敛情况加密；出现两个结构重排台阶时，引入中间态并分段计算。
NEB 起步阶段可先关闭 climbing image，让整条路径完成基本弛豫；曲线形状和最高能区域稳定后，再打开 `LCLIMB = .TRUE.` 精修鞍点。CI-NEB 适合精修最高能区域，端点身份、图像身份和电子态连续性仍由单独检查确认。带电体系中，某个图像收敛到不同自旋态或不同缺陷态分支时，需要回到该图像的初始磁矩和占据状态。
/// INCAR 参考示例：带电缺陷 NEB 路径
IBRION = 3
POTIM = 0
SPRING = -5
IMAGES = 5
LCLIMB = .FALSE. # 路径稳定后再开启
EDIFF = 1E-5
EDIFFG = -0.03
ISYM = 0
NELECT = same_as_endpoints
LWAVE = .FALSE.
LCHARG = .FALSE.
每张图像要看结构，也要看电子态
图3中的 DFT 与机器学习力场迁移路径说明，同一材料中的不同 q 态可以有不同峰位置和不同能量剖面。`neb.dat` 的最高点只给出能量序列，结构核对还需要逐帧打开 00 到终态的 CONTCAR，确认缺陷中心沿预期通道移动，周围配位连续变化。
对于带电缺陷，图像串的电子态检查和结构检查同等重要。若 VCd2- 的多余电子在某个图像转移到远端表面态，或 VI+ 的空穴在峰附近局域到不同原子，同一 q 态的路径比较才有能量意义。跨 q 态曲线需要同时给出端点身份和形成能背景；能量峰混入电子态切换时，势垒解释要同时包含几何鞍点和电子态分支变化。
图3 CsPbI3 中不同电荷态碘间隙和碘空位的 DFT/MLFF 迁移能量路径。曲线展示了电荷态、图像数和势垒高度之间的耦合关系。DOI: 10.1021/acs.jpclett.5c01139。
逐图像记录至少包含图像编号、相对能量、迁移原子到两个端点位点的距离、最近邻最短键长、总磁矩、缺陷原子局域磁矩、带隙内态能量位置。几何变量连续而磁矩突变时，优先检查电子态分支；磁矩连续而最短键长异常时，优先检查插值路径和图像数。
正向势垒和反向势垒都要报告
半导体缺陷的两个端点可能等价，也可能因为局域应变、表面电场、界面或有机阳离子取向而不等价。MAPbBr3 与 FAPbBr3 中溴空位、溴间隙的路径差异说明，晶格柔性会改变通道形状和峰位置。等价位点主要比较势垒高度和可用通道数；不等价位点还要报告端点能差和反向势垒。
图4 MAPbBr3 与 FAPbBr3 中溴空位和溴间隙迁移路径及能量剖面。不同晶格环境会改变通道形状、峰位置和迁移势垒。DOI: 10.1021/acsmaterialsau.4c00095。
式(2)把正向、反向和端点能差放在一起。ETSq 可来自 CI-NEB 的最高能图像或后续鞍点优化；不等价端点必须保留方向信息。Eiq 和 Efq 必须来自同一 q 态端点。
Em,fq = ETSq – Eiq；Em,rq = ETSq – Efq；ΔEq = Efq – Eiq 式(2)
带电缺陷还要说明修正口径。相同超胞、相同 q 态的 NEB 相对能量中，均匀背景电荷和部分静电误差可能相互抵消，但形成能、不同 q 态之间的比较、以及活化能 Ef+Em 仍要使用一致的电势对齐和有限尺寸修正方案。结论表需要逐列标明形成能和势垒采用的修正状态。
扩散趋势要区分浓度和迁移两件事
图5给出 CsPbI3 不同碘缺陷电荷态的 Arrhenius 扩散趋势。迁移势垒 Em 控制单个缺陷跳跃，缺陷形成能 Ef 控制热平衡浓度；扩散判断要分开浓度来源和迁移势垒。若缺陷来自非平衡制备或外源掺杂，浓度不一定由 Ef 决定。
图5 CsPbI3 中不同碘缺陷电荷态的扩散系数 Arrhenius 趋势。不同 q 态的迁移势垒和扩散行为需要在同一电荷态框架下解释。DOI: 10.1021/acs.jpclett.5c01139。
最后的结果表应包含缺陷类型、q 态、`NELECT`、端点能量、端点局域结构、图像数、最高能图像编号、Em,fq、Em,rq、ΔEq、是否做频率确认、形成能稳定区和修正方法。计算记录应追溯到每个端点和图像，后续补算也能准确定位到端点、图像或静态检查步骤。
鞍点确认可以从最高能图像出发做频率计算。单缺陷跳跃应出现主要沿迁移方向的一个虚频；若虚频混入大范围晶格软模，或模式方向与目标跳跃无关，当前峰值未必对应目标迁移通道的鞍点。此时需要回看图像串、端点定义和可能的中间态，过渡态表述以虚频模式和结构连续性为依据。
要点：半导体缺陷迁移计算要先确定缺陷类型、q 态、迁移方向和两个端点位点。
要点：端点和所有 NEB 图像必须使用相同 `NELECT`、相同超胞、相同原子顺序和相同计算口径。
要点：带电缺陷的图像检查同时包含结构连续性、总磁矩、局域磁矩和带隙内缺陷态连续性。
要点：势垒表应同时给出正向势垒、反向势垒、端点能差、修正方法和是否通过虚频方向确认。
