24小时热门版块排行榜    

查看: 4125  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +4 Tide man 2026-08-12 5/250 2026-08-19 07:41 by 201120518022
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 什么时候开奖? +4 CrisMessi 2026-08-18 4/200 2026-08-19 00:53 by gatelove
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 27/1350 2026-08-18 21:44 by Noways
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见