计算化学公社

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

作者
Author: Yaqi 时间: 2023-12-9 19:51
标题: CP2K进行晶体表面吸附小分子的结构优化过程不收敛的原因求助
计算模型的正视图和俯视图如下图所示，每个小分子带一个负电荷。
优化过程分别尝试了(1) GFN1-xTB，冻结全部晶体原子；（2）GFN1-xTB，冻结晶体底部3层原子；（3）PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH，冻结晶体底部3层原子。
输入输出文件以上传附件。

(, 下载次数 Times of downloads: 13) (, 下载次数 Times of downloads: 12)

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

附件：
(, 下载次数 Times of downloads: 29)

作者
Author: sobereva 时间: 2023-12-9 21:05
(, 下载次数 Times of downloads: 12)

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

作者
Author: Yaqi 时间: 2023-12-9 22:37

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

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

作者
Author: Eudaimonia 时间: 2023-12-10 08:01

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

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

作者
Author: sobereva 时间: 2023-12-10 19:28

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

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

作者
Author: Eudaimonia 时间: 2023-12-11 10:16

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

谢谢老师

作者
Author: Yaqi 时间: 2023-12-18 16:02

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，起初害怕优化轮数太多，但后来发现每一轮收敛需要的时间大大减小，整体速度可以接受。

作者
Author: faylovesnow 时间: 2024-8-26 09:28

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

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

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

作者
Author: HAC 时间: 2024-9-29 17:27

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

效果怎么样

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