求助：CP2K计算原子6000+大体系的电荷，SCF不收敛

John_Tao

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

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

sobereva

提供完整的输入文件。当前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跟你没关系

John_Tao

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

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

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

sobereva

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

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

干嘛跟DDEC6比，又没多理想

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

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

John_Tao

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

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

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