HSE06杂化泛函能带计算优化：精度提升与计算成本平衡策略

PBE算出的带隙偏小是行业共识。Si的实验带隙1.17 eV，PBE给出约0.55 eV——差了超过一倍。这不是某个体系的偶然，而是半局域交换关联泛函在描述电子-空穴相互作用时的系统性缺陷。HSE06（Heyd-Scuseria-Ernzerhof, 2006）通过混入一部分精确Hartree-Fock交换能，部分修复了这个问题。

代价也很直接：HSE06的计算时间是PBE的50-100倍。一个在PBE下1小时能跑完的Si原胞能带计算，用HSE06可能需要3-5天。所以问题的核心不是”要不要用HSE06″，而是”什么时候必须用，怎么用才能活下来”。

HSE06做了什么PBE做不了的事

PBE的交换关联泛函完全基于局域电子密度和它的梯度，天然忽略了一个关键物理：当电子从占据态跃迁到空态时，电子和留下的空穴之间存在库仑屏蔽。这个屏蔽减弱了电子-空穴对的束缚能——而在半导体中，束缚能就是带隙的”补丁”。

HSE06的处理方式很聪明：把交换能拆成短程和长程两部分。短程部分（r < ω⁻¹，ω是屏蔽参数，默认0.11 bohr⁻¹）混入25%的精确HF交换；长程部分仍然用PBE。混入的比例（α=0.25）和屏蔽长度（ω=0.11）是两个核心可调参数，不过绝大多数文献使用默认值，改动需要提供充分的物理理由。

这个设计不是随意选的。纯HF交换（100%混合）给出的带隙往往偏大（Si的HF预测约5 eV，实验1.17 eV），原因是HF完全缺失了对关联效应的描述。PBE给出的带隙偏小。HSE06 25%的混合比例是在大范围测试（包括约40种半导体）后确定的折中值，对sp键合半导体的平均带隙误差约0.2-0.3 eV，远优于PBE的0.5-1.0 eV。

对于含d/f电子的体系，HSE06和PBE+U的精度对比不是一边倒的。两种方法修复的是不同的物理缺陷：+U修复现场库仑排斥，HSE06修复交换关联的非局域性。对于强关联绝缘体（NiO、MnO），+U往往给出更好的磁性和带隙；对于sp半导体（GaAs、CdTe），HSE06明显优于+U。

VASP中HSE06的INCAR设置与计算策略

核心INCAR参数：

LHFCALC = .TRUE.     # 开启杂化泛函
HFSCREEN = 0.2       # 屏蔽参数（对应ω=0.11 bohr⁻¹）
ALGO = All           # 杂化泛函推荐All或Damped
TIME = 0.4           # All算法的阻尼参数
PRECFOCK = Normal    # Fock交换精度（Fast/Normal/Accurate）

ALGO的选择直接影响收敛速度。All（同时对角化所有KS轨道）比默认的Normal（逐带对角化）快2-3倍，但内存需求也大得多——大约需要O(N²)的内存。Damped是All的一个变体，在收敛边缘的体系中更稳定。实际选择：小体系（<100原子）用All，中等体系（100-300原子）用Damped，大体系（>300原子）重新考虑是否真的需要HSE06。

PRECFOCK控制Fock交换积分的精度。Fast版本对部分积分做了近似，对带隙的影响通常<0.05 eV，但计算量只有Normal的1/3-1/2。初期测试建议用Fast，关键结果用Normal做一次验证——这是文献审稿中最不敏感的精度差异。

四步接力策略：把HSE06的墙拆成台阶

HSE06直接怼上去且不说时间——收敛失败的概率也很高（Fock交换引入的非线性让SCF迭代比PBE敏感得多）。分步接力是降低失败率的核心策略：

第一步：PBE结构优化和SCF。 标准的PBE流程，输出WAVECAR。这是基础——HSE06继承PBE的电荷密度和波函数作为初猜。

第二步：PBE能带计算。 不是必需的，但可以作为基准对照。知道PBE的带隙偏多少，才能判断HSE06修正了多少。

第三步：HSE06 SCF + 粗K点。 用PBE的WAVECAR作为初猜，在比目标K点密度低一档的网格上做HSE06自洽。比如目标6×6×6，先用3×3×3收敛。这一步的目的是在低成本下得到HSE06的电荷密度，避免在细网格上从头迭代——后者可能不收敛。

第四步：HSE06能带计算 + 目标K点。 读取第三步的CHGCAR（ICHARG=11），在目标K点路径上做非自洽能带计算。这一步不再做自洽迭代，速度快得多（通常只需要SCF步骤时间的1/10）。

如果直接做第四步跳过第三步，HSE06在细网格上的SCF收敛可能需要几十甚至上百步——每一步的Fock交换计算都是靠时间堆出来的。第三步的粗网格收敛虽然也不快，但比细网格快一个数量级。

降低HSE06成本的其它技巧

K点缩减。 总能和能带对K点密度的敏感度不同。HSE06 SCF可以用比PBE低一档的K点密度（比如PBE验证收敛在6×6×6，HSE06 SCF用4×4×4），代价是带隙可能有0.02-0.05 eV的偏差——对于大多数横比分析来说，这个量级可以接受。

对称性利用。 HSE06对Fock交换的k点求和可以利用时间反演对称性减少约一半的k点数量。确保POSCAR中对称性被正确识别（ISYM=2），VASP会自动优化积分过程。

并行策略。 HSE06的并行效率比PBE好（因为Fock交换的k点循环天然可并行）。如果你的集群有多个节点，NCORE和KPAR需要针对杂化泛函重新调参——KPAR设成k点总数的约数（而不是NCORE的约数，这点和PBE相反）。

何时不需要HSE06

PBE+U已经够用的情况。 如果体系的带隙误差主要来自强关联d/f电子的局域化缺陷，+U比HSE06更对症。一个判断方法：PBE是金属、实验是绝缘体→+U优先；PBE是半导体但带隙偏小50%→HSE06优先。

相对趋势而非绝对值。 如果研究的是同系列材料中带隙的排序（比如ABO₃钙钛矿A位替换对带隙的影响），PBE的相对趋势通常已经正确——绝对值的系统偏移在横比中不构成误导。

mBJ势是一个轻量替代。 修改的Becke-Johnson（mBJ）势是一个meta-GGA级别的泛函，对带隙的预测精度接近HSE06（平均误差约0.2-0.3 eV），但计算成本只有PBE的1.2-1.5倍——几乎没有额外开销。mBJ不能做结构优化（它是势而非能量泛函），但做能带计算已经非常成熟。对于高通量筛选，mBJ是HSE06的最佳替代。

HSE06是当前带隙预测的精度的实用上限，但策略性地使用它比盲目全部切换要高效得多。科研学术网（https://www.keyanxueshu.com）的完整升级路径均有实战案例参考。