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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】vasp静态计算能量值错误
221/1返回列表
查看: 6083  |  回复: 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还有调剂机会吗 +25 xaw. 2026-04-11 25/1250 2026-04-12 14:39 by 豆酱蔬果沙拉
[考研] 272分材料子求调剂 +38 Loy0361 2026-04-10 49/2450 2026-04-12 12:46 by BruceLiu320
[考研] 调剂 +25 不逢春 2026-04-07 26/1300 2026-04-12 11:53 by 大力水手力大无�
[考研] 调剂 +10 月@163.com 2026-04-11 10/500 2026-04-12 09:14 by zhouyuwinner
[考研] 考研求调剂 +3 ban班小七 2026-04-11 3/150 2026-04-11 20:48 by may_新宇
[考研] 工科273调剂 +6 X1999 2026-04-09 7/350 2026-04-11 10:23 by zhq0425
[考研] 机械专硕270求调剂，接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
[材料工程] 材料调剂推荐 +8 蛋糕x2 2026-04-07 8/400 2026-04-10 23:13 by Ftglcn90
[考研] 本科西工大 324求调剂 +4 wysyjs25 2026-04-10 4/200 2026-04-10 20:00 by 来看流星雨10
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 一志愿211 0703化学 346分求调剂 +22 土豆er? 2026-04-09 23/1150 2026-04-10 10:58 by 高维春
[考研] 环境专硕调剂 +16 会说话的肘子 2026-04-06 16/800 2026-04-10 10:30 by asy1wn
[考研] 生物学求调剂 一志愿沪9，326分 +7 刘墨墨 2026-04-06 7/350 2026-04-10 08:11 by kangsm
[考研] 297求调剂 +27 GENJIOW 2026-04-07 30/1500 2026-04-09 23:20 by wolf97
[考研] 085404，334分，求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 一志愿南昌大学，085600，344分求调剂 +11 调剂上岸玘 2026-04-05 12/600 2026-04-08 16:17 by luoyongfeng
[考研] 化学0703-一志愿211-338分求调剂 +10 vants 2026-04-05 11/550 2026-04-08 16:02 by screening
[考研] 生物医药调剂｜SCI中科院三区一作+多项科研成果 +8 likangxing 2026-04-07 11/550 2026-04-08 00:02 by lys0704
[考研] 一志愿河北工业大学材料工程，初试344求专硕调剂 +6 15933906766 2026-04-05 6/300 2026-04-06 13:21 by 无际的草原
信息提示
关闭
请填处理意见
关闭