VASP用于进行自洽场(SCF)计算,INCAR文件配置

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计算所需的基本设置。可以根据具体的需求调整这些参数,例如增加精度、改变展宽参数或开启自旋极化等。

 

  • 17
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值