CP2K进行晶体表面吸附小分子的结构优化过程不收敛的原因求助

Yaqi · 发表于 Post on 2023-12-9 19:51:57

计算模型的正视图和俯视图如下图所示，每个小分子带一个负电荷。
优化过程分别尝试了(1) GFN1-xTB，冻结全部晶体原子；（2）GFN1-xTB，冻结晶体底部3层原子；（3）PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH，冻结晶体底部3层原子。
输入输出文件以上传附件。

疑问：
1. 目前遇到的是SCF不收敛的问题吗？
2. 为什么输出文件最后收敛条件没有达到四个YES，但是有satisfied .... run CONVERGED!的表述？
3. 以上情况应该如何进一步调整计算？

附件：

输入输出文件.rar (299.96 KB, 下载次数 Times of downloads: 29)

sobereva · 发表于 Post on 2023-12-9 21:05:27

我就看了PBE-D3(BJ)的计算，最后一步SCF明显没收敛。尝试用OT，或下文其它的做法帮助SCF收敛
解决CP2K中SCF不收敛的方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

Yaqi · 发表于 Post on 2023-12-9 22:37:34

sobereva 发表于 2023-12-9 21:05
我就看了PBE-D3(BJ)的计算，最后一步SCF明显没收敛。尝试用OT，或下文其它的做法帮助SCF收敛
解决CP2K ...

好的，我再尝试一下，谢谢老师解答

Eudaimonia · 发表于 Post on 2023-12-10 08:01:06

sobereva 发表于 2023-12-9 21:05
我就看了PBE-D3(BJ)的计算，最后一步SCF明显没收敛。尝试用OT，或下文其它的做法帮助SCF收敛
解决CP2K ...

借楼问一下sob老师，LBFGS这个特殊的判定收敛标准结果可以正常使用吗？还是需要保证该体系下所有计算结果都使用LBFGS进行优化才可以使用？或者是针对特殊情况才能使用？

sobereva · 发表于 Post on 2023-12-10 19:28:58

Eudaimonia 发表于 2023-12-10 08:01
借楼问一下sob老师，LBFGS这个特殊的判定收敛标准结果可以正常使用吗？还是需要保证该体系下所有计算结果 ...

如果受力偏离0较明显，则不能用

Eudaimonia · 发表于 Post on 2023-12-11 10:16:22

sobereva 发表于 2023-12-10 19:28
如果受力偏离0较明显，则不能用

谢谢老师

Yaqi · 发表于 Post on 2023-12-18 16:02:18

sobereva 发表于 2023-12-9 21:05
我就看了PBE-D3(BJ)的计算，最后一步SCF明显没收敛。尝试用OT，或下文其它的做法帮助SCF收敛
解决CP2K ...

跟进一下：

根据卢老师的建议，做了如下几点调整，调整后优化过程可以顺利结束。
（1）优化方法从LBFGS改为BFGS；
（2）CUTOFF和REL_CUTOFF从350，50提高到400，55；
（3）使用OT；
（4）TRUST_RADIUS从0.2降低到0.05，起初害怕优化轮数太多，但后来发现每一轮收敛需要的时间大大减小，整体速度可以接受。

faylovesnow · 发表于 Post on 2024-8-26 09:28:52

Yaqi 发表于 2023-12-18 16:02
跟进一下：

根据卢老师的建议，做了如下几点调整，调整后优化过程可以顺利结束。

感谢您的跟进，我现在在优化小分子吸附在Cu表面的结构，也是遇到了优化震荡问题，准备尝试减小置信半径的做法，看看效果。

HAC · 发表于 Post on 2024-9-29 17:27:54

faylovesnow 发表于 2024-8-26 09:28
感谢您的跟进，我现在在优化小分子吸附在Cu表面的结构，也是遇到了优化震荡问题，准备尝试减小置信半径的 ...

效果怎么样

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

[CP2K] CP2K进行晶体表面吸附小分子的结构优化过程不收敛的原因求助