Hi!
I tried many possible combinations of (N_k points, mesh cut_off, LDA/GGA, Temperature, Integral lower bound), al of which fail to avoid the disappearance of charge. As I describe before, only constrained SCF calculation for the real space density gives nonzero chage! You'd better tryo to fix this problem.