求助，VASP进行含Yb的周期性材料的优化，出现报错：新旧电荷不符

松鼠

在VASP中，对块体材料Yb2(SeO3)2(SO4)(H2O)2进行结构优化，下面是我的INCAR设置：

在之前的报错中，由于原胞的晶胞参数a，b，c相近，出现了Sub-Space-Matrix is not hermitian in DAV的报错，所以怀疑是k点设置的不合理的原因。添加了KPOINTS文件，用的4*4*4的mesh：
Regular 4 x 4 x 4 mesh centered at Gamma
0
Gamma
4 4 4
依旧会出现以下报错：
scaLAPACK will be used
LDA part: xc-table for Pade appr. of Perdew
POSCAR found type information on POSCAR  Yb S  O  Se H
POSCAR found :  5 types and      21 ions
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...
WAVECAR not read
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.339200265743E+05   -0.33920E+05   -0.53062E+04 24464   0.588E+03
DAV:   2     0.104879156101E+03    0.34025E+05   -0.64070E+04 23968   0.120E+03
DAV:   3    -0.120856618963E+06   -0.12096E+06   -0.40168E+03 23192   0.310E+02
DAV:   4    -0.184757105443E+08   -0.18355E+08   -0.35067E+05 21264   0.159E+02
DAV:   5    -0.302980425897E+09   -0.28450E+09   -0.68095E+08 20128   0.968E+02BRMIX: very serious problems
the old and the new charge density differ
old charge density:   142.00000 new  153.84313
    0.718E+02
WARNING: Sub-Space-Matrix is not hermitian in DAV

希望各位大神指点帮助一下，谢谢！

wsz

添加以下环境变量
export I_MPI_ADJUST_REDUCE=3
export MPIR_CVAR_COLL_ALIAS_CHECK=0

archer

把CHGCAR, WAVECAR删除试试

logzzz

出现这种“BRMIX: very serious problems the old and the new charge density differ” 有的时候不是因为你自己的输入文件问题。可以尝试，减小核数。或者换个其他版本计算。我一般用6.3.2出现这个问题，会用5.4.4算。一般都能解决。。这个跟软件编译，和运行环境也有关系。