什么是掺杂？从晶格身份、电子计数到电荷补偿，杂质能级与第一性原理判断范围详解

说明：本文华算科技主要介绍掺杂的晶格身份、电子计数、电荷补偿、杂质能级和第一性原理判断范围。

掺杂的基本概念怎样限定原子身份？

掺杂描述少量异质原子进入宿主材料后形成可辨认的占位，并改变局部原子环境或电子态。宿主仍保留原有晶体骨架，掺杂原子只占有限比例。掺杂量持续升高后，结构可能转入固溶体、合金或有序化合物区间，稀释杂质模型随之失去对应对象。

按原子收支可区分多种占位。外来原子 D 替换宿主 A 位时，模型含有替位缺陷 D_A；D 加入晶格空隙时形成间隙缺陷 D_i。表面 top、bridge 或 hollow 位上的 D 仍处于晶体外侧，许多催化研究把它纳入表面掺杂；体相缺陷研究通常称其为吸附原子。术语必须绑定原子位置，元素名称本身无法区分四类结构。

采用周期性超胞描述掺杂时，一个 D 会随晶格平移重复，4×4 单层替换一个原子与 2×2 单层替换一个原子具有不同名义浓度。掺杂原子之间的距离、应变场重叠和杂质能带宽度均会改变。超胞尺寸同时规定浓度和周期镜像间距，单个结构的结果只对应当前掺杂比例。

化学配方中的名义浓度与晶格中的有效浓度并非同一个量。掺杂原子可能进入目标位点，也可能停在晶界、表面或第二相中；进入晶格的部分还会呈现随机分布、短程有序或团簇。计算浓度由模型原子数精确定义，实验有效浓度则要由占位率、空间分布和相组成确定。

初始坐标只代表候选位。弛豫后，表面原子可能嵌入空位，间隙原子可能挤出宿主原子，替位原子附近还可能生成补偿空位。As 在单层 MoS₂ 的表面位、间隙位与 S 替位具有不同形成能；费米能级和富 Mo／富 S 化学势还会改变各电荷态的能量斜率与占优构型。

形成能较高的构型仍可能通过离子注入、快速退火或非平衡生长暂时保留，形成能较低的构型也可能受扩散势垒限制而难以达到。热力学占位与制备后的动力学分布对应不同时间尺度，掺杂定义本身不包含“易合成”或“长期稳定”的保证。

掺杂原子改变的电子数去了哪里？

等价掺杂常指掺杂离子与宿主位点具有相同名义价态，例如 Mo⁶⁺ 替换 WO₃ 中的 W⁶⁺。价电子计数不变，离子半径、电负性、轨道能量与键长仍会改变。异价掺杂产生有效电荷差，例如 La³⁺ 替换 ZrO₂ 中的 Zr⁴⁺ 后记作 La_Zr^′。

晶体整体满足电中性。La_Zr^′ 可由电子数调整形成的电子补偿承担，也可由两个 La_Zr^′ 配一个 V_O^•• 的离子补偿承担。空穴、电子、氧空位、阳离子空位、间隙离子和反位缺陷都可能参与补偿。名义价态差规定电荷账目，局部能量决定补偿物种。

平衡载流子与带电缺陷满足 Σc_jq_j + p − n = 0。c_j 与 q_j分别表示第 j 类缺陷浓度和有效电荷，p、n 分别为空穴与电子浓度。费米能级通过各缺陷形成能与载流子统计共同确定，施主浓度增加并不保证 n 以相同比例增加，受主缺陷可抵消一部分电子。

相同 6.25 cat.% La 浓度可以对应多种 La–La 间距与数十个氧空位位置。电子补偿模型的能量分布较窄，离子补偿模型因 V_O^•• 位置增加了构型离散。掺杂浓度只给出成分比例，并未指定掺杂原子的空间排列和补偿缺陷位置。

在离子步弛豫中，预设缺陷类型可能随原子位置重排而改变。La:ZrO₂ 中一个 O 原子离开晶格位并进入间隙后，产生 V_O^••–O_i^″ Frenkel 对；另加氧空位时，部分位置触发复合，部分位置保留 Frenkel 对。弛豫后的拓扑和原子收支才定义实际补偿结构。

掺杂会产生浅能级/深能级和局域磁矩？

施主释放电子后从中性态转入正电态，受主俘获电子后从中性态转入负电态。两种电荷态形成能相等时的 E_F 位置称为热力学电荷转变能级 ε(q/q′)。靠近导带底的施主能级和靠近价带顶的受主能级具有较低电离代价；位于带隙中部的深能级常表现为载流子俘获中心。

价电子多一个或少一个并不自动产生浅能级。掺杂原子附近的键长重排、介电屏蔽和电子–晶格耦合会使载流子局域化。单层 MoS₂ 中 Re_Mo 施主和 Nb_Mo 受主都出现局域电荷分布；PBEsol 与满足广义 Koopmans 条件的 HSE给出不同转变能级和局域形态。

经典浅杂质模型把束缚载流子视为受介电屏蔽的氢样态，其结合能随有效质量增大而升高、随介电常数增大而降低。原子尺度的中心势、强局域弛豫或小极化子形成会破坏氢样近似。浅能级描述低电离能，深能级描述强局域或较高电离能，两者不由元素族序单独决定。

过渡金属掺杂还引入 d 轨道占据、自旋交换和晶场分裂。同一种 dopant 位于体相、平坦晶面或台阶位时，连续配位数与近邻距离不同，d 态峰宽随之变化。Zr 和 Rh 嵌入 Cu、Ag 宿主的计算显示，宿主晶格常数与掺杂位配位约束会改变 d-PDOS 的展宽，孤立原子身份不足以规定轨道形态。

能带里的三种变化需要分别命名：费米能级移动描述电子化学势和占据变化，带边移动描述宿主 VBM／CBM 的能量变化，带隙内新峰描述杂质态或缺陷态。超胞能带还含折叠效应；原子与轨道投影 DOS、带展开、局域电荷密度及自旋密度可定位态的来源。

稀释浓度下的杂质波函数彼此分离，较高浓度会增加波函数重叠并形成具有有限带宽的杂质带。过渡金属掺杂中的未成对 d 电子还可能产生局域磁矩；磁矩大小受晶场、轨道占据、电荷态和掺杂原子间交换耦合控制。自旋极化计算要比较多组初始磁矩，单一初始自旋可能收敛到局域亚稳态。

第一性原理中的掺杂效应由哪些量限定

带电掺杂缺陷 D^q 的形成能通常写为 E_f(D^q) = E_tot(D^q) − E_tot(host) − Σn_iμ_i + q(E_F + E_VBM) + E_corr。n_iμ_i记录加减原子的化学势，q 项记录电子交换，E_corr修正周期带电超胞的有限尺寸误差。

不同位点的 E_f 比较必须使用一致的宿主超胞、泛函、原子化学势和电荷校正。带隙低估会压缩允许的 E_F 区间，并移动缺陷能级相对带边的位置；过渡金属氧化物还受U 值与磁序影响。化学势范围须满足宿主相稳定，并排除竞争相析出的区域。

形成能排序描述给定化学势与 E_F 下的热力学候选；ε(q/q′)描述电荷态切换；键长、配位数、ELF 与 COHP描述局域成键；Bader 电荷和差分电荷密度描述分区电子数与空间重排；能带、DOS 和有效质量描述能级与输运相关电子态。每个量对应不同物理对象，性能结论要由目标过程直接相关的能量或动力学量承担。

WO₃ 中 Mo⁶⁺ 属等价掺杂，V⁵⁺ 与 Nb⁵⁺ 属异价掺杂。DFT 结果把氢结合能与氧空位形成能按三种局域 H 位置分组，含空位和无空位构型呈现不同斜率。价态差、补偿空位和质子落位共同改变局域热化学，单个“掺杂后电荷增加”描述覆盖不了图中的分支。

当掺杂原子间距由稀释极限缩短到有限浓度时，离散杂质态可展宽为杂质带，局域磁矩可能耦合，补偿缺陷也会重新排布。计算记录必须固定掺杂位、浓度、化学势、电荷态、自旋态与补偿缺陷，并报告对应的 E_f、ε(q/q′)、局域键长和轨道投影信号。