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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 明明达到了收敛标准，却还在计算
51/1返回列表
查看: 1460  |  回复: 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的回帖

zyqh

铁杆木虫 (著名写手)
是不是 能量达到收敛标准，而力还没有达到收敛标准呢？
赞一下 回复此楼
雨后方能见彩虹
5楼2015-01-25 01:18:52
已阅   回复此楼   关注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的回帖

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的回帖

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的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085602化学工程求调剂。 +3 平乐乐乐 2026-03-26 3/150 2026-03-26 16:37 by 3Strings
[考研] 334分 一志愿武理 材料求调剂 +4 李李不服输 2026-03-26 4/200 2026-03-26 16:00 by 不吃魚的貓
[考研] 材料科学与工程 317求调剂 +4 JKSOIID 2026-03-26 4/200 2026-03-26 15:58 by 不吃魚的貓
[考研] 299求调剂 +4 15188958825 2026-03-25 4/200 2026-03-25 22:56 by 418490947
[考研] 07化学303求调剂 +5 睿08 2026-03-25 5/250 2026-03-25 22:46 by 418490947
[考研] 材料与化工 322求调剂 +6 然11 2026-03-19 6/300 2026-03-25 18:37 by haxia
[考研] 0854人工智能方向招收调剂 +4 章小鱼567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
[考研] 284求调剂 +15 Zhao anqi 2026-03-22 15/750 2026-03-25 12:51 by wht0531
[考研] 材料学硕，求调剂 6+4 糖葫芦888ll 2026-03-22 9/450 2026-03-25 11:19 by greychen00
[考研] 311求调剂 +3 冬十三 2026-03-24 3/150 2026-03-24 21:31 by peike
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 一志愿南航材料专317分求调剂 +5 炸呀炸呀炸薯条 2026-03-23 5/250 2026-03-24 16:52 by 星空星月
[考研] 291求调剂 +3 HanBeiNingZC 2026-03-24 3/150 2026-03-24 16:34 by barlinike
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[考研] 291求调剂 +8 hhhhxn.. 2026-03-23 8/400 2026-03-23 23:15 by peike
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
信息提示
关闭
请填处理意见
关闭