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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】vasp静态计算能量值错误
汕头大学海洋科学接受调剂
221/1返回列表
查看: 6102  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] RY：中国产出的科学垃圾论文，绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 生物学调剂 +7 纸扇zhishan 2026-04-13 7/350 2026-04-14 14:21 by jyl0317
[考研] 335求调剂 +19 想上岸呀！！ 2026-04-12 20/1000 2026-04-14 13:32 by 第一天好
[考研] 各位老师好，求调剂，本科211，一志愿天津大学生物与医药学硕，差两名录取。 +10 路六六jjj 2026-04-13 10/500 2026-04-14 07:08 by babysonlkd
[考研] 22专硕求调剂 +8 haoyun上岸 2026-04-11 10/500 2026-04-13 22:15 by zhq0425
[考研] B区0809 ，数一英一，290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 一志愿西交机械专硕求调剂 +9 求上岸的小王 2026-04-10 9/450 2026-04-13 16:08 by jiangguiquan11
[考研] 一志愿中国科学院上海有机所，有机化学356分找调剂 +12 Nadiums 2026-04-09 12/600 2026-04-13 08:50 by lhj2009
[考研] 。 +10 李多米lee. 2026-04-12 11/550 2026-04-12 22:58 by yuyin1233
[考研] 求调剂，262机械专硕 +8 嗯yyl 2026-04-08 8/400 2026-04-12 02:31 by 秋豆菜芽
[找工作] 山东高校教师考核超级无底线，员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 085400 328分 求调剂 +10 喂你一个大橙子 2026-04-09 14/700 2026-04-11 19:53 by lqspecial
[考研] 广东省 085601 329分求调剂 +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[考研] 材料类284调剂 +40 想换手机不想解� 2026-04-08 48/2400 2026-04-10 23:28 by 314126402
[考研] 一志愿中南大学物理学，英一66，求调剂 +4 长烟旖旎 2026-04-08 5/250 2026-04-10 10:31 by 颖果儿
[考研] 311求调剂 +6 surte 2026-04-08 13/650 2026-04-09 14:00 by surte
[考研] 0860004 求调剂 309分 +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by 啊李999
[考研] 334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 09:28 by wdyheheeh
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
信息提示
关闭
请填处理意见
关闭