the old and the new charge density differ

最近遇到了一个棘手的问题,有个体系能量不收敛;集群上提交多个任务,有的可以正常计算,有的不可以。如果中断了错误的任务,会出现如下报错,

DAV:  11    -0.185185649493E+05   -0.43364E+04   -0.24178E+03  1548   0.208E+02BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:   903.27265 new  897.41748
    0.502E+02

网上关于the old and the new charge density differ的大多数方法都用了,除了ALGO = FAST有效果,其余没有解决问题;但既然有的任务可以正常结束,且为了保持课题参数的一致性,就只能放弃这个方法了。

后来我尝试了1 重新编译VASP;2 重新编译Intel MKL MPI;3 修改/etc/profile等;未果。

NSC网页上发现了解决方法,只需于job.sh中添加,

export I_MPI_ADJUST_REDUCE=3

所有任务完美结束。