计算化学公社

标题: CP2K使用CINEB搜索过渡态scf不收敛 [打印本页]

作者
Author: Wednesday 时间: 2026-3-24 15:10
标题: CP2K使用CINEB搜索过渡态scf不收敛
本帖最后由 Wednesday 于 2026-3-24 18:13 编辑

(, 下载次数 Times of downloads: 2) 这是一个副本文件的截图，每一个电子步还是很快的，但是最大位移和最大受力一直震荡到了300步上限，这个scf不收敛想请教下各位老师是什么原因呢？

作者
Author: northern_mar 时间: 2026-3-24 19:57
本帖最后由 northern_mar 于 2026-3-24 20:00 编辑

你这个看起来已经快收敛了。在检查能量和轨迹后，你可以简单查看一下有没有收敛趋势。例：grep -e 'MAX FORCE' neb-BAND5.out
如果有收敛趋势，拿优化出来的结构接着做就行。如果出现震荡的情况，可以试着减小K_SPRING到0.02。

作者
Author: Wednesday 时间: 2026-3-24 20:39

northern_mar 发表于 2026-3-24 19:57
你这个看起来已经快收敛了。在检查能量和轨迹后，你可以简单查看一下有没有收敛趋势。例：grep -e 'MAX FOR ...

老师，这个我看了一下具体150步就已经开始震荡了，按理说早该收敛的，就这两个最大值一直震荡，调整k_spring就可以么？

作者
Author: northern_mar 时间: 2026-3-24 20:53

Wednesday 发表于 2026-3-24 20:39
老师，这个我看了一下具体150步就已经开始震荡了，按理说早该收敛的，就这两个最大值一直震荡，调整k_spr ...

这是指CINEB可以调整的办法，或者你也可以切换dimer方法。
其他的解决办法可以参考讲义中“结构优化难收敛的解决”部分，有不少办法可以尝试。

作者
Author: Wednesday 时间: 2026-3-25 14:33
@Stardust0831 星尘老师您好，受力和位移已经出现震荡，还需要麻烦您再看一下

作者
Author: Stardust0831 时间: 2026-3-25 14:35

Wednesday 发表于 2026-3-25 14:33
@Stardust0831 星尘老师您好，受力和位移已经出现震荡，还需要麻烦您再看一下

提供尽可能多的原始文件，打个压缩包后上传。

作者
Author: Wednesday 时间: 2026-3-25 14:46
本帖最后由 Wednesday 于 2026-3-25 15:32 编辑

Stardust0831 发表于 2026-3-25 14:35
提供尽可能多的原始文件，打个压缩包后上传。

做的是C6烯烃在ZSM-5上的脱氢反应，C6烯烃上的一个H原子脱到Al附近的O上，3.25的输出文件是根据您在群里的建议修改的，一个电子步需要10s左右，所以目前跑的时间不太长，但是从结果来看已经开始震荡。3.24的是昨天尝试过程的输出文件和轨迹，一直震荡到了300步上限。每个副本基本一致，就只上传了一个文件

作者
Author: northern_mar 时间: 2026-3-25 18:56
本帖最后由 northern_mar 于 2026-3-25 18:57 编辑

Wednesday 发表于 2026-3-25 14:46
做的是C6烯烃在ZSM-5上的脱氢反应，C6烯烃上的一个H原子脱到Al附近的O上，3.25的输出文件是根据您在群里 ...

先检查能量和结构。我看了下你3.24的能量，副本3能量比初态低，你可以取副本3优化后作为初态再做。

作者
Author: Wednesday 时间: 2026-3-25 18:59

northern_mar 发表于 2026-3-25 18:56
先检查能量和结构。我看了下你3.24的能量，副本3能量比初态低，你可以取副本3优化后作为初态再做。

好的老师，我再试试

作者
Author: Stardust0831 时间: 2026-3-25 19:12

Wednesday 发表于 2026-3-25 18:56
@sobereva sob老师，麻烦您帮我再看一下，没办法的话我只能是换软件了

cp2k几乎就是你的最优解了。
妥善使用GFN1-xTB这种半经验的级别，可以先拿个大概的结构/轨迹，然后再dft优化。

作者
Author: nuonuo123456 时间: 2026-3-25 20:23

Wednesday 发表于 2026-3-25 18:56
@sobereva sob老师，麻烦您帮我再看一下，没办法的话我只能是换软件了

你这个体系也不是很大，可以用vasp试一下，我感觉vasp在过渡态方面还是更稳定更优越一些。

作者
Author: Wednesday 时间: 2026-4-2 11:06
本帖最后由 Wednesday 于 2026-4-2 11:08 编辑

Stardust0831 发表于 2026-3-25 19:12
cp2k几乎就是你的最优解了。
妥善使用GFN1-xTB这种半经验的级别，可以先拿个大概的结构/轨迹，然后再dft ...

星尘老师，GFN1-xTB提任务报错我发现是因为cp2k自带的image有bug，改用sobneb和idpp后可以正常提交，但是CI-NEB还是同样的问题，跑到283步时max force依然震荡，而且收敛性似乎比用PBE的要差，看不到减小的趋势，想请教下您还有其他可用的方法么？

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