24小时热门版块排行榜    

查看: 3685  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601专硕,总分342求调剂,地区不限 +4 share_joy 2026-03-16 4/200 2026-03-17 21:32 by hmn_wj
[考研] 0703化学调剂 ,六级已过,有科研经历 +8 曦熙兮 2026-03-15 8/400 2026-03-17 20:31 by xilongliang
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 材料专硕306英一数二 +7 z1z2z3879 2026-03-16 9/450 2026-03-17 17:31 by ccjequ
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 5/250 2026-03-17 17:27 by Little-xue
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 312求调剂 +4 陌宸希 2026-03-16 5/250 2026-03-17 17:09 by ruiyingmiao
[考研] 274求调剂0856材料化工 +13 z2839474511 2026-03-11 14/700 2026-03-17 16:51 by share_joy
[考研] 267一志愿南京工业大学0817化工求调剂 +6 SUICHILD 2026-03-12 6/300 2026-03-17 09:24 by 雾散后相遇lc
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 0703化学调剂,求各位老师收留 +8 秋有木北 2026-03-14 8/400 2026-03-16 15:21 by 哦哦123
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
信息提示
请填处理意见