I have a system, where I get a not logical result in bader charge analysis. If I plot the CHGCARs for the same calculation, it seems to be campletly different, than what the bader analysis says.
My system is a Graphene layer with a Sodium adsorbed. The Sodium is in the center of a ring. Bader analysis says, that the electron from the sodium in transfered to only 3 of the carbons, but if I look at the electron density, it seems t be transferred to all 6 carbons.
The bader volumes seem to be almost the same size at every carbon.
You could compute the DDEC6 net atomic charges. These are the state-of-the-art method for charge assignment across diverse material types. The free Chargemol program (http://ddec.sourceforge.net) can compute these from the VASP files (AECCAR0, AECCAR2, CHGCAR, and POTCAR). See the following journal articles for details: Article Introducing DDEC6 atomic population analysis: Part 3. Compre...
,
Article Introducing DDEC6 Atomic Population Analysis: Part 1. Charge...
The step by step procedure to perform the Bader charge analysis using VASP and charge density difference plot using VESTA is clearly explained in a Youtube channel. You can find it here: