计算化学公社

标题: 求助使用cp2k结合multiwfn得到密度差 [打印本页]

作者
Author: mmxll 时间: 2025-3-5 22:59
标题: 求助使用cp2k结合multiwfn得到密度差
老师，您好。在按照使用CP2K结合Multiwfn绘制密度差图、平面平均密度差曲线和电荷位移曲线 - 第一性原理 (First Principle) - 计算化学公社帖子的内容计算结构的电子密度的cub文件，但是其中结构中铝13的部分scf不收敛，对角化及OT算法都已经计算过了，也已经按照CP2K中遇到SCF难收敛时的解决方法 - 第一性原理 (First Principle) - 计算化学公社帖子修改尝试过了，整体结构及两个片段生成cub文件的inp文件在下面，麻烦老师看一下是什么问题，谢谢老师。

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-3-5 23:14
能看出来这个模型是单层石墨烯吸附铝团簇（疑似Superatom，但是电荷和自旋多重度存疑），但是且不说有没有做过几何优化，让擅长处理周期性体系的软件用孤立簇模型加极大真空区做计算是个什么操作？多少有点委屈CP2K了……

作者
Author: sobereva 时间: 2025-3-6 02:39
这个体系用CP2K计算极其不合适。建议用Gaussian等主流量子化学程序，然后按下文得到密度差图
使用Multiwfn作电子密度差图
http://sobereva.com/113

如果你的目的就是计算石墨烯负载Au13，应当让CP2K以二维周期性模型来算（用石墨烯而不是氢饱和的石墨烯片段）

当前设置也莫名其妙。用ALPHA 0.001完全是开玩笑。PSOLVER MT也用的不是时候。开smearing或许可能帮助收敛

作者
Author: mmxll 时间: 2025-3-6 08:57
谢谢老师

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