科研避坑指南：VASP模拟STM时，你的INCAR文件可能设错了这3个参数

科研避坑指南：VASP模拟STM时INCAR文件三大关键参数解析

在表面科学和材料表征领域，扫描隧道显微镜(STM)模拟已成为理论验证实验的重要手段。许多研究者在使用VASP进行STM图像模拟时，往往能够顺利生成PARCHG文件，却在最终图像对比环节发现与实验数据存在显著偏差。问题的根源通常不在于计算流程本身，而在于几个关键参数的物理意义理解不足。本文将聚焦三个最易被忽视却直接影响模拟精度的INCAR参数，帮助您从"能计算"进阶到"会分析"。

1. EINT能量区间：扫描偏压的物理对应关系

EINT参数定义了STM模拟中能量积分的上下限，直接决定了模拟图像的电子态来源范围。许多研究者简单地将其设置为文献中的常见值（如-1.25 eV），却忽略了这一参数与实际STM实验中扫描偏压的对应关系。

典型错误配置：

EINT = -1.25 # 仅设置单边能量区间

优化建议配置：

EINT = -1.25 0.5 # 明确指定能量下限和上限

当仅设置单个数值时，VASP默认将另一端设为费米能级（0 eV）。这种简化处理可能导致：

1. 正负偏压模拟不对称：STM实验中正负偏压反映不同电子态，单边设置无法完整模拟
2. 表面态漏检：某些材料的表面态可能分布在特定能量区间
3. 与实验条件不匹配：实际STM扫描通常记录特定偏压范围内的信息

表：EINT参数设置对STM图像的影响对比

提示：建议先通过DOS计算确定感兴趣的能态分布范围，再设置EINT参数。对于未知体系，可尝试从±1 eV开始，逐步扩大范围观察图像变化。

2. NBMOD参数：-3与其他选项的深层差异

NBMOD参数控制着分波电荷密度的计算方式，其不同取值对应着截然不同的物理图像。许多用户直接套用教程中的-3设置，却不理解其背后的物理含义。

参数选项解析：

• NBMOD = -3：计算指定能量区间内所有本征态的叠加电荷密度（最接近真实STM图像）
• NBMOD = -2：计算单个能带的电荷密度（用于特定能带分析）
• NBMOD = -1：计算单个k点的电荷密度（特殊k点分析）
• NBMOD = 正整数n：计算前n个能带的电荷密度

常见误区：

1. 认为所有NBMOD负值设置效果相似
2. 不清楚-3与其他选项的适用场景差异
3. 忽略NBMOD与EINT的协同作用

表：NBMOD不同设置的应用场景

# p4vasp中查看能带范围的实用代码片段 from p4vasp import * system = System() system.read('vasprun.xml') print("能带范围：", min(system.eigenvalues), "to", max(system.eigenvalues))

3. p4vasp后处理：恒高与恒流模式的参数映射

获得PARCHG文件后，p4vasp中的处理方式同样关键。许多研究者对"恒高模式"(constant height)和"恒流模式"(constant current)的物理对应关系认识模糊，导致图像解读错误。

核心区别：

• 恒高模式：直接显示特定高度截面的局域态密度(LDOS)
  • 对应实验中的固定针尖高度扫描
  • 图像衬度反映电子态密度差异
• 恒流模式：模拟保持隧道电流恒定时的针尖高度变化
  • 对应大多数常规STM实验条件
  • 图像反映等电子态密度面形貌

关键参数设置：

1. 亮度/对比度调节：应参考实验图像的强度范围
2. 色阶范围：建议初始设置为±10% LDOS变化
3. 平滑处理：适度应用高斯平滑消除数值噪声

注意：恒流模式图像需要正确设置衰减长度(decay length)，典型值在0.5-1.0 Å之间，需根据体系功函数调整。

4. 综合案例：石墨烯表面空位缺陷的STM模拟

让我们通过一个具体案例展示参数设置如何影响最终结果。假设我们要模拟石墨烯表面单空位的STM图像。

初始错误配置：

EINT = -1.0 # 仅设置负偏压 NBMOD = -3 LPARD = .TRUE.

优化后配置：

EINT = -1.0 0.8 # 包含正负偏压范围 NBMOD = -3 LPARD = .TRUE. ISTART = 1 # 重用波函数 ICHARG = 1 # 重用电荷密度

模拟结果对比：

1. 单边EINT设置：仅显示π*态特征，缺失π态信息
2. 双边EINT设置：完整呈现空位周围的电子态重分布
3. 不同NBMOD设置：
   • -3：显示所有相关电子态叠加
   • -2：仅显示单个能带贡献（不完整）

操作建议流程：

1. 完成常规静态计算，确认电子结构合理性
2. 分析DOS确定感兴趣的能量窗口
3. 设置合适的EINT和NBMOD参数
4. 生成PARCHG后，在p4vasp中：
   • 先尝试恒高模式快速评估
   • 切换到恒流模式精细调节
   • 保存多组参数结果便于对比

5. 高级技巧与疑难排查

即使参数设置正确，STM模拟仍可能遇到各种问题。以下是几个常见情况的解决方案：

图像噪声过大：

• 检查KPOINTS网格密度，建议至少11×11×1
• 增加ENCUT（通常比默认值高20-30%）
• 在p4vasp中应用适度平滑

图像与实验偏差明显：

• 确认费米能级位置是否正确（检查OUTCAR）
• 尝试调整EINT范围（逐步扩大±0.2 eV）
• 检查赝势是否适合表面计算

p4vasp显示异常：

• 确保PARCHG文件完整（检查文件大小）
• 尝试重新生成PARCHG（设置LSEPB=.TRUE.获取分波文件）
• 更新p4vasp至最新版本

表：常见问题与解决措施