计算得到的气相分子能量是否可用？

zhaoshf901

用VASP计算N2H4的能量，在添加IDIPOL=4后，虽然收敛，但能量异常，这样得到的能量是否可用？还是有什么问题？

INCAR
SYSTEM = N2H4
ALGO = Normal
EDIFF = 0.0001
EDIFFG = -0.02
ENCUT = 400
GGA = PE
IBRION = 2
ICHARG = 2
ISIF = 2
#ISTART = 0
ISMEAR = 0
ISPIN = 2
IVDW = 12
LCHARG = .FALSE.
LDAU = .FALSE.
LORBIT = 11
LREAL = Auto
LWAVE = .TRUE.
NELM = 60
NELMIN = 8
NSW = 50
PREC = Accurate
SIGMA = 0.05
IDIPOL = 4


OSZICAR
   1 F= -.30247596E+02 E0= -.30247596E+02  d E =-.302476E+02  mag=    -0.0003
   2 F= -.30323435E+02 E0= -.30323435E+02  d E =-.758395E-01  mag=    -0.0000
   3 F= -.30315348E+02 E0= -.30315348E+02  d E =-.677518E-01  mag=     0.0000
   4 F= -.30323074E+02 E0= -.30323074E+02  d E =-.754787E-01  mag=     0.0000
   5 F= -.30007860E+02 E0= -.30007860E+02  d E =0.315215E+00  mag=     0.0000
   6 F= -.30323594E+02 E0= -.30323594E+02  d E =-.519599E-03  mag=     0.0000
   7 F= -.30323856E+02 E0= -.30323856E+02  d E =-.261617E-03  mag=     0.0000
   8 F= -.25362899E+02 E0= -.25362899E+02  d E =0.496070E+01  mag=     0.0000
   9 F= -.30272368E+02 E0= -.30272368E+02  d E =0.512264E-01  mag=     0.0000
  10 F= -.30080011E+02 E0= -.30080011E+02  d E =0.243583E+00  mag=    -0.0000
  11 F= -.30324424E+02 E0= -.30324424E+02  d E =-.244413E+00  mag=    -0.0000
  12 F= -.30159867E+02 E0= -.30159867E+02  d E =-.798558E-01  mag=    -0.0000
  13 F= -.30324698E+02 E0= -.30324698E+02  d E =-.164831E+00  mag=    -0.0000
  14 F= -.28072250E+02 E0= -.28072250E+02  d E =0.208762E+01  mag=     0.0000
  15 F= -.30218702E+02 E0= -.30218702E+02  d E =-.214645E+01  mag=    -0.0000

喵星大佬

你可以按你的盒子大小算算密度看看这玩意还是不是气体

zhaoshf901

喵星大佬 发表于 2021-3-16 23:13
你可以按你的盒子大小算算密度看看这玩意还是不是气体

6.6 kg m-3，还是气体。

wzkchem5

你是怎么知道能量异常的

卡开发发

大概率不能用，检查是否每一步SCF都收敛了，开idipol的情况很容易SCF难收敛。另外分子体系动能截断400eV可能会偏低。