计算化学公社

标题: VASP计算中小分子结构优化中能量收敛于不同值的问题 [打印本页]

作者
Author: liuxr 时间: 2021-11-26 10:23
标题: VASP计算中小分子结构优化中能量收敛于不同值的问题
各位老师好，本人是做VASP表面吸附的初学者。最进做CO分子吸附在金属Ni表面的论文复现，分别放进去两个盒子跑CO的结构优化，计算结果是二者皆收敛，但收敛于不同的能量之下。不知何故，特来请教。

作者
Author: sobereva 时间: 2021-11-26 11:25
“收敛于不同的能量之下”语义不明
说清楚是什么的能量

作者
Author: liuxr 时间: 2021-11-27 19:38

sobereva 发表于 2021-11-26 11:25
“收敛于不同的能量之下”语义不明
说清楚是什么的能量

我跑CO分子的结构优化，构形一样，就只有INCAR不一样，两个构型都能量收敛了，一个能量是-14.79565783eV,一个能量是-1.06510973eV,附上两个INCAR
第一个INCAR
SYSTEM = VWTi-PbO
ISTART = 1
ICHARG = 1
LWAVE = .TRUE.
LCHARG = .TRUE.
LVTOT = .FALSE.
LVHAR = .FALSE.
LELF = .FALSE.

ENCUT = 400
ISMEAR = 0
SIGMA = 0.05
EDIFF = 1E-6
NELMIN = 5
NELM = 300
GGA = PE
LREAL = Auto
LDIPOL = .TRUE.
IDIPOL = 3
ISYM = 0
ALGO = Normal
PREC = Normal
ADDGRID = .TRUE.

EDIFFG = -0.1
IBRION = 2
POTIM = 0.2
NSW = 500
ISIF = 2

第二个INCAR
SYSTEM = abx3

ISTART = 0
ICHARG = 2

ISMEAR = 0
SIGMA = 0.1

LREAL = Auto

IBRION = 2
ISIF=3
NSW = 500
POTIM = 0.5

EDIFFG =-0.01
EDIFF = 0.0001
ENCUT =500 eV
NELMIN=5

LCHARG = .F.
LWAVE = .F.

ALGO = Fast

IVDW = 11
VDW_RADIUS = 50.2
VDW_CNRADIUS = 20.0

作者
Author: sobereva 时间: 2021-11-28 14:13

liuxr 发表于 2021-11-27 19:38
我跑CO分子的结构优化，构形一样，就只有INCAR不一样，两个构型都能量收敛了，一个能量是-14.79565783eV, ...

我不用VASP。这种问题你就逐条对比测试找原因就完了
而且截断能都不一样有什么可比的

欢迎光临计算化学公社 (http://bbs.keinsci.com/)