计算化学公社

标题: cp2k使用sccs不收敛 [打印本页]

作者
Author: 古与春分 时间: 2025-7-21 18:24
标题: cp2k使用sccs不收敛
老师，我用cp2k模拟表面活性剂在方解石表面的吸附。在真空条件下这个任务很好的收敛，但是引入sccs后（用真空收敛的波函数当初猜，PBEsol/DZVP-MOLOPT-SR-GTH），不管是OT还是对角化，都不收敛。按照博客http://sobereva.com/665里方法尝试操作，也收敛不了。那我是不是应该换成GAPW结合全电子基组来试试这个了，比如B97M-rV泛函结合pob-TZVP-rev2，如果这能收敛，后面涉及sccs的模拟都得用B97M-rV泛函结合pob-TZVP-rev2

作者
Author: ultramanm87 时间: 2025-7-21 23:38
1. cp2k的SCCS收敛有不同于那篇博文的方法可以尝试， (, 下载次数 Times of downloads: 5)
2. cp2k2024版本用全电子基组+sccs会报错，要用2023版本。
3. 为啥要换泛函？
4. 如果以上1种方法都搞不定，可以尝试其他程序。QE的sccs比较稳健，但是体系一大就非常慢。ABACUS用的是原始版本的VASPSOL的方法，没法处理带电体系的溶解自由能，但是使用数值基速度快，收敛性还行。

作者
Author: sobereva 时间: 2025-7-22 05:46
修改源代码也可以让较新版本结合GAPW用SCCS，参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）幻灯片
(, 下载次数 Times of downloads: 4)

作者
Author: 古与春分 时间: 2025-7-23 11:06

ultramanm87 发表于 2025-7-21 23:38
1. cp2k的SCCS收敛有不同于那篇博文的方法可以尝试，
2. cp2k2024版本用全电子基组+sccs会报错，要用2023 ...

除了全电子基组外，其他全部都试过了，用cp2k就是不收敛。

作者
Author: 古与春分 时间: 2025-7-23 11:07

sobereva 发表于 2025-7-22 05:46
修改源代码也可以让较新版本结合GAPW用SCCS，参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.c ...

老师，我用的2025.1版本，按照要求删了那一行，但还是出现SCCS is not yet implemented with GAPW

作者
Author: sobereva 时间: 2025-7-25 15:32

古与春分发表于 2025-7-23 11:07
老师，我用的2025.1版本，按照要求删了那一行，但还是出现SCCS is not yet implemented with GAPW

随机应变。根据Fortran语法常识，总能找到引发报错的地方

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