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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】vasp静态计算能量值错误
汕头大学海洋科学接受调剂
221/1返回列表
查看: 6100  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 药学305求调剂 +4 玛卡巴卡boom 2026-04-11 4/200 2026-04-14 10:21 by dick_runner
[考研] 机械还有还有名额吗？太难了 +8 笑笑袁 2026-04-10 8/400 2026-04-14 08:44 by screening
[考研] 288求调剂，一志愿华南理工大学071005 +24 ioodiiij 2026-04-08 24/1200 2026-04-14 08:42 by nana_zhong
[考研] 各位老师好，求调剂，本科211，一志愿天津大学生物与医药学硕，差两名录取。 +10 路六六jjj 2026-04-13 10/500 2026-04-14 07:08 by babysonlkd
[考研] 297求调剂 +22 ORCHID1 2026-04-10 25/1250 2026-04-13 23:11 by pies112
[考研] 调剂 +11 月@163.com 2026-04-11 11/550 2026-04-13 22:52 by pies112
[考研] 339求调剂 +4 hanwudada 2026-04-12 4/200 2026-04-13 12:03 by 蓝云思雨
[考研] 266调剂 +10 daya sun 2026-04-07 11/550 2026-04-13 10:12 by fenglj492
[考研] 368化学求调剂 +14 wwwwabcde 2026-04-07 15/750 2026-04-13 08:36 by lhj2009
[考研] 电气专硕320求调剂 +6 小麻子111 2026-04-10 6/300 2026-04-12 10:54 by lemon6009
[考研] 求调剂 +11 翩翩一书生 2026-04-09 11/550 2026-04-11 19:57 by 逆水乘风
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[考研] 314求调剂 +18 xhhdjdjsjks 2026-04-09 19/950 2026-04-10 18:53 by HPUCZ
[考研] 一志愿矿大，材料工程专硕314分，0856可调都可以 +15 无懈可击的巨人 2026-04-09 15/750 2026-04-10 18:10 by hmn_wj
[考研] 282，电气工程专业，求调剂，不挑专业 +9 jggshjkkm 2026-04-10 9/450 2026-04-10 14:55 by 逆水乘风
[考研] 考研调剂-材料类-284 +28 想换手机不想解� 2026-04-08 28/1400 2026-04-09 20:08 by 倒数321?
[考研] 086000生物与医药调剂 +7 awwwwwooooo 2026-04-09 7/350 2026-04-09 13:31 by 北极159263
[考研] 085404，334分，求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 313求调剂 +3 十六拾陆 2026-04-07 3/150 2026-04-07 23:20 by lbsjt
信息提示
关闭
请填处理意见
关闭