Dimer方法计算过渡态相关问题

小螃蟹

我根据VASP手册中的6.61.2  给出的例子，将它应用到自己的体系中。首先猜测过渡态的初始构型，然后分别进行一下两步运算：       
第二步计算未能正常结束，出现了：Error EDDDAV: Call to ZHEGV failed. Returncode =  54 2  32的错误提示。请各位大神帮忙指点。

hakuna

1.那么简单的分子就不要+U了
2.换一下搜索TS的优化方法试试
3.参考例子参数设置

小螃蟹

hakuna 发表于 2019-4-27 11:02
1.那么简单的分子就不要+U了
2.换一下搜索TS的优化方法试试
3.参考例子参数设置

嗯，首先谢谢你的建议
1，我+U是因为自己的体系计算的是CeO2表面催化N2O，共有93个原子，对Ce进行了+U，
2. 我也查找了其它的过渡态方法，其中有CI-NEB，我也进行了计算，使用指令nebmake.pl系统生成了5个中间体，但是到最后中间过渡态构型都没收敛成功，下边有关NEB的输入：
然后计算完500圈后，几乎每个构型都是收敛到 500 F= -.75182661E+03 E0= -.75181080E+03  d E =0.325320E+00或500 F= -.75182661E+03 E0= -.75181080E+03  d E =0.938024E-01
然后OUTCAR里面没有收敛成功的显示
不知道您有遇到过此种情况吗？还请你指点。

hakuna

小螃蟹 发表于 2019-4-27 15:50
嗯，首先谢谢你的建议
1，我+U是因为自己的体系计算的是CeO2表面催化N2O，共有93个原子，对Ce进行了+U， ...

没看清楚,还以为你跑手册里面的例子.
搜TS需要有好的初猜,经验性比较强,特别是Dimer,初猜不对,是找不到TS的

CNEB的话,你这些信息不足以判断那里出了问题
反应物和产物设计不合理,找不到TS
插入的结构不合理也找不到TS
......
在保证IS和FS合理情况下,插入结构的合理性至关重要......
检查一下你的各个Image跑的情况, 看看达到指定迭代次数后, 结构是不是合理, 合理的话可以继续搜
如果不合理, 要调整初猜image,使之更合理

你可以降低收敛精度, 减少k点, 看看能不能粗收敛,然后在这个基础上, 再跑Dimer

小螃蟹

hakuna 发表于 2019-4-27 22:27
没看清楚,还以为你跑手册里面的例子.
搜TS需要有好的初猜,经验性比较强,特别是Dimer,初猜不对,是找不到T ...

好的，再次谢谢你的建议。