能带结构计算

能带结构计算是固体物理和材料计算的身份证。我第一次给导师交能带图，被一句话打回：你这路径怎么选的？Gamma到X之间断了。那一刻才意识到，能带结构计算最容易被忽视的，不是算不算得动，而是高对称路径选得对不对、连续不连续。一张断掉的能带图，在懂行的人眼里等于自报家门——路径都没搞明白。从那以后我养成习惯，能带图先查路径再查数据。

一、背景：能带回答了什么能带结构计算

能带结构描述电子能量随晶体动量k的变化，直接告诉你是金属、半导体还是绝缘体，带隙多大，有效质量多少。做光电材料看带边位置，做拓扑看能带反转，做输运看费米面。可以说，能带结构计算的质量，决定了你后续物性文章的立论根基。审稿人看能带图，第一眼就判断路径连不连续、带边标得清不清楚，这些细节直接关联可信度。我见过一篇稿子，数据漂亮但能带路径错俩点，直接被拒，可惜在细节。

二、核心原理：k路径与Bloch定理

布里渊区里只有高对称方向上的能带才有物理意义且便于对比。标准做法是先自洽得到电荷密度，再沿高对称k点连成的路径做非自洽计算。路径怎么取？用SeeK-path或vaspkit自动生成最稳，手动写容易在对称点衔接处漏一段，能带图就断开。我早年手敲KPOINTS，把一个六方体系漏了K到H段，导致文章里声称的直接带隙其实是间接带隙，差点闹笑话。布里渊区形状由空间群决定，不同空间群的高对称点和连线完全不同，不能套用立方系的模板。路径密度也要够，我一般每两点间插30-40个k点，否则能带看着断。

三、关键技术要点：带隙偏小怎么救

PBE泛函系统性低估带隙，这是教科书级的问题。金属体系无所谓，但半导体/绝缘体若带隙报出来比实验小一半，审稿人一眼就盯上。我的经验：先确认结构收敛，再上HSE06杂化泛函做单点修正，带隙能拉回实验值80%以上。HSE的计算量是PBE的数倍，所以策略是PBE跑全路径，HSE只在带边附近几个k点做修正。记得设LHFCALC=.TRUE.和HFSCREEN=0.2，别漏了PRECFOCK=Fast平衡精度与速度。强关联体系（如过渡金属氧化物）光HSE不够，还得加DFT+U，否则d带位置照样错。杂化比例也可调，但默认0.25最常用。

四、实操流程：从SCF到能带图

结构优化→SCF（ICHARG=2写CHGCAR，KPOINTS用常规密度）→复制目录，KPOINTS换高对称路径、ICHARG=11跑NSOC→读EIGENVAL和OUTCAR画图。我习惯在SCF步就把LCHARG=.TRUE.打开，后面NSOC直接吃CHGCAR，避免重复自洽。路径文件建议存成可复用模板，不同空间群对应不同路径，别每回都现找。画完图记得标出费米能级、高对称点字母和带边位置，三要素缺一个都不算完整。有效质量从带边曲率二阶导取，别靠眼睛估。能带结构计算

五、常见问题

Q1：能带图断开不连续？高对称路径端点没衔接，用工具生成连续路径。 Q2：带隙比文献小？上HSE或者检查结构是否真优化到位。 Q3：费米能级穿过能带？正常，金属就是这样，别误判成错误。 Q4：spin-polarized体系能带翻倍？设ISPIN=2，注意磁矩初值。 Q5：SOC怎么开？加LSORBIT=.TRUE.，重原子体系必开，计算量翻倍。 Q6：间接带隙误判成直接？核对价带顶和导带底所在k点，别只看最小值。 Q7：HSE不收敛？降ALGO或加NBANDS冗余，别硬扛。

讲一个我踩过的路径坑：六方体系的高对称路径里K到H那段，手动敲KPOINTS时最容易漏，我当年就漏过，导致文章里声称的直接带隙其实是间接带隙，差点闹笑话。现在一律用SeeK-path或vaspkit自动生成，再人工核对端点字母。路径密度也要够，我一般每两点间插30到40个k点，否则能带看着断、带边读不准。带隙系统性偏小是PBE的教科书级问题：金属无所谓，但半导体若带隙报出来比实验小一半，审稿人一眼盯上。我的策略是先确认结构收敛，再上HSE06只在带边附近几个k点做修正，带隙能拉回实验值八成以上，计算量可控。强关联体系光HSE不够，还得加DFT+U，否则d带位置照样错。画完图别忘标三要素：费米能级、高对称点字母、带边位置，缺一个都不算完整图。有效质量从带边曲率二阶导取，别靠眼睛估。遇到自旋轨道耦合体系要开LSORBIT，重原子必开，计算量翻倍但位置才对。

再说路径自动生成后的一个小陷阱：不同工具对空间群判定的容差不同，偶尔会把同一结构判成不同空间群，路径字母跟着变。我遇到过vaspkit和SeeK-path给的路径字母不一致，结果图上的标记对不上文献。现在的做法是生成路径后，把高对称点坐标手动和材料数据库（如Materials Project）比对一遍，确认字母与坐标一一对应再画图。带边位置的读取也有讲究：间接带隙要分别找价带顶和导带底所在的k点，别只看能量最小值；直接带隙则是同一k点。我习惯把这两点坐标都标在图上，审稿人一看就知道你区分清楚了。还有，HSE计算的NBANDS要比PBE多留余量，否则带边附近能带数不够，修正后的能带会缺带。

补充一个验证习惯：能带算完，我会抽一条已知能带（比如金属的sp带或绝缘体的价带顶）和文献或数据库比对形状，确认路径和方向没错。这一步几分钟，但能提前抓出路径端点写反、k点插错方向这类低级错误。能带这东西，画错一次可能埋雷到投稿才爆，前置半小时校验比事后撤稿划算太多。如果你手上的能带图老被质疑，先查路径端点和HSE设置，那两块最易出错。

补一句关于投影能带的：光看普通能带容易漏掉能带归属，我习惯用projwfc或VASP的投影能带（fatband）把每个原子/轨道的贡献画成带宽，能直接看出哪条带是d带、哪条是p带，对解释拓扑和简并特别有用。如果你手上的能带图被质疑，补一张投影能带往往能救回一城。还有，能带和DOS要互相印证：DOS上的峰应能在能带里找到对应色散，两者对不上说明路径或展宽有误。我交付前必把两者叠在一起看，这是最省事的互检。这套互检习惯让我避免过好几次把间接带隙写成直接的低级错误。

六、复盘总结

回过头看，能带结构计算最耗人的从来不是机器，而是路径和泛函选择上那点经验差。我把常用空间群的k路径都整理成了模板，新体系直接套，省下的时间够多跑三个项目。能带这东西，画错一次可能埋雷到投稿阶段才发现，所以前置多花半小时校验路径，比事后撤稿划算得多。需要现成的高对称路径模板和HSE输入示例，可以直接找我们要，少踩我当年踩过的坑。