差分电荷密度计算：揭示材料界面电荷重分布的可视化方法

差分电荷密度计算的核心思想是”做减法”——从复合体系的总电荷密度中减去各组分独立存在时的电荷密度，剩余的部分就是由于界面相互作用导致的电荷重分布。这个项目面对的是一个MoS₂/石墨烯范德华异质结，实验团队需要理解两层之间的是否存在电荷转移以及转移方向和量级。

参考态构建：计算成败的基石

差分电荷密度计算的首要难题不在减法本身，而在于参考态的精确构建。公式表达很简洁：Δρ = ρ(AB) – ρ(A) – ρ(B)，但三个ρ必须在相同的晶胞和原子位置下计算。这意味着孤立MoS₂层和孤立石墨烯层需要在异质结的晶胞中保持与复合体系完全一致的原子坐标，不能做结构弛豫。

项目组在初次尝试时犯了一个代价不菲的错误：对孤立层分别做了结构弛豫。弛豫后层内键长变化约0.02Å，这个微小位移导致差分图中出现了大量虚假的电荷重分布信号，淹没了真实的界面效应。发现这个问题耗费了将近一周时间。修正方法是冻结异质结中的原子位置，仅对孤立层计算自洽电子密度，确保三组计算共享同一套原子坐标。

计算参数与收敛控制

所有计算在VASP中完成，采用PBE泛函配合DFT-D3色散校正。色散校正在这里不是可选项——范德华异质结的层间距由弱相互作用主导，纯GGA计算会显著高估层间距，进而影响电荷密度的重叠区域。截断能设为520eV，k点网格取12×12×1，电子自洽收敛标准设为10⁻⁶ eV。

高收敛标准并非过度追求精度。差分电荷密度的信号量级通常在10⁻³至10⁻² e/ų范围内，如果自洽收敛不够严格，差分结果中会出现数值噪声，等值面呈现颗粒状的”盐椒”伪影。Physical Review B中相关研究的经验表明，10⁻⁶ eV收敛判据足以保证差分信号的可信度。

等值面可视化与数据解读

差分电荷密度的三维等值面是解读电荷重分布最直观的工具。正值区域表示复合体系比孤立组分电荷密度更高，电荷在此聚集；负值区域表示电荷耗散。本项目取±0.003 e/ų的等值面，正值用黄色显示，负值用蓝色显示。

结果显示，MoS₂层与石墨烯界面处出现明显的电荷重分布：石墨烯一侧呈现蓝色耗散区，MoS₂一侧呈现黄色聚集区。电荷转移方向从石墨烯指向MoS₂，转移量约0.012e每单胞。这个数值与Chemistry of Materials期刊报道的同类体系结果（0.008-0.015e）高度吻合。

面平均差分电荷密度

三维等值面直观但不便于定量分析。项目组进一步计算了沿z方向的面平均差分电荷密度，即将xy平面内的差分电荷密度积分，得到Δρ随z坐标的变化曲线。这种降维处理将复杂的三维分布压缩为一维曲线，使得不同体系之间的定量比较成为可能，也使得电荷转移的方向和量级一目了然。

曲线在界面区域出现一个明显的正峰和一个负峰，峰间距约3.4Å，与层间距一致。正峰积分面积减去负峰积分面积得到净转移电荷0.012e，与三维等值面的积分结果一致。两组数据相互验证，增强了结论的可信度。

方法局限与适用边界

差分电荷密度计算的局限性值得清醒认识。差分值是纯数学运算的结果，它不直接区分共价、离子或范德华贡献——所有类型的电荷重分布在差分图中都以相同的正负值呈现。区分成键性质需要结合电子局域函数（ELF）、Bader电荷分析等其他方法。

回过头看，差分电荷密度计算的价值在于它提供了一个直接、可视的电荷重分布图像，帮助研究者快速定位界面活性区域。当它与Bader分析、投影态密度等工具联合使用时，才能拼凑出完整的物理图景。值得警醒的是，差分电荷密度对参考态的构建极度敏感，任何坐标偏差或基组不一致都会引入系统性误差——这个看似简单的”减法”容错空间极小。