1. 课程导入与缺陷工程全景
在半导体、氧化物电子学、催化与能源材料研究中，缺陷往往同时扮演“性能杀手”和“功能来源”两种角色。氧空位既可能成为导电通道，也可能引入深能级复合中心；氮空位在GaN中既参与掺杂补偿，也可能形成量子发射中心。要判断一个缺陷究竟带来什么后果，仅凭“有缺陷”这一事实远远不够，必须回答三个定量问题：它有多容易形成？它稳定在哪个电荷态？它的电子能级落在带隙的什么位置？本课程围绕这三个问题，建立从DFT缺陷形成能计算到电荷态工程设计的完整工作流。
核心物理量：形成能与电荷态跃迁能级
缺陷形成能的标准定义是：在给定化学势和费米能级下，将缺陷引入完美晶体所需的能量。对电荷态为 q 的缺陷，其形成能写作：
E_f(D,q) = E_tot(D,q) - E_tot(bulk) - Σ n_i μ_i + q(E_F + E_VBM) + E_corr其中 E_tot(D,q) 是含缺陷超胞的总能，E_tot(bulk) 是完美超胞总能，n_i 和 μ_i 分别是增删原子的个数与化学势，E_F 是相对于价带顶 E_VBM 的费米能级，E_corr 是有限尺寸超胞的电荷修正项。这个公式的每一项都对应一个可操作的计算步骤，也是后续章节逐一展开的对象。
把同一缺陷不同电荷态的形成能画在同一张图上，就得到形成能随费米能级变化图。两条直线的交点即热力学跃迁能级，它决定了在什么掺杂条件下缺陷以何种电荷态存在。例如，若 (+/0) 跃迁能级靠近价带顶，说明该缺陷在n型材料中倾向带正电；若 (0/-) 能级位于带隙中央，则它很可能成为有效的复合中心。电荷态工程的目标，正是通过化学势、掺杂、应变或共掺杂，把这些跃迁能级移动到期望位置。
为什么必须做电荷修正与超胞收敛
带电缺陷在周期性超胞中会与其镜像发生库仑相互作用，导致形成能系统性偏移。常用的Freysoldt-Neugebauer-Van de Walle修正方案，通过求解缺陷电荷密度与介电张量，扣除镜像相互作用能。实际操作中，至少要做三件事：第一，检查超胞尺寸从2×2×2到3×3×3时跃迁能级是否漂移小于0.05 eV；第二，确认k点网格密度使总能收敛到1 meV/atom以内；第三，对带电缺陷检查是否出现“虚假的能带填充”或“缺陷态与导带混淆”。这些检验不是可选项，而是判断结果是否可信的前提。
学习路径与实战目标
本课程按“概念—计算—分析—设计”四层推进。前四章建立形成能、化学势、跃迁能级和有限尺寸修正的物理图像；中间章节进入VASP/Quantum ESPRESSO的实际输入文件编写、缺陷结构构建、自旋极化设置与后处理脚本；随后用氧化物、氮化物、二维材料等真实体系做案例拆解，展示如何从形成能图读出掺杂可行性与缺陷补偿机制；最后落到电荷态工程，讨论如何通过生长条件、共掺杂和应变把目标缺陷稳定在所需电荷态。
一个可执行的最小示例
以MgO中氧空位为例，假设使用3×3×3超胞（216原子），PBE泛函，平面波截断能500 eV，k点网格2×2×2。计算流程如下：
- 优化完美MgO超胞，得到 E_tot(bulk) 和价带顶 E_VBM。
- 移除一个氧原子，分别设置 q=0, +1, +2 三个电荷态，进行自旋极化弛豫。
- 对带电超胞施加背景电荷补偿，并计算Freysoldt修正能。
- 取氧的化学势 μ_O 为O₂分子总能的半值（O-rich极限），镁的化学势由热力学平衡条件确定。
- 代入形成能公式，绘制 E_f 随 E_F 从VBM到CBM的变化曲线。
若计算得到 ε(+/0) 位于VBM上方0.8 eV，ε(0/++) 位于1.6 eV，则说明在费米能级靠近VBM时氧空位倾向带+2价，在带隙中央区域以0价为主。这一结果直接告诉实验者：在p型MgO中氧空位可能充当补偿中心，而在弱n型条件下它更可能保持中性。这就是从DFT数字到材料设计判断的完整链条。
接下来的章节将逐一拆解上述每一步的物理含义、参数选择依据和代码实现细节。你最终要达成的目标不是“会跑一个缺陷计算”，而是能够独立判断：给定一个材料和一个缺陷，它在什么条件下形成、以什么电荷态存在、对电子结构产生什么影响，以及如何通过可控手段把它调整到想要的状态。