My Community

Dear vasp users

I have done a Band decomposed charge density (patial charge) calculation, and I think the data of grid points refer to the square of absolute values of wavefunctions.

It should all be positive,but I found that there are negative values in PARCHG.

Did I misunderstand something?

Any suggestions are welcomed!

Thanks

Yes you have understood correctly, a density should be positive definite...however (and this is the twist) the densities in the PARCHG file (just as for the CHGCAR file) come from pseudo wavefunctions, which can result in negative densities at the atomic centers (this is due to the construction of PAW). Another source of negative values is just simply noise. You will find in the vacuum reqions spaces where the density oscillates between positive and negative values...these are values that you should consider to be "zero" .

Danny

[quote="Danny"]Yes you have understood correctly, a density should be positive definite...however (and this is the twist) the densities in the PARCHG file (just as for the CHGCAR file) come from pseudo wavefunctions, which can result in negative densities at the atomic centers (this is due to the construction of PAW). Another source of negative values is just simply noise. You will find in the vacuum reqions spaces where the density oscillates between positive and negative values...these are values that you should consider to be "zero" .

Danny[/quote]

Thank you, Danny! It does help me very much!!