计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First Principle) › CP2K进行晶体表面吸附小分子的结构优化过程不收敛的原因 ...

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

查看数: 1406 | 评论数: 8 | 收藏 Add to favorites 2

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

Yaqi

发布时间: 2023-12-9 19:51

正文摘要:

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

回复 Reply

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

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

效果怎么样

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

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

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

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

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

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

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

谢谢老师

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

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

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

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进行优化才可以使用？或者是针对特殊情况才能使用？

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

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

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

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）

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