Castep 截断能与k-point计算报错怎么解决？显示negativecharge density encountered

柯--

各位大佬好，我接触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的时候，计算报错了，报错结果如图 于是我尝试720，这个时候计算结果正确，且720与700的差值小于0.001ev。由于已经收敛，我打算进一步计算k-point。由于在720时候达到收敛，我根据这个Use custom energy cutoff值，修改k-point的Grid parameters。相关设置如下.其他未改变，当计算到666的时候，出现了以下报错 Warning: ion_atom_descreen_psp: negative charge density encountered。这两个报错描述是一致的，但是我搜索了很久暂时没发现能够解决这个问题的方法，于是向各位请教。另外，关于Mn3O4的基础设置我是否正确？对于自旋和极化那一块我不是很能理解，所以设置也仅是照猫画虎的设置，是否在这里就已经出现错误了呢？