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

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

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

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

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

好的老师，我再试试

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

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

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

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

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

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

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

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

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

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

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

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