I feel the values of the partial charges in Job.out are somewhat strange ?
VASP evaluates charges only within a specific radius – the Wigner-Seitz (WS) radius. The total charge is the sum of all WS spheres and the contribution from the interstitial space. You can set the WS radius for each atom by specifying the parameter RWIGS in the Add to Input panel of the VASP user interfaces.