NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

引言

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

很多 VASP 初学者第一次遇到 NBANDS,是因为程序提示带数不够,或者光学性质算出来只有一段很短的谱。NBANDS 表面上只是“算多少条能带”,实际影响的是未占据态数量。对于普通结构优化,它常常被自动设置得够用;但一旦进入 DOS、光学吸收、介电函数、能带展开或 GW/HSE06 前处理,空带不足就会变成结果风险。

NBANDS 设少,本质上不是少输出几条曲线,而是把可参与电子跃迁和后处理分析的未占据态截掉了。这会让 DOS 的导带部分、光吸收高能区和空态相关计算同时受影响。

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

NBANDS 到底控制什么?

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

为什么它不只是输出带数?

VASP 在平面波基组中求解 Kohn-Sham 能级。占据态负责总能、电荷密度和力,未占据态虽然不直接贡献基态电荷,却会进入 DOS 高能区、光学跃迁、介电响应和许多后处理。NBANDS 设定的就是每个 k 点求多少条能带,其中一部分是占据带,另一部分是空带。

占据带决定基态是否能收敛,空带决定你能分析到多高的未占据能区。所以结构优化能正常结束,不代表 DOS 和光学结果已经足够。

Nbands = Nocc + Nempty

式(1)中 Nocc 是占据带数量,Nempty 是空带数量。NBANDS 过小,最先被压缩的通常就是 Nempty,也就是导带以上可分析的能级空间。

图1 给出了 Sr₃SbCl₃ 和 Sr₃PCl₃ 的能带与 DOS。对这种半导体或绝缘体,导带态的位置、成分和密度决定后续光学吸收边与高能峰,NBANDS 太少就会让导带窗口不完整。

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

图1 Sr₃MCl₃ 钙钛矿的能带结构与 DOS。图中展示了价带、导带和态密度分布,是判断空带窗口是否足够的基础。DOI: 10.1021/acsomega.6c02217。

什么时候默认值不够?

如果只做普通离子步优化,默认 NBANDS 往往足以完成收敛。但做静态 DOS、LOPTICS、介电函数、BSE/GW 前置计算、空带投影或高能导带分析时,默认值就可能偏少。尤其是包含重元素、大超胞、掺杂缺陷或低维材料的体系,未占据态分布很密,自动带数只保证“能算完”,不保证“能解释”。

NBANDS 的最低要求和计算目的有关。优化、DOS、光学性质和空带计算不能共用同一个“够不够”的标准。

实际判断时,可以先看 OUTCAR 里的电子数、NELECT、NBANDS 和每个自旋通道的占据情况。若最高几条能带仍有明显占据,说明带数至少不能再减;若要分析费米能级以上 5-10 eV 的 DOS 或吸收峰,还要确认输出能带覆盖到这个能量范围。对自旋极化体系,上下自旋通道的占据带数量可能不同,更不能只按非自旋体系的经验给一个小数值。

NBANDS 设定前要先看电子数和目标能量窗口。只按原子数乘一个固定倍数,容易在大超胞、磁性体系和缺陷体系里出错。

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

NBANDS 少了,DOS 会怎样变?

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

为什么导带 DOS 会先出问题?

DOS 通常包含价带、带隙附近态和导带态。价带主要由占据态构成,受 NBANDS 少量变化影响较小;导带和更高能区则依赖空带数量。若 NBANDS 太少,DOSCAR 里高能区会提前结束,或者导带部分看起来很薄,后续再讨论导带态、缺陷空态、吸附诱导空态就会缺证据。

NBANDS 不足会让 DOS 的未占据部分被截断。如果你要解释导带、反键态或空态投影,只看价带收敛是不够的。

图3 展示了 Sr₂TiGeO₆ 的能带和主要轨道分辨 DOS。可以看到,导带边主要来自特定原子轨道贡献;若空带数量不够,导带轨道贡献会被截短,分峰和高能态来源就不可靠。

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

图2 Sr₂TiGeO₆ 的能带结构与主要轨道分辨 DOS。图中展示了带边态和导带轨道来源,说明空态窗口会影响 DOS 解释。DOI: 10.1039/d6ra01242f。

NEDOS 能不能弥补 NBANDS?

NEDOS 只控制能量网格有多细,不能凭空补出没有求解的能带。把 NEDOS 从 2000 提到 6000,只会让已有能带的 DOS 曲线更细,不会恢复被 NBANDS 截掉的空态。很多人把曲线画平滑后以为导带信息够了,其实只是把稀疏数据插得更密。

NEDOS 负责分辨率,NBANDS 负责能级范围。一个控制横向采样细度,一个控制可用能带数量,不能互相替代。

还有一种常见误判,是只把 DOS 图画到费米能级附近 2 eV,然后认为 NBANDS 已经足够。这个做法对吸附态附近的低能导带可能勉强够用,但如果后面要讨论光生电子去向、缺陷能级以上的空态、反键轨道占据或宽能区光谱,2 eV 窗口就太窄。更稳的流程是先根据研究问题确定图中要展示的能量范围,再反推需要多少空带。

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

为什么光学性质对 NBANDS 更敏感?

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

光吸收为什么需要空带?

VASP 的独立粒子光学响应可以理解为从价带到导带的带间跃迁。价带给出初态,导带和更高空态给出末态。若 NBANDS 太小,高能末态没有被求出来,那么介电函数虚部、吸收系数、能量损失函数等都会缺少对应跃迁贡献。

ε2(ω) ∝ Σ |⟨ψc|p|ψv⟩|² δ(Ec – Ev – ℏω)

式(2)中 ψv 是价带态,ψc 是导带空态。NBANDS 太少时,可参与求和的 ψc 不够,高能吸收峰自然会缺失或偏弱。

图2 给出了 Sr₃MCl₃ 的介电函数、吸收系数、反射率和能量损失函数。它说明光学结果是一组由介电函数推导出的联动性质,空带不足不会只影响一条吸收曲线,而会影响整组光学量。

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

图3 Sr₃MCl₃ 的介电函数、吸收系数、反射率和能量损失函数。图中展示了光学响应对带间跃迁信息的依赖。DOI: 10.1021/acsomega.6c02217。

为什么低能区正常不代表全部正常?

很多时候,NBANDS 偏少时低能吸收边仍然看起来合理,因为最低几条导带已经包含在计算中。但随着光子能量升高,需要更高的空态作为跃迁末态,高能区就会开始失真。若研究对象是可见光吸收边,可以重点收敛 0-5 eV;若讨论紫外吸收、能量损失或宽能区介电函数,就要把空带数提高到目标能区稳定。

NBANDS 收敛必须绑定目标能量范围。只看吸收边,不等于 10 eV 以上的光谱也可信。

因此,光学性质的 NBANDS 测试最好不是只看一个最终图,而是把不同 NBANDS 下的 ε₂(ω) 和吸收系数叠在一起。若吸收边位置、主要峰位和目标能区内的峰强都基本不变,才说明空带数可以接受;如果高能峰随着 NBANDS 增大持续出现新峰,就说明此前的空带窗口还没有覆盖完整跃迁。

图4 展示了 Sr₂B′B″O₆ 的吸收系数、消光系数、折射率和反射率。对这类宽能区光学常数,如果空带不够,高能峰位和强度都会受到影响。

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

图4 Sr₂B′B″O₆ 的吸收系数、消光系数、折射率和反射率。图中显示宽能区光学常数需要足够的空带支撑。DOI: 10.1039/d6ra01242f。

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

NBANDS 应该怎样收敛?

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

先按占据带数估一个起点?

一个实用做法是先查看 OUTCAR 中实际电子数和带数,再按占据带数的倍数增加 NBANDS。普通 DOS 可从 1.2-1.5 倍占据带开始;光学性质常从 2 倍占据带起步;若要覆盖更高能区或做空态敏感后处理,可以继续测试 3 倍或更高。这里的倍数只是示例,最终要看目标能区的谱线是否稳定。

NBANDS 不是越大越好,而是要大到目标结果不再变化。空带过多会增加计算量,空带不足会损伤结果。

INCAR 参考示例可以怎样写?

# ===== NBANDS 收敛静态计算参考示例 =====

NSW = 0

ICHARG = 11

LORBIT = 11

NEDOS = 4000

NBANDS = 240 # 示例:按占据带数和目标能区测试

LOPTICS = .TRUE.

CSHIFT = 0.05

这组设置适合在结构优化之后做静态电子结构和光学输出。真正使用时,应保持同一结构、同一 KPOINTS、同一展宽和同一能量窗口,只改变 NBANDS 做对比。不要一边改 NBANDS,一边改 k 点或泛函,否则很难判断变化来自哪里。

NBANDS 收敛测试一次只改一个主要变量。否则吸收峰、DOS 峰和导带窗口的变化无法归因。

如果体系很大,NBANDS 成本会上升很快,可以先用较粗 k 点做空带范围预筛,再对最终 k 点做一组较少数量的确认测试。这个做法不能替代最终收敛,但能帮助判断大概需要多少空带。对于二维材料,还要同时留意真空层引入的高能自由电子态;这些态可能让高能 DOS 变得很密,进一步增加空带需求。

ΔA(E) = |ANB2(E) – ANB1(E)|

式(3)可以作为谱线差异的简单检查量。A(E) 可以是吸收系数、ε2(ω) 或 DOS。若提高 NBANDS 后目标能区内 ΔA(E) 很小,才说明该能区基本收敛。

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

哪些结果最容易被误判?

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

空态、反键态和光谱峰要特别小心?

在催化电子结构里,吸附物反键态、导带缺陷态、未占据 d 态和光诱导跃迁都可能落在费米能级以上。如果 NBANDS 太少,你可能看不到这些态,或者只看到一部分。结论就会从“没有空态参与”变成“没有算到空态参与”,这两句话差别很大。

没有看到未占据态,不等于未占据态不存在。先确认 NBANDS 和能量窗口足够,再讨论轨道参与或光学跃迁缺失。

图5 给出了 KSnI₃ 多晶型的介电函数和吸收系数。不同结构的吸收边与高能峰不同,说明比较多个体系时必须统一 NBANDS、KPOINTS、NEDOS 和展宽参数,否则差异可能来自参数而非材料本身。

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

图5 KSnI₃ 多晶型的介电函数与吸收系数。图中展示了不同结构对光学响应的影响,比较时需要统一空带和采样参数。DOI: 10.1021/acsomega.6c00784。

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

本章要点总结

NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析
NBANDS设少了会怎样?VASP 未占据态截断引发的系列计算误差剖析

要点:NBANDS 决定每个 k 点求解的能带总数,其中空带数量影响未占据态分析。

要点:DOS 的价带部分可能正常,但导带和高能空态会因 NBANDS 太少而被截断。

要点:光学性质依赖价带到导带空态的跃迁,NBANDS 不足会影响介电函数、吸收系数和能量损失函数。

要点:NBANDS 收敛要绑定目标能区,并保持结构、k 点、展宽和后处理设置一致。

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
赞 (0)
上一篇 2026年7月10日 下午3:02
下一篇 2026年7月13日 下午2:58

相关推荐