24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

nkleof

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 2093.2
帖子: 813
在线: 165.9小时
虫号: 152581

[交流] 【求助】vasp静态计算能量值错误

优化了一个结构后，进行静态计算，+U设置为零，发现电子步迭代时的能量值都非常大，如下：
  free energy TOTEN  = -386068.93443124 eV
  free energy TOTEN  = -468284.83293180 eV
  free energy TOTEN  = -483184.31035624 eV
  free energy TOTEN  = -488841.30556393 eV
  free energy TOTEN  = -491852.20379074 eV
  free energy TOTEN  = -278816.17117876 eV
  free energy TOTEN  = -10667.33547639 eV
  free energy TOTEN  =    36658.42027096 eV
  free energy TOTEN  =    3020.73440924 eV
  free energy TOTEN  = -86463.86284949 eV
  free energy TOTEN  = -120769.97367841 eV
  free energy TOTEN  = -142957.25416115 eV
  free energy TOTEN  = -175335.77020944 eV
  free energy TOTEN  = -171597.62982132 eV
  free energy TOTEN  = -165579.41016696 eV
  free energy TOTEN  = -170144.20096520 eV
  free energy TOTEN  = -144730.46667812 eV
  free energy TOTEN  = -200129.89254844 eV
  free energy TOTEN  = -136862.06349103 eV
优化的结构经过检查并作为初始结构再次优化，得到仍是这个结构，应该不会有问题。奇怪的是加上一定的U值后计算没有出现异常，而且再回头重新进行+U为零的计算又能得到正常的总能了。
请问可能造成这种情况的原因有哪些？

今天再检查发现从以上电子步以后迭代不再进行下去了，出现很多警告信息：
WARNING in EDDRMM: call to ZHEGV failed
求解！

发现以优化结果的CHGCAR作为初始电荷密度进行静态计算后，不会出现上述问题。

[ Last edited by nkleof on 2010-11-16 at 09:59 ]

回复此楼

» 收录本帖的淘帖专辑推荐

量子化学基础

» 猜你喜欢

27秋申博自荐已经有3人回复
郑州大学申博已经有7人回复
关于核磁制样已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

vasp5.2 openmpi+intel+mkl编译，能量正值问题已经有7人回复
vasp计算：优化能收敛，但同样的设置静态计算不收敛已经有3人回复
vasp计算静态出现错误！已经有5人回复
vasp计算在进行静态计算时，出错。已经有3人回复
【求助】vasp静态计算出现的WARNING 已经有18人回复
【求助】vasp做优化和静态计算的截断能一样吗？已经有4人回复
【求助】GaN的费米能是多少？要用VASP计算GaN，不知道算的对不对，请知道的给个参考值已经有4人回复
【求助】VASP怎样精确计算出费米能级处的态密度值？已经有22人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-11-14 21:06:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

nkleof

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 2093.2
帖子: 813
在线: 165.9小时
虫号: 152581

引用回帖:

Originally posted by ustbwgs at 2010-11-14 23:45:27:
可以提供一下你计算静态能和优化时的INCAR吗，你在哪上面设置+U呢，体系有多大啊？

在静态计算时+U，优化时没有+U。
具体设置如下：

PREC = High
LREAL = A
ALGO = Fast
NPAR = 4
NSIM = 8
ISTART = 0
ICHARG = 0
LMAXMIX = 4

NELMIN = 4
LORBIT = 12
EDIFF = 1E-5
EMAX = 20
EMIN = -20
NEDOS = 2000
ISMEAR = -5

ISPIN = 2
MAGMOM = 16*0 34*4 72*0

LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = 2 2 -1
LDAUU = 0 0 0
LDAUJ = 0 0 0
LDAUPRINT = 2

U值在0~2变化。