做VASP的NEB时,很多人前面不是不会插点,而是起算顺序没有收好,结果一跑就出现图像漂移、收敛发散,或者路径一开始就偏得很厉害。VASP Wiki对NEB的起算步骤写得很清楚,先优化初态和末态,再准备中间图像,再在父目录统一放INCAR、KPOINTS和POTCAR,最后从父目录启动整套计算;而且官方还特别提醒,中间图像不要贪多,先用尽量少的图像把路径跑顺,再决定要不要加密,会更稳。
一、VASP NEB计算怎么设置
NEB不是把几份POSCAR摆好就结束,真正影响后面能不能顺利收敛的,是目录结构、图像数量、弹簧参数和离子优化方法这几层有没有先站稳。VASP Wiki的标准做法是,在父目录下建立从00开始编号的子目录,00放初态,最后一个目录放末态,中间目录放插值后的中间图像;其中IMAGES填的是中间图像数量,不包含两端固定结构。
1、先把初态和末态单独优化干净
VASP Wiki把这一步放在NEB设置的第一步,而且说明得很直接,两端固定结构后面会作为整条弹性带的锚点,初态和末态本身如果没优化好,后面得到的路径就会跟着偏。
2、目录按两位编号排开
如果中间图像数是3,就建00、01、02、03、04这五个目录。00放初态,04放末态,01到03放中间图像。VASP Wiki还特别提醒,目录名要保持两位编号,这样多图像计算时最不容易出错。
3、IMAGES和SPRING先按标准起步
VASP Wiki说明,IMAGES设为中间图像数,SPRING负责控制图像之间的弹簧常数;NEB模式下常用负值,默认就是-5,而且官方明确提醒不要把SPRING设得过大,否则反而会拖慢收敛。起步阶段直接用默认的负五,通常比自己乱改更稳。
4、离子优化算法先避开共轭梯度
VASP Wiki专门提到,NEB的非线性约束会让IBRION等于2这一类共轭梯度优化更容易出收敛问题,因此更建议优先试IBRION等于1或IBRION等于3。要是前几步图像还比较乱,官方还建议先用更保守的预收敛方式把路径拉顺,再继续往下跑。
二、VASP NEB爬山算法怎么开启
爬山算法这件事,别和普通NEB混成一步。按公开文档来看,VASP Wiki的NEB页面本身把更高级的过渡态搜索直接指向了VTST工具,而VTST的官方说明写得很明确,NEB模式下如果要打开climbing image,就把LCLIMB设成TRUE;同时VTST的安装和使用说明也写到,NEB对应ICHAIN等于0。也就是说,普通NEB和爬山图像最好分成两段来做,不要一开始路径还没站稳就直接开爬山。
1、先把普通NEB跑到路径基本成形
VASP Wiki一直在强调中间图像要尽量接近真实过渡路径,这个提醒放到爬山算法这里也一样成立。更稳的做法通常是先把普通NEB跑到能看出明显最高能量图像,再考虑让最高点继续往鞍点爬。
2、开启爬山时把LCLIMB设成TRUE
VTST的NEB页面给出的正式写法就是LCLIMB等于TRUE,用它来把最高能量图像变成climbing image。这个设置的核心作用,是让最高图像不再只沿弹性带平衡,而是继续朝真正过渡态靠近。
3、NEB模式先确认ICHAIN等于0
VTST安装文档说明,ICHAIN用来指定链式方法的类型,其中ICHAIN等于0对应NEB。也就是说,如果你现在配的是VTST这套过渡态工具,NEB模式本身也要先放对,后面LCLIMB才是接在这条链上的下一步。
4、爬山阶段优化器不要乱换
如果你用的是VTST的力优化器路线,公开文档给出的常见做法是配合IOPT这类力基优化器来跑,像Quick Min就是常见起步选项之一;而VASP Wiki则提醒NEB本身不太适合用IBRION等于2这一类共轭梯度。把这两条放在一起看,更稳的做法就是在爬山阶段继续用适合NEB的优化器,不要为了图快又切回最容易发散的那类。
三、VASP NEB起算前哪些地方最容易出错
真正容易返工的地方,往往不在某一个参数,而在起算顺序做反了。明明初态和末态都还没优化干净,就先插很多图像;明明路径还没站稳,就急着开爬山;明明目录和输入文件放得不统一,后面又一直怀疑算法本身。这些问题单看都不大,凑在一起就最容易把整条路径带偏。VASP Wiki的标准流程其实已经把排错顺序给出来了,只要按它的节奏做,很多弯路本来就可以避开。
1、先查两端结构是不是已经独立收敛
只要初态和末态本身还不干净,后面的中间图像再怎么调,也只是建立在不稳的端点上。NEB最先该确认的永远是两端。
2、再查图像数量是不是一开始就放太多
VASP Wiki明确建议图像数从少起步,先用一张或少量中间图像把路径拉顺,再考虑加密。图像一开始就堆太多,收敛反而更慢,也更难看清真正的问题到底出在哪一段。
3、然后查普通NEB和爬山是不是混着开了
如果路径还没成形就直接开LCLIMB,最高图像经常会抢着往前跑,表面上看像是在找鞍点,实际却可能把原本还不稳的路径再拉歪一次。更稳的做法还是先普通NEB,后climbing image。这个判断是把VASP Wiki的标准起算步骤和VTST对LCLIMB的开启方式连起来以后,最容易落地的一条经验。
4、最后再查优化器和收敛方式
如果前面三层都没问题,结果还是总发散,再回头看IBRION、必要时的预收敛做法,以及你当前是不是用了更适合NEB的优化路线。这样排,通常比一上来就到处改INCAR更容易把问题收住。
总结
VASP的NEB真正难的,不是把几个目录摆出来,而是先把初态、末态、图像数量和优化方法这几层先收顺。路径还没站稳时,先把普通NEB跑清楚,确认最高图像落点基本合理,再去开climbing image,通常会比一开始就追求过渡态更稳。至于爬山算法本身,按公开文档最直接的开法就是在VTST路线下把LCLIMB设成TRUE,同时保持NEB模式和优化器设置一致。把这层顺序想明白,后面的路径和势垒结果通常会干净很多。