计算化学公社

标题: 求助：静态/非自洽计算时出现VERY BAD NEWS [打印本页]

作者
Author: shinkon 时间: 2023-12-3 03:30
标题: 求助：静态/非自洽计算时出现VERY BAD NEWS
我有一个超胞，其晶格常数如图所示。我之前是用的vaspkit的102功能-gamma方法-0.04来生成k点（如图所示，2x2x3）
(, 下载次数 Times of downloads: 13)
今天在检查slurm输出文件的时候（之前都是直接tail看能量和磁矩有没有太离谱），发现在做静态及非自洽计算的时候，有如下提示（结构优化没有）：
VERY BAD NEWS! internal error in subroutine IBZKPT:
Reciprocal lattice and k-lattice belong to different class of lattices. Often results are still useful... 168

但能量能正常收敛，磁矩什么的也没问题，画出来的DOS和能带图也是正常的。请问这种情况下，这个warning可以被忽略吗？计算结果可以采信吗？

PS：后续我用了更符合这个倒格子的KPOINTS（3x2x3）（其他参数不变）进行一次静态计算，得到的能量反而略高于原先用不太恰当的k点密度的静态计算获得的能量

作者
Author: 卡开发发 时间: 2023-12-3 09:02
这个warning可以被忽略，但有可能这个k点间距在三个轴向上有点差别，实际确实3:2:3会均匀一些。如果3x2x3能量会升高，可以关注一下这个能量升高是多少，如果对能量和力造成的影响不可被忽略，那么应该对k好好进行一下收敛性测试。

作者
Author: shinkon 时间: 2023-12-3 14:00

卡开发发发表于 2023-12-3 09:02
这个warning可以被忽略，但有可能这个k点间距在三个轴向上有点差别，实际确实3:2:3会均匀一些。如果3x2x3能 ...

好的，谢谢卡老师。

作者
Author: shinkon 时间: 2023-12-3 15:56

卡开发发发表于 2023-12-3 09:02
这个warning可以被忽略，但有可能这个k点间距在三个轴向上有点差别，实际确实3:2:3会均匀一些。如果3x2x3能 ...

我的问题，比错数据了。323的K点的静态自洽比223的K点的静态自洽的能量低了2x10^-6eV/atom，看来没必要专门再全部重算一遍了。。

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