差分电荷密度（Charge Density Difference, CDD）是计算材料科学中最常用的电子结构分析手段，通过定量展示成键 / 吸附过程中电荷的重新分布，直观揭示化学键本质、电荷转移方向与强度、催化活性位点作用机制等核心问题。本文系统讲解 VASP 平台下 CDD 的计算原理、标准流程、参数设置与结果解读。

一、基本原理与核心应用

1. 定义

差分电荷密度是总体系电荷密度减去各独立片段电荷密度的差值，数学表达式为：

Δρ = ρ_total​−∑ρ_{fragment​​}

其中，ρ_total​：优化后完整体系的电荷密度；ρ_fragmenti​​：保持总体系中原子位置不变的各个独立片段的电荷密度

2. 核心应用

• 化学键类型判断：离子键（电荷完全转移）、共价键（成键区电荷积累）、范德华作用（弱电荷转移）
• 催化反应分析：反应物与催化剂之间的电荷转移方向与数量，识别活性位点
• 异质结界面研究：界面处的电荷转移与内建电场形成
• 缺陷 / 掺杂效应：缺陷 / 掺杂原子对周围电子结构的影响

差分电荷密度物理意义示意图 [DOI: 10.1016/j.jechem.2026.04.027]（黄色区域表示电荷积累Δρ>0，青色区域表示电荷耗尽Δρ<0；直观展示了 A-B 成键过程中电子从 A 转移到 B 并在成键区积累的过程）

 

二、标准计算流程（三步法）

第一步：结构优化（必做）

首先对目标体系进行充分的结构优化，得到稳定的原子构型。

输入文件：INCAR、POSCAR、KPOINTS、POTCAR

关键参数：EDIFF=1E-6、EDIFFG=-0.01、ISIF=3（体相）或ISIF=2（表面 / 分子）

输出文件：CONTCAR（后续所有计算的结构输入）

 

第二步：总体系静态电荷密度计算

使用优化后的结构（CONTCAR 重命名为 POSCAR）进行高精度静态自洽计算，输出总电荷密度。

核心 INCAR 参数设置incarSYSTEM = Total system charge density

ENCUT = 520          # 比结构优化高10-20%，保证电荷密度精度

ISTART = 0           # 从头开始计算

ICHARG = 2           # 从原子电荷密度初始化

ISMEAR = 0           # 半导体/绝缘体用0，金属用1

SIGMA = 0.05

PREC = Accurate      # 高精度模式

LAECHG = .TRUE.      # 输出AECCAR0(芯电子)和AECCAR2(价电子)

LCHARG = .TRUE.      # 输出CHGCAR(总电荷密度)

LWAVE = .FALSE.      # 不输出波函数，节省空间

NSW = 0              # 不进行结构优化

 

输出文件：CHGCAR、AECCAR0、AECCAR2

 

第三步：各独立片段的静态电荷密度计算

保持总体系的晶格参数和原子位置不变，分别计算每个独立片段的电荷密度。

 

关键操作：

复制总体系的 INCAR、KPOINTS、POTCAR

修改 POSCAR，只保留对应片段的原子，其余原子删除

保持LAECHG=.TRUE.和LCHARG=.TRUE.参数不变

 

注意事项：所有片段的计算参数必须与总体系完全一致（ENCUT、K 点、赝势等）

 

三、差分电荷密度计算与可视化

1. 电荷密度相减（VASPKIT 一键完成）

VASPKIT 是处理 VASP 输出的常用工具，提供了差分电荷密度的自动计算功能：

 

将总体系的 CHGCAR 命名为CHGCAR_total

将各个片段的 CHGCAR 分别命名为CHGCAR_frag1、CHGCAR_frag2...

运行 VASPKIT，输入功能号113（Charge Density Difference）

按照提示输入总体系和片段的 CHGCAR 路径

输出文件：CHGDIFF.vasp（差分电荷密度文件）

 

2. 价电子差分电荷密度（推荐）

若需要排除芯电子的干扰，计算价电子差分电荷密度：

 

总体系价电子电荷密度：ρ_{total,valence}​ = AECCAR2_total​

各片段价电子电荷密度：ρ_{frag,valence​}=AECCAR2_frag​

价电子差分：Δρ_valence​=ρ_{total,valence}​−∑ρ_frag,valence​

 

 

VASPKIT 功能号114可直接计算价电子差分电荷密度

 

3. 可视化（VESTA）

 

用 VESTA 打开CHGDIFF.vasp

点击Edit → Properties → Isosurfaces

设置等值面数值：通常从0.01 e/ų开始调整，根据体系差异优化

设置颜色：正等值面（电荷积累）设为黄色，负等值面（电荷耗尽）设为青色

调整显示方式：可同时显示原子结构、键和等值面

 

四、关键注意事项与常见问题

1. 片段选取原则

必须保持片段中原子的位置与总体系完全一致，不能重新优化

对于表面吸附体系，通常分为 "催化剂基底" 和 "吸附分子" 两个片段

对于双原子催化剂，可分为 "载体" 和 "双原子活性中心" 两个片段

避免片段之间存在重叠或空隙

 

2. 计算精度控制

ENCUT 必须足够大，否则电荷密度会出现严重的数值噪声

K 点密度要高于结构优化，建议使用Gamma-centered网格

对于大体系，可适当降低 K 点密度，但需保证收敛性

 

3. 结果解读误区

 

不要仅凭等值面的大小判断电荷转移的多少，需结合 Bader 电荷分析定量计算

等值面数值的选择会显著影响图的外观，需在文中明确标注所用的等值面大小

差分电荷密度只能展示电荷的相对变化，不能给出绝对电荷密度

 

4. 常见错误排查

错误 1：片段计算时修改了晶格参数 → 导致电荷密度无法对齐

错误 2：总体系和片段使用了不同的赝势 → 电荷密度单位不一致

错误 3：LAECHG=.FALSE. → 无法输出 AECCAR 文件，无法计算价电子差分

错误 4：等值面数值过小 → 出现大量无意义的噪声点