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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】vasp静态计算能量值错误
221/1返回列表
查看: 5963  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 中国科学技术大学材料与化工281求调剂，有科研和获奖经历 +6 wsxw 2026-03-02 6/300 2026-03-02 23:49 by 半截的诗0927
[考研] 11408，学硕276求调剂 +3 崔wj 2026-03-02 5/250 2026-03-02 22:14 by 255671
[考研] 调剂 +5 13853210211 2026-03-02 6/300 2026-03-02 22:10 by 无际的草原
[考研] 290分材料工程085601求调剂 数二英一 +8 llx0610 2026-03-02 9/450 2026-03-02 22:09 by 无际的草原
[考研] 085600 英一数二272求调剂 5+6 vida_a 2026-03-01 18/900 2026-03-02 22:08 by sunny81
[考研] 理学，工学，农学调剂，少走弯路，这里欢迎您！ +5 likeihood 2026-03-02 8/400 2026-03-02 20:39 by ji493940045
[考研] 285求调剂 +9 满头大汗的学生 2026-02-28 9/450 2026-03-02 20:29 by hypershenger
[考研] 材料270求调剂 6+4 Eiiiio 2026-03-01 7/350 2026-03-02 20:17 by hypershenger
[考研] 0856求调剂285 +11 吕仔龙 2026-02-28 11/550 2026-03-02 20:15 by hypershenger
[考研] 302材料工程求调剂 +5 Doleres 2026-03-01 6/300 2026-03-02 19:53 by 张晓芳0105
[考研] 0805总分292，求调剂 +10 幻想之殇 2026-03-01 10/500 2026-03-02 19:51 by zhukairuo
[考研] 275求调剂 +7 明远求学 2026-03-01 7/350 2026-03-02 19:22 by zhukairuo
[考研] 化工京区271求调剂 +6 11ing 2026-03-02 6/300 2026-03-02 18:52 by caszguilin
[考研] 292求调剂 +7 yhk_819 2026-02-28 7/350 2026-03-02 12:43 by 无际的草原
[考研] 化工专硕342，一志愿大连理工大学，求调剂 +6 kyf化工 2026-02-28 7/350 2026-03-02 10:56 by 无际的草原
[基金申请] 成果系统访问量大，请一小时后再尝试。---NSFC啥时候好哦，已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 272求调剂 +6 田智友 2026-02-28 6/300 2026-03-01 21:40 by 公瑾逍遥
[考研] 304求调剂 +6 曼殊2266 2026-02-28 7/350 2026-03-01 15:14 by wjLi2017
[考研] 寻找调剂 +4 LYidhsjabdj 2026-02-28 4/200 2026-03-01 10:56 by sunny81
信息提示
关闭
请填处理意见
关闭