24小时热门版块排行榜    

查看: 6020  |  回复: 21
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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 ]
回复此楼

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

量子化学基础

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注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文件的最后面部分。
15楼2010-11-16 09:41:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 22 个回答

swx0789

铜虫 (小有名气)



nkleof(金币+2):谢谢参与
我不能回答你的问题。但我弱弱的问下,你的能量为什么取toten而不是取减掉熵的能量。
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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程085602 305分求调剂 +13 RichLi_ 2026-03-25 13/650 2026-03-26 14:35 by xiongkun222
[考研] 0856求调剂 +8 zhn03 2026-03-25 9/450 2026-03-26 13:42 by zzll406
[考研] 材料调剂 5+3 想要一壶桃花水 2026-03-25 8/400 2026-03-26 09:51 by 王小欠i
[考研] 308求调剂 +5 墨墨漠 2026-03-25 5/250 2026-03-25 22:19 by 544594351
[考研] 求调剂 +3 QiMing7 2026-03-25 3/150 2026-03-25 21:13 by 给你你注意休息
[考研] 333求调剂 +6 87639 2026-03-21 11/550 2026-03-25 16:17 by 87639
[考研] 材料学求调剂 +6 Stella_Yao 2026-03-20 6/300 2026-03-25 00:37 by baoball
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 一志愿华东理工大学081700,初试分数271 +5 kotoko_ik 2026-03-23 6/300 2026-03-24 10:29 by 学术搬砖er
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 298求调剂 +8 上岸6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 一志愿东华大学化学070300,求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
信息提示
请填处理意见