计算化学公社

标题: cp2k计算gap小的配合物能带后出现空轨道之间能量差比较大的问题 [打印本页]

作者
Author: tianmafei 时间: 2024-12-25 12:10
标题: cp2k计算gap小的配合物能带后出现空轨道之间能量差比较大的问题
本帖最后由 tianmafei 于 2024-12-25 18:18 编辑

各位老师好！

我使用cp2k 2024.1计算了一个二维配合物，开启了smearing才能让SCF迭代收敛，之后做了能带任务，结果显示带隙为0.059 ev。在能带图中LUCO似乎是一条直线，它与其他未占据轨道的能量差比较大，显得孤零零的。 (, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 0)

附件中是能带任务的输入输出文件和bs文件。我也做了单点任务，输出了轨道能级。老师们帮忙指导一二，这正常吗？是否计算的哪些因素出问题了？谢谢！

(, 下载次数 Times of downloads: 8)
(, 下载次数 Times of downloads: 4)

(, 下载次数 Times of downloads: 3)
(, 下载次数 Times of downloads: 1)

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

作者
Author: sobereva 时间: 2024-12-25 18:10
能带计算和单点计算输入文件明显不一样，前者SCF过程考虑了k点而后者没考虑，莫名其妙。要比就统一了再比。2D_AgS_sp.out对应的输入文件也必须上传

作者
Author: tianmafei 时间: 2024-12-25 18:16

sobereva 发表于 2024-12-25 18:10
此体系没有开smearing的必要，不开smearing也有办法弄收敛

能带计算和单点计算输入文件明显不一样，前者 ...

sob老师好！我曾试过能级移动，改变新密度矩阵的掺入比例等方法，SCF迭代了400多轮仍是震荡不收敛。您有推荐优先尝试的方法吗？先谢谢啦！

作者
Author: tianmafei 时间: 2024-12-25 18:23

sobereva 发表于 2024-12-25 18:10
能带计算和单点计算输入文件明显不一样，前者SCF过程考虑了k点而后者没考虑，莫名其妙。要比就统一了再比。 ...

抱歉！sob老师，单点的输入文件已补充上。我是小白一枚，做完能带任务后，觉得有问题但不知哪里出错了，就又做了单点任务，想看一下各轨道的能级。

作者
Author: sobereva 时间: 2024-12-25 18:42
试试PBE0，如果gap也很小，可能前线轨道对应连成一串的Ag的部分，可以看看相应轨道图像

作者
Author: tianmafei 时间: 2024-12-27 11:44
本帖最后由 tianmafei 于 2024-12-27 11:48 编辑

抱歉，给大家添麻烦了！这个结构中有一个孤立的水分子，但氧原子上没加氢，它处于晶体盒子的边缘，在进行几何优化时，把它和相邻晶胞中的镜像原子在一起优化成了氧气。把这个氧原子去掉，SCF迭代很容易就收敛了，用不上smearing。估计，拿着优化后结构做的能带任务，能带图中的那条直线与氧气相关。

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