Thanks!
I want to further ask about the density matrix. For a two-probe system with different electrodes, how is the total charge calculated? By the formula Q=Tr(DS) (where D is the density matrix and S is the basis overlap matrix), just like in Siesta? Or just by integrating the electron density in the real space?