求助：差分电荷密度组分分解以及与实验值不符的问题

jio

我在计算FeCoSNC双原子催化剂的差分电荷密度，引入了一个-OH预吸附，主要想看看Fe、Co与S之间的电荷转移。差分电荷密度公式：Δρ = ρ_system - ρ _atom1 - ρ_atom2 - ... - ρ_atomN，我应该如何分解组分比较好呢？初步打算是Δρ = ρ_system - ρ _without(FeCo) - ρ_without(S)，这样对不对呢？
以及，XPS测试结果表明，S向Fe提供电子，但是我的差分电荷密度和Bader电荷分析却表明Fe失电子，S得电子，这是不是算错了呢？应该如何修正呢？以下是我的INCAR文件：
SYSTEM =  qvasp

INCAR-basic
   PREC   = Normal
   ENCUT  = 520 eV
   NELMIN= 5
   LREAL = Auto

   ALGO = Fast
   EDIFF  = 1E-4
   ISMEAR = 0
   ISPIN = 2
   MAGMOM = 50*0 6*0 1*-1.0 1*5.0 1*0 1*0 1*0
   LORBIT = 11
   SIGMA = 0.1
   NPAR = 4     ## important to reduce calculate time, equal sqrt(cpu cores)

Ionic Relaxation
   ICHARG = 2
   EDIFFG = -0.03        force
   NSW    = 0          max steps for IOM
   IBRION = -1
   POTIM = 0.1
   ISIF   = 2            suface:2 crystal:3
   PSTRESS= 0
   LCHARG = .TRUE.
   LWAVE  = .FALSE.
   LAECHG = .TRUE.

DFT-D3
   IVDW = 12

DFT+U
   LDAU = .TRUE.
   LDAUTYPE = 2
   LDAUL = -1 -1 2 2 -1 -1 -1
   LDAUU = 0.0 0.0 4.0 3.0 0.0 0.0 0.0
   LDAUJ = 0.0 0.0 0.0 0.0 0.0 0.0 0.0
   LMAXMIX = 4
   LASPH = .TRUE.
   LDAUPRINT = 2

山里面的和尚

如果按照电负性和极性来看的，电荷密度差算出来Fe失电子，S得电是正常的