能带理论计算

能带理论计算，是每一个做凝聚态和材料模拟的人绕不开的基本功，却也是最容易被”算出来却读不懂”的一步。我当年第一次自己画能带图，高对称路径照着文献抄，画出来带隙是负的，还硬着头皮写了”金属性”的结论，后来被导师一句”你路径取错了”打回重算，才真正开始琢磨能带背后到底在算什么。能带理论计算难的从来不是跑完那一步，而是你选的高对称路径、k 点密度和泛函，到底能不能把材料的电子结构讲清楚。

一、背景：为什么能带是第一性原理的”门面”能带理论计算

能带直接决定材料是金属、半导体还是绝缘体，决定带隙多大、载流子有效质量多少、能不能被光激发。做催化剂的人看 d 带中心判断吸附强弱；做热电的人看能带简并度判断功率因子；做拓扑材料的人看能带翻转判断拓扑不变量。可以说，能带理论计算是你跟审稿人、跟合作者解释”这个材料为什么好用”的第一张图。很多新手把能带当成黑箱输出，画完就交差，等到别人问”为什么你的带隙跟实验差 1 eV”，才发现自己泛函选错了都没察觉。更现实的是，带隙和有效质量直接喂给输运和光学后续计算，源头错了后面全错。我见过一个合作项目，能带带隙算错导致光吸收波长预判偏了 100 nm，下游器件方案全推倒重来——能带这张图，从来不是装饰。我带过的几个学生里，最常见的翻车就是拿 PBE 的带隙直接往文章里写，等 reviewer 问”你这半导体带隙比实验小一半”才慌。带隙是能带图里最显眼的数，也是最容易错的那个，所以我现在要求组里出能带图必须先标泛函，PBE 的一律在图注里写”带隙未修正”，免得被当成定论拿去支撑结论。

二、核心原理：从 Bloch 波到能带折叠

第一性原理能带是在倒空间对 Bloch 波本征值 Eₙ(k) 的采样：对每个 k 点解 Kohn-Sham 方程得到一组能级数。关键认知是，能带图的横轴必须沿着布里渊区的高对称路径（如 Γ–X–W–K–Γ），否则你画出来的”带”在物理上没有连续意义。费米能级把占据态和未占据态切开，带隙就是价带顶到导带底的最小能量差；若这个差出现在同一 k 点就是直接带隙，不在同一 k 点就是间接带隙——这决定了光学跃迁要不要声子参与。有效质量则从能带曲率 1/m* = (1/ħ²)·d²E/dk² 读出，曲率越平质量越大、载流子越”重”。还有一个点：对超胞或表面体系，能带会沿高对称方向”折叠”，看起来带变密了，其实是小胞能带在新 BZ 里的重复投影，读图时要先还原折叠关系，否则有效质量和带隙都会读错。还有一点容易被忽略：拓扑非平庸的体系，能带在 Time-reversal 对称的 k 点（如 Γ）会发生带反转，看着像普通的能带交叉，其实是拓扑保护态的标志。这类体系如果只画普通能带不标拓扑不变量，很容易把拓扑特征当成普通简并漏掉。我算 Bi₂Se₃ 时就是先看到 Γ 点附近的带反转，才顺着去算 Z₂ 不变量确认拓扑绝缘体身份的。

三、关键技术要点：路径、k 点与泛函

高对称路径要用 SeeK-path 或 Materials Project 核对，别凭记忆手敲，路径断了带就断了。k 点密度要够，半导体一般 0.02–0.03 Å⁻¹ 间距，金属更密。泛函选择是分水岭：PBE 系统性低估带隙（常差 30%–50%），对强关联体系（如过渡金属氧化物）几乎失效；HSE06 杂化泛函能拉回大部分带隙但机时陡增；DFT+U 专治 d/f 电子强关联，U 值要按体系标定。自旋轨道耦合（SOC）对重元素（Pb、Bi、Pt）不可忽略，会让能带劈裂几十到上百 meV。我早年算 ZnO 用裸 PBE 给的带隙 0.8 eV，跟实验 3.4 eV 差出一个量级，加 HSE06 才救回来。另一个常被忽略的点：SOC 和 HSE06 同时开时计算量翻倍再翻倍，对含重元素的拓扑材料才值得，普通半导体不必。

四、实操流程：从结构到能带图能带理论计算

第一步，结构充分优化，晶格常数弛豫到应力收敛；第二步，静态自洽（ISMEAR 对金属用 -1 或 1，半导体用 -5），得到 CHGCAR；第三步，读取高对称路径，生成线性的 KPOINTS（或用 pymatgen 自动生成）；第四步，非自洽能带计算（ICHARG=11），提取 EIGENVAL；第五步，pymatgen/sumo 画图，按元素的 s/p/d/f 投影做 fat band，把特征能带按轨道着色。我习惯先画一条 PBE 能带确认路径和费米面没问题，再换 HSE06 修正带隙——直接上 HSE06 如果路径错了会白等几天机时。对磁性体系还要先定磁序，不同磁构型能带形状天差地别。对异质结或表面，要先用真空能级把两侧费米能级对齐，再谈能带偏移，否则”偏移”是假的。

五、常见问题

Q1：画出来带隙是负的？路径取错或结构没优化，先查高对称点。 Q2：带隙跟实验差很多？PBE 低估，换 HSE06 或加 GW。 Q3：能带在路径接头处断开？KPOINTS 路径不连续，核对 SeeK-path。 Q4：重元素体系带劈裂不明显？没开 SOC，加上再算。 Q5：过渡金属氧化物带隙几乎为零？强关联，上 DFT+U 或 HSE06。 Q6：fat band 看着没特征？投影轨道选错，换 s/p/d 组合看。 Q7：有效质量算出来离谱？曲率采样点太少，加密路径局部 k 点。

讲一个具体的坑：有次算一个钙钛矿，PBE 给的直接带隙，实验却是间接带隙，我一度怀疑结构。后来在路径上补了更多 k 点，发现导带底其实在另一个高对称点，是路径太稀把间接带隙的”鼓包”漏掉了。从那以后我看带隙类型必先确认路径覆盖了全部高对称点，而不是只看 Γ 附近。另一个常被忽略的点：费米能级对齐，异质结能带对齐前要统一真空能级参考，否则”偏移”是假的，我交付能带图一定先在图里标清楚费米能级位置。对高通量筛选，我不会每个材料都手动画能带。先用自动流水线跑 PBE 能带抓候选，再对前几名上 HSE06 精修带隙，这样机时分配最合理。pymatgen 的 Boltztrap 还能直接从能带插值出输运系数，把能带和有效质量直接喂给热电功率因子估算，一环扣一环。

六、复盘总结

回过头看，能带理论计算最容易被当成”画一条线”的体力活，但它背后每一步——路径、k 点、泛函、自旋、对齐——都在影响你最终敢不敢拿这张图去支撑结论。我现在的习惯是先 PBE 跑通结构确认路径，再上 HSE06 谈带隙，最后用 fat band 把物理讲透。能带不是装饰图，它是你整个电子结构论证的骨架，骨架歪了后面全歪。把路径和泛函这两根支柱钉死，带隙和有效质量才敢写进文章和方案。这些年我审过不少学生的能带图，最低级的错误永远是路径和泛函——路径抄错导致带隙类型读反，泛函选错导致带隙偏低。把这两件事钉死，至少能避开八成的返工，剩下的两成才是 fat band 和有效性质量这些精细活。