应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 明明达到了收敛标准，却还在计算
51/1返回列表
查看: 1455  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yysskk

银虫 (小有名气)

[求助] 明明达到了收敛标准，却还在计算 已有1人参与

今天在弛豫一个体系的时候，无意中发现明明早已达到收敛标准了，计算却还在继续。
这是INCAR文件：
ISTART = 0 ; ICHARG= 2
ENCUT = 400
ISMEAR = 2 ; SIGMA = 0.2
dynamic:
IBRION = 1
NSW = 60
POTIM = 0.2
EDIFFG = -0.05
LREAL = A
ALGO=Fast
NSIM=4
NCORE=24
收敛标准是力小于0.05eV
计算在第25步达到收敛标准终止，这是OUTCAR中第20步的结果：
step: 20
     9.46384       7.01935       4.17889         0.00326      -0.00108      -0.00008
    -2.71429       4.70215       4.17596         0.00494      -0.00182      -0.00012
    12.16736       2.34078       4.18754         0.00222      -0.00008      -0.00075
    -1.36790       7.01668       4.18510         0.00350      -0.00327       0.00066
    -0.05056       4.70705       4.16075        -0.00173       0.00096      -0.00206
    -0.01194       0.02207       4.18084        -0.00111       0.00130       0.00214
     1.31064       7.06307       4.17596         0.00213      -0.00165      -0.00580
     1.33486       2.31176       4.17722         0.00019      -0.00219       0.00119
     2.69068       0.00669       4.17737         0.00416      -0.00086       0.00015
     2.65993       4.66485       4.30412         0.00050      -0.00111       0.00028
     4.06083       2.32909       4.18447        -0.00073      -0.00346       0.00002
     4.06516       7.03286       4.18795        -0.00056       0.00212       0.00197
     5.43410       4.68419       4.17705         0.00031       0.00112       0.00459
     5.40923       0.01108       4.17625         0.00113      -0.00065       0.00248
     6.76245       2.33950       4.17716         0.00157      -0.00248       0.00280
     6.76425       7.02348       4.18846         0.00104       0.00034       0.00590
     8.10727       0.01938       4.17374        -0.00038      -0.00301      -0.00017
     8.11014       4.69639       4.17638         0.00089      -0.00236       0.00479
     9.47190       2.33668       4.17668        -0.00075       0.00036       0.00102
    10.81960       0.00401       4.17358        -0.00140       0.00045       0.00189
     8.10158       7.81695       2.08881        -0.00076       0.00045       0.00001
     9.45834       5.47559       2.08584         0.00150      -0.00020       0.00088
    -2.71358       7.82110       2.08338         0.00161       0.00036       0.00296
    10.81672       3.13254       2.08894         0.00151      -0.00058       0.00362
    -1.36258       5.47446       2.09050         0.00265      -0.00028       0.00764
    12.16647       0.78915       2.08912         0.00102      -0.00024       0.00154
    -0.00961       3.12356       2.07565         0.00035       0.00058       0.00685
    -0.01004       7.82529       2.09491        -0.00048       0.00152       0.00579
     1.35284       5.48351       2.08345         0.00013      -0.00046       0.00267
     1.34166       0.79285       2.09358         0.00064       0.00051       0.00303
     2.69559       3.15225       2.09265         0.00113      -0.00282       0.00124
     2.70090       7.82003       2.08123         0.00289       0.00210       0.00194
     4.05866       5.47888       2.07023         0.00368       0.00065      -0.00246
     4.04886       0.79908       2.09502         0.00125      -0.00149      -0.00024
     5.40228       7.81472       2.09643         0.00218       0.00043      -0.00091
     5.40089       3.13419       2.08366         0.00255      -0.00142      -0.00177
     6.75627       0.79039       2.08640         0.00132      -0.00163       0.00093
     6.75654       5.47721       2.10119         0.00249      -0.00157      -0.00670
     8.11037       3.13243       2.08893         0.00064      -0.00070      -0.00254
     9.46172       0.79076       2.08680         0.00058      -0.00227       0.00135
    -4.05870       7.02990       0.00000        -0.05658       0.07477      -0.00860
    -2.70580       4.68660       0.00000        -0.03048      -0.06867      -0.04895
    -1.35290       7.02990       0.00000        -0.05607       0.09248      -0.03840
    -1.35290       2.34330       0.00000        -0.04192       0.09206      -0.09107
     0.00000       0.00000       0.00000        -0.01884      -0.05779      -0.00963
     0.00000       4.68660       0.00000        -0.01839      -0.04547       0.04743
     1.35290       2.34330       0.00000        -0.01700      -0.02258      -0.02606
     1.35290       7.02990       0.00000         0.03824      -0.03998      -0.06070
     2.70580       4.68660       0.00000         0.05156       0.02615       0.31155
     2.70580       0.00000       0.00000         0.04280      -0.06079       0.01474
     4.05870       7.02990       0.00000        -0.04736       0.04942      -0.02004
     4.05870       2.34330       0.00000        -0.00997       0.02566      -0.03565
     5.41160       4.68660       0.00000         0.00947       0.01142      -0.02814
     5.41160       0.00000       0.00000        -0.06625       0.00504      -0.13090
     6.76450       2.34330       0.00000        -0.00196       0.08162       0.01979
     6.76450       7.02990       0.00000        -0.00243      -0.00834      -0.09579
     8.11740       0.00000       0.00000         0.03133      -0.04949      -0.03659
     8.11740       4.68660       0.00000         0.01062      -0.05799       0.02902
     9.47030       2.34330       0.00000        -0.05510       0.04946      -0.01725
    10.82320       0.00000       0.00000        -0.03716      -0.05629       0.07876
     1.27989       5.11303       6.94810         0.03791       0.01757       0.00890
     3.39901       4.65296       8.68506         0.01873      -0.00449       0.00277
     3.27976       3.95499       7.32438         0.01314      -0.00630       0.01211
     2.39192       4.65821       6.33101         0.02605       0.00594       0.01073
     4.03176       4.05718       9.35500         0.01971       0.00588       0.00035
     3.85154       5.64880       8.58061         0.00708      -0.00725      -0.00048
     2.41455       4.77243       9.15292         0.02689      -0.00164       0.00806
     2.84688       2.94737       7.45616         0.00482      -0.00218       0.01336
     4.26642       3.80759       6.85854         0.00994      -0.01724       0.01615
     0.69760       5.58606       6.27939         0.02434       0.01690       0.00337
     3.05540       5.70801       5.72041         0.04079      -0.02348       0.02443
  free  energy   TOTEN  =      -567.724581 eV
  energy  without entropy=     -567.781882  energy(sigma->0) =     -567.738906
只有中间Z轴坐标为0的原子是非弛豫原子，其余所有原子早已达到收敛标准，为了防止我眼花看错了，我还特意写程序帮我查找，都没有发现弛豫原子上的力超过0.05ev的。
OUTCAR文件开头的参数里也写明了   EDIFFG = -.5E-01   stopping-criterion for IOM
POSCAR里面弛豫和非弛豫原子的标记我也检查过了没问题。
想问一下是我对收敛标准理解错了么？还有其它什么限制条件？
回复此楼

» 猜你喜欢
1楼 2015-01-24 20:26:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yysskk

银虫 (小有名气)
引用回帖:
2楼: Originally posted by ggdh at 2015-01-24 21:24:21
5.41160       0.00000       0.00000        -0.06625       0.00504      -0.1309
    -4.05870       7.02990       0.00000        -0.05658       0.07477      -0.00860
    -1.35290       7.02990        ...

那是非弛豫原子，不会计入收敛标准
赞一下 回复此楼
3楼2015-01-24 21:37:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

ggdh

专家顾问 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
ljw4010: 金币+2, 谢谢指导解答！ 2015-01-26 21:47:00
5.41160       0.00000       0.00000        -0.06625       0.00504      -0.1309
    -4.05870       7.02990       0.00000        -0.05658       0.07477      -0.00860
    -1.35290       7.02990       0.00000        -0.05607       0.09248      -0.03840
9.47030       2.34330       0.00000        -0.05510       0.04946      -0.01725

这有一堆超过0.05的
赞一下 回复此楼
2楼2015-01-24 21:24:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ggdh

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
yysskk: 金币+50, ★★★★★最佳答案 2015-01-24 21:53:36
ljw4010: 金币+2, 谢谢指导解答！ 2015-01-26 21:47:16
引用回帖:
2楼: Originally posted by ggdh at 2015-01-24 21:24:21
5.41160       0.00000       0.00000        -0.06625       0.00504      -0.1309
    -4.05870       7.02990       0.00000        -0.05658       0.07477      -0.00860
    -1.35290       7.02990        ...

我没有仔细审题
吧所有Z！=0的行的x y z 方向上受力平方和开根号后可得
0.00343517102922
0.00526596619814
0.00234463216731
0.00483513184929
0.00285623878554
0.00273892314606
0.00639526387259
0.00249965997688
0.00425061172068
0.00124919974384
0.00353622680268
0.00294769401397
0.00473482840238
0.00280174945347
0.00405651328113
0.00600059997
0.00303865101649
0.00541348316706
0.00131624465811
0.00239470248674
0.000883289307079
0.00175054277297
0.00338870181633
0.00396495901618
0.00809138430678
0.0018626862323
0.0068834148502
0.00600540589802
0.00271245276457
0.00313856655179
0.00328129547588
0.00406518142277
0.00447398033076
0.0019596428246
0.00240112473645
0.00341347330442
0.00229438444904
0.00731812817598
0.00271130964665
0.00270403402345
0.0427209901571
0.0194588257611
0.0189473401827
0.0287926900445
0.0205713635912
0.0101449149824
0.028119838193
0.0143692170977
0.0256289621327
0.0298228519763
0.0530279115184
我想这应该是判断标准
赞一下 回复此楼
4楼2015-01-24 21:39:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyqh

铁杆木虫 (著名写手)
是不是 能量达到收敛标准，而力还没有达到收敛标准呢？
赞一下 回复此楼
雨后方能见彩虹
5楼2015-01-25 01:18:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 275求调剂 +3 Micky11223 2026-03-25 4/200 2026-03-25 19:51 by 莫xiao
[考研] 07化学303求调剂 +4 睿08 2026-03-25 4/200 2026-03-25 19:15 by qingfeng258
[考研] 材料与化工304求B区调剂 +3 邱gl 2026-03-25 3/150 2026-03-25 19:03 by Ainin_
[考研] 考研调剂 +5 呼呼？~+123456 2026-03-24 5/250 2026-03-25 18:15 by xcjcqu
[考研] 求b区院校调剂 +4 周56 2026-03-24 5/250 2026-03-25 17:12 by yishunmin
[考研] 318求调剂 +3 plum李子 2026-03-23 3/150 2026-03-25 09:42 by 雾散后相遇lc
[考研] 318求调剂 +5 plum李子 2026-03-21 8/400 2026-03-25 09:26 by aa331100
[考研] 0854 考研调剂 招生了！AI 方向 +5 pk3725069 2026-03-19 17/850 2026-03-24 17:30 by zhouxuan..
[考研] 一志愿211 初试270分 求调剂 +5 谷雨上岸 2026-03-23 6/300 2026-03-24 16:32 by laoshidan
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗？ +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
[考研] 求调剂材料学硕080500，总分289分 5+3 @taotao 2026-03-19 21/1050 2026-03-23 10:17 by 冠c哥
[考研] 352求调剂 +3 大米饭！ 2026-03-22 3/150 2026-03-22 23:28 by king123！
[考研] 一志愿东华大学化学070300，求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 317求调剂 +12 申子申申 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[考研] 一志愿华中农业071010，总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 考研调剂 +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
信息提示
关闭
请填处理意见
关闭