能带模拟计算：三大主流软件的赝势路线与收敛策略

同一个体系用不同软件算能带，结果差出0.3 eV的带隙——这种事情在计算材料圈子里不算新闻。能带模拟计算的软件选择，看起来是个操作偏好的问题，背后其实是赝势方案、基组质量和自洽策略三条技术路线的分歧。一个项目在初期如果没把这个选型问题搞清楚，后面所有参数的微调都可能是在给方向偏差填坑。

赝势策略的分野：USPP vs PAW vs NCPP

VASP走的是PAW路线，这也是目前公认精度最高的赝势方案。PAW在原子核附近保留了一个”增强区域”，通过投影算符还原全电子波函数的信息，对磁性体系（过渡金属、稀土元素）和核区敏感量（EFG、超精细耦合）的精算价值无可替代。代价是PAW的数据集制备复杂，用户基本只能依赖VASP官方发布的版本。

CASTEP默认走超软赝势（USPP），可选模守恒（NCPP）。USPP的优势在于截断能要求低——氧化物体系（ZnO、TiO₂）在PAW需要400 eV才能收敛到1 meV/atom的地方，USPP可能300 eV就够了。但对于含3d电子的过渡金属体系，USPP容易在芯-价交换关联中出现伪像，这是CASTEP用户最常踩的坑。

Quantum ESPRESSO（QE）走的是模守恒赝势加SSP（标准固态赝势）路线，核心数据来自Pslibrary和SSSP效率库。QE的优势是赝势完全开源可审计——用户能追溯每个赝势的生成参数和验证标准。但在稀土和锕系元素的精算中，QE的模守恒赝势对d/f电子的处理需要更严格的验证，平面波截断能往往推到极高值（>100 Ry，即~1360 eV）才收敛。

k点收敛策略的差异：Gamma-centered vs MP

VASP默认Γ-centered网格，QE默认Monkhorst-Pack。这个差异在六方晶系中产生的影响最明显——Γ-centered在Γ-A路径的带边精度优于MP，因为六方晶体的对称性使得带边往往落在Γ相关的k路径上，MP可能恰好避开这些关键采样点。

CASTEP的k点采样同时支持两种方案，但需要用户手动指定。默认设置下CASTEP使用MP栅格，切换为Γ-centered后六方GaN的带隙可以从2.85 eV跳到2.88 eV（PBE级别）——30 meV看起来小，但能带对齐和异质结建模里，这个偏移量足以改变界面态的能级排布。

实际项目中，真正花时间的不在k点网格本身，而在不同k点密度对总能和Fermi能的收敛性测试。每增加一个k点维度（如4×4×4→6×6×6），总的SCF迭代步数可能增加20-30%，因为更大的k点集需要更精细的电荷密度混合。这是能带模拟计算中时间成本的隐性大头。

并行效率与体系规模的匹配

VASP在CPU上的并行效率已经到了相当成熟的阶段，对20-100原子的体系，8-16核基本线性加速。但在GPU加速下，VASP的FFT变换会成为新瓶颈——GPU的显存带宽限制了散射项的计算吞吐。对含稀土元素的重体系（如CeO₂的4f电子），VASP的CPU版本仍是最稳的选择。

QE对大规模体系的并行扩展性优于VASP——QE的线性标度DFT算法（通过在适当截止半径内限制轨道重叠实现）在几百到上千原子的超胞计算中，墙钟时间是VASP标准对角化方法的1/3到1/5。这个优势在超胞能带展开（band unfolding）和缺陷超胞分析中特别有用。

CASTEP的并行方案跟MS平台深度绑定，优势不在纯算力而在前后处理一体化——从建晶胞到出能带图全在一个窗口完成。对需要快速产出能带结果的对比分析场景（如多相能带对齐），操作效率高于算力效率。

选型结论

氧化物、简单半导体、金属体系：VASP的PAW方案精度最高，是首选。但入门门槛和许可证费用不低。

大规模超胞、缺陷、非晶：QE的赝势透明性和并行扩展性更优，适合需要大量试算和参数扫描的项目。

教学、快速预研、Materials Studio生态用户：CASTEP的集成工作流和图形界面降低了配置出错概率，但对高精度要求的磁性体系需要额外的赝势验证。

三者没有绝对优劣——选什么取决于体系对近核区精度的要求、超胞规模、以及计算结果后续用途。能带模拟计算的终极标准，是用选定的工具算出来的物理预测，跟实验落差在可控范围内。

更多内容请访问 https://www.keyanxueshu.com/