the old and the new charge density differ

问题

最近遇到了一个棘手的问题，有个体系能量不收敛。即集群上提交多个任务，有的可以正常计算，有的不可以。

如果中断了错误的任务，会出现如下报错，

1
2
3
4

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中添加，

1

export I_MPI_ADJUST_REDUCE=3

所有任务完美结束。