计算化学公社

标题: Castep 截断能与k-point计算报错怎么解决?显示negativecharge density encountered [打印本页]

作者
Author:
柯--    时间: 2025-10-15 21:31
标题: Castep 截断能与k-point计算报错怎么解决?显示negativecharge density encountered
本帖最后由 柯-- 于 2025-10-15 21:34 编辑

各位大佬好,我接触MS才两周时间,学的不多不精很多地方不太理解,请多包涵。下面是我遇到的问题:
首先,我建立了一个Mn3O4的晶胞,通过网站直接给出的cif以及Findit对比,确认建立的模型是正确的。
接着,我通过cestep--calculation做了如下设置:
虽说列出了这么多,但是实际修改的主要是basis中的Use custom energy cutoff以及SCF cycles到500,还有k-points的Grid parameters修改为9 9 9.此时,确保了Properties中未勾选任何选项,打开Job control开始运行。在我检验截断能的过程中,energy cutoff从200,250,...,700,750.结果表明,在100~700的时候都没出现任何问题。但是当我想去计算750的时候,计算报错了,报错结果如图 (, 下载次数 Times of downloads: 10) 于是我尝试720,这个时候计算结果正确,且720与700的差值小于0.001ev。由于已经收敛,我打算进一步计算k-point。由于在720时候达到收敛,我根据这个Use custom energy cutoff值,修改k-point的Grid parameters。相关设置如下.其他未改变,当计算到666的时候,出现了以下报错 (, 下载次数 Times of downloads: 12) Warning: ion_atom_descreen_psp: negative charge density encountered。这两个报错描述是一致的,但是我搜索了很久暂时没发现能够解决这个问题的方法,于是向各位请教。另外,关于Mn3O4的基础设置我是否正确?对于自旋和极化那一块我不是很能理解,所以设置也仅是照猫画虎的设置,是否在这里就已经出现错误了呢?




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3