计算化学公社

标题: 用CI-NEB搜索过渡态时难以收敛怎么解决？ [打印本页]

作者
Author: 科研好难 时间: 2023-12-14 09:53
标题: 用CI-NEB搜索过渡态时难以收敛怎么解决？
大家好，我在cp2k中用CI-NEB搜索过渡态，始末结构是提前优化好的（分子吸附在晶面上，共331个原子），REPLICA =7，K_SPRING =0.08，计算了将近100步还是没有收敛趋势，且中间插点的能量比始末结构的能量更低，请问我这些情况怎么解决？望懂的老师不吝赐教，非常感谢！输入和部分输出文件在附件：
(, 下载次数 Times of downloads: 12) (, 下载次数 Times of downloads: 12) (, 下载次数 Times of downloads: 9) (, 下载次数 Times of downloads: 12)

作者
Author: sobereva 时间: 2023-12-14 13:39
拿中间能量更低的结构做几何优化看得到什么
可能本来就不是基元过程
并且最好再多插几个点
CUTOFF 300对于GTH赝势基组算含Mg的体系来说太小了

作者
Author: 科研好难 时间: 2023-12-14 15:07

sobereva 发表于 2023-12-14 13:39
拿中间能量更低的结构做几何优化看得到什么
可能本来就不是基元过程
并且最好再多插几个点

好的，sob老师，我再试试

作者
Author: reid 时间: 2023-12-14 19:55
我做过离子迁移的Cl-NEB，发现虽然插点多算的慢，但是结果好、收敛快

作者
Author: 科研好难 时间: 2023-12-14 20:18

reid 发表于 2023-12-14 19:55
我做过离子迁移的Cl-NEB，发现虽然插点多算的慢，但是结果好、收敛快

谢谢解答，我把插点加多再试试，非常感谢

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