DFT计算能带结构：高对称路径选取与投影能带分析的实战经验

能带图是固体物理论文里的”身份证”，审稿人第一眼看能带形状判断材料是金属还是半导体、是直接带隙还是间接带隙。我做DFT计算能带结构十多年，见过太多因为K点路径选错导致能带图完全没有物理意义的情况——一条本该平滑的导带底出现锯齿状跳变，价带顶在不该出现的地方突起，都是路径选取和参数设置出了问题。

能带计算本身不复杂，VASP的标准流程就三步：结构优化→自洽计算（SCF）→能带计算（band）。但每一步都有容易踩的坑，而且这些坑往往要到画图阶段才暴露出来，那时候已经浪费了几天机时。

Kohn-Sham本征值与能带色散：为什么SCF和band必须分两步跑

DFT计算能带结构的物理基础是Kohn-Sham方程的本征值。在自洽计算中，VASP在布里渊区均匀K点网格上求解Kohn-Sham方程，得到收敛的电荷密度。能带计算则是在高对称路径上做非自洽计算——从SCF的电荷密度出发，沿路径取密集K点，计算每个K点的本征值，画出E(k)色散曲线。

为什么不能在自洽计算中直接沿高对称路径取K点？因为电荷密度的收敛需要布里渊区均匀采样。高对称路径上的K点分布不均匀，直接用它们做自洽计算，电荷密度收敛不到物理基态，后续的能带就是建立在错误的电荷密度之上。这就是SCF和band必须分两步跑的根本原因。

SCF的K点网格选择遵循收敛性原则。体相材料用Gamma中心网格，密度取决于体系类型：金属需要更密的K点（如15×15×15），半导体可以稍稀（如9×9×9）。KSPACING参数是一个更直观的控制方式，设为0.25-0.30Å⁻¹对应比较密的采样。表面slab模型在面内方向用和体相类似的密度，真空方向取1个K点。

band计算的K点路径取决于晶系。立方晶系标准路径Γ-X-W-K-Γ-L-U-W-L-K，六方晶系Γ-M-K-Γ-A-L-H-A，四方晶系Γ-X-M-Γ-Z-R-A-Z。这些路径不是随便定的，是Setyawan-Curtarolo或Bradley-Cracknell给出的空间群标准高对称点序列。用SeeK-path工具输入空间群号，自动生成路径和K点坐标，是最可靠的做法。

高对称K点路径的选取方法：从查表到自动化工具

能带计算的K点路径选取有一个原则：路径必须经过布里渊区所有高对称点，且高对称点之间的连线方向对应物理上有意义的色散方向。不同晶系的高对称点坐标不同，路径也不同，不能混用。

我遇到过最典型的错误是拿立方晶系的路径套到四方晶系上。立方晶系的X点在(0, 1/2, 0)，四方晶系的X点也在(0, 1/2, 0)，但四方晶系多了一个M点(1/2, 1/2, 0)，路径是Γ-X-M-Γ-Z-A-M-X-R-Z，比立方多了M点和对应的线段。如果用立方路径跑四方晶系，能带图会缺少Γ-M方向的色散信息，间接带隙可能被误判为直接带隙。

自动化工具推荐SeeK-path（materialscloud.org上在线版）和sumo（Python包）。SeeK-path输入CIF文件，输出标准路径的K点坐标和线段名，直接可以写成VASP的KPOINTS文件。sumo不仅能生成路径，还能直接画能带图，自动标注高对称点位置。

K点路径上每段取多少个点也有讲究。太少（<20）能带曲线不平滑，会出现折线感；太多（>80）计算量增加但物理信息没有增多。我的经验是每段40-50个点，对于高对称路径总共有8-10段的体系，路径上总K点数约300-500个。

ISMEAR设置对能带形状的决定性影响

在DFT计算能带结构的实践中，ISMEAR是最容易被忽视但对能带形状影响最大的参数之一。

金属体系用Methfessel-Paxton方案（ISMEAR=1或2），SIGMA取0.2eV。这个方案在费米面附近用高斯展宽处理电子占据，避免费米面的数值震荡。金属体系如果误用ISMEAR=-5（tetrahedron方法），费米面附近会出现锯齿状震荡，能带在费米能级附近来回跳变，看起来像打开了一个假带隙。

半导体和绝缘体用tetrahedron方法（ISMEAR=-5），不设SIGMA。tetrahedron方法不做展宽，能给出最精确的总能量和态密度。半导体如果误用ISMEAR=1，费米面附近出现非物理电子占据，导带底被部分占据，带隙被压小甚至消失——这是我见过的带隙为零最常见的非物理原因。

有一个特殊情况：结构优化阶段半导体体系不能用ISMEAR=-5，因为tetrahedron方法不能计算力和应力。结构优化时用ISMEAR=0（Gaussian展宽），SIGMA=0.05eV，优化完再切到ISMEAR=-5做静态SCF。这个切换不是可选项，是VASP的硬性要求。

投影能带（Fat Band）的绘制与轨道分析

普通能带图只显示每个K点的本征值，看不出各轨道的贡献。投影能带（fat band）通过LORBIT=11参数输出每个原子、每个轨道（s、p_x、p_y、p_z、d_xy等）在每条能带上的投影权重，用颜色深浅或线条粗细展示轨道成分。

投影能带的核心价值在于分析轨道杂化。做TiO₂的DFT计算能带结构时，价带顶主要是O 2p轨道贡献，导带底主要是Ti 3d轨道贡献——这个信息从总能带图看不出来，从投影能带上一目了然。如果要做更细致的d轨道分析（区分t₂g和e_g），需要LORBIT=11配合LMMIX=.TRUE.，让VASP输出分轨道投影矩阵。

用sumo或pymatgen读取PROCAR文件画投影能带图，可以指定要投影的元素和轨道。画图时有一个细节：投影权重的归一化方式。pymatgen默认对每条能带做归一化（每个K点每条能带的投影权重之和为1），这种方式适合分析轨道成分比例。如果要看绝对投影量（某轨道在能带上的绝对贡献），需要关掉归一化。

HSE06能带计算也能输出投影信息，但需要额外设置LORBIT=11。HSE06的PROCAR格式和PBE略有不同，某些版本的pymatgen解析时可能出错。我的做法是HSE06能带计算后用vaspkit解析PROCAR，兼容性更好。

能带计算常见异常的排查与修复

能带曲线出现锯齿状震荡——最常见原因是SCF电荷密度没有充分收敛。检查SCF的OUTCAR最后几步能量变化，如果还有0.1meV量级波动，说明EDIFF不够紧。把EDIFF从1E-6调到1E-7或1E-8重跑SCF。另一个原因是K点路径上取点太少，增加到每段50-60个点。

带隙为负值——对金属体系正常，但如果体系预期是半导体，首先查ISMEAR。ISMEAR=0或1会压小带隙，改成ISMEAR=-5重跑SCF。其次查结构是否充分优化，晶格参数偏差大于1%会导致能带拓扑变化，间接带隙可能变成直接带隙或反过来。

能带图在Γ点出现非物理简并——查K点路径是否正确。某些高对称点在特定空间群下有额外的简并性，如果路径经过的不是真正的高对称点，简并会被错误解除。用SeeK-path重新生成路径。

HSE06能带和PBE能带形状差异大——查HSE06的SCF是否收敛。HSE06的SCF收敛比PBE困难，ALGO=All或ALGO=Damped可能需要更多电子步。如果HSE06的SCF没收敛就去做band计算，能带形状会明显异常。检查OUTCAR中”reached required accuracy”是否出现。

在DFT计算能带结构的实际项目中，能带计算还有一个不常被提及的坑：POTCAR中价电子构型的选择。Ti的PAW势有Ti（3p⁶3d²4s²，10个价电子）和Ti_sv（3p⁶3d²4s²，实际上是把3p半芯态纳入价态，12个价电子）两种。Ti_sv算出的能带在费米面附近的Ti 3p-3d杂化更完整，带隙通常比Ti高0.1-0.3eV。同一个项目里必须统一POTCAR版本，混用会导致能带不可比。

能带计算的经验沉淀：三步走流程的可靠性

做了这么多能带计算，我总结出一套三步走的标准流程，在95%以上的体系中都能给出可靠结果。

第一步PBE结构优化+静态SCF，拿到能带骨架。PBE的能带形状——导带底位置、价带顶位置、直接/间接带隙类型——是足够准的，带隙绝对值偏小是已知系统误差。先确认PBE能带形状合理，没有锯齿、没有非物理简并解除，再进入下一步。

第二步根据体系特征选择修正方案。含d电子的氧化物用DFT+U，sp半导体用HSE06，层状材料加范德华修正。修正方案选定后，从PBE的CONTCAR和CHGCAR出发做修正计算，不要从头重新优化结构。

第三步画投影能带做轨道分析。LORBIT=11输出PROCAR，用sumo或pymatgen画fat band图，分析价带顶和导带底的轨道成分。这一步是能带计算最有物理意义的部分——带隙数值本身只是一个数字，轨道杂化信息才是理解材料电子结构的钥匙。

能带计算的可靠性建立在每一步的严谨之上。结构没收敛就跑SCF，SCF没收敛就跑band，K点路径选错，ISMEAR设错——任何一个环节出问题，最终的能带图都会以某种方式”看起来不对”。养成跑计算前过一遍参数清单的习惯，比出问题后再排查高效得多。