CP2K使用CINEB搜索过渡态scf不收敛

Wednesday · 发表于 Post on 2026-3-24 15:10:44

本帖最后由 Wednesday 于 2026-3-24 18:13 编辑

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

northern_mar · 发表于 Post on 2026-3-24 19:57:15

本帖最后由 northern_mar 于 2026-3-24 20:00 编辑

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

Wednesday · 发表于 Post on 2026-3-24 20:39:18

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

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

northern_mar · 发表于 Post on 2026-3-24 20:53:01

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

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

Wednesday · 发表于 Post on 2026-3-25 14:33:38

@Stardust0831 星尘老师您好，受力和位移已经出现震荡，还需要麻烦您再看一下

Stardust0831 · 发表于 Post on 2026-3-25 14:35:33

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

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

Wednesday · 发表于 Post on 2026-3-25 14:46:52

本帖最后由 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步上限。每个副本基本一致，就只上传了一个文件

northern_mar · 发表于 Post on 2026-3-25 18:56:12

本帖最后由 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优化后作为初态再做。

Wednesday · 发表于 Post on 2026-3-25 18:59:13

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

好的老师，我再试试

Stardust0831 · 发表于 Post on 2026-3-25 19:12:12

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

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

nuonuo123456 · 发表于 Post on 2026-3-25 20:23:42

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

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

Wednesday · 发表于 Post on 2026-4-2 11:06:51

本帖最后由 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的要差，看不到减小的趋势，想请教下您还有其他可用的方法么？

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] CP2K使用CINEB搜索过渡态scf不收敛

浏览过的版块