求助 VASP过渡态计算鞍点不是插入的整数点怎么查看对应的结构

lyzhang

VASP计算过渡态时插入3个点，计算结果收敛，但是极值点没有出现在插入的1,2等整数点，而是出现在2.5，那么怎么查看2.5这个点对应的结构呢？虚频验证的时候就用2.5这个结构吗

含光君

寻找过渡态时就应该用CI-NEB方法去做。
已经用NEB算了也没关系，neb2dim.pl脚本会自动计算极值点附近的两个结构取中点自动生成DIMER计算的POSCAR。
参见刘锦程老师的博客 https://blog.shishiruqi.com/2019/05/01/scf/

lyzhang

含光君 发表于 2023-2-17 00:39
寻找过渡态时就应该用CI-NEB方法去做。
已经用NEB算了也没关系，neb2dim.pl脚本会自动计算极值点附近的两 ...

您好，谢谢您的指点。我计算的时候设置了CINEB的方法，LCLIME=.TRUE.，但是给出的极值点不是插入点，这是计算出错了吗？下面是输入的INCAR
SYSTEM=W
ISTART=0
ICHARGE=2
ENCUT=500
ISMEAR=1
SIGMA=0.15
IBRION=3
POTIM=0
IOPT=1
ICHAIN=0
LCLIMB=.TRUE.
IMAGES=3
SPRING=-5
ISIF=2
NSW=100
EDIFF=1E-3
EDIFFG=-0.01
ISYM=0
PREC=Normal
LREAL=Auto
LWAVE=.F.
LCHARGE=.F.
NCORE=2

ahxb

lyzhang 发表于 2023-2-17 10:47
您好，谢谢您的指点。我计算的时候设置了CINEB的方法，LCLIME=.TRUE.，但是给出的极值点不是插入点，这是 ...

mep.eps里绘制的曲线是通过最后各image中结构的能量插值得到的，有时不代表真实的能量变化情况。
如果只是为了确定过渡态结构，在路径合理时，可以去对应文件夹里查OUTCAR，看是否收敛，而不必过于依赖nebresults.pl给出的信息。

lyzhang

ahxb 发表于 2023-2-17 14:51
mep.eps里绘制的曲线是通过最后各image中结构的能量插值得到的，有时不代表真实的能量变化情况。
如果只 ...

我是计算原子迁移能垒，插入了三个点，计算收敛，曲线显示鞍点出现在大概2.5的位置，此时算出来的能垒是可靠的吗，还是需要调整插点位置，使计算出来的插入点处于曲线最高点才行呢

ahxb

lyzhang 发表于 2023-2-17 20:21
我是计算原子迁移能垒，插入了三个点，计算收敛，曲线显示鞍点出现在大概2.5的位置，此时算出来的能垒是 ...

”鞍点“一词对一维曲线不是合适的用词，我猜你说的是极大值。
mep.eps中的曲线是通过插值算法算出的，除插值点外的其余位置的值不由vasp计算得到，一般不作为绝对有效的判据。如果各点的能量变化趋势合理，且能量最高的结构达到了收敛限，可以不理mep.eps中的曲线的问题。

lyzhang

ahxb 发表于 2023-2-18 13:45
”鞍点“一词对一维曲线不是合适的用词，我猜你说的是极大值。
mep.eps中的曲线是通过插值算法算出的， ...

好的，清楚了，非常感谢您