My Community

If I am calculating band- (and k-) decomposed charge density within certain limits with respect to the Fermi energy, at a single k-point (Gamma), and band number N minority spin is within these limits, but the band number N majority spin is not, will my PARCHG file contain spin-up+spin-down components of charge anyway? Because within those limits it should only contain minority spin charge. But ...

the second set in the PARCHG file ((spin-up)-(spin-down) charge) has both negative and positive values, indicating magnetization density, for a single band.

Any help is greatly appreciated.

yes, the current implementation of the partial band-decomposed charge density adds up the density of the n-th band of up and down spin, irrespective of whether both or just one of the bands are in the specified energy range or not.

but if you used the specific bands, the parchg file only contain the spinup component, do not have the spindown . the two components in parchg file was same as each other ,thy are both the spinup component ,why ?

@ xinhuanwen,

You need to get the fixed parsden.F and then recompile VASP. It was a bug.

Cheers!

We're seeing a similar problem (both grids printed out in PARCHG.00NM have the same information for band projected densities), but only for hybrid functionals. For PBE, it seems to behave sensibly (first grid is consistent with total density, and second with magnetization). It seems to be related to the fact that when you're doing a PAW calculation without exact exchange ISYM is set to 2 by default, but for hybrid functionals it's set to 3.

Has anyone else seen this problem? Would an example input showing the issue help?

thanks,
Noam

P.S. I'm using the latest 5.2.11 version from the ftp site (date Jan 18 2011). The web portal version is dated Nov 2010.