计算化学公社

标题: cp2k精度设置以及加弥散后“Cholesky decompose failed”错误求助 [打印本页]

作者
Author: hhhhhhhhh 时间: 2024-1-23 11:29
标题: cp2k精度设置以及加弥散后“Cholesky decompose failed”错误求助
本帖最后由 hhhhhhhhh 于 2024-1-23 15:13 编辑

各位老师好！我计算的目标是针对周期性系统和非周期性系统，分别使用cp2k和Gaussian进行dft计算得到能量和受力，并使得这两者的计算结果在相同的理论水平（或者说具有可比性）。Gaussian中计算设置的关键词是“ub3lyp/6-311+g(d,p) em=GD3BJ guess=read nosymm force scf=noincore”，由于使用Gaussian计算周期性结构无法加色散，以及已经使用Gaussian计算并得到了大量结果不方便更改计算软件，因此想要利用cp2k进行周期性体系的计算，对此我有以下几个问题
1. 如何在cp2k里面对dft的积分格点进行设置，使得cp2k计算的能量和受力的结果具有和Gaussian计算结果的可比性？
2. 在阅读了社长的博文后：CP2K做杂化泛函计算的关键要点和简单例子 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)，计划先做“blyp-d3(bj)/6-311+G**”的计算得到wfn文件，再做“b3lyp-d3(bj)/6-311+G**”的计算得到最后的结果，但是我在cp2k中设置的泛函和基组为“blyp/6-311+G**”，就会出现 “Cholesky decompose failed: the matrix is not positive definite or ill-conditioned”错误，而相同的结构我在不加弥散时，不管是blyp或是b3lyp都能计算且scf最后收敛，我尝试加大结构和更换基组为6-311++G**（6-311+G**的参数是我自己加进去的），仍然出现相同的错误，请问各位老师，这个错误应该如何解决呢？附件给出了我使用的“6-311+G**”基组参数，我把这些参数添加到文件EMSL_BASIS_SETS中去了。
3. 我计算的结构是石墨烯，一般说来，在计算的时候应该开启smearing，但是我发现在开启smearing时，耗时很长，且有很有可能scf不收敛，即使此时已经读取纯泛函的wfn文件，而开启ot计算却能节省大量时间，因此我的问题是由于我只需要结构的能量和受力，能否直接使用ot的计算结果？附件给出了分别使用smearing和ot计算的out文件，由于弥散的问题我还没有解决，因此计算的泛函和基组为“b3lyp(d3bj)/6-311G**”，开启smearing计算时读取的wfn文件是ot计算得到的（直接读纯泛函的wfn文件时scf不收敛）。

作者
Author: sobereva 时间: 2024-1-23 19:35
1 CUTOFF越大积分格点越精细，还可以&DFT / &XC / &XC_GRID里写USE_FINER_GRID T使用更精细的格点（四倍的CUTOFF）算交换-相关泛函部分。这只是影响可比性的一部分，参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里的讲解：

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

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

2 计算凝聚态周期性体系绝对不要用常规的带弥散函数的基组，这是常识知识，否则线性相关问题会非常厉害。必须要弥散函数就用MOLOPT系列基组，带了指数很小的高斯函数（尤其是不带SR的版本），而由于它是广义收缩又不至于出现线性依赖问题导致数值不稳定性。这点我之前在论坛里说过多遍

3 smearing并不会影响每一轮SCF的耗时。算导体体系不开smearing通常明显更难收敛。石墨烯开smearing时SCF收敛不了一定是自己使用不当，没必要试图靠切换成OT来回避。

作者
Author: hhhhhhhhh 时间: 2024-1-23 20:23
本帖最后由 hhhhhhhhh 于 2024-1-23 20:47 编辑

sobereva 发表于 2024-1-23 19:35
1 CUTOFF越大积分格点越精细，还可以&DFT / &XC / &XC_GRID里写USE_FINER_GRID T使用更精细的格点（四倍的C ...

好的，谢谢社长的解答。

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