计算化学公社

标题: 利用CP2K计算石墨时的问题求助 [打印本页]

作者
Author: Yjc 时间: 2023-7-18 12:44
标题: 利用CP2K计算石墨时的问题求助
本帖最后由 Yjc 于 2023-7-18 13:33 编辑

各位老师好，我在利用CP2K计算石墨能量时输出文件出现报错，而且是在还未开始第一轮SCF迭代的时候就出现的：

Cholesky decompose failed: matrix is not positive definite or ill-conditioned

这次任务的k点数目是经过收敛性测试后得到的（13,13,8）而且考虑到石墨是导体，也开启了smearing；烦请各位老师帮忙看下计算过程中存在哪些问题，如何解决这一报错，输入文件和输出文件放在附件里了，谢谢各位老师！

(, 下载次数 Times of downloads: 8)

(, 下载次数 Times of downloads: 1)

PS:这次的报错貌似是基组问题，当换用DZVP-MOLOPT-SR-GTH或TZVP-MOLOPT-SR-GTH之后很快就能算完；之所以用TZV2P-GTH是因为想尝试下精度和耗时介于DZVP-MOLOPT-SR-GTH和TZVP-MOLOPT-SR-GTH之间的泛函，没想到发生了报错；这可能是因为TZV2P恰好不适合算这个体系嘛？

作者
Author: sobereva 时间: 2023-7-20 13:17
若无必要甭用TZV2P-GTH这种东西，老的GTH赝势基组的比较大尺寸的版本很容易因为线性相关问题导致崩溃。MOLOPT系列线性相关问题比老GTH赝势基组好得多

作者
Author: Yjc 时间: 2023-7-20 16:20

sobereva 发表于 2023-7-20 13:17
若无必要甭用TZV2P-GTH这种东西，老的GTH赝势基组的比较大尺寸的版本很容易因为线性相关问题导致崩溃。MOLO ...

谢谢sob老师解答
那若是TZVP-MOLOPT-SR-GTH算不动而又感觉DZVP-MOLOPT-SR-GTH档次稍低，那应该如何合理选择基组呢？

作者
Author: sobereva 时间: 2023-7-20 19:25

Yjc 发表于 2023-7-20 16:20
谢谢sob老师解答
那若是TZVP-MOLOPT-SR-GTH算不动而又感觉DZVP-MOLOPT-SR-GTH档次稍低，那应该如何合理 ...

不至于DZVP-MOLOPT-SR-GTH算得轻轻松松而TZVP-MOLOPT-SR-GTH完全算不动。何况你的体系非常小，像样机子下TZVP-MOLOPT-SR-GTH算不动是不可能的

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