计算化学公社

标题: 计算得到的气相分子能量是否可用？ [打印本页]

作者
Author: zhaoshf901 时间: 2021-3-16 22:34
标题: 计算得到的气相分子能量是否可用？
用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

作者
Author: 喵星大佬 时间: 2021-3-16 23:13
你可以按你的盒子大小算算密度看看这玩意还是不是气体

作者
Author: zhaoshf901 时间: 2021-3-17 17:29

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

6.6 kg m-3，还是气体。

作者
Author: wzkchem5 时间: 2021-3-17 18:12
你是怎么知道能量异常的

作者
Author: 卡开发发 时间: 2021-3-17 20:19
本帖最后由卡开发发于 2021-3-17 20:21 编辑

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

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