VASP用于进行标准的自洽场(SCF)计算,INCAR文件配置如下:
INCAR文件
# INCAR 自洽场计算
# 初始化选项
ISTART = 0 # 从头开始计算
ICHARG = 1 # 从POTCAR初始化电荷密度
# 电子占据状态展宽
ISMEAR = 0 # 费米-狄拉克分布
SIGMA = 0.05 # 展宽参数 (eV)
# 平面波截断能
ENCUT = 500 # 平面波截断能 (eV)
# 自洽循环收敛准则
EDIFF = 1E-6 # 能量收敛阈值 (eV)
NELM = 50 # 最大自洽循环次数
# 实空间加速度
LREAL = AUTO # VASP自动选择实空间加速度
# 自旋极化
ISPIN = 1 # 不考虑自旋极化
# 局部轨道信息输出
LORBIT = 11 # 输出每个原子的投影态密度 (PDOS)
# 算法选择
ALGO = Fast # 使用快速算法
# 精度设置
PREC = Normal # 使用标准精度
# 结构优化
NSW = 0 # 不进行结构优化
ISIF = 2 # 优化原子位置 (仅当NSW > 0时有效)
# 结构优化收敛准则
EDIFFG = -0.02 # 原子力收敛阈值 (eV/Å)
# 波函数和电荷密度保存
LWAVE = .TRUE. # 保存波函数
LCHARG = .TRUE. # 保存电荷密度
# 离子步长控制
POTIM = 0.5 # 离子步长 (Å)
# 电荷守恒
NELECT = 0 # 自动计算电子数
NELMIN = 4 # 最小电子步数
# 其他选项
NELMDL = -4 # 最小电子步数 (用于混合)
NWRITE = 1 # 写入波函数的频率
参数解释
1. ISTART:
设置VASP从头开始计算还是从波函数文件(WAVECAR)继续计算。
0: 从头开始计算。
1: 从WAVECAR文件继续计算。
2. ICHARG:
设置电荷密度的初始猜测。
1: 从POTCAR初始化电荷密度。
11: 从WAVECAR文件读取波函数。
3. ISMEAR:
斯梅尔函数的类型,用于电子占据状态的布罗伊特-温纳函数的展宽。
0: 费米-狄拉克分布(默认)。
5: 高斯展宽,适用于金属系统。
4. SIGMA:
斯梅尔函数的宽度。
0.05 eV: 适用于大多数情况。
对于绝缘体,可以设置为较小的值,例如0.01 eV。
5. ENCUT:
平面波截断能,定义了基组的最大能量。
推荐设置为至少比默认值大20%。
默认值由POTCAR文件给出,可以通过ENCUT = x * POTCAR来确定,其中x通常取1.3到1.5之间。
6. EDIFF:
自洽循环收敛准则的能量差阈值。
1E-6 eV: 适用于大多数情况。
对于高精度计算,可以设置为更小的值,例如`1E-8 eV`。
7. NELM:
最大自洽循环次数。
50: 适用于大多数情况。
可以根据实际需要调整。
8. LREAL:
控制是否使用实空间加速度。
.FALSE.: 关闭实空间加速度。
AUTO: VASP自动选择。
9. ISPIN:
自旋极化选项。
1: 不考虑自旋极化。
2: 考虑自旋极化。
10. LORBIT:
控制局部轨道信息的输出。
11: 输出每个原子的投影态密度(PDOS)。
11. ALGO:
算法选择。
Fast: 使用快速算法。
All: 使用所有可用算法。
12. PREC:
精度设置。
Normal: 使用标准精度。
Accurate: 使用高精度。
13. NSW:
最大步骤数,用于结构优化。
0: 不进行结构优化,只做单点能计算。
14. ISIF:
结构优化的自由度。
2: 优化原子位置。
3: 优化原子位置和体积。
15. EDIFFG:
结构优化的收敛准则。
-0.02 eV/Å: 适用于大多数情况。
对于更严格的优化,可以设置为更小的值。
16. LWAVE 和 LCHARG:
决定是否保存波函数和电荷密度。
.TRUE.: 保存波函数。
.TRUE.: 保存电荷密度。
17. POTIM:
离子步长控制,仅在进行结构优化时有效。
0.5: 适用于大多数情况。
18. NELECT:
电子数,设置为0时自动计算电子数。
19. NELMIN:
最小电子步数。
20. NELMDL:
最小电子步数(用于混合)。
21. NWRITE:
写入波函数的频率。
无解释简洁版:
ISTART = 0
ICHARG = 1
ISMEAR = 0
SIGMA = 0.05
ENCUT = 500
EDIFF = 1E-6
NELM = 50
LREAL = AUTO
ISPIN = 1
LORBIT = 11
ALGO = Fast
PREC = Normal
NSW = 0
ISIF = 2
EDIFFG = -0.02
LWAVE = .TRUE.
LCHARG = .TRUE.
POTIM = 0.5
NELECT = 0
NELMIN = 4
NELMDL = -4
NWRITE = 1
总结
以上就是一份典型的用于自洽场计算的INCAR文件配置。这些参数覆盖了VASP中进行SCF计算所需的基本设置。可以根据具体的需求调整这些参数,例如增加精度、改变展宽参数或开启自旋极化等。