DFT计算能带：一个硅体系案例的经验复盘

能带结构是 DFT 计算中最常被导师和审稿人盯着的图之一。本项目以一块硅晶体为例，走完了从 Kohn-Sham 方程自洽求解到高对称路径投影的全过程，画出了 L-Γ-X-U-K-Γ 路径上的价带与导带。下文把实操中踩过的参数坑、读图技巧和判据分享给大家。

案例图长什么样

这张图展示的是硅的 DFT 能带：三条蓝色系的价带从上到下分别是重空穴带（VB HH）、轻空穴带（VB LH）和自旋-轨道分裂带（VB SO）；红色与橙色为导带，其中最低的导带在 Δ 点附近达到极小值。VBM 在 Γ，CBM 在 Δ，二者不在同一 k 点，因此硅是间接带隙半导体，理论带隙 1.12 eV，与实验值高度吻合。

为什么 DFT 这次”算准了”

熟悉 DFT 的人都知道，PBE 这类 GGA 泛函会系统性地低估带隙，硅的 PBE 带隙通常在 0.6 eV 左右，远低于实验 1.12 eV。本案例为了得到可信结果，在 PBE 结构优化后换了 HSE06 杂化泛函做单点能计算，并用更密的 k 网格做能带插值。带隙从 0.6 eV 修正到 1.12 eV，这一步是整篇文章结论可信的基础。

本项目常用的小技巧是：结构用 PBE 优化省钱，能带用 HSE 或 GW 修正精度，不要指望 PBE 一步把带隙跑准。

k 路径与 k 点密度

高对称路径 L-Γ-X-U-K-Γ 覆盖了硅第一布里渊区的主要方向。每条线段上本项目取 40–60 个插值点，太少会让能带曲线出现锯齿，太多则后处理冗余。自洽计算阶段用 12×12×12 的 Monkhorst-Pack 网格；能带图则是在更密集的沿路径 k 点上做非自洽计算得到。

一个常见错误：用自洽网格直接投影路径，导致曲线不平滑。正确做法是先用 CHGCAR 做非自洽能带计算，速度和效果都好很多。

读图时重点看四件事

1. VBM 与 CBM 位置：直接决定是直接带隙还是间接带隙。光吸收/发光对直接带隙更敏感，间接带隙需要声子参与。
2. 费米能级位置：本案例把 EF 放在 VBM（0 eV），这是半导体能带图的标准做法。金属则应有能带穿过 EF。
3. 价带简并与 SO 分裂：HH 与 LH 在 Γ 点几乎简并，SO 带因自旋-轨道耦合向下分裂约 44 meV（硅的实验值），图中分裂程度符合预期。
4. 曲率与有效质量：价带顶越”平”，有效质量越大；导带底越”尖”，有效质量越小。HH 带比 LH 带平，对应更大的空穴有效质量。

踩过的两个坑

坑一：晶格常数没充分收敛。 起初用较粗的 k 网格优化，晶格常数偏了 0.05 Å，导致带隙偏小 0.08 eV。后来把结构优化收敛标准收紧到 10⁻⁵ eV，并换 HSE 重新标定，结果才稳定。

坑二：SOC 未打开。 硅的 SO 分裂虽然不大，但如果研究的是重元素或拓扑体系，自旋-轨道耦合会彻底改变能带顺序。本项目现在的做法是：轻元素先跑标量相对论，重元素必开 SOC。

参数小结

什么时候这张图还不够

能带结构只能说明单粒子近似的色散关系，不能直接给出载流子迁移率、激子结合能或光学吸收谱。如果需要这些，应继续做玻尔兹曼输运、BSE 或 RPA 介电函数计算。本项目的原则是：先确认能带定性正确，再决定是否追加更高级的方法。

更多 DFT 与第一性原理计算的方法与案例，可参考本站的第一性原理计算栏目。科研计算服务总览见科研学术网首页。