计算化学公社

标题: Dimer方法计算过渡态相关问题 [打印本页]

作者
Author: 小螃蟹 时间: 2019-4-26 21:13
标题: Dimer方法计算过渡态相关问题
我根据VASP手册中的6.61.2 给出的例子，将它应用到自己的体系中。首先猜测过渡态的初始构型，然后分别进行一下两步运算： (, 下载次数 Times of downloads: 52) (, 下载次数 Times of downloads: 55)
第二步计算未能正常结束，出现了：Error EDDDAV: Call to ZHEGV failed. Returncode = 54 2 32的错误提示。请各位大神帮忙指点。

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

作者
Author: 小螃蟹 时间: 2019-4-27 15:50

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的输入： (, 下载次数 Times of downloads: 59)
然后计算完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里面没有收敛成功的显示
不知道您有遇到过此种情况吗？还请你指点。

作者
Author: hakuna 时间: 2019-4-27 22:27

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

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

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

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

作者
Author: 小螃蟹 时间: 2019-4-28 08:30

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)