24小时热门版块排行榜    

查看: 3686  |  回复: 9

dengban

新虫 (初入文坛)

[求助] 菜鸟急求,判断VASP是否已经计算结束

用VASP计算应力,不知道为什么在OUTCAR里出现了多次“ reached required accuracy - stopping structural energy minimisation”,这样正常么?还有在OSZICAR里为什么只有一个离子步啊,我设置的EDIFFG=default(EDIFF×10),也就是说以两离子步能量差为判据判断收敛啊。一个离子步怎么判断的出来啊
这种情况到底算不算计算完成了呀,我看别人的OSZICAR最后一行好像都会提示“reached required accuracy”的
求高手~
INCAR:
LREAL=.FALSE.
LWAVE= .False.
#ICHARG= 11
ISMEAR = 1
SIGMA= 0.2
#EDIFFG= -0.0001
POTIM = 0.05
#SYMPREC = 1e-3
ISIF=2
IBRION= 2
NSW=20
#NFREE = 2
#ISPIN= 2
#NPAR     = 8
#IALGO    = 48
#LPLANE  = .True.
LCHARG=.TRUE.
ENCUT= 600 eV
PREC=HIGH
grep 'reached' OUTCAR:
------------------------ aborting loop because EDIFF is reached ----------------------------------------
reached required accuracy - stopping structural energy minimisation
------------------------ aborting loop because EDIFF is reached ----------------------------------------
reached required accuracy - stopping structural energy minimisation
------------------------ aborting loop because EDIFF is reached ----------------------------------------
reached required accuracy - stopping structural energy minimisation
------------------------ aborting loop because EDIFF is reached ----------------------------------------
reached required accuracy - stopping structural energy minimisation
------------------------ aborting loop because EDIFF is reached ----------------------------------------
reached required accuracy - stopping structural energy minimisation
------------------------ aborting loop because EDIFF is reached ----------------------------------------
reached required accuracy - stopping structural energy minimisation
OSZICAR:
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.146960014218E+03    0.14696E+03   -0.60121E+03   780   0.107E+03
DAV:   2    -0.395589104059E+01   -0.15092E+03   -0.14257E+03  1393   0.178E+02
DAV:   3    -0.229532452180E+02   -0.18997E+02   -0.18836E+02  1015   0.830E+01
DAV:   4    -0.233993104206E+02   -0.44607E+00   -0.44188E+00  1049   0.156E+01
DAV:   5    -0.234127203622E+02   -0.13410E-01   -0.13388E-01   986   0.266E+00    0.241E+01
DAV:   6    -0.222253427034E+02    0.11874E+01   -0.55098E+00   936   0.201E+01    0.875E+00
DAV:   7    -0.220454482390E+02    0.17989E+00   -0.30426E+00   936   0.172E+01    0.708E+00
DAV:   8    -0.217970488728E+02    0.24840E+00   -0.38976E-01   925   0.647E+00    0.151E+00
DAV:   9    -0.217969189812E+02    0.12989E-03   -0.71349E-02   870   0.209E+00    0.581E-01
DAV:  10    -0.217929921273E+02    0.39269E-02   -0.10763E-02   930   0.107E+00    0.139E-01
DAV:  11    -0.217924641875E+02    0.52794E-03   -0.10912E-03   895   0.326E-01    0.101E-01
DAV:  12    -0.217923696827E+02    0.94505E-04   -0.17703E-04   644   0.114E-01
   1 F= -.21792370E+02 E0= -.21793217E+02  d E =-.217924E+02
~
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Lathander

新虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2012-11-03 13:20:20
体系高度对称,因为对称合力严格为0

[ 发自手机版 http://muchong.com/3g ]
2楼2012-11-03 07:28:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dengban

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by Lathander at 2012-11-03 07:28:20
体系高度对称,因为对称合力严格为0

所以这个到底算计算完成了么
还有我的EDIFFG不是负值,应该是按照两步迭代的能量差来判断收敛的啊,和合力为零没什么关系吧,而且OSZICAR里只写了一个离子步
3楼2012-11-03 08:59:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

T_liyi_T

铜虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2012-11-03 13:20:29
精度达到要求了吧,计算完成的话会有个文件(output.xxxxxxxxx)描述任务信息的。或者看看还在不在计算就好了.
4楼2012-11-03 09:30:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dengban

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by T_liyi_T at 2012-11-03 09:30:55
精度达到要求了吧,计算完成的话会有个文件(output.xxxxxxxxx)描述任务信息的。或者看看还在不在计算就好了.

你好,谢谢你的回答,我查过了,没有output.XXX这个文件,不过确实没在计算了
我再总结下我的问题好了
1、 为何采用能量判据,OSZICAR里只显示了一个离子步;在OSZICAR里没有显示“reached required accuracy”,是否表示计算未完成
2、OUTCAR里多次出现“reached required accuracy”是否正常
5楼2012-11-03 10:49:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

T_liyi_T

铜虫 (小有名气)

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2012-11-04 19:47:39
我也菜鸟一枚。不客气
我也看了看OSZICAR这个文件里面没有reached required accuracy的信息,它记录了每次计算的过程,一般不会只有一个离子步的。有可能是计算的时候因机器出现错误,退出来了。可以重新在算一下。
POTIM = 0.05这个参数,是不是有些小了,默认是0.5,这个小了10倍,在结构优化的时候可能每一步之间的变化也不大,就多次出现了达到精读。
不过我算的时候每次任务完了,会有output.XXXXX的文件,这个里面和OSZICAAR的内容差不多。
6楼2012-11-04 13:34:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dengban

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by T_liyi_T at 2012-11-04 13:34:03
我也菜鸟一枚。不客气
我也看了看OSZICAR这个文件里面没有reached required accuracy的信息,它记录了每次计算的过程,一般不会只有一个离子步的。有可能是计算的时候因机器出现错误,退出来了。可以重新在算一下。 ...

设potim这个参数主要是怕在用CG算法的时候不收敛,不过这个INCAR里我把算法改成quasi-newton了,没有再把它改大,也是希望搜索的时候能精细点。
我现在纠结的就是就OSZICAR的信息来看,只有一个离子步。。。真不知道这个结果算是收敛没有。。。
7楼2012-11-05 01:15:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

T_liyi_T

铜虫 (小有名气)

OSZICAR一步的情况是正常的。
8楼2012-11-05 14:05:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dengban

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by T_liyi_T at 2012-11-05 14:05:19
OSZICAR一步的情况是正常的。

但是一步的话怎么会有能量差呢?EDIFFG的收敛条件是怎么达到的呢?
9楼2012-11-06 07:25:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
本帖内容被屏蔽

10楼2022-06-19 16:22:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dengban 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302求调剂 +7 呼呼呼。。。。 2026-03-17 7/350 2026-03-18 08:00 by ms629
[考研] 303求调剂 +3 睿08 2026-03-17 5/250 2026-03-17 22:01 by 睿08
[考研] 299求调剂 +4 △小透明* 2026-03-17 4/200 2026-03-17 20:09 by peike
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 274求调剂0856材料化工 +13 z2839474511 2026-03-11 14/700 2026-03-17 16:51 by share_joy
[考研] 求调剂,总分315,考的生物医药,一志愿湖南师范大学。调剂到任何专业都可以 +4 小丁想进步 2026-03-11 5/250 2026-03-17 16:05 by 外星文明
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
信息提示
请填处理意见