第7章｜2026新版VASP吸附与催化教程：vasp差分电荷密度计算解析！

差分电荷密度的物理意义与理论基础

差分电荷密度是材料电子结构分析中最重要的可视化工具之一，它能够清晰地揭示化学键形成、电荷转移、吸附相互作用、界面耦合等电子层面的物理过程。

差分电荷密度的定义

差分电荷密度定义为复合体系电荷密度与组成组分电荷密度之和的差值：

其中：

为复合体系（如吸附体系、界面、缺陷等）的总电荷密度；

为各独立组分（原子、分子、片段）在复合体系构型下的电荷密度；

 n为组分数目

该定义源自VASP官方文档和文献研究，反映了化学键合过程中电荷的重新分布。

差分电荷密度的物理图像

差分电荷密度的主要应用场景

VASP差分电荷密度计算工作流程

标准三步计算法

差分电荷密度计算需要完成三次独立的VASP计算：

# ============ 计算1：复合体系 ============ 

# INCAR配置 

SYSTEM = Composite System 

PREC = Accurate 

ENCUT = 450 

ISPIN = 2          # 如需自旋极化 

LCHARG = .TRUE.    # 输出CHGCAR 

LWAVE = .FALSE. 

EDIFF = 1E-6 

# 结构文件：POSCAR_system 

# ============ 计算2：孤立组分1 ============ 

# INCAR配置（与计算1相同参数） 

SYSTEM = Component 1 

PREC = Accurate 

ENCUT = 450 

ISPIN = 2 

LCHARG = .TRUE. 

ISMEAR = 0 

SIGMA = 0.05 

EDIFF = 1E-6 

# 结构文件：POSCAR_component1 

# 注意：保持与复合体系相同的晶胞参数 

# ============ 计算3：孤立组分2（如需） ============ 

# INCAR配置同上 

# 结构文件：POSCAR_component2 

INCAR关键参数设置

晶胞与K点设置要求

关键注意事项：

1.晶胞尺寸必须一致：所有三种计算的晶胞参数必须完全相同

2.K点网格必须一致：确保相同的布里渊区采样

3.赝势版本一致：POTCAR文件不能混用

4.ENCUT一致：平面波截断能必须相同

5.对称性设置一致：ISYM参数应相同

# K点设置示例（三种计算需完全相同） 

KPOINTS

Auto generation 

0

Monkhorst-Pack

9 9 9 

0 0 0 

差分电荷密度可视化工具

VESTA可视化操作指南

步骤1：加载差分电荷密度文件

# 启动VESTA 

vesta CHGCAR_diff 

步骤2：设置等值面显示

Edit → Edit Data → Register 

Data category: Charge density 

# 设置等值面值 

Value: 0.005 (e/ų)     # 典型值范围0.001-0.01 e/ų 

Isosurface: 显示等值面 

Slicing: 生成切片图 

步骤3：颜色映射

Edit → Properties → Isosurface 

Color:  

– 正差分：红色/蓝色（电荷积累） 

– 负差分：蓝色/黄色（电荷消耗） 

Contour Level: 15 (等值面层级数) 

可视化技巧：

等值面：显示电荷密度的三维分布

切片图：观察特定平面的二维分布

透明度调节：叠加原子位置理解空间关系

截面分析：沿特定方向生成剖面图

本章要点总结

VASP差分电荷密度计算知识点梳理

差分电荷密度的物理意义与理论基础：差分电荷密度的定义、物理图像、主要应用场景

VASP差分电荷密度计算工作流程：标准三步计算法、INCAR关键参数设置、晶胞与K点设置要求

差分电荷密度可视化工具：加载差分电荷密度文件、设置等值面显示、颜色映射

下一步学习建议

下一章将正式引入本次教程的核心—VASP结合能计算。我们将从结合能计算的基本原理、结合能计算的前期准备、结合能计算完整流程方面详细介绍VASP结合能计算。