24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 3757  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291分调剂 +8 上岸小莹加油 2026-04-09 9/450 2026-04-15 00:41 by zzzggc
[考研] 调剂求收留 +33 果然有我 2026-04-10 34/1700 2026-04-15 00:37 by zzzggc
[考研] 290调剂生物0860 +35 哇哈哈,。 2026-04-11 41/2050 2026-04-14 23:47 by Xurambo2014
[考研] 材料工程281还有调剂机会吗 +41 xaw. 2026-04-11 42/2100 2026-04-14 20:05 by honglizhao
[考研] 求调剂 +9 小聂爱学习 2026-04-11 13/650 2026-04-14 19:18 by Art1977
[考研] 通信工程求调剂!!! +4 zlb770521 2026-04-14 4/200 2026-04-14 18:19 by lbsjt
[考研] 材料工程085601,270求调剂 +42 @ASDF1234 2026-04-08 46/2300 2026-04-14 18:12 by 逍遥三郎
[考研] 一志愿沪9,326求生物学调剂 +10 刘墨墨 2026-04-13 10/500 2026-04-14 15:16 by zs92450
[考研] 材料085601调剂 +32 何润采123 2026-04-10 34/1700 2026-04-14 08:47 by 木木mumu~
[考研] 考研英一数一338分 +9 长江大学东校区 2026-04-13 10/500 2026-04-14 00:41 by 王珺璞
[考研] B区0809 ,数一英一,290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 一志愿0807 数一英一 313 有没有二轮调剂 +12 emokidd 2026-04-08 13/650 2026-04-13 08:32 by lhj2009
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 346分,工科0854求调剂,专硕 +6 moser233 2026-04-12 7/350 2026-04-12 22:11 by fqwang
[考研] 0854调剂 +12 长弓傲 2026-04-09 13/650 2026-04-12 09:56 by 逆水乘风
[考研] 求调剂,262机械专硕 +8 嗯yyl 2026-04-08 8/400 2026-04-12 02:31 by 秋豆菜芽
[考研] 调剂求助 +6 果然有我 2026-04-11 7/350 2026-04-11 16:22 by 明月此时有
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
[论文投稿] mdpi小修rvr时间四五天了 20+3 哈哈high 2026-04-08 5/250 2026-04-10 16:02 by 北京莱茵润色
[考研] 296求调剂 +6 汪!?! 2026-04-08 6/300 2026-04-10 11:02 by mattzhming
信息提示
请填处理意见