计算化学公社

标题: 求助，PyridineN_Fe_S吸附OOH，计算能量不收敛 [打印本页]

作者
Author: zhaoquanyou 时间: 2025-3-9 18:39
标题: 求助，PyridineN_Fe_S吸附OOH，计算能量不收敛
求助大家，为什么我用cp2k计算Pyridine_N_S的结构对OOH分子吸附，发现几何优化前边不收敛，后边慢慢收敛的，但是我用收敛后的restart文件创建能量计算，但是发现收敛性特别差，我搞不明白为什么，按照道理来说既然几何优化的后边每一轮能量都是可以收敛的，就代表这个结构是合理的，但是为什么单独计算能量时，SCF就不收敛，而且完全没有收敛的趋势呢？

作者
Author: sobereva 时间: 2025-3-9 22:38
因为几何优化默认会用之前几步的波函数外推出的波函数当初猜，下文里搜“外推”，自然比起以自动产生的初猜更容易SCF收敛。如果之前优化任务留下了波函数文件，读取之当初猜
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

作者
Author: zhaoquanyou 时间: 2025-3-9 22:54
好的好的，我明白了，用之前的几何优化的波函数文件作为能量计算的初猜，就可以解决这个问题，谢谢老师

作者
Author: zhaoquanyou 时间: 2025-3-10 00:46
老师，我不明白的是为什么几何优化最终能够收敛，但是能量计算用几何优化的波函数进行初猜，还是不收敛

作者
Author: sobereva 时间: 2025-3-10 02:06

zhaoquanyou 发表于 2025-3-10 00:46
老师，我不明白的是为什么几何优化最终能够收敛，但是能量计算用几何优化的波函数进行初猜，还是不收敛

且不说别的，PyridineN_Fe_S_OOH_energy5 (1).inp里的晶胞信息都和PyridineN_Fe_S_OOH_opt2-2.inp里的不一样。能弄一致的都弄得严格一致再说

要回复我，点击帖子下面的回复按钮，否则>99%概率我看不到你后续发的东西

作者
Author: zhaoquanyou 时间: 2025-3-10 08:35

sobereva 发表于 2025-3-10 02:06
且不说别的，PyridineN_Fe_S_OOH_energy5 (1).inp里的晶胞信息都和PyridineN_Fe_S_OOH_opt2-2.inp里的不 ...

谢谢老师，算的太多了，脑子昏掉了，现在已经成功收敛

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