My Community

If I understand correctly, you specify in array IBAND the number of the bands whose charge density is to be printed. However, how do you select among spin up or spin down bands in spin-polarized calculations?

The resulting PARCHG* file contains charge density grid for spin_up_+_spin_down and spin_up_-_spin_down . You need to add these two and divide by 2 to get spin_up or subtract to get spin_down...

The resulting PARCHG* file contains charge density grid for spin_up_+_spin_down and spin_up_-_spin_down . You need to add these two and divide by 2 to get spin_up or subtract to get spin_down...

Does it mean that the data written to the PARCHG file correspond to the sum of the charge density of the spin up band number n and the charge density of the spin down band number n and to their difference?

Actually, the second part of the density written to the PARCHG file is the same as the first part so the substraction gives zero.

They are not. You should check your read statement