静电势是描述电荷分布产生的电势能分布的物理量，其数学表达式为：
其中：
为电子密度和原子核电荷密度的总和；
为距离
即为空间任意点r 处的静电势
在VASP计算中，静电势主要包含以下成分：
|
成分 |
说明 |
|
Hartree势 |
电子–电子相互作用产生的势场 |
|
离子势 |
原子核（赝势）产生的势场 |
|
交换关联势 |
电子交换关联效应产生的势场 |
|
总静电势 |
上述所有成分的综合 |
LOCPOT文件是VASP输出静电势信息的标准格式文件，包含三维空间的静电势分布数据：
# LOCPOT格式说明
第一行：系统标题
第二行：缩放因子
第三至五行：晶格基矢
第六行：原子数量
第七行至原子信息行：原子位置
后续行：静电势数据（3D网格数据）
在INCAR中设置以下参数以生成LOCPOT文件：
|
参数 |
推荐设置 |
说明 |
|
LVTOT |
.TRUE. |
输出总静电势到LOCPOT |
|
LVHAR |
.TRUE. |
输出Hartree势到LOCPOT |
|
LOCPOT |
.TRUE. (旧版本) |
输出LOCPOT文件 |
|
LAECHG |
.TRUE. |
输出赝势局域电荷密度 |
推荐参数组合：
LVTOT = .TRUE. # 输出总静电势
LVHAR = .TRUE. # 输出Hartree势
LOCPOT = .TRUE. # 生成LOCPOT文件（VASP 5.x及以下版本）
VASP 5.3.3及以下：使用LOCPOT参数生成文件
VASP 5.4及以上：默认生成LOCPOT文件，可通过LVTOT/LVHAR控制具体内容
平面平均电势是静电势分析最常用的方法，将三维静电势投影到特定方向形成一维曲线：
处理流程：
1.从LOCPOT文件中提取三维静电势数据
2.沿选定方向（通常为z轴）进行平面平均
3.得到沿该方向的电势分布曲线
工具一：Multiwfn
Multiwfn是功能强大的分子波函数分析软件，可处理LOCPOT文件：
# 启动Multiwfn
Multiwfn LOCPOT
# 选择静电势分析模块
选择静电势相关功能
选择平面平均或等值面分析
工具三：vasp2plot / VESTA
VESTA软件处理流程：
1.导入LOCPOT文件
2.使用Volume/Isosurface功能可视化
3.可生成二维切片或三维等值面图
工具三：自定义Fortran/C程序
原VASP官方提供vtotav.f程序用于平面平均电势计算，也可使用Fortran 95重写版本：
program vtotav
! 读取LOCPOT文件
! 计算平面平均电势
! 输出到potential.dat
end program
功函数计算是静电势分析最重要的应用之一：
计算步骤：
1.构建slab模型，设置足够厚的真空层（≥15 Å）
2.计算总静电势（LVTOT=.TRUE.）
3.平面平均电势在真空区域形成平台
4.真空平台值即为
5.从OUTCAR提取费米能级
6.两者相减得到功函数
偶极修正的重要性：
对于不对称slab模型，必须设置：
LDIPOL = .TRUE.
IDIPOL = 3
DIPOL = 0.0 0.0 0.5 # 偶极中心位置
在异质结、界面体系中，静电势分析可揭示界面势垒高度：
分别计算两侧材料的真空能级
分析界面区域的电势变化
确定能带偏移（band alignment）
界面势垒 = E_vacuum(interface) – E_vacuum(bulk)
吸附分子后，局部电势变化可指示：
吸附位点稳定性：电势低谷对应稳定吸附位
电荷转移方向：通过电势梯度判断
键合类型判断：离子型vs共价型
基本物理概念：静电势定义、物理意义
后处理方法：平面平均电势、常用后处理工具
具体应用场景：功函数计算、界面势垒分析、表面吸附位点识别
下章将正式引入本次教程的核心—VASP态密度计算。我们将从态密度物理意义与定义、结构优化、INCAR参数设置详解方面详细介绍VASP态密度计算。
