|
各位老师同学,请教个问题,我在用vasp计算FeO(111)面(O暴露面)吸附氢气时总是第一个离子步电子都不收敛 然后直接断掉, 但是计算纯FeO时虽然第一个离子步的电子不收敛 但后续能正常收敛,磁矩也较为正常,Fe原子磁矩大概为4左右,为什么表面放个氢气就这么难收敛呢,主要的参数都一样的?INCAR如下:
SYSTEM = VB2
ISTART = 0
ICHARGE = 2
PREC = Normal
LREAL = Auto
ALGO = Normal
IBRION = 2
ISIF = 2
NSW = 500
POTIM = 0.1
EDIFFG = -0.02
ISPIN = 2
LDAU= .TRUE.
LDAUTYPE= 2
LDAUL= -1 -1 2
LDAUU= 0 0 4
LDAUJ= 0 0 0
MAGMOM = 2*0 20*0 16*5
LDIPOL = .TRUE.
IDIPOL = 3
LREAL = .TRUE.
LMAXMIX = 4
AMIX = 0.2
BMIX = 0.0001
AMIX_MAG = 0.8
BMIX_MAG = 0.0001
NELM = 150
ENCUT = 500
EDIFF = 1E-6
ISMEAR = 0
SIGMA = 0.1
LCHARG = .F.
LWAVE = .F.
NELMIN = 4
NPAR = 8
IVDW = 11
|
|