为什么掺杂会改变 HOMO-LUMO 能级？计算对象、影响因素与核验方法

说明：本文华算科技主要介绍掺杂前后 HOMO、LUMO 的计算对象，原子势、电子数、几何弛豫与库仑环境对前线能级的作用，以及轨道重排和计算参照的核验方法。

掺杂前后比较的HOMO和LUMO是什么

HOMO 是体系基态中能量最高的占据轨道，LUMO 是能量最低的空轨道，二者的单粒子本征值差写作 Δε_HL=ε_LUMO−ε_HOMO。轨道名称由电子占据和能量排序共同确定；掺杂计算完成自洽迭代后，程序会按新的本征值与占据重新命名 HOMO、LUMO，原来的 HOMO 未必保留原编号。

“掺杂”在计算中既可指替位或间隙原子，也可指载流子注入。前一种操作改变原子种类、局部配位和价电子数；后一种操作通过分子受体、分子给体、栅压或带电模型改变载流子数，并由反离子或均匀补偿电荷维持电中性。替位或间隙原子主要改动 V_ion 与价电子构成，载流子注入主要改动占据、V_H[n] 和 V_xc[n]。

有限分子、团簇和充分隔离的量子点适合使用 HOMO/LUMO。周期晶体中的态随波矢 k 色散，术语通常改为 VBM、CBM 与 E_F。图1所用 C₃₀H₁₄ 量子点在三个方向设置 20 Å 真空层并只取 Γ 点，HOMO/LUMO 对应零维模型的离散轨道；把同样的异原子放进无限周期石墨烯后，讨论对象会变成杂质带、狄拉克点和费米能级。

掺杂前后轨道比较包含三项：本征值、占据数和波函数形状。只列 ε_HOMO、ε_LUMO 会遗漏轨道是否换序；只画等值面又会遗漏能量参照。轨道重叠、原子投影和节点位置能追踪同一电子态，轨道编号只记录当前排序，两者承担的物理含义不同。

原子势/电子数/几何弛豫怎样移动前线能级？

Kohn–Sham 方程中的离子势 V_ion 由元素和原子位置给定，Hartree 势 V_H[n] 与交换相关势 V_xc[n] 又随自洽电子密度 n(r) 更新。异原子替换 C 原子后，核电荷、赝势和价电子数一起改变；局域 p 轨道能量、C—X 耦合矩阵元及电子占据随之变化。给体更常抬高占据态或引入额外电子，受体更常形成低能空态，但实际方向还受位点和周围碳骨架控制。

图2中的紫色等值面仍大量分布在锯齿边缘，掺杂原子及邻近 C 原子的权重却随元素和位点改变。局域势扰动会选择性移动具有该区域权重的轨道：某个态在掺杂位附近的振幅越大，它感受到的势能修正通常越明显；远离掺杂位的边缘态可能保持相近形状和能量。HOMO 与 LUMO 对局域扰动的空间重叠不同，二者很少等量平移。

几何弛豫进一步改写轨道耦合。该量子点本征 C—C 键长约 1.4 Å；内部 N—C 键仍约 1.41 Å，内部 P—C 和 S—C 则增至约 1.65 Å 与 1.77 Å。键长控制轨道重叠，键角和翘曲控制对称匹配，较大的 P、S 原子会削弱原有 sp² 共轭并产生新的局域组合。固定未掺杂几何与充分弛豫后的掺杂几何会给出不同 ε_HOMO 和 ε_LUMO。

电荷转移掺杂不替换骨架原子，却改变电子总数、局域极化和静电势。三芳胺–芴共聚物经氧化形成空穴，SbF₆⁻ 或 BArF⁻ 负责电荷补偿。空穴—反离子距离产生 Madelung 势，反离子尺寸差异可使强掺杂薄膜的价带与功函数移动 0.3–0.4 eV；前线轨道的绝对位置由骨架电荷、反离子和介电极化共同设定。

替位掺杂中的 V_ion 改变、结构弛豫中的耦合项改变、载流子掺杂中的占据改变、反离子引起的静电势移动会在同一次自洽计算中叠加。HOMO 上移、LUMO 下移或两者同向移动都可能出现，掺杂元素的“给体/受体”标签只描述电子交换倾向，并未固定 Δε_HL 的变化方向。

能级移动与前线轨道换身份为什么会同时发生？

闭壳层中每个分子轨道容纳两个反平行自旋电子。体系失去一个电子后，原 HOMO 变成单占轨道 SOMO；自旋极化计算会给出占据与空的自旋分量。电子数改变会重新定义前线轨道，此时继续把旧 HOMO 与新 HOMO 按编号配对，常会配到不同波函数。

mTFF 模型中，空穴掺杂使原 HOMO 分裂为 SOMO 与空 SOMO*，CAM-B3LYP 计算得到的 Hubbard 间隔为 3.9 eV。原 HOMO−1 随后成为新的 HOMO′。“HOMO 能量改变”在这里包含轨道换序：新 HOMO′ 的空间骨架更接近旧 HOMO−1，SOMO 才保留原 HOMO 的主要三芳胺特征。库仑充电又把 SOMO 压到更低能量，反离子的吸引部分抵消该位移。

轨道重排也出现在实验价带谱中。掺杂程度 DL 从 0 增至接近 1 h⁺/重复单元时，原 HOMO 强度衰减，新谱重进入原能带间谷区，HOMO−1 相关谱带展宽。谱带中心、轨道波函数和电子占据三项相互对应，可区分整体静电位移、局域杂质态和轨道顺序交换。

当掺杂中心从单个缺陷扩展为有限浓度，局域轨道之间出现耦合，离散能级逐渐展宽成 DOS 分布。费米能级由占据态和空态尾部的热平衡位置确定，并非固定停在一条半占据轨道上。TAF 共聚物在高掺杂区的功函数随掺杂浓度每十倍约移动 95 meV，DOS 宽度、SOMO* 与 HOMO′ 的相对位置共同控制该斜率。

替位掺杂若引入未成对电子，也会出现 α、β 两套前线轨道。SOMO、磁矩与自旋密度应保持同一占据描述，用闭壳层 HOMO/LUMO gap 表示开壳层体系会抹去自旋分裂。氧化态、总电荷、多重度或 MAGMOM 初值不同，都可能得到另一组占据顺序。

怎样确认能级变化来自掺杂本身？

冻结几何与弛豫几何构成一组对照。本征结构、本征结构中直接替换元素、掺杂后充分弛豫三套计算使用同一泛函、基组或赝势、总电荷和自旋设置。第一处能量差主要含元素势与电子数贡献，第二处差值记录键长、键角和翘曲贡献。再把掺杂原子放到边缘、内部与不同配位位点，位点响应便与元素响应分开。

绝对能级比较要使用统一能量零点。孤立量子点可用真空势对齐 ε_HOMO、ε_LUMO，带电超胞还要评估有限尺寸与背景电荷修正；周期固体改报 VBM、CBM、E_F 和投影权重。同一模型内的 Δε_HL 不受整体常数势平移影响，单条 HOMO 或 LUMO 的绝对数值则会随真空参照、反离子距离和表面偶极移动。

Kohn–Sham gap、准粒子 gap 和光学 gap 属于不同计算对象。C₃₀H₁₄ 本征量子点的 PBE gap 为 1.32 eV，GW 修正后为 4.28 eV；N-IN 构型对应 0.198 eV 与 3.114 eV。方法改变量级可大于掺杂差值，所以一组掺杂序列应保持同一数值口径。光吸收还包含电子—空穴相互作用，GW+BSE 谱峰承担的物理含义与 ε_LUMO−ε_HOMO 不同。

轨道追踪可用相邻计算间的波函数重叠或最大投影，避免能级交叉时沿用错误编号。ε_HOMO 与 ε_LUMO 对应的原子/轨道权重应与 PDOS 的近前线峰一致；局域在掺杂位的空态还应在该位点投影中出现。冻结几何差、弛豫差、轨道重叠矩阵和真空对齐分别回答元素势、结构、轨道身份与绝对能量四个问题。

电荷分区和自旋密度负责核对占据变化，ΔSCF 电离能与电子亲和能负责核对加减电子代价，形成能负责限定掺杂构型出现的热力学条件。新 LUMO 若局域在掺杂位附近，对应 PDOS 应出现同一元素与轨道权重，Bader 分区或自旋密度也应在相同邻域产生响应。