计算化学公社

标题: 石墨烯计算不收敛及电荷设置等问题 [打印本页]

作者
Author: yoghurt117 时间: yesterday 19:04
标题: 石墨烯计算不收敛及电荷设置等问题
本帖最后由 yoghurt117 于 2026-1-22 19:14 编辑

(, 下载次数 Times of downloads: 5) (, 下载次数 Times of downloads: 0)

在做单空位缺陷石墨烯相关计算的时候遇到点问题，总结了一下
1. 做了一下单空位石墨烯，单空位石墨烯和Na+离子的结构，开始用cp2k晶胞优化得到的pbe泛函下晶胞边长大概2.6212埃(与实验值差别较大，大了8%左右)，然后做了一个32原子的表面结构，接着删除了一个碳原子，然后用multiwfn对附件中的xyz文件制作了cp2k的输入文件。单独缺陷石墨烯的文件做出来，没有设置电荷，波函数优化不收敛；

2. 缺陷石墨烯和钠离子的结构设置了电荷+1，用multiwfn制作的输入文件，里面有了uks，电荷和自旋多重度，生成的multiplicity值为2，然后cp2k运行报错，我注释掉了自旋多重度，能计算，但是还是不收敛。

我也不太清楚哪里出了问题，网上查了一些资料，也试着改了参数，不知道为什么不收敛，做的这个正交的石墨烯结构周期性扩展以后是没有问题的，我现在怀疑是结构太小了，不知道还有没有其他原因，谢谢大家。

作者
Author: sobereva 时间: 12 hour ago
这么小的晶胞连k点都没考虑，完全都是胡算瞎算。十分不推荐用PSOLVER ANALYTIC。

PS：十分建议通过北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）完整系统学学再计算，免得像这样因为缺乏基础和常识踩坑。Multiwfn完全无法把CP2K变成黑箱程序，只不过是能作为创建输入文件的很方便的工具。

作者
Author: yoghurt117 时间: 4 hour ago
好的，谢谢社长

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