说明:本文华算科技主要介绍能带的形成来源、第一性原理计算中的能带构造、占据态与禁带的物理含义,以及用能带判断材料性质时需要保留的计算前提。
一、能带这个概念本身是什么?
孤立原子的电子能级是离散的,电子只能占据某些特定能量。固体晶体中有大量同类原子按照周期排列,外层轨道发生空间重叠,原来相同的原子能级会分裂成大量彼此很近的能级。数量达到宏观尺度后,这些密集能级形成近似连续的能量区间,这个允许电子存在的能量区间就是能带。相邻能带之间若没有可用电子态,就形成禁带或带隙。
能带的“带”并不是一条人为画出的线,而是周期势场中电子波函数的允许解集合。晶体平移对称性使电子态可用 Bloch 波描述,每个电子态带有能带编号 n 和波矢 k。同一个能带在不同 k 点有不同能量,曲线弯曲程度来自轨道重叠、键合反键合关系和晶体对称性。轨道相互作用越强,能级随 k 的变化通常越明显;局域化轨道贡献较多时,能带可能更平。

图1. 紧束缚框架将原子轨道、Bloch 轨道和晶体能带联系起来。DOI:10.1039/D4MH01038H
在第一性原理文章中,能带常沿布里渊区的高对称路径呈现。横轴不是实空间距离,而是倒空间中的 k 路径;纵轴是电子态能量,常用 eV 表示。能带概念的对象是允许能量和 k 空间波函数,不是某个原子静止不动时的单个电子能级。
高对称点的选择会改变图中横轴标签,却不改变能带来自周期体系本征态这一事实。对低对称晶体,路径更长、分段更多,能带曲线也更密集。同一张能带图同时包含能量、对称性和轨道重叠的信息,读它时要先确认晶胞、布里渊区和能量零点。

图2. Cs3Sb、CsK2Sb 和 Cs2Te 的晶胞与布里渊区高对称路径。DOI:10.3390/mi12091002
从原子能级过渡到晶体能带时,内层轨道和价电子轨道的表现并不相同。内层轨道与邻近原子的重叠弱,形成的带通常很窄;价电子轨道参与成键、反键和轨道杂化,能带宽度更容易随键长、配位和晶格对称性变化。能带只规定电子态允许出现的能量范围,具体填充情况还由电子数、温度近似和费米能级位置决定。
二、第一性原理计算里能带从哪些量来?
DFT 能带来自自洽电子密度对应的 Kohn-Sham 方程。计算先在给定晶体结构、赝势、交换关联泛函和 k 点取样下得到电子密度,再沿指定高对称路径求解一组本征值 εn,k。能带曲线实际记录的是 εn,k 随 k 的变化,其中 n 标记第几条能带,k 标记倒空间位置。
能量零点常设在价带顶、费米能级 EF 或真空能级。若不同论文的零点不同,同一材料的曲线位置看起来会移动。真正可比较的是带隙、带边相对位置、曲线色散和轨道贡献,不是图上某条线距离页面底部有多高。
Kohn-Sham 能级常用于分析电子结构趋势,但它不完全等同于实验光电子能级。PBE、SCAN、HSE、GW、SOC 和 U 值都会改变带隙、带边次序或带分裂。图中同一材料用 SCAN 与 G0W0 给出的能带可以接近,也可能在导带细节上有偏移,方法差异会进入能带判断本身,特别是半导体带隙、重元素 SOC 和强关联 d/f 电子体系。

图3. Cs3Sb、CsK2Sb 和 Cs2Te 的 SCAN 与 G0W0 能带对比。DOI:10.3390/mi12091002
晶格常数、应变、层数、缺陷、电荷态和磁序会重新塑造周期势场,使能带计算对结构模型高度敏感。结构优化后得到的原子位置一旦改变,εn,k 也会随之移动。同一化学组成在不同结构模型下可能对应不同能带,二维材料、表面 slab、异质结和掺杂超胞便不能直接互相套用能带结论。
当能带和 DOS、PDOS 同时出现时,二者承担不同任务。能带保留 k 方向上的色散,DOS 把所有 k 点的态按能量汇总。能带负责告诉电子态怎样随晶体动量变化,DOS 更适合判断某个能量范围内有多少态、由哪些原子轨道贡献。
能带线条很多并不代表判断对象变多,关键仍是能量窗口和轨道来源。费米能级附近、VBM 附近、CBM 附近通常承担不同物理含义:前者与低能激发和金属性有关,后两者与光吸收、空穴或电子传输有关。能带分析要先锁定目标能区,再比较色散、带边位置和投影轨道。
三、如何划分材料占据态、未占据态与禁带
在零温近似下,电子按能量从低到高填充。已占据能带中最高的一组常称为价带,未占据能带中最低的一组常称为导带。半导体或绝缘体的价带顶 VBM 与导带底 CBM 之间没有电子态,能量差就是带隙 Eg。带隙描述的是已占据态与未占据态之间的能量间隔,不是单独某条曲线的高度。VBM 和 CBM 指的是带边极值点,它们的轨道组成会影响光学跃迁和载流子有效质量。
若 VBM 和 CBM 位于同一个 k 点,材料具有直接带隙;若二者位于不同 k 点,则为间接带隙。直接带隙通常更利于动量守恒的光吸收或辐射复合,间接带隙材料的跃迁还涉及声子参与。直接或间接是带边位置关系,不能只凭带隙数值大小决定。

图4. Ca3Y2Ge3O12 在零压力下沿高对称方向的能带结构。DOI:10.3390/cryst13010029
在金属体系中,费米能级附近是否存在可用电子态比带隙概念本身更直接。若有能带穿过 EF,或价带和导带发生重叠,费米能级附近存在可占据的空态和已占据态,电子在微小外场下就能改变占据分布。金属性首先看 EF 附近是否有能带穿越或有限态密度,再讨论有效质量、散射和输运各向异性。

图7. In3Sc 立方相的能带和 DOS 中,多个能带跨过费米能级。DOI:10.3390/cryst15110946
平带提供另一类信息。曲线很平时,能量对 k 的变化小,电子群速度和有效质量会呈现不同于强色散能带的特征,常与局域化轨道、几何阻挫或特定晶格连通方式有关。平带不是“没有电子态”,它常对应能量集中、态密度增大的电子结构特征。

图5. Mielke 模型与 Pb2Nb2O7、Tl2Nb2O7 中的准平带特征。DOI:10.3390/nano9060876
四、能带分析材料性质需保留哪些计算前提
能带能直接支持的判断主要包括金属/半导体/绝缘体分类、直接或间接带隙、带边所在 k 点、带的色散强弱、缺陷态是否进入禁带以及自旋或 SOC 造成的分裂。这些判断都绑定具体模型和计算方法,离开晶胞、泛函、自旋态、U 值、SOC、应变和缺陷浓度后,数值不宜跨体系搬用。
缺陷和掺杂会把新的局域态引入禁带,或让原有价带、导带发生轨道混合。TiO2 掺杂模型中的能带变化说明,杂质态的位置、是否孤立、是否与宿主 O 2p 或 Ti 3d 态混合,都会影响吸收边与载流子复合倾向。禁带中的孤立态与带边混合态对应不同物理后果,这类差异通常要和 DOS/PDOS、电荷密度分布、光吸收或实验谱学共同校验。

图6. TiO2 及 In 掺杂、氧空位相关模型的能带结构对比。DOI:10.3390/ma11060952
带隙大小也不能单独承担器件性能判断。光催化材料还涉及带边相对氧化还原电位的位置,热电材料要比较有效质量、谷简并和散射,电催化表面则更关注费米能级附近态密度、吸附态轨道杂化和界面电荷重新分布。能带给出电子态的能量骨架,具体性能还需对应到吸收、输运、反应路径或界面电场。
计算能带时,k 点收敛、截断能、真空层、偶极修正、磁矩初值和结构弛豫都会影响曲线细节。二维材料和 slab 模型还要检查真空能级参考,强关联氧化物要说明 U 值,重元素体系要比较 SOC。一张能带图适合回答“电子态允许在哪里、怎样随 k 变化”,当问题转向载流子寿命、反应能垒或真实电导率时,仍需相应的动力学、输运或实验信号给出约束。
在异质结、表面吸附和电化学界面中,局域电场、真空层厚度、偶极修正和电荷补偿方式会改变能级参照。若把体相能带、表面态和吸附态放在同一个能量标尺下比较,必须先统一真空能级或电势参照。能带位置的比较依赖同一参照系,否则带边移动可能只是零点选择差异,而非真实电子结构改变。