很多做界面计算的人，跑完结构优化之后顺手就做了一个差分电荷密度。LEV00 = .TRUE.，开个 VESTA 看看等值面，红蓝两色一标，图片放进文章里，好像就算分析过了。
但差分电荷密度真正能回答的问题，远比”红色区域有电荷聚集”要多。把这个分析做透，需要回到计算流程里多问几层。
差分电荷密度的定义本质上很简单：总体系的电荷密度，减去各部分单独存在时的电荷密度。问题在于”各部分”怎么定义。
如果你算的是 MoS₂ 吸附 Li 原子的体系，总电荷密度是 MoS₂+Li 的，参考体系是独立的 MoS₂ 层和独立的 Li 原子。这一步很多人会做，但容易忽略一个关键点：参考体系必须使用与总体系相同的晶格和截断能。 如果总体系用了 15Å 真空层的超胞，参考体系的盒子大小必须和总体系完全一致。
不止一次看到有人把 Li 原子单独算在一个小盒子里，然后用 VESTA 做减法，结果等值面上全是周期边界造成的虚假条纹。这不是数据有问题，是参考体系选错了。
更隐蔽的陷阱在于表面和吸附质的”几何冻结”。总体系里的 MoS₂ 在吸附 Li 之后有弛豫，但差分电荷密度减掉的是”未吸附时的独立 MoS₂”——这个独立 MoS₂ 用的是弛豫前的构型还是弛豫后的构型？如果用了弛豫前的构型，差分出来的电荷重排就混入了”晶格变形引起的电荷变化”和”界面电荷转移引起的电荷变化”两种效应。这两种效应在物理上完全不同，但等值面不会告诉你哪个是哪个。正确的做法是：参考体系中的各组分，使用总体系弛豫后冻结的几何构型进行单点计算。
计算差分电荷密度，很多人只在总体系里加上 LCHARG = .TRUE. 和 LAECHG = .TRUE.。但如果你用的是 PBE+U、范德华修正或偶极修正，参考体系的计算必须使用完全一致的设置。
LDAU 参数是第一个容易出问题的地方。MoS₂ 里的 Mo 加了 U 修正，参考体系里单独算 MoS₂ 层时也必须加同样的 U。参数不一致带来的电荷密度差异，会被 VESTA 忠实地渲染成红蓝两色，但你无法判断这个差异是物理上的电荷转移还是参数差异造成的。
偶极修正（LDIPOL = .TRUE.）是第二个容易忽略的点。总体系在 z 方向加了偶极修正，参考体系如果是独立的 MoS₂ 层，在有限真空层下它的电荷密度分布也会受到虚假偶极的影响。去掉偶极修正、重新跑参考体系的单点计算，看似多了一步，实际上避免了等值面上虚假的 z 方向电荷堆积。
NGXF/NGYF/NGZF 的 grid 大小必须在三个体系里完全一致。很多人让 VASP 自动决定 FFT grid，结果三个体系的 grid 尺寸不一样，用程序做减法之前还要插值。插值本身就引入了数值误差，在等值面阈值较低（如 0.002 e/ų）时尤为明显。手动设置 NGX/NGY/NGZ 或用 PREC = Accurate 迫使用相同 grid，能让后续分析少一层不确定性。
VESTA 打开差分电荷密度文件，默认等值面阈值是 0.05 e/ų。这个值在很多体系里偏大——界面电荷重排的典型量级在 0.01-0.03 e/ų 之间，用 0.05 的阈值几乎看不到明显的界面电荷特征。
靠谱的做法是从低阈值（0.005）开始逐步增加，观察等值面的形态变化。如果等值面在 0.01 时呈连续片状、到 0.03 时收缩为团簇、到 0.05 时消失——那说明电荷重排的强度大概在 0.02-0.03 量级，最终的展示阈值选在 0.01 左右更合理。
另外，颜色映射也是个容易被低估的决策。默认的红-蓝配色（红=电荷积累，蓝=电荷耗散）在化学上没有任何问题，但对色弱读者（男性约 8%）不友好。换成蓝-白-红 diverge colormap，白色作为零基准，是一个更包容的选择。
差分电荷密度最终的目的是回答物理问题：电荷从哪里来、去了哪里、为什么会这样转移。等值面只是可视化手段，不是答案本身。
吸附 Li 原子之后 MoS₂ 的等值面在 Mo-S 之间出现了蓝色区域——这可能意味着 Mo-S 键的电荷被抽取了，也可能只是几何弛豫导致的电子云形变。不做 Bader 电荷分析或者 PDOS 对照，仅凭一个等值面就下结论说”电荷从 Mo-S 转移到 Li”，是不够严谨的。
从十几年的实操回头看，差分电荷密度是一个很好的”发现问题”的工具——它在等值面上呈现异常分布时，强迫你回到电子结构和化学键层面去追问原因。但它本身不提供”解决问题”的答案。用好它的关键，是把参考体系选对、参数设置统一、阈值取合理，然后带着这个等值面去和 Bader charge / ELF / PDOS 一起看。
