VASP计算中小分子结构优化中能量收敛于不同值的问题

liuxr

各位老师好，本人是做VASP表面吸附的初学者。最进做CO分子吸附在金属Ni表面的论文复现，分别放进去两个盒子跑CO的结构优化，计算结果是二者皆收敛，但收敛于不同的能量之下。不知何故，特来请教。

sobereva

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

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

liuxr

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

sobereva

“收敛于不同的能量之下”语义不明
说清楚是什么的能量