计算化学公社

标题: 求助：CP2K计算原子6000+大体系的电荷，SCF不收敛 [打印本页]

作者
Author: John_Tao 时间: 2023-1-24 14:39
标题: 求助：CP2K计算原子6000+大体系的电荷，SCF不收敛
本帖最后由 John_Tao 于 2023-4-12 16:31 编辑

各位老师好，
我想用cp2k跑一个COF的单点。体系很大，6000多原子，元素是C、H、O、N。
我用Multiwfn在默认设置下产生的输入文件（PBE/DZVP-MOLOPT-SR-GTH）试了一下，但完全收敛不了。请问这样是可行的吗？
我想算一下体系的Hirshfeld-I或者REPEAT电荷，不知道用什么方法和参数的设置可以实现，请各位老师指点！

另外，我在用cp2k计算电荷时，输出文件中有DDAP FULL DENSITY charges，想请问各位老师这个电荷的含义是什么？

作者
Author: sobereva 时间: 2023-1-24 15:59
提供完整的输入文件。当前inp里根本不完整，连你算的结构是什么样都不知道，根本没法回答
&COORD
   N       30.84532117 57.11336679 14.25069218
   O       26.76704066 52.58831040 7.13555737
   C       37.34053486 52.03104744 11.25538338
   ...
   ...  #the number of atoms are more than 6000
   ...
   H       38.00375564 51.61512909 10.75337133
&END COORD

DDAP FULL DENSITY charges跟你没关系

作者
Author: John_Tao 时间: 2023-1-24 16:33
本帖最后由 John_Tao 于 2023-1-24 17:00 编辑

sobereva 发表于 2023-1-24 15:59
提供完整的输入文件。当前inp里根本不完整，连你算的结构是什么样都不知道，根本没法回答
&COORD
...

感谢老师回复，已更新好完整的。

因为我是周期性多孔材料，原先我是用VASP+DDEC6算电荷。我尝试cp2k以后，发现Hirshfeld-I和REPEAT算出来的都和用DDEC6差的很多，但是无意中却发现DDAP FULL DENSITY charge和DDEC6的结果很接近。

作者
Author: sobereva 时间: 2023-1-25 09:05

John_Tao 发表于 2023-1-24 16:33
感谢老师回复，已更新好完整的。

因为我是周期性多孔材料，原先我是用VASP+DDEC6算电荷。我尝试cp2k以 ...

干嘛跟DDEC6比，又没多理想

老老实实用REPEAT，对此类体系原理上最理想

这么大体系应该开OT，不仅更快，往往也容易收敛得多

作者
Author: John_Tao 时间: 2023-1-25 09:56

sobereva 发表于 2023-1-25 09:05
干嘛跟DDEC6比，又没多理想

老老实实用REPEAT，对此类体系原理上最理想

谢谢sob老师，我开了OT，收敛了！

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