应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】vasp静态计算能量值错误
汕头大学海洋科学接受调剂
221/1返回列表
查看: 6099  |  回复: 21
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

nkleof

银虫 (正式写手)

[交流] 【求助】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 ]
回复此楼

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

量子化学基础

» 猜你喜欢

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

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

查看全部散金贴
1楼 2010-11-14 21:06:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

shadowlisong

铁杆木虫 (正式写手)
★
nkleof(金币+2): 谢谢参与
你可以试试在结构优化的时候就加U，看看是否有区别。出现ZHEGV 可能会是LAGO的影响 切换到38 也就是Normal试试，我出这个错误的时候更多的是结构不太合理。
赞一下(1人) 回复此楼
21楼2015-08-20 16:43:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

swx0789

铜虫 (小有名气)
★
nkleof(金币+2):谢谢参与
我不能回答你的问题。但我弱弱的问下，你的能量为什么取toten而不是取减掉熵的能量。
赞一下(1人) 回复此楼
2楼2010-11-14 21:22:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)
引用回帖:
Originally posted by swx0789 at 2010-11-14 21:22:37:
我不能回答你的问题。但我弱弱的问下，你的能量为什么取toten而不是取减掉熵的能量。

帖子里列出的是电子步迭代的能量变化，在达到收敛后，TOTEN值和减去entropy之后的能量值是一样的，关于这点我也不了解具体原理，不过如果TOTEN值和减去entropy后的能量值不一样应该是说明计算有问题吧？希望多给建议。
赞一下 回复此楼
3楼2010-11-14 21:28:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swx0789

铜虫 (小有名气)
我也是个初学者。祝你好运！
赞一下 回复此楼
4楼2010-11-14 21:58:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)
★
nkleof(金币+2):谢谢参与
我也路过，友情帮顶！希望尽早得到解决，祝愿你！！
赞一下(1人) 回复此楼
5楼2010-11-14 22:13:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustbwgs

木虫 (正式写手)

可以把问题写得更详细点吗？？

★ ★
nkleof(金币+2):谢谢参与
zzy870720z(金币+1):谢谢回复，呵呵 2010-11-15 08:05:37
可以提供一下你计算静态能和优化时的INCAR吗，你在哪上面设置+U呢，体系有多大啊？
赞一下(2人) 回复此楼
6楼2010-11-14 23:45:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qpf398632826

木虫 (小有名气)
★
nkleof(金币+2):谢谢参与
帮你顶下吧，希望能让更多的人看到来帮你解决
赞一下(1人) 回复此楼
7楼2010-11-15 00:20:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)
引用回帖:
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变化。
赞一下 回复此楼
8楼2010-11-15 08:50:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

真实的自我

金虫 (小有名气)
★ ★
nkleof(金币+2):谢谢参与
sunyang1988(金币+1):谢谢交流 2010-11-17 21:19:17
引用回帖:
Originally posted by nkleof at 2010-11-14 21:06:19:
优化了一个结构后，进行静态计算，+U设置为零，发现电子步迭代时的能量值都非常大，如下：
  free energy    TOTEN  =   -386068.93443124 eV
  free energy    TOTEN  =   -468284.83293180 eV
  free energy ...

你是不是你的收敛标准太高了？
赞一下(2人) 回复此楼
9楼2010-11-15 09:08:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

harrymia

禁虫 (小有名气)
★ ★
nkleof(金币+2):谢谢参与
sunyang1988(金币+1):谢谢交流 2010-11-17 21:19:30
本帖内容被屏蔽
10楼2010-11-15 11:13:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)
引用回帖:
Originally posted by harrymia at 2010-11-15 11:13:32:
我也是初学者，想问下静态计算能量不是看OSZICAR 最后一行的那个能量吗？这里列出的OUTCAR 中的free energy TOTEN是做什么用的？

默认的EDIFF是1E-4吧，你的设置好高啊，会不会是这里的问题，我不确定啊，初学 ...

应该不是精度设置的问题，我设置的精度也只是比默认精度高一个数量级而已，不算高的。
赞一下 回复此楼
11楼2010-11-15 19:24:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

harrymia

禁虫 (小有名气)
本帖内容被屏蔽
12楼2010-11-16 08:24:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinhua1983

金虫 (小有名气)
★
nkleof(金币+2):谢谢参与
看能量看的不是自由能，而是能量个差值！
赞一下(1人) 回复此楼
13楼2010-11-16 08:56:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skyscorpioa

新虫 (初入文坛)
★ ★
nkleof(金币+2):谢谢参与
zzy870720z(金币+1):谢谢建议 2010-11-17 21:19:39
经常遇到这种问题，一直不知道原因，可能和用的编译器和lapack库有关系
如果你用的是pbe，可以换成pw91试试，有时候就可以的
赞一下(2人) 回复此楼
14楼2010-11-16 08:59:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)
★ ★
cenwanglai(金币+2):谢谢回复！ 2010-11-17 16:26:45
引用回帖:
Originally posted by harrymia at 2010-11-16 08:24:49:


还是想请问静态计算能量是不是看OSZICAR 最后一行的那个能量？这里列出的OUTCAR 中的free energy TOTEN是做什么用的？多谢啦

一般要查看能量我都是找OUTCAR的，OSZICAR的数值和OUTCAR应该是一样，可能有小差别，我觉得还是应该以OUCAR为准。
OUTCAR里面会记录每一个电子步优化的能量结果，就是我前面列出的free energy TOTEN，达到收敛后会有另外一个free  energy TOTEN，以及energy without entropy，是最终的基态能量，在OUTCAR文件的最后面部分。
赞一下(1人) 回复此楼
15楼2010-11-16 09:41:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)
引用回帖:
Originally posted by skyscorpioa at 2010-11-16 08:59:12:
经常遇到这种问题，一直不知道原因，可能和用的编译器和lapack库有关系
如果你用的是pbe，可以换成pw91试试，有时候就可以的

多谢建议，vasp官网论坛上也提到过可能的原因包括lapack的问题。
换赝势可能也是一个办法，不过我的大部分计算都用的是PBE，只这一个换赝势的话就不自洽了。
赞一下 回复此楼
16楼2010-11-16 09:43:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustbwgs

木虫 (正式写手)
★
cenwanglai(金币+1):谢谢回复！ 2010-11-17 16:27:04
收敛标准高可能使体系在要求的精度内收敛不了，但不会因为收敛精度高而出错
赞一下(1人) 回复此楼
17楼2010-11-17 16:04:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy870720z

荣誉版主 (文坛精英)
★
nkleof(金币+2):谢谢参与
我也遇到过一次
建议你看一下最后的结构
也有可能是结构不合理
赞一下(1人) 回复此楼
18楼2010-11-17 21:22:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weio_china

铁杆木虫 (正式写手)
★
nkleof(金币+2): 谢谢参与
新手学习
回复此楼
19楼2012-05-04 09:50:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zuocuiping

木虫 (职业作家)
★
nkleof(金币+2): 谢谢参与
楼主问题解决了吗？
我也遇到这个问题了谢谢
赞一下 回复此楼
20楼2015-08-20 14:46:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

din5g

铜虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
21楼: Originally posted by shadowlisong at 2015-08-20 16:43:06
你可以试试在结构优化的时候就加U，看看是否有区别。出现ZHEGV 可能会是LAGO的影响 切换到38 也就是Normal试试，我出这个错误的时候更多的是结构不太合理。

怎样看出是结构不合理啊？以及怎样解决？ 多谢！希望不吝赐教！
赞一下 回复此楼
22楼2018-04-15 21:07:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nkleof 的主题更新
221/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290求调剂 +18 luoziheng 2026-04-10 18/900 2026-04-14 09:29 by 带壳的鸡蛋
[考研] 材料相关专业344求调剂双非工科学校或课题组 +19 hualkop 2026-04-12 20/1000 2026-04-14 07:02 by laoshidan
[考研] 化工学硕294分，求导师收留 +27 yzyzx 2026-04-12 30/1500 2026-04-14 05:42 by 西北望—风沙
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 材料299专硕求调剂 +13 +21 2026-04-09 13/650 2026-04-13 14:16 by 张zhihao
[考研] 335求调剂 +18 想上岸呀！！ 2026-04-12 19/950 2026-04-13 14:15 by 张zhihao
[考研] 290调剂生物0860 +28 哇哈哈，。 2026-04-11 31/1550 2026-04-13 01:16 by 幸免 ..
[考研] 280求调剂 +13 wzzz王 2026-04-09 13/650 2026-04-12 00:31 by 勇攀高峰0126
[考研] 296求调剂 +8 汪！？！ 2026-04-09 8/400 2026-04-11 21:02 by 逆水乘风
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 269求调剂 +11 啊啊我我 2026-04-07 11/550 2026-04-11 16:45 by vgtyfty
[考研] 求调剂，一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 求调剂 +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[考研] 调剂 +4 电气300求调剂不 2026-04-08 7/350 2026-04-11 10:44 by 紫曦紫棋
[考研] 085500求调剂材料 +10 易11122 2026-04-09 10/500 2026-04-11 10:39 by maddjdld
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 293求调剂 +6 勇远库爱314 2026-04-08 6/300 2026-04-11 10:08 by zhq0425
[考研] 考研调剂 +13 冰冰，，， 2026-04-07 13/650 2026-04-09 17:01 by Lilly_Li
[考研] 生物学学硕，初试351分，求调剂 +4 …～、王…～ 2026-04-08 5/250 2026-04-08 21:49 by limeifeng
[考研] 专硕085403，291分，有两篇专利，一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
信息提示
关闭
请填处理意见
关闭