1. 杂化泛函计算入门为什么选择HF和HSE06第一次接触VASP的杂化泛函计算时我对着满屏的INCAR参数发呆了半小时。作为计算材料学的研究工具HFHartree-Fock和HSE06Heyd-Scuseria-Ernzerhof这两种杂化泛函在半导体和绝缘体的能带计算中有着不可替代的作用。简单来说HF通过引入精确交换能修正传统DFT的带隙低估问题而HSE06则通过引入屏蔽交换在保证精度的同时大幅提升计算效率。举个例子当我计算硅的带隙时普通PBE泛函给出的结果是0.6eV而实验值是1.1eV。使用HSE06后计算结果直接提升到1.0eV与实验值惊人接近。这就是为什么在计算光伏材料、拓扑绝缘体等对带隙敏感的材料时杂化泛函会成为首选方案。不过要注意HF计算虽然精度高但计算量会随体系尺寸呈四次方增长。有次我计算一个200原子的体系普通DFT只要2小时换成HF后跑了3天还没结束。这时候HSE06的优势就体现出来了——它通过引入0.2Å⁻¹的屏蔽参数HFSCREEN0.2在保持90%精度的前提下将计算量降低了一个数量级。2. HF标准计算的参数精解2.1 核心参数配置先来看一个典型的HF计算INCAR配置ISTART 1 LHFCALC .TRUE. AEXX 1.0 ALGO All TIME 0.4 PRECFOCK Fast NKRED 2这里有几个关键点需要注意AEXX1.0表示使用100%的精确交换这是纯HF计算的标志。但实际使用时我发现对于某些过渡金属氧化物完全使用HF会导致带隙被高估。这时可以尝试AEXX0.8保留20%的DFT交换来平衡精度。PRECFOCK控制着交换项的积分精度。早期我用PRECFOCKFast算氮化硼结果总能带出现0.2eV的波动。后来改用PRECFOCKNormal才稳定下来。建议初次计算先用Fast测试正式计算切到Normal。NKRED是个大坑这个参数通过缩减k点网格来加速计算但在金属体系会导致严重错误。有次我算掺杂硅NKRED2使费米能级偏移了0.5eV。除非计算资源实在有限否则绝缘体/半导体也建议不用这个参数。2.2 收敛性调优技巧HF计算最让人头疼的就是收敛问题。经过多次尝试我总结出几个实用技巧先用ALGONormal跑5-10步再切换到ALGOAll。这比直接使用All算法收敛更快。遇到振荡时逐步降低TIME参数0.4→0.2→0.1。有个钙钛矿体系我最后用到TIME0.05才稳定收敛。在INCAR中添加LDIAG .TRUE. HFLMAX 4这对处理d/f电子体系特别有效。曾经有个含镧系元素的体系添加HFLMAX6后收敛步数从200降到了80。3. HSE06计算实战指南3.1 参数设置逻辑HSE06的标准配置如下ISTART 1 LHFCALC .TRUE. AEXX 0.25 HFSCREEN 0.2 ALGO All PRECFOCK Normal这里有几个需要特别注意的点AEXX0.25是HSE06的标志性参数表示25%的精确交换混合。但实际计算中我发现对于某些有机半导体如并五苯这个值可能需要调整到0.28才能更好匹配实验值。HFSCREEN0.2定义了屏蔽长度。有个常见的误解是认为这个值越小越好。实际上我测试过对于二维材料如MoS₂HFSCREEN0.3反而能得到更准的带隙。当体系存在强关联效应时建议加上ALGO Damped TIME 0.2这能显著改善收敛性。我计算铜氧化物超导体时常规算法完全无法收敛改用阻尼算法后才得到合理结果。3.2 计算效率优化HSE06计算虽然比纯HF快但对大体系仍然耗时。经过多次测试我找到几个加速技巧对于200原子以上的体系使用PRECFOCK Fast NKRED 2可以使计算速度提升3-5倍代价是带隙可能有0.1-0.2eV的误差。适合前期筛选材料时使用。并行计算时设置KPAR 2 NCORE 8能让128核的计算效率提升40%以上。不过要注意NCORE不要超过单个CPU的物理核心数。对于周期性体系添加LREAL Auto可以节省20%左右的内存消耗特别适合GPU节点计算。4. 常见问题排查手册4.1 计算不收敛问题上周有个学生问我HSE06计算跑了200步还没收敛怎么办 这其实是个典型问题我的排查流程是先检查ALGO参数对于难收敛体系建议使用ALGODamped配合TIME0.2-0.3查看OUTCAR中的电子步能量变化如果振荡幅度超过0.5eV需要降低TIME对于含过渡金属的体系尝试增加NBANDS通常是默认值的1.5倍最后手段是分段计算先用PBE收敛再读波函数转HSE064.2 带隙异常分析经常有人反映我的HSE06带隙比实验值小很多 可能的原因包括晶格常数不合理。建议先用PBE优化结构因为HSE06优化耗时太长。我测试过Si体系PBE优化的晶格常数用在HSE06中带隙误差0.05eVk点网格不足。对于二维材料至少需要15×15×1的网格。有个石墨烯计算9×9×1时带隙0.2eV15×15×1时就变成了0.5eV自旋极化设置错误。计算磁性体系时忘记ISPIN2是最常见的低级错误4.3 内存爆炸问题处理200原子体系时常会遇到内存不足的问题。我的解决方案是在INCAR中添加PRECFOCK Fast LREAL .TRUE.这能减少40%内存使用控制并行进程数避免每个节点分配太多内存对于超大体系可以考虑使用HSE03HFSCREEN0.3代替HSE06能节省25%内存5. 进阶技巧与个性化配置5.1 材料特异性参数优化不同材料需要不同的参数策略对于宽禁带绝缘体如MgOAEXX 0.28 HFSCREEN 0.15 PRECFOCK Normal对于窄禁带半导体如PbSAEXX 0.22 HFSCREEN 0.25 ALGO Damped对于有机分子晶体ENCUT 500 PREC Accurate5.2 精度与效率的平衡经过上百次计算测试我总结出几个经验法则筛选材料阶段PRECFOCK Fast NKRED 2 EDIFF 1E-5计算速度快3倍精度损失约0.1eV最终精度计算PRECFOCK Normal NKRED (不设置) EDIFF 1E-6特别敏感体系如拓扑材料PRECFOCK Accurate ADDGRID .TRUE.5.3 杂化泛函的特殊应用除了常规能带计算杂化泛函还有一些高阶用法激子效应研究LOPTICS .TRUE. CSHIFT 0.1配合HSE06可以准确预测激子结合能缺陷能级计算ISMEAR 0 SIGMA 0.01需要特别小的SIGMA值来分辨浅能级非线性光学性质LEPSILON .TRUE. LRPA .FALSE.使用HSE06计算二阶极化率更准确